/Users/jack/Code/basso_dev/inc/SuperLUSolverCSR.h

Go to the documentation of this file.

#ifndef _SUPERLU_SOLVER_CSR_BASSO_H_
#define _SUPERLU_SOLVER_CSR_BASSO_H_

#include "basso.h"
#include "Solver.h"
#include "gmm_superlu_interface.h"

namespace basso {


        class SuperLUSolverCSR : public Solver<DynamicCSRMatrix>
        {
        public:
                SuperLUSolverCSR( DynamicCSRMatrix &K, nArray &f ) 
                        : Solver<DynamicCSRMatrix>(K,f) { ordering=1; }

                virtual int Solve( nArray &d );

                void SetOrdering( int ord ) { if ( ord>=0 && ord<3 ) ordering=ord; }
                
                Numeric ConditionEstimate() const 
                        { 
                                if ( Solved() ) 
                                        return condest; 
                                WarningMessage("SuperLUSolverCSR::ConditionEstimate","sytem not yet solved");
                                return 0;
                        }

        protected:
                int ordering;
                Numeric condest;

        };


        int SuperLUSolverCSR::Solve( nArray &d )
        {
                CSRMatrix Kcsr(NumEquations(),NumEquations());
                clean(*Kptr,1e-12);
                copy(*Kptr,Kcsr);
                SuperLU_solve( Kcsr, d, *fptr, condest, ordering );
                solved=true;
                return 0;
        }

} // end namespace

#endif

---