hypre-compute-errors.hpp

Go to the documentation of this file.

#ifndef HYPRE_COMPUTE_ERRORS_HPP
#define HYPRE_COMPUTE_ERRORS_HPP
 
#include "hypre.hpp"
 
namespace HArDCore2D {  
  
  namespace DSL {
 
    struct ComputeErrorsParameters {
      bool use_threads = true;
      std::ostream & output = std::cout;
    };
    
    std::map<std::string, double>
    compute_discrete_errors(
                            const HYPRE * scheme,
                            const Eigen::VectorXd & uph,
                            const Eigen::VectorXd & upI,
                            const Eigen::VectorXd & brh,
                            const Eigen::VectorXd & brI,
                            const double & viscosity,
                            const double & magnetic_diffusivity,
                            const ComputeErrorsParameters & parameters = {}
                            );
 
    std::map<std::string, double>
    compute_errors(const HYPRE * scheme,
                   const HArDCore2D::NavierStokesSolutions::IExactSolution * isolution,
                   const Eigen::VectorXd & uh,
                   const Eigen::VectorXd & ph,
                   const Eigen::VectorXd & bh,
                   const Eigen::VectorXd & rh,
                   const ComputeErrorsParameters & parameters = {},
           const double & t = 0.
                   );
 
    std::map<std::string, double>
    compute_l2_dof_errors(
                          const HYPRE * scheme,
                          const Eigen::VectorXd & uph,
                          const Eigen::VectorXd & upI,
                          const Eigen::VectorXd & brh,
                          const Eigen::VectorXd & brI,
                          const ComputeErrorsParameters & parameters
                          );
 
 
    std::map<std::string, double>
    compute_energy_errors(
                          const HYPRE * scheme,
                          const HArDCore2D::NavierStokesSolutions::IExactSolution * isolution,
                          const Eigen::VectorXd & uph,
                          const Eigen::VectorXd & upI,
                          const Eigen::VectorXd & brh,
                          const Eigen::VectorXd & brI,
                          const double & viscosity,
                          const double & magnetic_diffusivity,
                          const double & current_time,
                          const ComputeErrorsParameters & parameters = {}
                          );                      
    
    std::map<std::string, double>
    compute_energy_errors(
                          const HYPRE * scheme,
                          const HArDCore2D::NavierStokesSolutions::IExactSolution * isolution,
                          const Eigen::VectorXd & uph_half,
                          const Eigen::VectorXd & upI_half,
                          const Eigen::VectorXd & uph,
                          const Eigen::VectorXd & upI,
                          const Eigen::VectorXd & brh_half,
                          const Eigen::VectorXd & brI_half,
                          const Eigen::VectorXd & brh,
                          const Eigen::VectorXd & brI,
                          const double & viscosity,
                          const double & magnetic_diffusivity,
                          const double & time_step,
                          const double & current_time,
                          const ComputeErrorsParameters & parameters = {}
                          );
 
  } // namespace DSL
 
} // namespace HArDCore2D
 
#endif