solver_info.hh

#ifndef SOLVER_INFO_HH
#define SOLVER_INFO_HH
 
namespace Dune{
    namespace Composites{
 
        class SolverInfo{
 
            public:
                int verb;
                bool UG;
                int MaxIt;
                double KrylovTol;
                double eigenvalue_threshold;
                bool widlund_part_unity;
                int nev;
                int nev_arpack;
                double eigen_shift;
                std::vector<double> times;
                Dune::InverseOperatorResult results;
                std::string solver, preconditioner, subdomainSolver, coarseSolver;
                int overlap;
                bool coarseSpaceActive;
 
                SolverInfo(int verb_ = 0): verb(verb_){
 
                    solver = "CG";
                    preconditioner = "GenEO";
                    subdomainSolver = "UMFPack";
                    coarseSolver = "UMFPack";
                    UG = false;
 
                    // Default Solver Options for Domain Decomposition Solver
                    overlap = 2;
 
                    // Default Krylov Solver Options
                    MaxIt = 1000;
                    KrylovTol = 1e-5;
 
                    // Default GenEO Preconditioner Options
                    verb = 10;
                    eigenvalue_threshold = 0.2;
                    widlund_part_unity = true;
                    nev = 10; // Number of eigenvalues used 
                    nev_arpack = 10; // Use ArPack's inbuilt error estimator
                    eigen_shift = 0.001;
 
                }
 
                void inline setTimes(std::vector<double>& t){
                    times.resize(t.size());
                    t = times;
                }
 
                void inline recordResult(Dune::InverseOperatorResult& res){
                    results.clear();
                    results.iterations = res.iterations;
                    results.reduction = res.reduction;
                    results.converged = res.converged;
                    results.conv_rate = res.conv_rate;
                    results.elapsed = res.elapsed;
                }
 
                bool inline converged(){return results.converged;}
        };
    }
}
 
#endif