lininterp.h

Go to the documentation of this file.

 
#pragma once
 
#include <petscksp.h>
#include <petscsys.h>
 
#include <petibm/mesh.h>
#include <petibm/type.h>
 
namespace petibm
{
namespace misc
{
class LinInterpBase
{
public:
    LinInterpBase() = default;
 
    ~LinInterpBase();
 
    PetscErrorCode destroy();
 
    virtual PetscErrorCode interpolate(const DM &da, const Vec &vec, PetscReal &v) = 0;
 
protected:
    type::RealVec1D target;
 
    type::RealVec1D bl;
 
    type::RealVec1D tr;
 
    type::IntVec1D idxDirs;
 
    MPI_Comm comm;
 
    PetscMPIInt commSize;
 
    PetscMPIInt commRank;
 
    PetscErrorCode init(const MPI_Comm &comm,
                        const type::RealVec1D &point,
                        const type::Mesh &mesh,
                        const type::Field &field);
 
    PetscErrorCode getBLGridlineIndices(const type::Mesh &mesh,
                                        const type::Field &field);
 
    PetscErrorCode getBoxCoords(const type::Mesh &mesh,
                                const type::Field &field);
 
};  // LinInterpBase
 
class TriLinInterp: public LinInterpBase
{
public:
    TriLinInterp(const MPI_Comm &comm,
                 const type::RealVec1D &point,
                 const type::Mesh &mesh,
                 const type::Field &field);
    
    ~TriLinInterp() = default;
 
    PetscErrorCode interpolate(const DM &da, const Vec &vec, PetscReal &v);
 
};  // TriLinInterp
 
class BiLinInterp: public LinInterpBase
{
public:
    BiLinInterp(const MPI_Comm &comm,
                const type::RealVec1D &point,
                const type::Mesh &mesh,
                const type::Field &field);
    
    ~BiLinInterp() = default;
 
    PetscErrorCode interpolate(const DM &da, const Vec &vec, PetscReal &v);
 
};  // BiLinInterp
 
}  // end of namespace misc
 
namespace type
{
typedef std::shared_ptr<misc::LinInterpBase> LinInterp;
 
}  // end of namespace type
 
namespace misc
{
PetscErrorCode createLinInterp(const MPI_Comm &comm,
                               const type::RealVec1D &point,
                               const type::Mesh &mesh,
                               const type::Field &field,
                               type::LinInterp &interp);
 
}  // end of namespace misc
 
}  // end of namespace petibm