voidlinearsolver.h

Go to the documentation of this file.

#ifndef VOID_LINEAR_SOLVER_H_
#define VOID_LINEAR_SOLVER_H_

#include <deal.II/lac/vector.h>
#include <deal.II/lac/block_sparsity_pattern.h>
#include <deal.II/lac/block_sparse_matrix.h>
#include <deal.II/lac/compressed_simple_sparsity_pattern.h>
#include <deal.II/lac/solver_cg.h>
#include <deal.II/lac/precondition.h>
#include <deal.II/lac/full_matrix.h>
#include <deal.II/lac/sparse_direct.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/numerics/vector_tools.h>

#include <vector>

namespace DOpE
{
  template <typename VECTOR>
    class VoidLinearSolver
  {
  public:
    VoidLinearSolver(ParameterReader &param_reader);
    ~VoidLinearSolver();
   
    static void declare_params(ParameterReader &param_reader);
 
    template<typename PROBLEM>
      void ReInit(PROBLEM& pde);

    template<typename PROBLEM, typename INTEGRATOR>
      void Solve(PROBLEM& pde, INTEGRATOR& integr, VECTOR &rhs, VECTOR &solution, bool force_matrix_build=false);
     
  protected:

  private:
  
  };

/*********************************Implementation************************************************/
 
  template <typename VECTOR>
    void VoidLinearSolver<VECTOR>::declare_params(ParameterReader &/*param_reader*/) 
  {
  }

  /******************************************************/

  template <typename VECTOR>
    VoidLinearSolver<VECTOR>::VoidLinearSolver(ParameterReader &/*param_reader*/) 
  {
  }

/******************************************************/

template <typename VECTOR>
 VoidLinearSolver<VECTOR>::~VoidLinearSolver()
{
}

/******************************************************/

template <typename VECTOR>
  template<typename PROBLEM>
  void  VoidLinearSolver<VECTOR>::ReInit(PROBLEM& /*pde*/)
{
 
}

/******************************************************/

template <typename VECTOR>
  template<typename PROBLEM, typename INTEGRATOR>
  void VoidLinearSolver<VECTOR>::Solve(PROBLEM& /*pde*/, INTEGRATOR& /*integr*/, VECTOR &rhs, VECTOR &solution, bool /*force_matrix_build*/)
{
  solution = rhs;
}


}
#endif