t8_geometry_base.hxx

Go to the documentation of this file.

/*
  This file is part of t8code.
  t8code is a C library to manage a collection (a forest) of multiple
  connected adaptive space-trees of general element classes in parallel.
 
  Copyright (C) 2015 the developers
 
  t8code is free software; you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation; either version 2 of the License, or
  (at your option) any later version.
 
  t8code is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.
 
  You should have received a copy of the GNU General Public License
  along with t8code; if not, write to the Free Software Foundation, Inc.,
  51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
 
#ifndef T8_GEOMETRY_BASE_HXX
#define T8_GEOMETRY_BASE_HXX
 
#include <t8.h>
#include <t8_cmesh.h>
 
T8_EXTERN_C_BEGIN ();
 
struct t8_geometry
{
public:
 
  /* Basic constructor that sets the dimension, the name, and the name for the attribute. */
  t8_geometry (int dimension, const char *name, const char *attribute_name =
               NULL)
:  dimension (dimension), name (name) {
  }
 
  /* Base constructor with no arguments. We need this since it
   * is called from derived class constructors.
   * Sets dimension and name to invalid values. */
  t8_geometry         ():t8_geometry (-1, "Invalid")
  {
  }
 
  virtual ~ t8_geometry () {
  }
 
  virtual void        t8_geom_evaluate (t8_cmesh_t cmesh,
                                        t8_gloidx_t gtreeid,
                                        const double *ref_coords,
                                        double out_coords[3]) const = 0;
 
  virtual void        t8_geom_evalute_jacobian (t8_cmesh_t cmesh,
                                                t8_gloidx_t gtreeid,
                                                const double *ref_coords,
                                                double *jacobian) const = 0;
 
  virtual void        t8_geom_load_tree_data (t8_cmesh_t cmesh,
                                              t8_gloidx_t gtreeid) = 0;
 
  inline int          t8_geom_get_dimension () const
  {
    return dimension;
  }
 
  inline const char  *t8_geom_get_name () const
  {
    return name;
  }
 
protected:
 
  int                 dimension;
  const char         *name;
};
 
class               t8_geometry_w_vertices:public t8_geometry
{
public:
  /* Basic constructor that sets the dimension, the name, and the name for the attribute. */
  t8_geometry_w_vertices (int dimension, const char *name,
                          const char *attribute_name = NULL)
:  t8_geometry (dimension, name, attribute_name) {
    active_tree_vertices = NULL;
    active_tree = -1;
  }
 
  /* Base constructor with no arguments. We need this since it
   * is called from derived class constructors.
   * Sets dimension and name to invalid values. */
  t8_geometry_w_vertices ():t8_geometry_w_vertices (-1, "Invalid")
  {
    active_tree_vertices = NULL;
    active_tree = -1;
  }
 
  virtual ~ t8_geometry_w_vertices () {
  }
 
  virtual void        t8_geom_load_tree_data (t8_cmesh_t cmesh,
                                              t8_gloidx_t gtreeid);
 
protected:
  t8_gloidx_t         active_tree;      /*< The tree of which currently vertices are loaded. */
  t8_eclass_t         active_tree_class;        /*< The class of the currently active tree. */
  const double       *active_tree_vertices;     /*< The vertices of the currently active tree. */
};
 
T8_EXTERN_C_END ();
 
#endif /* !T8_GEOMETRY_BASE_HXX! */