wmstructuredgeometry.h

Go to the documentation of this file.

#ifndef WMSTRUCTUREDGEOMETRY_H_
#define WMSTRUCTUREDGEOMETRY_H_
 
// warpxm includes
#include "warpxm/warpxm_config.h"
#include "basis/wmbasisarrayset.h"
 
// STL includes
#include <vector>
 
class WmStructuredGeometry
{
public:
    WmStructuredGeometry(const WmBasisArraySet& PhaseSpaceBasisArraySet);
 
    virtual ~WmStructuredGeometry();
 
    void setup();
 
    virtual void generateNormals() = 0;
 
    virtual void generatePositionDirectedFaceNumbers() = 0;
 
    virtual void generateVelocityDirectedFaceNumbers() = 0;
 
    virtual int getOppositeVelocityDirectedFaceNode(const int globalNodeNumber,
                                                    const int faceNumber) = 0;
 
    const std::vector<real>& getRealNormal(const int face)
    {
        return _real_normals[face];
    }
 
    const std::vector<real>& getIsoparametricNormal(const int face)
    {
        return _isoparametric_normals[face];
    }
 
    const std::vector<int>& getPositionDirectedFaceNumbers()
    {
        return _position_directed_faces;
    }
 
    const std::vector<int>& getVelocityDirectedFaceNumbers()
    {
        return _velocity_directed_faces;
    }
 
    int getNodeOrientation()
    {
        return _node_orientation;
    }
 
    const std::vector<int>& getVelocityDirectedOppositeFaceNumbers()
    {
        return _velocity_directed_opposite_faces;
    }
 
protected:
    std::vector<int> _position_directed_faces; // faces with some component of normal in
                                               // physical space direction (x)
 
    int _num_position_directed_faces;
 
    std::vector<int> _velocity_directed_faces; // faces with no copoonent of normal in
                                               // physical space direction
    int _num_velocity_directed_faces;
 
    std::vector<int>
        _velocity_directed_opposite_faces; // the opposite faces with no copoonent of
                                           // normal in physical space direction
 
    std::vector<std::vector<real>>
        _real_normals; // normals real space [N_f_phase_space] x [N_d_phase_space]
 
    std::vector<std::vector<real>>
        _isoparametric_normals; // isoparametric real space [N_fp_phase_space] x
                                // [N_d_phase_space]
 
    int _node_orientation;
 
    WmBasisArraySet _phase_space_basis_array_set;
};
 
#endif