9 #ifndef __IPMA27TSOLVERINTERFACE_HPP__
10 #define __IPMA27TSOLVERINTERFACE_HPP__
33 const std::string& prefix);
63 Index numberOfNegEVals);
224 Index numberOfNegEVals);
FORTRAN_INTEGER_TYPE ipfint
Interface to the symmetric linear solver MA27, derived from SparseSymLinearSolverInterface.
bool InitializeImpl(const OptionsList &options, const std::string &prefix)
overloaded from AlgorithmStrategyObject
Number pivtol_
Pivol tolerance.
virtual Index NumberOfNegEVals() const
Number of negative eigenvalues detected during last factorization.
ipfint icntl_[30]
integer control values
bool initialized_
Flag indicating if internal data is initialized.
virtual ~Ma27TSolverInterface()
Destructor.
ipfint nsteps_
MA27's NSTEPS.
ipfint * ikeep_
MA27's IKEEP.
Ma27TSolverInterface(const Ma27TSolverInterface &)
Copy Constructor.
bool refactorize_
Flag that is true if we just requested the values of the matrix again (SYMSOLVER_CALL_AGAIN) and have...
virtual bool IncreaseQuality()
Request to increase quality of solution for next solve.
double * a_
factor A of matrix
Number la_init_factor_
Factor for estimating initial value of la.
static void RegisterOptions(SmartPtr< RegisteredOptions > roptions)
Methods for IpoptType.
bool ignore_singularity_
Flag indicating if MA27 should continue if a singular matrix is detected, but right hands sides are s...
bool warm_start_same_structure_
Flag indicating whether the TNLP with identical structure has already been solved before.
virtual ESymSolverStatus InitializeStructure(Index dim, Index nonzeros, const Index *airn, const Index *ajcn)
Method for initializing internal stuctures.
Number liw_init_factor_
Factor for estimating initial value of liw.
ESymSolverStatus SymbolicFactorization(const Index *airn, const Index *ajcn)
Call MA27AD and reserve memory for MA27 data.
double cntl_[5]
real control values
bool skip_inertia_check_
Flag indicating if the interia is always assumed to be correct.
void operator=(const Ma27TSolverInterface &)
Overloaded Equals Operator.
bool liw_increase_
flag indicating that liw should be increased before next factorization
EMatrixFormat MatrixFormat() const
Query of requested matrix type that the linear solver understands.
ipfint * iw_
integer work space
bool pivtol_changed_
Flag indicating if the matrix has to be refactorized because the pivot tolerance has been changed.
virtual bool ProvidesInertia() const
Query whether inertia is computed by linear solver.
Index dim_
Number of rows and columns of the matrix.
ESymSolverStatus Factorization(const Index *airn, const Index *ajcn, bool check_NegEVals, Index numberOfNegEVals)
Call MA27BD to factorize the Matrix.
Number pivtolmax_
Maximal pivot tolerance.
bool la_increase_
flag indicating that la should be increased before next factorization
Index negevals_
Number of negative eigenvalues.
Ma27TSolverInterface()
Constructor.
ipfint maxfrt_
MA27's MAXFRT.
ipfint liw_
length of integer work space
Index nonzeros_
Number of nonzeros of the matrix.
virtual double * GetValuesArrayPtr()
Method returing an internal array into which the nonzero elements (in the same order as airn and ajcn...
virtual ESymSolverStatus MultiSolve(bool new_matrix, const Index *airn, const Index *ajcn, Index nrhs, double *rhs_vals, bool check_NegEVals, Index numberOfNegEVals)
Solve operation for multiple right hand sides.
ESymSolverStatus Backsolve(Index nrhs, double *rhs_vals)
Call MA27CD to do the backsolve.
Number meminc_factor_
Factor for increaseing memory.
This class stores a list of user set options.
Template class for Smart Pointers.
Base class for interfaces to symmetric indefinite linear solvers for sparse matrices.
EMatrixFormat
Enum to specify sparse matrix format.
@ Triplet_Format
Triplet (MA27) format.
ESymSolverStatus
Enum to report outcome of a linear solve.
int Index
Type of all indices of vectors, matrices etc.
double Number
Type of all numbers.