t8_geometry.h

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_H
#define T8_GEOMETRY_H
 
#include <t8.h>
#include <t8_refcount.h>
 
typedef enum t8_geometry_type {
  T8_GEOMETRY_TYPE_ZERO = 0,
  T8_GEOMETRY_TYPE_LINEAR,
  T8_GEOMETRY_TYPE_LINEAR_AXIS_ALIGNED,
  T8_GEOMETRY_TYPE_ANALYTIC,
  T8_GEOMETRY_TYPE_OCC,
  T8_GEOMETRY_TYPE_COUNT,
  T8_GEOMETRY_TYPE_UNDEFINED
} t8_geometry_type_t;
 
typedef struct t8_geometry t8_geometry_c;
 
/* The t8_geometry_c type must be know to cmesh.h, thus we
 * include it after the typedef. */
#include <t8_cmesh.h>
 
typedef struct t8_geometry_handler
{
  sc_array_t registered_geometries;
  t8_geometry_c *active_geometry;
  t8_gloidx_t active_tree;
  int is_committed;
  t8_refcount_t rc;
} t8_geometry_handler_t;
 
T8_EXTERN_C_BEGIN ();
 
void
t8_geom_handler_init (t8_geometry_handler_t **pgeom_handler);
 
void
t8_geom_handler_ref (t8_geometry_handler_t *geom_handler);
 
void
t8_geom_handler_unref (t8_geometry_handler_t **pgeom_handler);
 
void
t8_geom_handler_destroy (t8_geometry_handler_t **pgeom_handler);
 
void
t8_geom_handler_register_geometry (t8_geometry_handler_t *geom_handler, const t8_geometry_c *geometry);
 
void
t8_geom_handler_commit (t8_geometry_handler_t *geom_handler);
 
/* Check if a geometry handler was committed. */
int
t8_geom_handler_is_committed (const t8_geometry_handler_t *geom_handler);
 
size_t
t8_geom_handler_get_num_geometries (const t8_geometry_handler_t *geom_handler);
 
const t8_geometry_c *
t8_geom_handler_get_unique_geometry (const t8_geometry_handler_t *geom_handler);
 
t8_geometry_c *
t8_geom_handler_find_geometry (const t8_geometry_handler_t *geom_handler, const char *name);
 
void
t8_geometry_evaluate (t8_cmesh_t cmesh, t8_gloidx_t gtreeid, const double *ref_coords, const size_t num_coords,
                      double *out_coords);
 
void
t8_geometry_jacobian (t8_cmesh_t cmesh, t8_gloidx_t gtreeid, const double *ref_coords, const size_t num_coords,
                      double *jacobian);
 
t8_geometry_type_t
t8_geometry_get_type (t8_cmesh_t cmesh, t8_gloidx_t gtreeid);
 
T8_EXTERN_C_END ();
 
#endif /* !T8_GEOMETRY_H! */