Mesh.h

//
// Class Mesh
//   The Mesh base class. Right now, this mainly acts as a standard base
//   class for all meshes so that other objects can register as users of
//   the mesh and can be notified if the mesh changes (e.g., it is rescaled
//   or restructured entirely).
//
#ifndef IPPL_MESH_H
#define IPPL_MESH_H
 
#include "Types/Vector.h"
 
namespace ippl {
    template <typename T, unsigned Dim>
    class Mesh {
    public:
        typedef T value_type;
        enum {
            Dimension = Dim
        };
 
        typedef Vector<T, Dim> vector_type;
        typedef Vector<vector_type, Dim> matrix_type;
 
        KOKKOS_INLINE_FUNCTION Mesh(){};
 
        KOKKOS_INLINE_FUNCTION virtual ~Mesh(){};
 
        // Get the origin of mesh vertex positions
        KOKKOS_INLINE_FUNCTION vector_type getOrigin() const;
 
        // Set the origin of mesh vertex positions
        KOKKOS_INLINE_FUNCTION void setOrigin(const vector_type& origin);
 
        KOKKOS_INLINE_FUNCTION const vector_type& getGridsize() const;
 
        KOKKOS_INLINE_FUNCTION virtual T getCellVolume() const = 0;
 
        KOKKOS_INLINE_FUNCTION virtual T getMeshVolume() const = 0;
 
        KOKKOS_INLINE_FUNCTION T getGridsize(size_t dim) const;
 
        KOKKOS_INLINE_FUNCTION virtual vector_type getVertexPosition(
            const NDIndex<Dim>& ndi) const = 0;
 
        KOKKOS_INLINE_FUNCTION virtual const vector_type& getMeshSpacing() const = 0;
 
        KOKKOS_INLINE_FUNCTION virtual vector_type getDeltaVertex(
            const NDIndex<Dim>& ndi) const = 0;
 
    protected:
        vector_type origin_m;     // Origin of mesh coordinates (vertices)
        vector_type gridSizes_m;  // Sizes (number of vertices)
    };
}  // namespace ippl
 
#include "Meshes/Mesh.hpp"
 
#endif