desdeo 2.0.0__py3-none-any.whl → 2.1.0__py3-none-any.whl

desdeo/api/tests/test_models.py

@@ -1,20 +1,29 @@
 """Tests related to the SQLModels."""
 
+import numpy as np
 from sqlmodel import Session, select
 
 from desdeo.api.models import (
-    ArchiveEntryBase,
-    ArchiveEntryDB,
     Bounds,
     ConstantDB,
     ConstraintDB,
     DiscreteRepresentationDB,
+    ENautilusState,
     ExtraFunctionDB,
+    ForestProblemMetaData,
+    Group,
+    GroupIteration,
     InteractiveSessionDB,
+    NIMBUSClassificationState,
+    NIMBUSFinalState,
+    NIMBUSInitializationState,
+    NIMBUSSaveState,
     ObjectiveDB,
     PreferenceDB,
     ProblemDB,
+    ProblemMetaDataDB,
     ReferencePoint,
+    RepresentativeNonDominatedSolutions,
     RPMState,
     ScalarizationFunctionDB,
     SimulatorDB,
@@ -22,9 +31,13 @@ from desdeo.api.models import (
     TensorConstantDB,
     TensorVariableDB,
     User,
+    UserSavedSolutionDB,
     VariableDB,
 )
-from desdeo.mcdm import rpm_solve_solutions
+from desdeo.api.models.gdm.gnimbus import (
+    OptimizationPreference,
+)
+from desdeo.mcdm import enautilus_step, generate_starting_point, rpm_solve_solutions, solve_sub_problems
 from desdeo.problem.schema import (
     Constant,
     Constraint,
@@ -46,6 +59,7 @@ from desdeo.problem.testproblems import (
     dtlz2,
     momip_ti2,
     momip_ti7,
+    multi_valued_constraint_problem,
     nimbus_test_problem,
     pareto_navigator_test_problem,
     re21,
@@ -63,7 +77,7 @@ from desdeo.problem.testproblems import (
     spanish_sustainability_problem,
     zdt1,
 )
-from desdeo.tools import available_solvers
+from desdeo.tools import PyomoBonminSolver, available_solvers
 
 
 def compare_models(
@@ -135,7 +149,7 @@ def test_tensor_constant(session_and_user: dict[str, Session | list[User]]):
     # check that original added TensorConstant and fetched match
     assert db_tensor == from_db_tensor
 
-    from_db_tensor_dump = from_db_tensor.model_dump()
+    from_db_tensor_dump = from_db_tensor.model_dump(exclude={"id", "problem_id"})
     t_tensor_validated = TensorConstant.model_validate(from_db_tensor_dump)
 
     assert t_tensor_validated == t_tensor
@@ -159,7 +173,7 @@ def test_constant(session_and_user: dict[str, Session | list[User]]):
 
     assert db_constant == from_db_constant
 
-    from_db_constant_dump = from_db_constant.model_dump()
+    from_db_constant_dump = from_db_constant.model_dump(exclude={"id", "problem_id"})
     constant_validated = Constant.model_validate(from_db_constant_dump)
 
     assert constant_validated == constant
@@ -191,7 +205,7 @@ def test_variable(session_and_user: dict[str, Session | list[User]]):
 
     assert db_variable == from_db_variable
 
-    from_db_variable_dump = from_db_variable.model_dump()
+    from_db_variable_dump = from_db_variable.model_dump(exclude={"id", "problem_id"})
     variable_validated = Variable.model_validate(from_db_variable_dump)
 
     assert variable_validated == variable
@@ -224,7 +238,7 @@ def test_tensor_variable(session_and_user: dict[str, Session | list[User]]):
 
     assert db_t_variable == from_db_t_variable
 
-    from_db_t_variable_dump = from_db_t_variable.model_dump()
+    from_db_t_variable_dump = from_db_t_variable.model_dump(exclude={"id", "problem_id"})
     t_variable_validated = TensorVariable.model_validate(from_db_t_variable_dump)
 
     assert t_variable_validated == t_variable
@@ -264,7 +278,7 @@ def test_objective(session_and_user: dict[str, Session | list[User]]):
 
     assert db_objective == from_db_objective
 
-    from_db_objective_dump = from_db_objective.model_dump()
+    from_db_objective_dump = from_db_objective.model_dump(exclude={"id", "problem_id"})
     objective_validated = Objective.model_validate(from_db_objective_dump)
 
     assert objective_validated == objective
@@ -300,7 +314,7 @@ def test_constraint(session_and_user: dict[str, Session | list[User]]):
 
     assert db_constraint == from_db_constraint
 
-    from_db_constraint_dump = from_db_constraint.model_dump()
+    from_db_constraint_dump = from_db_constraint.model_dump(exclude={"id", "problem_id"})
     constraint_validated = Constraint.model_validate(from_db_constraint_dump)
 
     assert constraint_validated == constraint
@@ -333,7 +347,7 @@ def test_scalarization_function(session_and_user: dict[str, Session | list[User]
 
     assert db_scalarization == from_db_scalarization
 
-    from_db_scalarization_dump = from_db_scalarization.model_dump()
+    from_db_scalarization_dump = from_db_scalarization.model_dump(exclude={"id", "problem_id"})
     scalarization_validated = ScalarizationFunction.model_validate(from_db_scalarization_dump)
 
     assert scalarization_validated == scalarization
@@ -366,7 +380,7 @@ def test_extra_function(session_and_user: dict[str, Session | list[User]]):
 
     assert db_extra == from_db_extra
 
-    from_db_extra_dump = from_db_extra.model_dump()
+    from_db_extra_dump = from_db_extra.model_dump(exclude={"id", "problem_id"})
     extra_validated = ExtraFunction.model_validate(from_db_extra_dump)
 
     assert extra_validated == extra
@@ -395,7 +409,7 @@ def test_discrete_representation(session_and_user: dict[str, Session | list[User
 
     assert db_discrete == from_db_discrete
 
-    from_db_discrete_dump = from_db_discrete.model_dump()
+    from_db_discrete_dump = from_db_discrete.model_dump(exclude={"id", "problem_id"})
     discrete_validated = DiscreteRepresentation.model_validate(from_db_discrete_dump)
 
     assert discrete_validated == discrete
@@ -425,7 +439,7 @@ def test_simulator(session_and_user: dict[str, Session | list[User]]):
 
     assert db_simulator == from_db_simulator
 
-    from_db_simulator_dump = from_db_simulator.model_dump()
+    from_db_simulator_dump = from_db_simulator.model_dump(exclude={"id", "problem_id"})
     simulator_validated = Simulator.model_validate(from_db_simulator_dump)
 
     assert simulator_validated == simulator
@@ -464,6 +478,7 @@ def test_from_problem_to_d_and_back(session_and_user: dict[str, Session | list[U
         spanish_sustainability_problem(),
         zdt1(10),
         dtlz2(5, 3),
+        multi_valued_constraint_problem(),
         momip_ti2(),
         momip_ti7(),
         nimbus_test_problem(),
@@ -493,49 +508,84 @@ def test_from_problem_to_d_and_back(session_and_user: dict[str, Session | list[U
         assert compare_models(problem, problem_db)
 
 
-def test_archive_entry(session_and_user: dict[str, Session | list[User]]):
-    """Test that the archive works as intended."""
+def test_user_save_solutions(session_and_user: dict[str, Session | list[User]]):
+    """Test that user_save_solutions correctly saves solutions to the usersavedsolutiondb in the database."""
     session = session_and_user["session"]
     user = session_and_user["user"]
 
-    problem = dtlz2(n_variables=5, n_objectives=3)
-    problem_db = ProblemDB.from_problem(problem, user)
-
-    session.add(problem_db)
-    session.commit()
-    session.refresh(problem_db)
+    # Create test solutions with proper dictionary values
+    objective_values = {"f_1": 1.2, "f_2": 0.9}
+    variable_values = {"x_1": 5.2, "x_2": 8.0, "x_3": -4.2}
+
+    user_id = user.id
+    problem_id = 1
+
+    test_solutions = [
+        UserSavedSolutionDB(
+            name="Solution 1",
+            objective_values=objective_values,
+            variable_values=variable_values,
+            user_id=user_id,
+            problem_id=problem_id,
+            origin_state_id=1,
+        ),
+        UserSavedSolutionDB(
+            name="Solution 2",
+            objective_values=objective_values,
+            variable_values=variable_values,
+            solution_index=2,
+            user_id=user_id,
+            problem_id=problem_id,
+            origin_state_id=2,
+        ),
+    ]
 
-    variable_values = {"x_1": 0.3, "x_2": 0.8, "x_3": 0.1, "x_4": 0.6, "x_5": 0.9}
-    objective_values = {"f_1": 1.2, "f_2": 0.9, "f_3": 1.5}
-    constraint_values = {"g_1": 0.1}
-    extra_func_values = {"extra_1": 5000, "extra_2": 600000}
+    num_test_solutions = len(test_solutions)
 
-    archive_entry = ArchiveEntryBase(variable_values=variable_values, objective_values=objective_values)
+    # Create NIMBUSSaveState
+    save_state = NIMBUSSaveState(solutions=test_solutions)
 
-    archive_entry_db = ArchiveEntryDB.model_validate(
-        archive_entry,
-        update={
-            "user_id": user.id,
-            "problem_id": problem_db.id,
-            "constraint_values": constraint_values,
-            "extra_func_values": extra_func_values,
-        },
-    )
+    # Create StateDB with NIMBUSSaveState
+    state_db = StateDB.create(session, problem_id=problem_id, state=save_state)
 
-    session.add(archive_entry_db)
+    session.add(state_db)
     session.commit()
-    session.refresh(archive_entry_db)
-
-    from_db = session.get(ArchiveEntryDB, archive_entry_db.id)
-
-    assert from_db.user_id == user.id
-    assert from_db.problem_id == problem_db.id
-    assert from_db == user.archive[0]
-    assert compare_models(from_db.problem, problem_db)
-    assert from_db.variable_values == variable_values
-    assert from_db.objective_values == objective_values
-    assert from_db.constraint_values == constraint_values
-    assert from_db.extra_func_values == extra_func_values
+    session.refresh(state_db)
+
+    # Verify the solutions were saved
+    saved_solutions = session.exec(select(UserSavedSolutionDB)).all()
+    assert len(saved_solutions) == num_test_solutions
+
+    # Verify the content of the first solution
+    first_solution = saved_solutions[0]
+    assert first_solution.name == "Solution 1"
+    assert first_solution.objective_values == objective_values
+    assert first_solution.variable_values == variable_values
+    assert first_solution.origin_state_id == 1
+    assert first_solution.solution_index is None
+    assert first_solution.user_id == user.id
+    assert first_solution.problem_id == problem_id
+
+    # Verify the content of the second solution
+    second_solution = saved_solutions[1]
+    assert second_solution.name == "Solution 2"
+    assert second_solution.objective_values == objective_values
+    assert second_solution.variable_values == variable_values
+    assert second_solution.origin_state_id == 2
+    assert second_solution.solution_index == 2
+    assert second_solution.user_id == user.id
+    assert second_solution.problem_id == problem_id
+
+    # Verify state relationship
+    saved_state = session.exec(select(StateDB).where(StateDB.id == state_db.id)).first()
+    assert isinstance(saved_state.state, NIMBUSSaveState)
+    assert len(saved_state.state.solutions) == num_test_solutions
+
+    # Check that relationships are formed
+    session.refresh(user)
+
+    assert len(user.archive) == len(test_solutions)
+    assert len(session.get(ProblemDB, problem_id).solutions) == len(test_solutions)
 
 
 def test_preference_models(session_and_user: dict[str, Session | list[User]]):
@@ -579,9 +629,6 @@ def test_preference_models(session_and_user: dict[str, Session | list[User]]):
     assert from_db_bounds.user == user
     assert from_db_ref_point.user == user
 
-    assert from_db_bounds.solutions == []
-    assert from_db_ref_point.solutions == []
-
 
 def test_rpm_state(session_and_user: dict[str, Session | list[User]]):
     """Test the RPM state that it works correctly."""
@@ -609,30 +656,26 @@ def test_rpm_state(session_and_user: dict[str, Session | list[User]]):
         problem,
         asp_levels_1,
         scalarization_options=scalarization_options,
-        solver=available_solvers[solver],
+        solver=available_solvers[solver]["constructor"],
         solver_options=solver_options,
     )
 
     # create preferences
 
     rp_1 = ReferencePoint(aspiration_levels=asp_levels_1)
-    preferences = PreferenceDB(user_id=user.id, problem_id=problem_db.id, preference=rp_1)
-
-    session.add(preferences)
-    session.commit()
-    session.refresh(preferences)
 
     # create state
 
     rpm_state = RPMState(
+        preferences=rp_1,
         scalarization_options=scalarization_options,
         solver=solver,
         solver_options=solver_options,
         solver_results=results,
     )
 
-    state_1 = StateDB(
-        problem_id=problem_db.id, preference_id=preferences.id, session_id=isession.id, parent_id=None, state=rpm_state
+    state_1 = StateDB.create(
+        database_session=session, problem_id=problem_db.id, session_id=isession.id, state=rpm_state
     )
 
     session.add(state_1)
@@ -641,8 +684,6 @@ def test_rpm_state(session_and_user: dict[str, Session | list[User]]):
 
     asp_levels_2 = {"f_1": 0.6, "f_2": 0.4, "f_3": 0.5}
 
-    problem = Problem.from_problemdb(problem_db)
-
     scalarization_options = None
     solver = "pyomo_bonmin"
     solver_options = None
@@ -651,31 +692,36 @@ def test_rpm_state(session_and_user: dict[str, Session | list[User]]):
         problem,
         asp_levels_2,
         scalarization_options=scalarization_options,
-        solver=available_solvers[solver],
+        solver=available_solvers[solver]["constructor"],
+        solver_options=solver_options,
+    )
+
+    # create state
+
+    rpm_state = RPMState(
+        scalarization_options=scalarization_options,
+        solver=solver,
         solver_options=solver_options,
+        solver_results=results,
     )
 
     # create preferences
 
     rp_2 = ReferencePoint(aspiration_levels=asp_levels_2)
-    preferences = PreferenceDB(user_id=user.id, problem_id=problem_db.id, preference=rp_2)
-
-    session.add(preferences)
-    session.commit()
-    session.refresh(preferences)
 
     # create state
 
     rpm_state = RPMState(
+        preferences=rp_2,
         scalarization_options=scalarization_options,
         solver=solver,
         solver_options=solver_options,
         solver_results=results,
     )
 
-    state_2 = StateDB(
+    state_2 = StateDB.create(
+        database_session=session,
         problem_id=problem_db.id,
-        preference_id=preferences.id,
         session_id=isession.id,
         parent_id=state_1.id,
         state=rpm_state,
@@ -690,8 +736,8 @@ def test_rpm_state(session_and_user: dict[str, Session | list[User]]):
     assert len(state_1.children) == 1
     assert state_1.children[0] == state_2
 
-    assert state_1.preference.preference == rp_1
-    assert state_2.preference.preference == rp_2
+    assert state_1.state.preferences == rp_1
+    assert state_2.state.preferences == rp_2
 
     assert state_2.problem == problem_db
     assert state_2.session.user == user
@@ -699,3 +745,435 @@ def test_rpm_state(session_and_user: dict[str, Session | list[User]]):
     assert state_2.children == []
     assert state_2.parent.problem == problem_db
     assert state_2.parent.session.user == user
+
+
+def test_problem_metadata(session_and_user: dict[str, Session | list[User]]):
+    """Test that the problem metadata can be put into database and brought back."""
+    session = session_and_user["session"]
+    user = session_and_user["user"]
+
+    # Just some test problem to attach the metadata to
+    problem = ProblemDB.from_problem(dtlz2(5, 3), user=user)
+
+    session.add(problem)
+    session.commit()
+    session.refresh(problem)
+
+    representative_name = "Test solutions"
+    representative_description = "These solutions are used for testing"
+    representative_variables = {"x_1": [1.1, 2.2, 3.3], "x_2": [-1.1, -2.2, -3.3]}
+    representative_objectives = {"f_1": [0.1, 0.5, 0.9], "f_2": [-0.1, 0.2, 199.2], "f_1_min": [], "f_2_min": []}
+    solution_data = representative_variables | representative_objectives
+    representative_ideal = {"f_1": 0.1, "f_2": -0.1}
+    representative_nadir = {"f_1": 0.9, "f_2": 199.2}
+
+    metadata = ProblemMetaDataDB(
+        problem_id=problem.id,
+    )
+
+    session.add(metadata)
+    session.commit()
+    session.refresh(metadata)
+
+    forest_metadata = ForestProblemMetaData(
+        metadata_id=metadata.id,
+        map_json="type: string",
+        schedule_dict={"type": "dict"},
+        years=["type:", "list", "of", "strings"],
+        stand_id_field="type: string",
+    )
+
+    repr_metadata = RepresentativeNonDominatedSolutions(
+        metadata_id=metadata.id,
+        name=representative_name,
+        description=representative_description,
+        solution_data=solution_data,
+        ideal=representative_ideal,
+        nadir=representative_nadir,
+    )
+
+    session.add(forest_metadata)
+    session.add(repr_metadata)
+    session.commit()
+    session.refresh(forest_metadata)
+    session.refresh(repr_metadata)
+
+    statement = select(ProblemMetaDataDB).where(ProblemMetaDataDB.problem_id == problem.id)
+    from_db_metadata = session.exec(statement).first()
+
+    assert from_db_metadata.id is not None
+    assert from_db_metadata.problem_id == problem.id
+
+    metadata_forest = from_db_metadata.forest_metadata[0]
+
+    assert isinstance(metadata_forest, ForestProblemMetaData)
+    assert metadata_forest.map_json == "type: string"
+    assert metadata_forest.schedule_dict == {"type": "dict"}
+    assert metadata_forest.years == ["type:", "list", "of", "strings"]
+    assert metadata_forest.stand_id_field == "type: string"
+
+    metadata_representative = from_db_metadata.representative_nd_metadata[0]
+
+    assert isinstance(metadata_representative, RepresentativeNonDominatedSolutions)
+    assert metadata_representative.name == representative_name
+    assert metadata_representative.solution_data == solution_data
+    assert metadata_representative.ideal == representative_ideal
+    assert metadata_representative.nadir == representative_nadir
+
+    assert problem.problem_metadata == from_db_metadata
+
+
+def test_group(session_and_user: dict[str, Session | list[User]]):
+    """なに？！ちょっとまって。。。ドクメンタはどこですか？？？"""
+    session: Session = session_and_user["session"]
+    user: User = session_and_user["user"]
+
+    group = Group(
+        name="TestGroup",
+        owner_id=user.id,
+        user_ids=[user.id],
+        problem_id=1,
+    )
+
+    session.add(group)
+    session.commit()
+    session.refresh(group)
+
+    assert group.id == 1
+    assert group.user_ids[0] == user.id
+    assert group.name == "TestGroup"
+
+
+def test_gnimbus_datas(session_and_user: dict[str, Session | list[User]]):
+    session: Session = session_and_user["session"]
+    user: User = session_and_user["user"]
+
+    group = Group(
+        name="TestGroup",
+        owner_id=user.id,
+        user_ids=[user.id],
+        problem_id=1,
+    )
+
+    session.add(group)
+    session.commit()
+    session.refresh(group)
+
+    giter = GroupIteration(
+        problem_id=1,
+        group_id=group.id,
+        info_container=OptimizationPreference(
+            set_preferences={},
+        ),
+        notified={},
+        parent_id=None,
+        parent=None,
+        child=None,
+    )
+    session.add(giter)
+    session.commit()
+    session.refresh(giter)
+
+    assert type(giter.info_container) is OptimizationPreference
+    assert giter.problem_id == 1
+    assert giter.group_id == group.id
+
+
+def test_enautilus_state(session_and_user: dict[str, Session | list[User]]):
+    """Test the E-NAUTILUS state that it works correctly."""
+    session = session_and_user["session"]
+    user = session_and_user["user"]
+
+    # create interactive session
+    isession = InteractiveSessionDB(user_id=user.id)
+
+    session.add(isession)
+    session.commit()
+    session.refresh(isession)
+
+    # use dummy problem
+    dummy_problem = Problem(
+        name="Synthetic-4D",
+        description="Unit-test Problem for E-NAUTILUS",
+        variables=[Variable(name="x", symbol="x", variable_type=VariableTypeEnum.real)],
+        objectives=[
+            Objective(name="f1", symbol="f1", maximize=False),
+            Objective(name="f2", symbol="f2", maximize=True),
+            Objective(name="f3", symbol="f3", maximize=False),
+            Objective(name="f4", symbol="f4", maximize=True),
+        ],
+    )
+
+    x = np.array([1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0])
+    f1 = np.array([0.40, 0.60, 0.50, 0.70, 0.45, 0.55, 0.65, 0.48])
+    f2 = np.array([4.00, 3.80, 4.10, 3.70, 4.05, 3.90, 3.60, 4.20])
+    f3 = np.array([1.00, 1.30, 1.10, 1.40, 1.05, 1.20, 1.35, 1.15])
+    f4 = np.array([2.50, 2.30, 2.60, 2.20, 2.55, 2.40, 2.10, 2.65])
+
+    nadir = {"f1": np.max(f1), "f2": np.min(f2), "f3": np.max(f3), "f4": np.min(f4)}
+    ideal = {"f1": np.min(f1), "f2": np.max(f2), "f3": np.min(f3), "f4": np.max(f4)}
+
+    non_dom_data = {
+        "x": x.tolist(),
+        "f1": f1.tolist(),
+        "f1_min": f1.tolist(),
+        "f2": f2.tolist(),
+        "f2_min": (-f2).tolist(),
+        "f3": f3.tolist(),
+        "f3_min": f3.tolist(),
+        "f4": f4.tolist(),
+        "f4_min": (-f4).tolist(),
+    }
+
+    # add problem to DB and refresh it
+    problemdb = ProblemDB.from_problem(dummy_problem, user)
+
+    session.add(problemdb)
+    session.commit()
+    session.refresh(problemdb)
+
+    metadata = ProblemMetaDataDB(
+        problem_id=problemdb.id,
+    )
+
+    # add metadata to DB
+    session.add(metadata)
+    session.commit()
+    session.refresh(metadata)
+
+    reprdata = RepresentativeNonDominatedSolutions(
+        metadata_id=metadata.id,
+        name="Dummy data",
+        description="Dummy data for a problem",
+        solution_data=non_dom_data,
+        ideal=ideal,
+        nadir=nadir,
+    )
+
+    # add reprdata to DB
+    session.add(reprdata)
+    session.commit()
+    session.refresh(reprdata)
+
+    # test the nautilus step state
+    # first iteration
+    selected_point = nadir
+    reachable_indices = list(range(len(x)))  # entire front reachable
+
+    total_iters = 2  # DM first thinks 2 iterations are enough
+    n_points = 3  # DM wants to see 3 points at first
+    current = 0
+
+    # First iteration
+    res = enautilus_step(
+        problem=dummy_problem,
+        non_dominated_points=non_dom_data,
+        current_iteration=current,
+        iterations_left=total_iters - current,
+        selected_point=selected_point,
+        reachable_point_indices=reachable_indices,
+        number_of_intermediate_points=n_points,
+    )
+
+    enautilus_state = ENautilusState(
+        non_dominated_solutions_id=reprdata.id,
+        current_iteration=current,
+        iterations_left=total_iters - current,
+        selected_point=selected_point,
+        reachable_point_indices=reachable_indices,
+        number_of_intermediate_points=n_points,
+        enautilus_results=res,
+    )
+
+    state_1 = StateDB.create(
+        database_session=session,
+        problem_id=problemdb.id,
+        session_id=isession.id,
+        parent_id=None,
+        state=enautilus_state,
+    )
+
+    session.add(state_1)
+    session.commit()
+    session.refresh(state_1)
+
+    # Second iteration
+    res_2 = enautilus_step(
+        problem=dummy_problem,
+        non_dominated_points=non_dom_data,
+        current_iteration=res.current_iteration,
+        iterations_left=res.iterations_left,
+        selected_point=res.intermediate_points[0],
+        reachable_point_indices=res.reachable_point_indices[0],
+        number_of_intermediate_points=n_points,
+    )
+
+    enautilus_state_2 = ENautilusState(
+        non_dominated_solutions_id=reprdata.id,
+        current_iteration=res.current_iteration,
+        iterations_left=res.iterations_left,
+        selected_point=res.intermediate_points[0],
+        reachable_point_indices=res.reachable_point_indices[0],
+        number_of_intermediate_points=n_points,
+        enautilus_results=res_2,
+    )
+
+    state_2 = StateDB.create(
+        database_session=session,
+        problem_id=problemdb.id,
+        session_id=isession.id,
+        parent_id=enautilus_state.id,
+        state=enautilus_state_2,
+    )
+
+    session.add(state_2)
+    session.commit()
+    session.refresh(state_2)
+
+    assert state_1.problem_id == problemdb.id
+    assert state_2.problem_id == problemdb.id
+
+    assert state_1.session_id == isession.id
+    assert state_2.session_id == isession.id
+
+    assert state_1.parent is None
+    assert state_2.parent == state_1
+
+    assert len(state_1.children) == 1
+    assert state_1.children[0] == state_2
+    assert state_2.children == []
+
+
+def test_nimbus_models(session_and_user: dict[str, Session | list[User]]):
+    """Test that the NIMBUS models are in working order."""
+    session: Session = session_and_user["session"]
+    user: User = session_and_user["user"]
+    problem_db = user.problems[0]
+    problem = Problem.from_problemdb(problem_db)
+
+    isession = InteractiveSessionDB(user_id=user.id)
+
+    session.add(isession)
+    session.commit()
+    session.refresh(isession)
+
+    # 1. Initialize the NIMBUS problem (NIMBUSInitializationState)
+    results_1 = generate_starting_point(
+        problem=problem,
+    )
+    nimbus_init_state = NIMBUSInitializationState(solver_results=results_1)
+
+    state_1 = StateDB.create(
+        database_session=session, problem_id=problem_db.id, session_id=isession.id, state=nimbus_init_state
+    )
+
+    session.add(state_1)
+    session.commit()
+    session.refresh(state_1)
+
+    actual_state_1: NIMBUSInitializationState = state_1.state
+
+    assert type(actual_state_1) is NIMBUSInitializationState
+    assert np.allclose(
+        [x for _, x in actual_state_1.solver_results.optimal_objectives.items()],
+        [x for _, x in results_1.optimal_objectives.items()],
+        0.001,
+    )
+
+    # 2. Solve sub problems (NIMBUSClassificationState)
+    aspirations = {"f_1": 0.1, "f_2": 0.9, "f_3": 0.6}
+
+    results_2 = solve_sub_problems(
+        problem=problem, current_objectives=results_1.optimal_objectives, reference_point=aspirations, num_desired=4
+    )
+    nimbus_classification_state = NIMBUSClassificationState(
+        preferences=ReferencePoint(aspiration_levels=aspirations),
+        current_objectives=results_1.optimal_objectives,
+        previous_preferences=ReferencePoint(aspiration_levels=aspirations),
+        solver_results=results_2,
+    )
+
+    state_2 = StateDB.create(
+        database_session=session, problem_id=problem_db.id, session_id=isession.id, state=nimbus_classification_state
+    )
+
+    session.add(state_2)
+    session.commit()
+    session.refresh(state_2)
+
+    actual_state_2: NIMBUSClassificationState = state_2.state
+
+    assert type(actual_state_2) is NIMBUSClassificationState
+    assert np.allclose(
+        [x for _, x in actual_state_2.preferences.aspiration_levels.items()], [x for _, x in aspirations.items()], 0.001
+    )
+    assert np.allclose(
+        [x for _, x in actual_state_2.solver_results[0].optimal_objectives.items()],
+        [x for _, x in results_2[0].optimal_objectives.items()],
+        0.001,
+    )
+
+    # 3. (TODO) Save a found solution (NIMBUSSaveState)
+    # 4. Finalize the NIMBUS process (NIMBUSFinalState)
+    nimbus_final_state = NIMBUSFinalState(
+        solution_origin_state_id=state_2.state.id,
+        solution_result_index=0,
+        solver_results=results_2[0]
+    )
+
+    state_3 = StateDB.create(
+        database_session=session, problem_id=problem_db.id, session_id=isession.id, state=nimbus_final_state
+    )
+
+    session.add(state_3)
+    session.commit()
+    session.refresh(state_3)
+
+    actual_state_3: NIMBUSFinalState = state_3.state
+    assert type(actual_state_3) is NIMBUSFinalState
+    assert np.allclose(
+        [x for _, x in actual_state_3.solver_results.optimal_objectives.items()],
+        [x for _, x in results_2[0].optimal_objectives.items()],
+        0.01,
+    )
+    assert np.allclose(
+        [x for _, x in actual_state_3.solver_results.optimal_variables.items()],
+        [x for _, x in results_2[0].optimal_variables.items()],
+        0.01,
+    )
+
+
+def test_nimbus_initialize_w_multidimensional_constraints(session_and_user: dict[str, Session | list[User]]):
+    """Test that the NIMBUS initialization model works with multidimensional constraints."""
+    session: Session = session_and_user["session"]
+    user: User = session_and_user["user"]
+    problem = multi_valued_constraint_problem()
+
+    problem_db = ProblemDB.from_problem(problem, user)
+    isession = InteractiveSessionDB(user_id=user.id)
+
+    session.add(isession)
+    session.add(problem_db)
+    session.commit()
+    session.refresh(isession)
+    session.refresh(problem_db)
+
+    results_1 = generate_starting_point(problem=problem, solver=PyomoBonminSolver)
+    nimbus_init_state = NIMBUSInitializationState(solver_results=results_1)
+
+    state_1 = StateDB.create(
+        database_session=session, problem_id=problem_db.id, session_id=isession.id, state=nimbus_init_state
+    )
+
+    session.add(state_1)
+    session.commit()
+    session.refresh(state_1)
+
+    actual_state_1: NIMBUSInitializationState = state_1.state
+
+    assert type(actual_state_1) is NIMBUSInitializationState
+    assert np.allclose(
+        [x for _, x in actual_state_1.solver_results.optimal_objectives.items()],
+        [x for _, x in results_1.optimal_objectives.items()],
+        0.001,
+    )