Standalone scheme

src/CMakeLists.txt

@@ -156,6 +156,7 @@ target_sources( T8 PRIVATE
     t8_schemes/t8_default/t8_default_tri/t8_dtri_connectivity.c
     t8_schemes/t8_default/t8_default_vertex/t8_default_vertex.cxx
     t8_schemes/t8_default/t8_default_vertex/t8_dvertex_bits.c
+    t8_schemes/t8_standalone/t8_standalone.cxx
     t8_vtk/t8_vtk_polydata.cxx 
     t8_vtk/t8_vtk_unstructured.cxx 
     t8_vtk/t8_vtk_parallel.cxx 

src/t8_schemes/t8_scheme.hxx

@@ -41,6 +41,8 @@
 #include <t8_schemes/t8_default/t8_default_tet/t8_default_tet.hxx>
 #include <t8_schemes/t8_default/t8_default_prism/t8_default_prism.hxx>
 #include <t8_schemes/t8_default/t8_default_pyramid/t8_default_pyramid.hxx>
+#include <t8_schemes/t8_standalone/t8_standalone.hxx>
+#include <t8_schemes/t8_standalone/t8_standalone_implementation.hxx>
 #include <string>
 #if T8_ENABLE_DEBUG
 // Only needed for t8_debug_print_type
@@ -92,7 +94,11 @@ class t8_scheme {
                                 t8_default_scheme_hex,
                                 t8_default_scheme_tet,
                                 t8_default_scheme_prism,
-                                t8_default_scheme_pyramid
+                                t8_default_scheme_pyramid,
+                                t8_standalone_scheme<T8_ECLASS_VERTEX>,
+                                t8_standalone_scheme<T8_ECLASS_LINE>,
+                                t8_standalone_scheme<T8_ECLASS_QUAD>,
+                                t8_standalone_scheme<T8_ECLASS_HEX>
                                 >;
   /* clang-format on */
 

src/t8_schemes/t8_standalone/t8_standalone.cxx

@@ -0,0 +1,72 @@
+/*
+  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) 2025 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.
+*/
+
+#include <t8_schemes/t8_standalone/t8_standalone.hxx>
+#include <t8_schemes/t8_scheme_builder.hxx>
+
+/* We want to export the whole implementation to be callable from "C" */
+T8_EXTERN_C_BEGIN ();
+
+t8_scheme *
+t8_scheme_new_standalone (void)
+{
+  t8_scheme_builder builder;
+
+  builder.add_eclass_scheme<t8_standalone_scheme<T8_ECLASS_VERTEX>> ();
+  builder.add_eclass_scheme<t8_standalone_scheme<T8_ECLASS_LINE>> ();
+  builder.add_eclass_scheme<t8_standalone_scheme<T8_ECLASS_QUAD>> ();
+  builder.add_eclass_scheme<t8_default_scheme_tri> ();
+  builder.add_eclass_scheme<t8_standalone_scheme<T8_ECLASS_HEX>> ();
+  builder.add_eclass_scheme<t8_default_scheme_tet> ();
+  builder.add_eclass_scheme<t8_default_scheme_prism> ();
+  builder.add_eclass_scheme<t8_default_scheme_pyramid> ();
+  // builder.add_eclass_scheme<t8_standalone_scheme<T8_ECLASS_PYRAMID>> (); //ADD if PYRAMID in "element_get_shape"
+  return builder.build_scheme ();
+}
+
+int
+t8_eclass_scheme_is_standalone (const t8_scheme *scheme, const t8_eclass_t eclass)
+{
+  switch (eclass) {
+  case T8_ECLASS_VERTEX:
+    return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_VERTEX>> (T8_ECLASS_VERTEX);
+  case T8_ECLASS_LINE:
+    return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_LINE>> (T8_ECLASS_LINE);
+  case T8_ECLASS_QUAD:
+    return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_QUAD>> (T8_ECLASS_QUAD);
+  // case T8_ECLASS_TRIANGLE:
+  //   return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_TRIANGLE>> (T8_ECLASS_TRIANGLE);
+  case T8_ECLASS_HEX:
+    return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_HEX>> (T8_ECLASS_HEX);
+  // case T8_ECLASS_TET:
+  //   return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_PRISM>> (T8_ECLASS_TET);
+  // case T8_ECLASS_PRISM:
+  //   return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_PRISM>> (T8_ECLASS_PRISM);
+  // case T8_ECLASS_PYRAMID:
+  //   return scheme->check_eclass_scheme_type<t8_standalone_scheme<T8_ECLASS_PYRAMID>> (T8_ECLASS_PYRAMID);
+  default:
+    SC_ABORT_NOT_REACHED ();
+  }
+  return 0; /* Default return value false */
+}
+
+T8_EXTERN_C_END ();

src/t8_schemes/t8_standalone/t8_standalone.hxx

@@ -0,0 +1,46 @@
+/*
+  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) 2025 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_STANDALONE_HXX
+#define T8_STANDALONE_HXX
+
+#include <t8_schemes/t8_scheme.hxx>
+#include <t8_schemes/t8_standalone/t8_standalone_implementation.hxx>
+
+T8_EXTERN_C_BEGIN ();
+
+/** Return the standalone element implementation of t8code. */
+t8_scheme *
+t8_scheme_new_standalone (void);
+
+/** Check whether a given eclass_scheme is one of the standalone schemes.
+ * \param [in] scheme   A (pointer to a) scheme
+ * \param [in] eclass   The eclass to check
+ * \return              True (non-zero) if \a scheme is one of the default schemes,
+ *                      false (zero) otherwise.
+ */
+int
+t8_eclass_scheme_is_standalone (const t8_scheme *scheme, const t8_eclass_t eclass);
+
+T8_EXTERN_C_END ();
+
+#endif /* !T8_STANDALONE_HXX */

src/t8_schemes/t8_standalone/t8_standalone_elements.hxx

@@ -0,0 +1,84 @@
+/*
+  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) 2025 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_STANDALONE_ELEMENTS_HXX
+#define T8_STANDALONE_ELEMENTS_HXX
+
+#include <t8.h>
+#include <array>
+#include <bitset>
+
+#define t8_standalone_element t8_standalone
+
+constexpr uint8_t T8_ELEMENT_DIM[T8_ECLASS_COUNT] = { 0, 1, 2, 2, 3, 3, 3, 3 };
+constexpr uint8_t T8_ELEMENT_MAXLEVEL[T8_ECLASS_COUNT] = { 255, 30, 29, 29, 21, 21, 21, 18 };
+constexpr uint8_t T8_ELEMENT_MAX_NUM_FACES[T8_ECLASS_COUNT] = { 1, 2, 4, 3, 6, 4, 5, 5 };
+constexpr uint8_t T8_ELEMENT_NUM_CHILDREN[T8_ECLASS_COUNT] = { 1, 2, 4, 4, 8, 8, 8, 10 };
+constexpr uint8_t T8_ELEMENT_NUM_CORNERS[T8_ECLASS_COUNT] = { 1, 2, 4, 3, 8, 4, 6, 5 };
+constexpr uint8_t T8_ELEMENT_NUM_FACES[T8_ECLASS_COUNT] = { 0, 2, 4, 3, 6, 4, 5, 5 };
+constexpr uint8_t T8_ELEMENT_MAX_NUM_FACECHILDREN[T8_ECLASS_COUNT] = { 0, 1, 2, 2, 4, 4, 4, 4 };
+
+constexpr uint8_t T8_ELEMENT_NUM_EQUATIONS[T8_ECLASS_COUNT] = { 0, 0, 0, 1, 0, 3, 1, 2 };
+
+/**PARENT CHILD BIJECTION*/
+template <t8_eclass_t TEclass>
+constexpr int8_t t8_element_type_Iloc_to_childtype[1 << T8_ELEMENT_NUM_EQUATIONS[TEclass]]
+                                                  [T8_ELEMENT_NUM_CHILDREN[TEclass]];
+template <t8_eclass_t TEclass>
+constexpr int8_t t8_element_type_Iloc_to_childcubeid[1 << T8_ELEMENT_NUM_EQUATIONS[TEclass]]
+                                                    [T8_ELEMENT_NUM_CHILDREN[TEclass]];
+template <t8_eclass_t TEclass>
+constexpr int8_t t8_element_type_cubeid_to_parenttype[1 << T8_ELEMENT_NUM_EQUATIONS[TEclass]]
+                                                     [1 << T8_ELEMENT_DIM[TEclass]];
+template <t8_eclass_t TEclass>
+constexpr int8_t t8_element_type_cubeid_to_Iloc[1 << T8_ELEMENT_NUM_EQUATIONS[TEclass]][1 << T8_ELEMENT_DIM[TEclass]];
+
+/**TYPE EQUATIONS*/
+template <t8_eclass_t TEclass>
+constexpr int8_t t8_type_edge_equations[T8_ELEMENT_NUM_EQUATIONS[TEclass]][2];
+
+/**VERTEX*/
+template <t8_eclass_t TEclass>
+constexpr int8_t t8_type_vertex_dim_to_binary[1 << T8_ELEMENT_NUM_EQUATIONS[TEclass]][T8_ELEMENT_NUM_CORNERS[TEclass]]
+                                             [T8_ELEMENT_DIM[TEclass]];
+
+typedef uint32_t t8_element_coord;
+typedef uint8_t t8_element_level;
+typedef int8_t t8_cube_id;
+
+template <t8_eclass_t TEclass>
+using t8_element_type = std::bitset<T8_ELEMENT_NUM_EQUATIONS[TEclass]>;
+
+template <t8_eclass_t TEclass>
+using t8_element_coords = std::array<t8_element_coord, T8_ELEMENT_DIM[TEclass]>;
+
+template <t8_eclass_t TEclass>
+struct t8_standalone_element
+{
+  /** The refinement level of the element relative to the root at level 0. */
+  t8_element_level level;
+
+  /** The coordinates of the anchor vertex of the element. */
+  t8_element_coords<TEclass> coords;
+};
+
+#endif /* T8_STANDALONE_ELEMENTS_HXX */