domainDecomp.hh

Go to the documentation of this file.

 
#ifndef solverDomainDecomp_hh
#define solverDomainDecomp_hh
 
#include "basics/typedefs.hh"
#include "operator/compositions.hh"
#include "operator/solverFabric.hh"
#include "operator/submatrix.hh"
#include "space/domainDecomp.hh"
#include "toolbox/sequence.hh"
#include "toolbox/set.hh"
 
namespace concepts {
 
 // ************************************************************** DDSolver **
 
  template<class F, class G = typename Realtype<F>::type>
  class DDSolver : public Operator<F> {
  public:
    DDSolver(DDSpace<G>& spc, Sequence<SparseMatrix<F>*>& A, bool explSchur,
             SolverFabric<F>& innerSolver,
             SolverFabric<F>& schurComplSolver,
             SolverFabric<F>* schurComplPrecond = 0);
    virtual ~DDSolver();
    virtual void operator()(const Function<F>& fncY, Function<F>& fncX);
    void operator()(const Vector<F>& fncY, Vector<F>& fncX);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    DDSpace<G>& spc_;
    Set<IndexRange> indicesB_;
    uint n_;
    Set<uint> domains_;
    DummySpace<G>* spc_B_;
    Sequence<const SparseMatrix<F>*> A_;
    Sequence<SubMatrixN<SparseMatrix<F> >*> A_II_, A_IB_, A_BI_, A_BB_;
    Sequence<Operator<F>*> A_II_inv_;
    Sequence<Operator<F>*> S_;
    Operator<F>* S_0_inv_;
    Operator<F>* S_inv_;
  };
 
 
} // namespace concepts
 
#endif // solverDomainDecomp_hh