diff --git a/uno/ingredients/inequality_handling_methods/interior_point_methods/PrimalDualInteriorPointMethod.cpp b/uno/ingredients/inequality_handling_methods/interior_point_methods/PrimalDualInteriorPointMethod.cpp index 4e62dc21..eb0dc264 100644 --- a/uno/ingredients/inequality_handling_methods/interior_point_methods/PrimalDualInteriorPointMethod.cpp +++ b/uno/ingredients/inequality_handling_methods/interior_point_methods/PrimalDualInteriorPointMethod.cpp @@ -212,7 +212,8 @@ namespace uno { // check the inertia [[maybe_unused]] auto [number_pos_eigenvalues, number_neg_eigenvalues, number_zero_eigenvalues] = this->linear_solver->get_inertia(); - assert(number_pos_eigenvalues == number_variables && number_neg_eigenvalues == number_constraints && number_zero_eigenvalues == 0); + assert(number_pos_eigenvalues == problem.number_variables && number_neg_eigenvalues == problem.number_constraints && + number_zero_eigenvalues == 0); // rhs this->assemble_augmented_rhs(current_multipliers, problem.number_variables, problem.number_constraints); diff --git a/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.cpp b/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.cpp index 970243d6..aa7c3f78 100644 --- a/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.cpp +++ b/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.cpp @@ -45,7 +45,7 @@ namespace uno { // preallocate a bunch of stuff BQPDSolver::BQPDSolver(size_t number_variables, size_t number_constraints, size_t number_objective_gradient_nonzeros, size_t number_jacobian_nonzeros, size_t number_hessian_nonzeros, BQPDProblemType problem_type, const Options& options): - QPSolver(), number_hessian_nonzeros(number_hessian_nonzeros), + QPSolver(), lower_bounds(number_variables + number_constraints), upper_bounds(number_variables + number_constraints), constraints(number_constraints), diff --git a/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.hpp b/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.hpp index 85cd6a1a..74e87d6f 100644 --- a/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.hpp +++ b/uno/ingredients/subproblem_solvers/BQPD/BQPDSolver.hpp @@ -61,8 +61,6 @@ namespace uno { [[nodiscard]] double hessian_quadratic_product(const Vector& primal_direction) const override; private: - const size_t number_hessian_nonzeros; - std::vector lower_bounds{}, upper_bounds{}; // lower and upper bounds of variables and constraints std::vector constraints; SparseVector linear_objective;