44 #ifndef EPETRA_SERIALDENSESOLVER_H 45 #define EPETRA_SERIALDENSESOLVER_H 169 void SolveWithTranspose(
bool Flag) {Transpose_ = Flag;
if (Flag) TRANS_ =
'T';
else TRANS_ =
'N';
return;};
178 void EstimateSolutionErrors(
bool Flag) ;
188 virtual int Factor(
void);
194 virtual int Solve(
void);
200 virtual int Invert(
void);
206 virtual int ComputeEquilibrateScaling(
void);
212 virtual int EquilibrateMatrix(
void);
218 int EquilibrateRHS(
void);
225 virtual int ApplyRefinement(
void);
231 int UnequilibrateLHS(
void);
240 virtual int ReciprocalConditionEstimate(
double & Value);
293 int M()
const {
return(M_);};
296 int N()
const {
return(N_);};
299 double *
A()
const {
return(A_);};
302 int LDA()
const {
return(LDA_);};
305 double *
B()
const {
return(B_);};
308 int LDB()
const {
return(LDB_);};
311 int NRHS()
const {
return(NRHS_);};
314 double *
X()
const {
return(X_);};
317 int LDX()
const {
return(LDX_);};
320 double *
AF()
const {
return(AF_);};
323 int LDAF()
const {
return(LDAF_);};
326 int *
IPIV()
const {
return(IPIV_);};
329 double ANORM()
const {
return(ANORM_);};
332 double RCOND()
const {
return(RCOND_);};
337 double ROWCND()
const {
return(ROWCND_);};
342 double COLCND()
const {
return(COLCND_);};
345 double AMAX()
const {
return(AMAX_);};
348 double *
FERR()
const {
return(FERR_);};
351 double *
BERR()
const {
return(BERR_);};
354 double *
R()
const {
return(R_);};
357 double *
C()
const {
return(C_);};
362 virtual void Print(std::ostream& os)
const;
367 void AllocateWORK() {
if (WORK_==0) {LWORK_ = 4*N_; WORK_ =
new double[LWORK_];}
return;};
double * FERR() const
Returns a pointer to the forward error estimates computed by LAPACK.
bool Transpose()
Returns true if transpose of this matrix has and will be used.
int LDAF() const
Returns the leading dimension of the factored matrix.
bool B_Equilibrated()
Returns true if RHS is equilibrated (RHS available via B() and LDB()).
Epetra_SerialDenseMatrix * FactoredMatrix() const
Returns pointer to factored matrix (assuming factorization has been performed).
bool SolutionRefined()
Returns true if the current set of vectors has been refined.
double ROWCND() const
Ratio of smallest to largest row scale factors for the this matrix (returns -1 if not yet computed)...
int M() const
Returns row dimension of system.
void SolveWithTranspose(bool Flag)
If Flag is true, causes all subsequent function calls to work with the transpose of this matrix...
bool SolutionErrorsEstimated()
Returns true if forward and backward error estimated have been computed (available via FERR() and BER...
double RCOND() const
Returns the reciprocal of the condition number of the this matrix (returns -1 if not yet computed)...
double * R() const
Returns a pointer to the row scaling vector used for equilibration.
Epetra_CompObject & operator=(const Epetra_CompObject &src)
Epetra_SerialDenseMatrix: A class for constructing and using real double precision general dense matr...
Epetra_SerialDenseMatrix * RHS_
double * C() const
Returns a pointer to the column scale vector used for equilibration.
double ANORM() const
Returns the 1-Norm of the this matrix (returns -1 if not yet computed).
double * A() const
Returns pointer to the this matrix.
virtual void Print(std::ostream &os) const
Print object to an output stream Print method.
Epetra_BLAS: The Epetra BLAS Wrapper Class.
double * X() const
Returns pointer to current solution.
double COLCND() const
Ratio of smallest to largest column scale factors for the this matrix (returns -1 if not yet computed...
Epetra_SerialDenseMatrix * LHS() const
Returns pointer to current LHS.
Epetra_Object: The base Epetra class.
virtual bool ShouldEquilibrate()
Returns true if the LAPACK general rules for equilibration suggest you should equilibrate the system...
int LDB() const
Returns the leading dimension of the RHS.
Epetra_SerialDenseMatrix * Matrix_
Epetra_CompObject: Functionality and data that is common to all computational classes.
Epetra_SerialDenseMatrix * Matrix() const
Returns pointer to current matrix.
double * AF() const
Returns pointer to the factored matrix (may be the same as A() if factorization done in place)...
bool Inverted()
Returns true if matrix inverse has been computed (inverse available via AF() and LDAF()).
bool SolutionErrorsEstimated_
bool ReciprocalConditionEstimated()
Returns true if the condition number of the this matrix has been computed (value available via Recipr...
Epetra_LAPACK: The Epetra LAPACK Wrapper Class.
int LDA() const
Returns the leading dimension of the this matrix.
void SolveToRefinedSolution(bool Flag)
Causes all solves to compute solution to best ability using iterative refinement. ...
double * B() const
Returns pointer to current RHS.
bool ReciprocalConditionEstimated_
bool A_Equilibrated()
Returns true if factor is equilibrated (factor available via AF() and LDAF()).
Epetra_SerialDenseMatrix * Factor_
int NRHS() const
Returns the number of current right hand sides and solution vectors.
void FactorWithEquilibration(bool Flag)
Causes equilibration to be called just before the matrix factorization as part of the call to Factor...
bool Factored()
Returns true if matrix is factored (factor available via AF() and LDAF()).
bool Solved()
Returns true if the current set of vectors has been solved.
int Solve(int, TYPE *, TYPE *, TYPE *)
int LDX() const
Returns the leading dimension of the solution.
Epetra_SerialDenseMatrix * LHS_
double AMAX() const
Returns the absolute value of the largest entry of the this matrix (returns -1 if not yet computed)...
double * BERR() const
Returns a pointer to the backward error estimates computed by LAPACK.
Epetra_SerialDenseMatrix * RHS() const
Returns pointer to current RHS.
int * IPIV() const
Returns pointer to pivot vector (if factorization has been computed), zero otherwise.
Epetra_SerialDenseSolver: A class for solving dense linear problems.
bool EstimateSolutionErrors_
int N() const
Returns column dimension of system.