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 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,
                                          double *out_coords);
 
void                t8_geometry_jacobian (t8_cmesh_t cmesh,
                                          t8_gloidx_t gtreeid,
                                          const double *ref_coords,
                                          double *jacobian);
 
T8_EXTERN_C_END ();
 
#endif /* !T8_GEOMETRY_H! */