mkl_dss.hpp

// FEAT3: Finite Element Analysis Toolbox, Version 3
// Copyright (C) 2010 by Stefan Turek & the FEAT group
// FEAT3 is released under the GNU General Public License version 3,
// see the file 'copyright.txt' in the top level directory for details.
 
#pragma once
 
// includes, FEAT
#include <kernel/base_header.hpp>
#include <kernel/solver/base.hpp>
#include <kernel/lafem/saddle_point_matrix.hpp>
#include <kernel/lafem/sparse_matrix_bcsr.hpp>
#include <kernel/lafem/sparse_matrix_csr.hpp>
#include <kernel/lafem/mean_filter.hpp>
#include <kernel/lafem/tuple_vector.hpp>
#include <kernel/lafem/dense_vector.hpp>
#include <kernel/lafem/dense_vector_blocked.hpp>
 
namespace FEAT
{
  namespace Solver
  {
#if defined(FEAT_HAVE_MKL) || defined(DOXYGEN)
    typedef double mkldss_native_data_type;
#ifdef MKL_ILP64
    typedef std::uint64_t mkldss_native_index_type;
#else
    typedef std::uint32_t mkldss_native_index_type;
#endif
 
    class MKLDSS :
      public SolverBase<LAFEM::DenseVector<double, Index>>
    {
    public:
      typedef LAFEM::SparseMatrixCSR<double, Index> MatrixType;
      typedef LAFEM::DenseVector<double, Index> VectorType;
 
      typedef SolverBase<VectorType> BaseClass;
 
    private:
      const MatrixType& _system_matrix;
      void* _mkl_dss_handle;
 
    public:
      explicit MKLDSS(const MatrixType& system_matrix);
 
      virtual ~MKLDSS();
 
      virtual String name() const override
      {
        return "MKLDSS";
      }
 
      virtual void init_symbolic() override;
      virtual void done_symbolic() override;
      virtual void init_numeric() override;
      virtual void done_numeric() override;
 
      virtual Status apply(VectorType& vec_sol, const VectorType& vec_rhs) override;
    }; // class MKLDSS
 
    inline std::shared_ptr<MKLDSS> new_mkl_dss(const LAFEM::SparseMatrixCSR<double, Index>& matrix)
    {
      return std::make_shared<MKLDSS>(matrix);
    }
#endif // defined(FEAT_HAVE_MKL) || defined(DOXYGEN)
  } // namespace Solver
} // namespace FEAT