solution.h

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
#include <string>
 
#include <petscsys.h>
 
#include <petibm/mesh.h>
 
namespace petibm
{
namespace solution
{
class SolutionBase
{
public:
    PetscInt dim;
 
    Vec UGlobal;
 
    Vec pGlobal;
 
    std::string info;
 
    SolutionBase() = default;
 
    SolutionBase(const type::Mesh &mesh){};
 
    virtual ~SolutionBase();
 
    virtual PetscErrorCode destroy();
 
    PetscErrorCode printInfo() const;
 
    virtual PetscErrorCode setInitialConditions(const YAML::Node &node) = 0;
 
    virtual PetscErrorCode convert2Velocity(const Mat &Rinv) = 0;
 
    virtual PetscErrorCode convert2Flux(const Mat &R) = 0;
 
    virtual PetscErrorCode write(const std::string &filePath) const = 0;
 
    virtual PetscErrorCode read(const std::string &filePath) = 0;
 
protected:
    virtual PetscErrorCode init(const type::Mesh &mesh) = 0;
 
    MPI_Comm comm;
 
    PetscMPIInt mpiSize;
 
    PetscMPIInt mpiRank;
 
    type::Mesh mesh;
 
};  // SolutionBase
}  // end of namespace solution
 
namespace type
{
typedef std::shared_ptr<solution::SolutionBase> Solution;
}  // end of namespace type
 
namespace solution
{
PetscErrorCode createSolution(const type::Mesh &mesh, type::Solution &solution);
}  // end of namespace solution
 
}  // end of namespace petibm