simplegrid.hxx

Go to the documentation of this file.

#ifndef NUMCXX_SIMPLEGRID_HXX
#define NUMCXX_SIMPLEGRID_HXX


#include "tarray1.hxx"
#include "tarray2.hxx"
#include "geometry.hxx"

namespace numcxx
{
  class SimpleGrid
  {
  public:

    SimpleGrid(){};
    SimpleGrid(const SimpleGrid& g)=delete;

    SimpleGrid(const Geometry &geometry, const char *triangle_flags);

    static std::shared_ptr<SimpleGrid> create(const Geometry &geometry, const std::string triangle_flags) { return std::make_shared<SimpleGrid>(geometry,triangle_flags.c_str());}

    static std::shared_ptr<SimpleGrid> create(std::shared_ptr<Geometry> geometry, const char * triangle_flags) { return std::make_shared<SimpleGrid>(*geometry,triangle_flags);}


    const TArray2<double>& get_points() const { return *points;}; 
    
    const TArray2<int>& get_cells() const { return *cells;};

    const TArray1<int>& get_cellregions() const {return  *cellregions;};

    const TArray2<int>& get_bfaces() const { return *bfaces;};

    const TArray1<int>& get_bfaceregions() const { return *bfaceregions;};

    const int spacedim()  const { return points->shape(1);}

    const int griddim() const { return bfaces->shape(1);}

    const int ncells() const {return cells->shape(0);}

    const int npoints() const {return points->shape(0);}

    const int nbfaces() const {return bfaces->shape(0);}

    void calc_hminmax(double& hmin, double& hmax) const;
    
  private:
    std::shared_ptr<TArray2<double>> points=nullptr;
    std::shared_ptr<TArray2<int>> cells=nullptr;
    std::shared_ptr<TArray1<int>> cellregions=nullptr;
    std::shared_ptr<TArray2<int>> bfaces=nullptr;
    std::shared_ptr<TArray1<int>> bfaceregions=nullptr;


  };
}


#endif