MLMGOptions.H

Go to the documentation of this file.

#ifndef MLMGOPTIONS_H
#define MLMGOPTIONS_H
 
#include <string>
 
#include "AMReX_REAL.H"
#include "AMReX_ParmParse.H"
#include "AMReX_MLLinOp.H"
#include "AMReX_MLMG.H"
 
using namespace amrex::literals;
 
namespace Hydro {
class NodalProjector;
class MacProjector;
} // namespace Hydro
 
namespace amr_wind {

struct MLMGOptions
{
    explicit MLMGOptions(const std::string& prefix);

    MLMGOptions(
        const std::string& default_prefix, const std::string& custom_prefix);
 
    void operator()(amrex::MLMG& /*mlmg*/);
    void operator()(Hydro::NodalProjector& /*nodal_proj*/);
    void operator()(Hydro::MacProjector& /*mac_proj*/);

    amrex::LPInfo& lpinfo() { return m_lpinfo; }
 
    // Linear operator options
    int max_order{2};
 
#ifdef BL_USE_FLOAT
    amrex::Real rel_tol{std::numeric_limits<amrex::Real>::epsilon() * 1.0e2_rt};

    amrex::Real abs_tol{std::numeric_limits<amrex::Real>::epsilon() * 1.0e1_rt};
#else
    amrex::Real rel_tol{std::numeric_limits<amrex::Real>::epsilon() * 1.0e5_rt};

    amrex::Real abs_tol{std::numeric_limits<amrex::Real>::epsilon() * 1.0e2_rt};
#endif
 
private:
    void parse_options(const std::string& /*prefix*/);

    amrex::LPInfo m_lpinfo;
 
    // MLMG options

    std::string m_bottom_solver_type{"bicgstab"};

    std::string m_hypre_namespace{"hypre"};

    std::string m_hypre_interface{"ij"};
 
#ifdef BL_USE_FLOAT
    amrex::Real m_bottom_rel_tol{
        std::numeric_limits<amrex::Real>::epsilon() * 1.0e3_rt};

    amrex::Real m_bottom_abs_tol{
        std::numeric_limits<amrex::Real>::epsilon() * 1.0e2_rt};
#else
    amrex::Real m_bottom_rel_tol{
        std::numeric_limits<amrex::Real>::epsilon() * 1.0e12_rt};

    amrex::Real m_bottom_abs_tol{
        std::numeric_limits<amrex::Real>::epsilon() * 1.0e4_rt};
#endif

    int m_verbose{0};
    int m_max_iter{200};
    int m_max_fmg_iters{0};
    int m_num_pre_smooth{2};
    int m_num_post_smooth{2};
    int m_num_final_smooth{8};
    int m_num_bottom_smooth{0};

    int m_bottom_verbose{0};

    int m_bottom_max_iter{200};

    bool m_do_fixed_iters{true};

    bool m_do_nsolve{false};

#ifdef AMREX_USE_GPU
    int m_nsolve_grid_size{64};
#else
    int m_nsolve_grid_size{16};
#endif
};
 
} // namespace amr_wind
 
#endif /* MLMGOPTIONS_H */