Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions src/ert/storage/local_ensemble.py
Original file line number Diff line number Diff line change
Expand Up @@ -678,6 +678,9 @@ def save_parameters_numpy(
if complete_df is None:
complete_df = ds
else:
complete_df = complete_df.drop(
[c for c in ds.columns if c != "realization"], strict=False
)
complete_df = (
complete_df.join(ds, on="realization", how="left")
.unique(subset=["realization"], keep="first")
Expand Down
31 changes: 31 additions & 0 deletions tests/ert/unit_tests/storage/test_local_storage.py
Original file line number Diff line number Diff line change
Expand Up @@ -1224,6 +1224,37 @@ def test_that_permission_error_is_logged_in_load_index(snake_oil_storage, caplog
index_path.chmod(0o744) # restore permissions


def test_that_multiple_save_parameters_numpy_calls_overwrite_previous_values(tmp_path):
"""This test is to ensure that updating parameters in local ensemble does not
raise duplicate column errors when saving with numpy arrays multiple times, and
instead overwrites the previous values as expected.
"""
writer = open_storage(tmp_path, mode="w")
gen_kw_parameter = GenKwConfig(
name="some_param", distribution={"name": "normal", "mean": 10, "std": 0.1}
)
exp = writer.create_experiment(parameters=[gen_kw_parameter])
num_reals = 5
parameter_data_0 = np.array([0.0] * num_reals)
parameter_data_1 = np.array([1.1] * num_reals)
parameter_data_2 = np.array([2.2] * num_reals)
ens: LocalEnsemble = exp.create_ensemble(ensemble_size=num_reals, name="uniq_ens")
iens_active_index = np.array(range(num_reals))
ens.save_parameters_numpy(
parameter_data_0, gen_kw_parameter.name, iens_active_index
)
ens.save_parameters_numpy(
parameter_data_1, gen_kw_parameter.name, iens_active_index
)
ens.save_parameters_numpy(
parameter_data_2, gen_kw_parameter.name, iens_active_index
)
saved_ens = ens.load_scalar_keys([gen_kw_parameter.name]).to_dict()[
gen_kw_parameter.name
]
assert saved_ens.to_list() == [2.2 for _ in range(num_reals)]


@dataclass
class Ensemble:
uuid: UUID
Expand Down
Loading