matrix_based_solver.h

Go to the documentation of this file.

 
#ifndef HYBRIDADRSOLVER_MATRIX_BASED_SOLVER_H
#define HYBRIDADRSOLVER_MATRIX_BASED_SOLVER_H
 
#include "core/problem_definition.h"
#include "core/solver.h"
#include "core/types.h"
 
namespace HybridADRSolver {
using namespace dealii;
 
// Forward declarations
template <int dim> struct ScratchData;
struct CopyData;

template <int dim> class MatrixBasedSolver : public ParallelSolverBase<dim> {
public:
    MatrixBasedSolver(const ProblemInterface<dim>& problem, unsigned int degree,
                      MPI_Comm comm,
                      const SolverParameters& params = SolverParameters());
 
    void run(unsigned int n_refinements) override;
 
    SolverType get_solver_type() const override {
        return SolverType::MatrixBased;
    }
 
    std::string get_name() const override { return "Matrix-Based Solver"; }

    const LADistributed::MPI::SparseMatrix& get_system_matrix() const {
        return system_matrix;
    }

    double compute_memory_usage() const;

    double compute_l2_error();
 
protected:
    void setup_dofs() override;
    void assemble_system() override;
    void solve() override;
    void output_results(unsigned int cycle) const override;
 
private:
    const ProblemInterface<dim>& problem;
    unsigned int fe_degree;
 
    LADistributed::MPI::SparseMatrix system_matrix;
    LADistributed::MPI::Vector system_rhs;
    LADistributed::MPI::Vector solution;
};
 
} // namespace HybridADRSolver
 
#endif // HYBRIDADRSOLVER_MATRIX_BASED_SOLVER_H