9 #ifndef __IP_AUGRESTOSYSTEMSOLVER_HPP__
10 #define __IP_AUGRESTOSYSTEMSOLVER_HPP__
35 bool skip_orig_aug_solver_init=
true);
43 const std::string& prefix);
71 Index numberOfNegEVals);
138 const Vector& any_vec_in_c);
146 const Vector& any_vec_in_d);
152 const Vector& any_vec_in_n_c);
157 const Vector& any_vec_in_p_c);
162 const Vector& any_vec_in_n_d);
167 const Vector& any_vec_in_p_d);
ESymSolverStatus
Enum to report outcome of a linear solve.
CachedResults< SmartPtr< Vector > > rhs_dR_cache_
bool skip_orig_aug_solver_init_
CachedResults< SmartPtr< Vector > > sigma_tilde_n_c_inv_cache_
Class that converts the an augmented system with compound restoration pieces into a smaller "pivoted"...
SmartPtr< const Vector > Sigma_tilde_n_c_inv(const SmartPtr< const Vector > &sigma_tilde_n_c, Number delta_x, const Vector &any_vec_in_n_c)
Sigma tilde is the sum of Sigma and delta_x times the identity.
virtual Index NumberOfNegEVals() const
Returns the number of negative eigenvalues from the original augmented system call.
SmartPtr< const Vector > Sigma_tilde_p_d_inv(const SmartPtr< const Vector > &sigma_tilde_p_d, Number delta_x, const Vector &any_vec_in_p_d)
AugRestoSystemSolver()
Default Constructor.
double Number
Type of all numbers.
CachedResults< SmartPtr< Vector > > d_x_plus_wr_d_cache_
SmartPtr< const Vector > Neg_Omega_d_plus_D_d(const Matrix &Pd_L, const SmartPtr< const Vector > &sigma_tilde_n_d_inv, const Matrix &neg_Pd_U, const SmartPtr< const Vector > &sigma_tilde_p_d_inv, const Vector *D_d, const Vector &any_vec_in_d)
CachedResults< SmartPtr< Vector > > neg_omega_c_plus_D_c_cache_
SmartPtr< const Vector > Sigma_tilde_p_c_inv(const SmartPtr< const Vector > &sigma_tilde_p_c, Number delta_x, const Vector &any_vec_in_p_c)
SmartPtr< const Vector > Neg_Omega_c_plus_D_c(const SmartPtr< const Vector > &sigma_tilde_n_c_inv, const SmartPtr< const Vector > &sigma_tilde_p_c_inv, const Vector *D_c, const Vector &any_vec_in_c)
virtual bool ProvidesInertia() const
Query whether inertia is computed by linear solver.
CachedResults< SmartPtr< Vector > > sigma_tilde_p_d_inv_cache_
This is the base class for all derived symmetric matrix types.
SmartPtr< AugSystemSolver > orig_aug_solver_
Template class for Smart Pointers.
This class stores a list of user set options.
Base class for Solver for the augmented system.
SmartPtr< const Vector > D_x_plus_wr_d(const SmartPtr< const Vector > &CD_x0, Number factor, const Vector &wr_d)
void operator=(const AugRestoSystemSolver &)
Overloaded Equals Operator.
SmartPtr< const Vector > Sigma_tilde_n_d_inv(const SmartPtr< const Vector > &sigma_tilde_n_d, Number delta_x, const Vector &any_vec_in_n_d)
SmartPtr< const Vector > Rhs_dR(const Vector &rhs_d, const SmartPtr< const Vector > &sigma_tilde_n_d_inv, const Vector &rhs_n_d, const Matrix &pd_L, const SmartPtr< const Vector > &sigma_tilde_p_d_inv, const Vector &rhs_p_d, const Matrix &pd_U)
int Index
Type of all indices of vectors, matrices etc.
bool InitializeImpl(const OptionsList &options, const std::string &prefix)
overloaded from AlgorithmStrategyObject
virtual bool IncreaseQuality()
Request to increase quality of solution for next solve.
CachedResults< SmartPtr< Vector > > sigma_tilde_n_d_inv_cache_
virtual ~AugRestoSystemSolver()
Default destructor.
virtual ESymSolverStatus Solve(const SymMatrix *W, double W_factor, const Vector *D_x, double delta_x, const Vector *D_s, double delta_s, const Matrix *J_c, const Vector *D_c, double delta_c, const Matrix *J_d, const Vector *D_d, double delta_d, const Vector &rhs_x, const Vector &rhs_s, const Vector &rhs_c, const Vector &rhs_d, Vector &sol_x, Vector &sol_s, Vector &sol_c, Vector &sol_d, bool check_NegEVals, Index numberOfNegEVals)
Translate the augmented system (in the full space of the restoration variables) into the smaller spac...
SmartPtr< const Vector > Rhs_cR(const Vector &rhs_c, const SmartPtr< const Vector > &sigma_tilde_n_c_inv, const Vector &rhs_n_c, const SmartPtr< const Vector > &sigma_tilde_p_c_inv, const Vector &rhs_p_c)
CachedResults< SmartPtr< Vector > > neg_omega_d_plus_D_d_cache_
CachedResults< SmartPtr< Vector > > rhs_cR_cache_
CachedResults< SmartPtr< Vector > > sigma_tilde_p_c_inv_cache_