singleboundary.h

Go to the documentation of this file.

 
#pragma once
 
// C++ STL
#include <memory>
 
// here goes PETSc headers
#include <petscsys.h>
#include <petscvec.h>
 
// here goes headers from our PetIBM
#include <petibm/mesh.h>
#include <petibm/type.h>
 
namespace petibm
{
namespace boundary
{
class SingleBoundaryBase
{
public:
    PetscInt dim;
 
    type::BCLoc loc;
 
    type::Field field;
 
    type::BCType type;
 
    PetscReal value;
 
    PetscReal normal;
 
    type::GhostPointsList points;
 
    PetscBool onThisProc;
 
    SingleBoundaryBase(const type::Mesh &mesh, const type::BCLoc &loc,
                       const type::Field &field, const type::BCType &type,
                       const PetscReal &value);
 
    SingleBoundaryBase() = default;
 
    virtual ~SingleBoundaryBase();
 
    virtual PetscErrorCode destroy();
 
    PetscErrorCode setGhostICs(const Vec &vec);
 
    PetscErrorCode updateEqs(const Vec &vec, const PetscReal &dt);
 
    PetscErrorCode updateGhostValues(const Vec &vec);
 
    PetscErrorCode copyValues2LocalVec(Vec &lclVec);
 
protected:
    PetscErrorCode init(const type::Mesh &mesh, const type::BCLoc &loc,
                        const type::Field &field, const type::BCType &type,
                        const PetscReal &bcValue);
 
    virtual PetscErrorCode setGhostICsKernel(const PetscReal &targetValue,
                                             type::GhostPointInfo &p) = 0;
 
    virtual PetscErrorCode updateEqsKernel(const PetscReal &targetValue,
                                           const PetscReal &dt,
                                           type::GhostPointInfo &p) = 0;
 
    MPI_Comm comm;
 
    PetscMPIInt mpiSize;
 
    PetscMPIInt mpiRank;
 
    type::Mesh mesh;
 
};  // SingleBoundaryBase
 
}  // end of namespace boundary
 
namespace type
{
typedef std::shared_ptr<boundary::SingleBoundaryBase> SingleBoundary;
}  // end of namespace type
 
namespace boundary
{
PetscErrorCode createSingleBoundary(const type::Mesh &mesh,
                                    const type::BCLoc &loc,
                                    const type::Field &field,
                                    const PetscReal &value,
                                    const type::BCType &bcType,
                                    type::SingleBoundary &singleBd);
}  // end of namespace boundary
 
}  // end of namespace petibm