HHO_DiffAdvecReac.hpp

Go to the documentation of this file.

// Implementation of the HHO scheme in 2D for the diffusion equation, with K piecewise constant, \beta \in H^1(\Omega)^d and \mu piecewise continuous
//
//   { -div(K \grad(u) + \beta u) + \mu u = f,       inside Omega
//   {                   K \grad(u) . nTF = g,       on GammaN
//   {                                  u = g,       on GammaD
//
//  Pure Neumann not yet working for nonzero reaction
 
#include "HHO-general/HHO2D.hpp"
#include "boost/program_options.hpp"
#include <unsupported/Eigen/SparseExtra>
 
const std::string mesh_dir = "../../typ2_meshes/";
 
bool program_options(
    const int,          
    const char *[],     
    std::string &,      
    std::string &,      
    std::vector<int> &, 
    size_t &,           
    size_t &,           
    std::string &,      
    bool &,             
    bool &              
);
 
// Parameters that get passed to program options
std::string mesh_name, bc_id, plot_file;
std::vector<int> id_tcase;
size_t L, K;
bool use_threads, export_matrix;
 
// @}