00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #ifndef SUNDANCE_INTERNALSOLVEMANAGER_HPP
00032 #define SUNDANCE_INTERNALSOLVEMANAGER_HPP
00033
00034 #include "SundanceDefs.hpp"
00035 #include "SundanceMesh.hpp"
00036 #include "SundanceExpr.hpp"
00037 #include "SundanceBlock.hpp"
00038 #include "SundanceDiscreteFunction.hpp"
00039 #include "PlayaLinearOperatorDecl.hpp"
00040 #include "PlayaVectorDecl.hpp"
00041 #include "PlayaVectorType.hpp"
00042 #include "PlayaSolverState.hpp"
00043
00044 namespace Sundance
00045 {
00046
00047 class LinearSolveDriver
00048 {
00049 public:
00050
00051 LinearSolveDriver() {}
00052
00053
00054
00055 Expr formSolutionExpr(const Array<Vector<double> >& solnVector,
00056 const Array<RCP<DiscreteSpace> >& solutionSpace,
00057 const Array<Array<string> >& names,
00058 int verb) const ;
00059
00060
00061 SolverState<double> solve(const LinearSolver<double>& solver,
00062 const LinearOperator<double>& A,
00063 const Array<Vector<double> >& rhs,
00064 const Array<RCP<DiscreteSpace> >& solutionSpace,
00065 const Array<Array<string> >& names,
00066 int verb,
00067 Expr& soln) const ;
00068
00069
00070
00071 void writeIntoSolutionExpr(
00072 const Array<Vector<double> >& solnVector,
00073 Expr soln, int verb) const ;
00074
00075
00076 static std::string& badMatrixFilename()
00077 {static std::string rtn = "badMatrix.dat"; return rtn;}
00078
00079
00080 static std::string& badVectorFilename()
00081 {static std::string rtn = "badVector.dat"; return rtn;}
00082
00083
00084 static bool& solveFailureIsFatal()
00085 {static bool rtn=true; return rtn;}
00086
00087
00088 static bool& dumpBadMatrix()
00089 {static bool rtn=false; return rtn;}
00090
00091
00092
00093
00094 private:
00095 };
00096
00097 }
00098
00099 #endif