PreconConfig.hpp

#pragma once
 
#include <string>
 
#include <Eigen/IterativeLinearSolvers>
 
#include <amdis/Initfile.hpp>
 
namespace AMDiS
{
  template <class Precon>
  struct PreconConfig
  {
    static void init(std::string const& prefix, Precon& solver) {}
  };
 
  template <class T, class I>
  struct PreconConfig<Eigen::IncompleteLUT<T,I>>
  {
    static void init(std::string const& prefix, Eigen::IncompleteLUT<T,I>& solver)
    {
      auto dropTol = Parameters::get<double>(prefix + "->drop tolerance");
      if (dropTol)
        solver.setDroptol(dropTol.value());
 
      auto fillFactor = Parameters::get<int>(prefix + "->fill factor");
      if (fillFactor)
        solver.setFillfactor(fillFactor.value());
    }
  };
 
} // end namespace AMDiS