\n Classes
Implement iterators over GridGlue intersections.
Model of the Intersection concept provided by GridGlue.
Describes the parallel communication interface class for Dune::GridGlue.
Model of the Intersection concept provided by GridGlue.
Describes the parallel communication interface class for Dune::GridGlue.
Definition: gridglue.hh:37
sequential adapter to couple two grids at specified close together boundaries
Definition: gridglue.hh:67
Dune::FieldVector< ctype, dimworld > Coords
The type used for coordinate vectors.
Definition: gridglue.hh:174
Intersection getIntersection(int i) const
Definition: gridglue.hh:388
GridVertex< 1 > Grid1Vertex
The type of the Grid1 vertices.
Definition: gridglue.hh:194
Dune::GridGlue::IntersectionIndexSet< P0, P1 > IndexSet
Type of remote intersection indexSet.
Definition: gridglue.hh:206
\n \n- \n+ \n
\n gridglue.cc File Reference
rangegenerators.hh File Reference
\n-
\n-
\n-
\n-
\n-
\n-
\n-\n+\n+
\n
\n
\n \n \n
\n Generated by![\"doxygen\"/](\"doxygen.svg\")
1.9.4\n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,28 +5,28 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Namespaces | Functions\n-gridglue.cc File Reference\n-#include \"intersection.hh\"\n-#include \n-#include \n-#include \"../gridglue.hh\"\n-#include \"../common/ringcomm.hh\"\n-#include \n+Classes | Namespaces | Functions\n+rangegenerators.hh File Reference\n+#include \n+Go_to_the_source_code_of_this_file.\n+ Classes\n+struct \u00a0Dune::GridGlue::Reverse<_reverse_>\n+\u00a0\n Namespaces\n namespace \u00a0Dune\n \u00a0\n namespace \u00a0Dune::GridGlue\n \u00a0\n Functions\n-template\n-void\u00a0Dune::GridGlue::printVector (const std::vector< T > &v, std::string name,\n- int rank)\n+template<... >\n+IteratorRange<... >\u00a0Dune::GridGlue::intersections (const GridGlue<... > &glue,\n+ const Reverse<... > &reverse=!reversed)\n+\u00a0 Iterate over all intersections of a GridGlue. More...\n \u00a0\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00008.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00008.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: gridgluevtkwriter.hh File Reference \n+dune-grid-glue: gridglue.cc File Reference \n \n \n \n \n \n \n \n@@ -63,48 +63,39 @@\n \n \n \n+ \n \n
#include "intersection.hh"\n-
#include <vector>\n-
#include <iterator>\n-
#include "../gridglue.hh"\n-
#include "../common/ringcomm.hh"\n-
#include <dune/common/unused.hh>\n-
\n+Classes | |
| struct | Dune::GridGlue::Reverse< reverse > |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
\n Functions | |
| template<typename T > | |
| void | Dune::GridGlue::printVector (const std::vector< T > &v, std::string name, int rank) |
| template<... > | |
| IteratorRange<... > | Dune::GridGlue::intersections (const GridGlue<... > &glue, const Reverse<... > &reverse=!reversed) |
| Iterate over all intersections of a GridGlue. More... | |
\n Generated by
1.9.4\n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,28 +5,28 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Namespaces | Functions\n-gridglue.cc File Reference\n-#include \"intersection.hh\"\n-#include gridglue.cc File Reference
\n-\n-
\n-
\n-
\n+
\n
\n-
\n-
\n-
\n-
\n-\n-\n-
\n+\n
\n+
\n-\n+\n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,34 +5,28 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces\n-gridgluevtkwriter.hh File Reference\n-Write all remote intersections to a vtk file for debugging. More...\n-#include \n-#include \n-#include \n+Namespaces | Functions\n+gridglue.cc File Reference\n+#include \"intersection.hh\"\n #include \n-#include \n-#include \n-#include \n-#include \n-Go_to_the_source_code_of_this_file.\n- Classes\n-class \u00a0Dune::GridGlue::GridGlueVtkWriter\n-\u00a0 Write remote intersections to a vtk file for debugging purposes.\n- More...\n-\u00a0\n+#include \n+#include \"../gridglue.hh\"\n+#include \"../common/ringcomm.hh\"\n+#include \n Namespaces\n namespace \u00a0Dune\n \u00a0\n namespace \u00a0Dune::GridGlue\n \u00a0\n-***** Detailed Description *****\n-Write all remote intersections to a vtk file for debugging.\n+ Functions\n+template\n+void\u00a0Dune::GridGlue::printVector (const std::vector< T > &v, std::string name,\n+ int rank)\n+\u00a0\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00014.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00014.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: intersectionindexset.hh File Reference \n+dune-grid-glue: intersection.hh File Reference \n \n \n \n \n \n \n \n@@ -64,35 +64,74 @@\n \n \n \n \n+ \n \n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,47 +5,62 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces | Macros\n-intersection.hh File Reference\n-Model of the Intersection concept provided by GridGlue. More...\n-#include \n-#include \n-#include \n-#include \n-#include \n-#include \n-#include \n-#include \n+Classes | Namespaces | Typedefs | Functions\n+gridgluecommunicate.hh File Reference\n+Describes the parallel communication interface class for Dune::GridGlue.\n+More...\n+#include \n+#include \n+#include \n+#include \n+#include \n Go_to_the_source_code_of_this_file.\n Classes\n- class \u00a0Dune::GridGlue::IntersectionData<_P0,_P1_>\n-\u00a0 storage class for Dune::GridGlue::Intersection related data More...\n+struct \u00a0Dune::GridGlue::GlobalId\n \u00a0\n-struct \u00a0Dune::GridGlue::IntersectionTraits<_P0,_P1,_inside,_outside_>\n+ class \u00a0Dune::GridGlue::CommDataHandle<_DataHandleImp,_DataTypeImp_>\n+\u00a0 describes the features of a data handle for communication in parallel\n+ runs using the GridGlue::communicate methods. More...\n \u00a0\n- class \u00a0Dune::GridGlue::Intersection<_P0,_P1,_I,_O_>\n-\u00a0 The intersection of two entities of the two patches of a GridGlue.\n- More...\n+ class \u00a0Dune::GridGlue::StreamingMessageBuffer<_DT_>\n+\u00a0\n+ class \u00a0Dune::GridGlue::CommunicationOperator<_dir_>\n+\u00a0 forward gather scatter to user defined CommInfo class More...\n+\u00a0\n+struct \u00a0Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n+\u00a0 collects all GridGlue data requried for communication More...\n+\u00a0\n+struct \u00a0Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,\n+ DataTypeImp_>_>\n+\u00a0 specialization of the CommPolicy struct, required for the\n+ ParallelIndexsets More...\n \u00a0\n Namespaces\n namespace \u00a0Dune\n \u00a0\n namespace \u00a0Dune::GridGlue\n \u00a0\n- Macros\n-#define\u00a0ONLY_SIMPLEX_INTERSECTIONS\n+ Typedefs\n+ typedef std::pair< int, int >\u00a0Dune::GridGlue::RankPair\n+\u00a0\n+typedef CommunicationOperator< Dune::\n+ ForwardCommunication >\u00a0Dune::GridGlue::ForwardOperator\n+\u00a0\n+typedef CommunicationOperator< Dune::\n+ BackwardCommunication >\u00a0Dune::GridGlue::BackwardOperator\n+\u00a0\n+ Functions\n+std::ostream &\u00a0Dune::GridGlue::operator<< (std::ostream &os, const GlobalId\n+ &id)\n \u00a0\n ***** Detailed Description *****\n-Model of the Intersection concept provided by GridGlue.\n+Describes the parallel communication interface class for Dune::GridGlue.\n Author\n Christian Engwer\n-***** Macro Definition Documentation *****\n-***** \u25c6\u00a0ONLY_SIMPLEX_INTERSECTIONS *****\n-#define ONLY_SIMPLEX_INTERSECTIONS\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00017_source.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00017_source.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: intersection.hh Source File \n+dune-grid-glue: gridgluecommunicate.hh Source File \n \n \n \n \n \n \n \n@@ -62,525 +62,323 @@\n \n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,758 +5,443 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-intersection.hh\n+gridgluecommunicate.hh\n Go_to_the_documentation_of_this_file.\n 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-\n 2// vi: set et ts=4 sw=2 sts=2:\n 3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 11#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n- 12#define DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n- 13\n- 14#include \n- 15#include \n- 16#include \n- 17\n- 18#include \n- 19#include \n- 20#include \n- 21#include \n- 22#include \n+ 5#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUECOMMUNICATE_HH\n+ 6#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUECOMMUNICATE_HH\n+ 7\n+ 13#include \n+ 14\n+ 15#include \n+ 16#include \n+ 17#include \n+ 18#include \n+ 19\n+ 20\n+ 21namespace Dune {\n+ 22 namespace GridGlue {\n 23\n-24#define ONLY_SIMPLEX_INTERSECTIONS\n- 25\n- 26namespace Dune {\n- 27 namespace GridGlue {\n- 28\n- 29 // forward declaration\n- 30 template\n- 31 class IntersectionIndexSet;\n- 32\n- 36 template\n-37 class IntersectionData\n- 38 {\n- 39 public:\n-40 typedef ::Dune::GridGlue::GridGlue GridGlue;\n- 41\n-42 typedef typename GridGlue::IndexType IndexType;\n- 43\n-45 static constexpr int coorddim = GridGlue::dimworld;\n- 46\n- 47 private:\n- 48 // intermediate quantities\n- 49 template\n- 50 static constexpr int dim() { return GridGlue::template GridView::\n-Grid::dimension - GridGlue::template GridPatch::codim; }\n- 51\n- 52 public:\n-54 static constexpr int mydim = dim<0>() < dim<1>() ? dim<0>() : dim<1>();\n- 55\n- 56 template\n-57 using GridLocalGeometry = AffineGeometry<\n- 58 typename GridGlue::template GridView::ctype, mydim, GridGlue::\n-template GridView::dimension>;\n- 59\n-60 using Grid0LocalGeometry [[deprecated(\"please use GridLocalGeometry<0>\n-instead\")]] = GridLocalGeometry<0>;\n-61 using Grid1LocalGeometry [[deprecated(\"please use GridLocalGeometry<1>\n-instead\")]] = GridLocalGeometry<1>;\n+24 typedef std::pair RankPair;\n+25 struct GlobalId : public std::pair\n+ 26 {\n+30 GlobalId() {\n+ 31 this->first.first = 0;\n+ 32 this->first.second = 0;\n+ 33 this->second = 0;\n+ 34 }\n+38 GlobalId(int i) {\n+ 39 this->first.first = i;\n+ 40 this->first.second = i;\n+ 41 this->second = 0;\n+ 42 }\n+48 GlobalId(int i, int j, unsigned int n) {\n+ 49 this->first.first = std::min(i,j);\n+ 50 this->first.second = std::max(i,j);\n+ 51 this->second = n;\n+ 52 }\n+ 53 };\n+ 54\n+55 inline std::ostream& operator<<(std::ostream& os, const GlobalId & id)\n+ 56 {\n+ 57 os << \"(\"\n+ 58 << id.first.first << \",\" << id.first.second << \",\"\n+ 59 << id.second << \")\";\n+ 60 return os;\n+ 61 }\n 62\n- 63 template\n-64 using GridGeometry = AffineGeometry<\n- 65 typename GridGlue::template GridView::ctype, mydim, GridGlue::\n-template GridView::dimensionworld>;\n- 66\n-67 using Grid0Geometry [[deprecated(\"please use GridGeometry<0> instead\")]] =\n-GridGeometry<0>;\n-68 using Grid1Geometry [[deprecated(\"please use GridGeometry<1> instead\")]] =\n-GridGeometry<1>;\n- 69\n- 70 template\n-71 using GridIndexType = typename GridGlue::template GridView::IndexSet::\n-IndexType;\n- 72\n-73 using Grid0IndexType [[deprecated(\"please use GridIndexType<0> instead\")]] =\n-GridIndexType<0>;\n-74 using Grid1IndexType [[deprecated(\"please use GridIndexType<1> instead\")]] =\n-GridIndexType<1>;\n- 75\n- 77 IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned\n-int offset, bool grid0local, bool grid1local);\n- 78\n-80 IntersectionData() = default;\n+ 75 template \n+76 class CommDataHandle\n+ 77 {\n+ 78 public:\n+80 typedef DataTypeImp DataType;\n 81\n- 82 /* Accessor Functions */\n- 83\n- 84 template\n-85 const GridLocalGeometry& localGeometry(unsigned int parentId = 0)\n-const\n- 86 { return *std::get(sideData_).gridlocalgeom[parentId]; }\n+ 82 protected:\n+ 83 // one should not create an explicit instance of this inteface object\n+84 CommDataHandle() {}\n+ 85\n+ 86 public:\n 87\n- 88 template\n-89 const GridGeometry& geometry() const\n- 90 { return *std::get(sideData_).gridgeom; }\n- 91\n- 92 template\n-93 bool local() const\n- 94 { return std::get(sideData_).gridlocal; }\n- 95\n- 96 template\n-97 IndexType index(unsigned int parentId = 0) const\n- 98 { return std::get(sideData_).gridindices[parentId]; }\n- 99\n- 100 template\n-101 IndexType parents() const\n- 102 { return std::get(sideData_).gridindices.size(); }\n- 103\n- 104 private:\n- 105 template\n- 106 void initializeGeometry(const GridGlue& glue, unsigned mergeindex);\n- 107\n- 108 /* M E M B E R V A R I A B L E S */\n+ 91 template\n+92 size_t size (RISType& i) const\n+ 93 {\n+ 94 CHECK_INTERFACE_IMPLEMENTATION((asImp().size(i)));\n+ 95 return asImp().size(i);\n+ 96 }\n+ 97\n+ 103 template\n+104 void gather (MessageBufferImp& buff, const EntityType& e, const RISType &\n+i) const\n+ 105 {\n+ 106 MessageBufferIF buffIF(buff);\n+ 107 CHECK_AND_CALL_INTERFACE_IMPLEMENTATION((asImp().gather(buffIF,e,i)));\n+ 108 }\n 109\n- 110 public:\n-112 IndexType index_;\n- 113\n- 114 private:\n- 115 template\n- 116 struct SideData {\n- 118 bool gridlocal = false;\n- 119\n- 121 std::vector< GridIndexType > gridindices;\n- 122\n- 124 std::vector< std::optional< GridLocalGeometry > > gridlocalgeom;\n- 125\n- 133 std::optional< GridGeometry > gridgeom;\n- 134 };\n+ 117 template\n+118 void scatter (MessageBufferImp& buff, const EntityType& e, const RISType &\n+i, size_t n)\n+ 119 {\n+ 120 MessageBufferIF buffIF(buff);\n+ 121 CHECK_AND_CALL_INTERFACE_IMPLEMENTATION((asImp().scatter(buffIF,e,i,n)));\n+ 122 }\n+ 123\n+ 124 private:\n+ 126 DataHandleImp& asImp () {\n+ 127 return static_cast (*this);\n+ 128 }\n+ 130 const DataHandleImp& asImp () const\n+ 131 {\n+ 132 return static_cast(*this);\n+ 133 }\n+ 134 }; // end class CommDataHandleIF\n 135\n- 136 std::tuple< SideData<0>, SideData<1> > sideData_;\n- 137 };\n- 138\n- 139 template\n- 140 template\n- 141 void IntersectionData::initializeGeometry(const GridGlue& glue,\n-unsigned mergeindex)\n- 142 {\n- 143 auto& data = std::get(sideData_);\n+ 140 template\n+141 class StreamingMessageBuffer {\n+ 142 public:\n+143 typedef DT value_type;\n 144\n- 145 const unsigned n_parents = glue.merger_->template parents\n-(mergeindex);\n- 146\n- 147 // init containers\n- 148 data.gridindices.resize(n_parents);\n- 149 data.gridlocalgeom.resize(n_parents);\n- 150\n- 151 // default values\n- 152 data.gridindices[0] = 0;\n- 153\n- 154 static constexpr int nSimplexCorners = mydim + 1;\n- 155 using ctype = typename GridGlue::ctype;\n- 156\n- 157 // initialize the local and the global geometries of grid `side`\n- 158\n- 159 // compute the coordinates of the subface's corners in codim 0 entity\n-local coordinates\n- 160 static constexpr int elementdim = GridGlue::template GridView::\n-template Codim<0>::Geometry::mydimension;\n- 161\n- 162 // coordinates within the subentity that contains the remote intersection\n- 163 std::array() >, nSimplexCorners>\n-corners_subEntity_local;\n- 164\n- 165 for (unsigned int par = 0; par < n_parents; ++par) {\n- 166 for (int i = 0; i < nSimplexCorners; ++i)\n- 167 corners_subEntity_local[i] = glue.merger_->template parentLocal\n-(mergeindex, i, par);\n+ 145 // Constructor\n+146 StreamingMessageBuffer (DT *p)\n+ 147 {\n+ 148 a=p;\n+ 149 i=0;\n+ 150 j=0;\n+ 151 }\n+ 152\n+ 153 // write data to message buffer, acts like a stream !\n+ 154 template\n+155 void write (const Y& data)\n+ 156 {\n+ 157 static_assert(std::is_same::value, \"DataType mismatch\");\n+ 158 a[i++] = data;\n+ 159 }\n+ 160\n+ 161 // read data from message buffer, acts like a stream !\n+ 162 template\n+163 void read (Y& data) const\n+ 164 {\n+ 165 static_assert(std::is_same::value, \"DataType mismatch\");\n+ 166 data = a[j++];\n+ 167 }\n 168\n- 169 // Coordinates of the remote intersection corners wrt the element\n-coordinate system\n- 170 std::array, nSimplexCorners>\n-corners_element_local;\n- 171\n- 172 if (data.gridlocal) {\n- 173 data.gridindices[par] = glue.merger_->template parent\n-(mergeindex,par);\n- 174\n- 175 typename GridGlue::template GridPatch::LocalGeometry\n- 176 gridLocalGeometry = glue.template patch().geometryLocal\n-(data.gridindices[par]);\n- 177 for (std::size_t i=0; i::Geometry\n- 196 gridWorldGeometry = glue.template patch().geometry(data.gridindices\n-[par]);\n- 197\n- 198 // world coordinates of the remote intersection corners\n- 199 std::array::\n-dimensionworld>, nSimplexCorners> corners_global;\n- 200\n- 201 for (std::size_t i=0; i\n-213 IntersectionData::IntersectionData(const GridGlue& glue, unsigned\n-int mergeindex, unsigned int offset,\n- 214 bool grid0local, bool grid1local)\n- 215 : index_(mergeindex+offset)\n- 216 {\n- 217 // if an invalid index is given do not proceed!\n- 218 // (happens when the parent GridGlue initializes the \"end\"-Intersection)\n- 219 assert (0 <= mergeindex || mergeindex < glue.index__sz);\n+169 size_t counter() const { return i; }\n+ 170\n+171 void clear()\n+ 172 {\n+ 173 i = 0;\n+ 174 j = 0;\n+ 175 }\n+ 176\n+ 177 // we need access to these variables in an assertion\n+ 178#ifdef NDEBUG\n+ 179 private:\n+ 180#endif\n+181 DT *a;\n+182 size_t i;\n+183 mutable size_t j;\n+ 184 };\n+ 185\n+ 192 template\n+193 class CommunicationOperator\n+ 194 {\n+ 195 public:\n+ 196 template\n+197 static const typename CommInfo::DataType& gather(const CommInfo& commInfo,\n+size_t i, size_t j = 0)\n+ 198 {\n+ 199 // get Intersection\n+ 200 typedef typename CommInfo::GridGlue::Intersection Intersection;\n+ 201 Intersection ris(commInfo.gridglue->getIntersection(i));\n+ 202\n+ 203 // fill buffer if we have a new intersection\n+ 204 if (j == 0)\n+ 205 {\n+ 206 commInfo.mbuffer.clear();\n+ 207 if (dir == Dune::ForwardCommunication)\n+ 208 {\n+ 209 // read from grid0\n+ 210 if(ris.self())\n+ 211 commInfo.data->gather(commInfo.mbuffer, ris.inside(), ris);\n+ 212 }\n+ 213 else // (dir == Dune::BackwardCommunication)\n+ 214 {\n+ 215 // read from grid1\n+ 216 if(ris.neighbor())\n+ 217 commInfo.data->gather(commInfo.mbuffer, ris.outside(), ris.flip());\n+ 218 }\n+ 219 }\n 220\n- 221 std::get<0>(sideData_).gridlocal = grid0local;\n- 222 std::get<1>(sideData_).gridlocal = grid1local;\n- 223\n- 224 initializeGeometry<0>(glue, mergeindex);\n- 225 initializeGeometry<1>(glue, mergeindex);\n- 226 }\n- 227\n- 232 template\n-233 struct IntersectionTraits\n- 234 {\n-235 using GridGlue = ::Dune::GridGlue::GridGlue;\n-236 using IntersectionData = Dune::GridGlue::IntersectionData;\n- 237\n-238 using InsideGridView = typename GridGlue::template GridView;\n-239 using OutsideGridView = typename GridGlue::template GridView;\n- 240\n-241 using InsideLocalGeometry = typename IntersectionData::template\n-GridLocalGeometry;\n-242 using OutsideLocalGeometry = typename IntersectionData::template\n-GridLocalGeometry;\n- 243\n-244 using Geometry = typename IntersectionData::template GridGeometry;\n-245 using OutsideGeometry = typename IntersectionData::template\n-GridGeometry;\n- 246\n-247 static constexpr auto coorddim = IntersectionData::coorddim;\n-248 static constexpr auto mydim = IntersectionData::mydim;\n-249 static constexpr int insidePatch = inside;\n-250 static constexpr int outsidePatch = outside;\n- 251\n-252 using ctype = typename GridGlue::ctype;\n-253 using LocalCoordinate = Dune::FieldVector;\n-254 using GlobalCoordinate = Dune::FieldVector;\n- 255 };\n- 256\n- 259 template\n-260 class Intersection\n- 261 {\n+ 221 // return the j'th value in the buffer\n+ 222 assert(j < commInfo.mbuffer.i);\n+ 223 return commInfo.buffer[j];\n+ 224 }\n+ 225\n+ 226 template\n+227 static void scatter(CommInfo& commInfo, const typename CommInfo::DataType&\n+v, std::size_t i, std::size_t j = 0)\n+ 228 {\n+ 229 // extract GridGlue objects...\n+ 230 typedef typename CommInfo::GridGlue::Intersection Intersection;\n+ 231 Intersection ris(commInfo.gridglue->getIntersection(i));\n+ 232\n+ 233 // get size if we have a new intersection\n+ 234 if (j == 0)\n+ 235 {\n+ 236 commInfo.mbuffer.clear();\n+ 237 commInfo.currentsize = commInfo.data->size(ris);\n+ 238 }\n+ 239\n+ 240 // write entry to buffer\n+ 241 commInfo.buffer[j] = v;\n+ 242\n+ 243 // write back the buffer if we are at the end of this intersection\n+ 244 if (j == commInfo.currentsize-1)\n+ 245 {\n+ 246 if (dir == Dune::ForwardCommunication)\n+ 247 {\n+ 248 // write to grid1\n+ 249 if(ris.neighbor())\n+ 250 commInfo.data->scatter(commInfo.mbuffer, ris.outside(), ris.flip(),\n+commInfo.currentsize);\n+ 251 }\n+ 252 else // (dir == Dune::BackwardCommunication)\n+ 253 {\n+ 254 // write to grid0\n+ 255 if(ris.self())\n+ 256 commInfo.data->scatter(commInfo.mbuffer, ris.inside(), ris,\n+commInfo.currentsize);\n+ 257 }\n+ 258 assert(commInfo.mbuffer.j <= commInfo.currentsize);\n+ 259 }\n+ 260 }\n+ 261 };\n 262\n- 263 public:\n- 264\n-265 typedef IntersectionTraits Traits;\n- 266\n-267 typedef typename Traits::GridGlue GridGlue;\n-268 typedef typename Traits::IntersectionData IntersectionData;\n- 269\n- 270\n-271 typedef typename Traits::InsideGridView InsideGridView;\n-272 typedef typename Traits::InsideLocalGeometry InsideLocalGeometry;\n- 273\n-274 typedef typename Traits::OutsideGridView OutsideGridView;\n-275 typedef typename Traits::OutsideLocalGeometry OutsideLocalGeometry;\n-276 typedef typename Traits::OutsideGeometry OutsideGeometry;\n- 277\n-278 typedef typename Traits::Geometry Geometry;\n-279 typedef typename Traits::ctype ctype;\n- 280\n-281 typedef typename InsideGridView::Traits::template Codim<0>::Entity\n-InsideEntity;\n-282 typedef typename OutsideGridView::Traits::template Codim<0>::Entity\n-OutsideEntity;\n+263 typedef CommunicationOperator ForwardOperator;\n+264 typedef CommunicationOperator\n+BackwardOperator;\n+ 265\n+ 270 template \n+271 struct CommInfo\n+ 272 {\n+273 typedef DataTypeImp value_type;\n+274 typedef GG GridGlue;\n+275 typedef DataTypeImp DataType;\n+ 276\n+277 CommInfo() : buffer(100), mbuffer(&buffer[0])\n+ 278 {}\n+ 279\n+ 280 // tunnel information to the policy and the operators\n+281 const GridGlue * gridglue;\n+282 ::Dune::GridGlue::CommDataHandle * data;\n 283\n-284 typedef typename Traits::LocalCoordinate LocalCoordinate;\n-285 typedef typename Traits::GlobalCoordinate GlobalCoordinate;\n- 286\n-288 static constexpr auto coorddim = Traits::coorddim;\n- 289\n-291 static constexpr auto mydim = Traits::mydim;\n+ 284 // state variables\n+285 std::vector buffer;\n+286 mutable ::Dune::GridGlue::StreamingMessageBuffer mbuffer;\n+287 size_t currentsize;\n+288 Dune::CommunicationDirection dir;\n+ 289 };\n+ 290\n+ 291 } // end namespace GridGlue\n 292\n-294 static constexpr int insidePatch = Traits::insidePatch;\n- 295\n-297 static constexpr int outsidePatch = Traits::outsidePatch;\n- 298\n- 299 // typedef unsigned int IndexType;\n- 300\n- 301 private:\n- 305 static constexpr int codimensionWorld = coorddim - mydim;\n- 306\n- 307 public:\n- 308 /* C O N S T R U C T O R S */\n- 309\n-311 Intersection(const GridGlue* glue, const IntersectionData* i) :\n- 312 glue_(glue), i_(i) {}\n- 313\n- 314 /* F U N C T I O N A L I T Y */\n- 315\n- 318 InsideEntity\n-319 inside(unsigned int parentId = 0) const\n- 320 {\n- 321 assert(self());\n- 322 return glue_->template patch().element(i_->template index\n-(parentId));\n- 323 }\n- 324\n- 327 OutsideEntity\n-328 outside(unsigned int parentId = 0) const\n- 329 {\n- 330 assert(neighbor());\n- 331 return glue_->template patch().element(i_->template index\n-(parentId));\n- 332 }\n- 333\n-335 bool conforming() const\n- 336 {\n- 337 throw Dune::NotImplemented();\n- 338 }\n- 339\n-342 const InsideLocalGeometry& geometryInInside(unsigned int parentId = 0)\n-const\n- 343 {\n- 344 return i_->template localGeometry(parentId);\n- 345 }\n- 346\n-349 const OutsideLocalGeometry& geometryInOutside(unsigned int parentId = 0)\n-const\n- 350 {\n- 351 return i_->template localGeometry(parentId);\n- 352 }\n- 353\n-360 const Geometry& geometry() const\n- 361 {\n- 362 return i_->template geometry();\n- 363 }\n- 364\n-371 const OutsideGeometry& geometryOutside() const // DUNE_DEPRECATED\n- 372 {\n- 373 return i_->template geometry();\n- 374 }\n- 375\n-377 Dune::GeometryType type() const\n- 378 {\n- 379 #ifdef ONLY_SIMPLEX_INTERSECTIONS\n- 380 # if DUNE_VERSION_NEWER(DUNE_GEOMETRY, 2, 6)\n- 381 return Dune::GeometryTypes::simplex(mydim);\n- 382 # else\n- 383 static const Dune::GeometryType type(Dune::GeometryType::simplex, mydim);\n- 384 return type;\n- 385 # endif\n- 386 #else\n- 387 #error Not Implemented\n- 388 #endif\n- 389 }\n- 390\n- 391\n-393 bool self() const\n- 394 {\n- 395 return i_->template local();\n- 396 }\n- 397\n-399 size_t neighbor(unsigned int g = 0) const\n- 400 {\n- 401 if (g == 0 && i_->template local()) {\n- 402 return i_->template parents();\n- 403 } else if (g == 1 && i_->template local()) {\n- 404 return i_->template parents();\n- 405 }\n- 406 return 0;\n- 407 }\n- 408\n-410 int indexInInside(unsigned int parentId = 0) const\n- 411 {\n- 412 assert(self());\n- 413 return glue_->template patch().indexInInside(i_->template index\n-(parentId));\n- 414 }\n- 415\n-417 int indexInOutside(unsigned int parentId = 0) const\n- 418 {\n- 419 assert(neighbor());\n- 420 return glue_->template patch().indexInInside(i_->template index\n-(parentId));\n- 421 }\n- 422\n-427 GlobalCoordinate outerNormal(const LocalCoordinate &local) const\n- 428 {\n- 429 GlobalCoordinate normal;\n- 430\n- 431 if (codimensionWorld == 0)\n- 432 DUNE_THROW(Dune::Exception, \"There is no normal vector to a full-\n-dimensional intersection\");\n- 433 else if (codimensionWorld == 1) {\n- 434 /* TODO: Implement the general n-ary cross product here */\n- 435 const auto jacobianTransposed = geometry().jacobianTransposed(local);\n- 436 if (mydim==1) {\n- 437 normal[0] = - jacobianTransposed[0][1];\n- 438 normal[1] = jacobianTransposed[0][0];\n- 439 } else if (mydim==2) {\n- 440 normal[0] = (jacobianTransposed[0][1] * jacobianTransposed[1][2] -\n-jacobianTransposed[0][2] * jacobianTransposed[1][1]);\n- 441 normal[1] = - (jacobianTransposed[0][0] * jacobianTransposed[1][2] -\n-jacobianTransposed[0][2] * jacobianTransposed[1][0]);\n- 442 normal[2] = (jacobianTransposed[0][0] * jacobianTransposed[1][1] -\n-jacobianTransposed[0][1] * jacobianTransposed[1][0]);\n- 443 } else\n- 444 DUNE_THROW(Dune::NotImplemented, \"Remote intersections don't implement the\n-'outerNormal' method for \" << mydim << \"-dimensional intersections yet\");\n- 445 } else\n- 446 DUNE_THROW(Dune::NotImplemented, \"Remote intersections don't implement the\n-'outerNormal' method for intersections with codim >= 2 yet\");\n- 447\n- 448 return normal;\n- 449 }\n- 450\n-455 GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const\n- 456 {\n- 457 GlobalCoordinate normal = outerNormal(local);\n- 458 normal /= normal.two_norm();\n- 459 return normal;\n- 460 }\n- 461\n-466 GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const\n- 467 {\n- 468 return (unitOuterNormal(local) *= geometry().integrationElement(local));\n- 469 }\n- 470\n-475 GlobalCoordinate centerUnitOuterNormal () const\n- 476 {\n- 477 return unitOuterNormal(ReferenceElements::general(type\n-()).position(0,0));\n- 478 }\n- 479\n-483 Intersection flip() const\n- 484 {\n- 485 return Intersection(glue_,i_);\n- 486 }\n- 487\n- 488#ifdef QUICKHACK_INDEX\n- 489 typedef typename GridGlue::IndexType IndexType;\n- 490\n- 491 IndexType index() const\n- 492 {\n- 493 return i_->index_;\n- 494 }\n- 495\n- 496#endif\n- 497\n- 498 private:\n- 499\n- 500 friend class IntersectionIndexSet;\n- 501\n- 502 /* M E M B E R V A R I A B L E S */\n- 503\n- 505 const GridGlue* glue_;\n- 506\n- 508 const IntersectionData* i_;\n- 509 };\n- 510\n- 511\n- 512 } // end namespace GridGlue\n- 513} // end namespace Dune\n- 514\n- 515#endif // DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n-gridglue.hh\n-Central component of the module implementing the coupling of two grids.\n+ 293#if HAVE_MPI\n+ 298 template\n+299 struct CommPolicy< ::Dune::GridGlue::CommInfo >\n+ 300 {\n+304 typedef ::Dune::GridGlue::CommInfo Type;\n+ 305\n+309 typedef DataTypeImp IndexedType;\n+ 310\n+ 314 // typedef SizeOne IndexedTypeFlag;\n+315 typedef VariableSize IndexedTypeFlag;\n+ 316\n+320 static size_t getSize(const Type& commInfo, size_t i)\n+ 321 {\n+ 322 // get Intersection\n+ 323 typedef typename Type::GridGlue::Intersection Intersection;\n+ 324 Intersection ris(commInfo.gridglue->getIntersection(i));\n+ 325\n+ 326 // ask data handle for size\n+ 327 return commInfo.data->size(ris);\n+ 328 }\n+ 329 };\n+ 330#endif\n+ 331\n+ 332} // end namespace Dune\n+ 333#endif\n Dune\n Definition: gridglue.hh:37\n-Dune::GridGlue::ProjectionImplementation::inside\n-bool inside(const Coordinate &x, const Field &epsilon)\n-Definition: projection_impl.hh:111\n-Dune::GridGlue::GridGlue\n-sequential adapter to couple two grids at specified close together boundaries\n-Definition: gridglue.hh:67\n-Dune::GridGlue::GridGlue::IndexType\n-unsigned int IndexType\n-Definition: gridglue.hh:147\n-Dune::GridGlue::GridGlue::ctype\n-PromotionTraits< typenameGridView< 0 >::ctype, typenameGridView< 1 >::ctype >::\n-PromotedType ctype\n-The type used for coordinates.\n-Definition: gridglue.hh:171\n-Dune::GridGlue::GridGlue::dimworld\n-static constexpr int dimworld\n-export the world dimension This is the maximum of the extractors' world\n-dimensions.\n-Definition: gridglue.hh:166\n-Dune::GridGlue::IntersectionData\n-storage class for Dune::GridGlue::Intersection related data\n-Definition: intersection.hh:38\n-Dune::GridGlue::IntersectionData::IndexType\n-GridGlue::IndexType IndexType\n-Definition: intersection.hh:42\n-Dune::GridGlue::IntersectionData::GridGeometry\n-AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim,\n-GridGlue::template GridView< side >::dimensionworld > GridGeometry\n-Definition: intersection.hh:65\n-Dune::GridGlue::IntersectionData::localGeometry\n-const GridLocalGeometry< side > & localGeometry(unsigned int parentId=0) const\n-Definition: intersection.hh:85\n-Dune::GridGlue::IntersectionData::coorddim\n-static constexpr int coorddim\n-Dimension of the world space of the intersection.\n-Definition: intersection.hh:45\n-Dune::GridGlue::IntersectionData::Grid0Geometry\n-GridGeometry< 0 > Grid0Geometry\n-Definition: intersection.hh:67\n-Dune::GridGlue::IntersectionData::GridIndexType\n-typename GridGlue::template GridView< side >::IndexSet::IndexType GridIndexType\n-Definition: intersection.hh:71\n-Dune::GridGlue::IntersectionData::index\n-IndexType index(unsigned int parentId=0) const\n-Definition: intersection.hh:97\n-Dune::GridGlue::IntersectionData::parents\n-IndexType parents() const\n-Definition: intersection.hh:101\n-Dune::GridGlue::IntersectionData::Grid1LocalGeometry\n-GridLocalGeometry< 1 > Grid1LocalGeometry\n-Definition: intersection.hh:61\n-Dune::GridGlue::IntersectionData::geometry\n-const GridGeometry< side > & geometry() const\n-Definition: intersection.hh:89\n-Dune::GridGlue::IntersectionData::local\n-bool local() const\n-Definition: intersection.hh:93\n-Dune::GridGlue::IntersectionData::GridGlue\n-::Dune::GridGlue::GridGlue< P0, P1 > GridGlue\n-Definition: intersection.hh:40\n-Dune::GridGlue::IntersectionData::mydim\n-static constexpr int mydim\n-Dimension of the intersection.\n-Definition: intersection.hh:54\n-Dune::GridGlue::IntersectionData::GridLocalGeometry\n-AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim,\n-GridGlue::template GridView< side >::dimension > GridLocalGeometry\n-Definition: intersection.hh:58\n-Dune::GridGlue::IntersectionData::Grid0LocalGeometry\n-GridLocalGeometry< 0 > Grid0LocalGeometry\n-Definition: intersection.hh:60\n-Dune::GridGlue::IntersectionData::index_\n-IndexType index_\n-index of this intersection after GridGlue interface\n-Definition: intersection.hh:112\n-Dune::GridGlue::IntersectionData::Grid1IndexType\n-GridIndexType< 1 > Grid1IndexType\n-Definition: intersection.hh:74\n-Dune::GridGlue::IntersectionData::Grid1Geometry\n-GridGeometry< 1 > Grid1Geometry\n-Definition: intersection.hh:68\n-Dune::GridGlue::IntersectionData::Grid0IndexType\n-GridIndexType< 0 > Grid0IndexType\n-Definition: intersection.hh:73\n-Dune::GridGlue::IntersectionData::IntersectionData\n-IntersectionData()=default\n-Default Constructor.\n+Dune::GridGlue::BackwardOperator\n+CommunicationOperator< Dune::BackwardCommunication > BackwardOperator\n+Definition: gridgluecommunicate.hh:264\n+Dune::GridGlue::ForwardOperator\n+CommunicationOperator< Dune::ForwardCommunication > ForwardOperator\n+Definition: gridgluecommunicate.hh:263\n+Dune::GridGlue::RankPair\n+std::pair< int, int > RankPair\n+Definition: gridgluecommunicate.hh:24\n+Dune::GridGlue::operator<<\n+std::ostream & operator<<(std::ostream &os, const GlobalId &id)\n+Definition: gridgluecommunicate.hh:55\n Dune::GridGlue::Intersection\n The intersection of two entities of the two patches of a GridGlue.\n Definition: intersection.hh:261\n-Dune::GridGlue::Intersection::OutsideLocalGeometry\n-Traits::OutsideLocalGeometry OutsideLocalGeometry\n-Definition: intersection.hh:275\n Dune::GridGlue::Intersection::flip\n Intersection< P0, P1, O, I > flip() const\n Return a copy of the intersection with inside and outside switched.\n Definition: intersection.hh:483\n-Dune::GridGlue::Intersection::conforming\n-bool conforming() const\n-Return true if intersection is conforming.\n-Definition: intersection.hh:335\n-Dune::GridGlue::Intersection::InsideEntity\n-InsideGridView::Traits::template Codim< 0 >::Entity InsideEntity\n-Definition: intersection.hh:281\n-Dune::GridGlue::Intersection::GridGlue\n-Traits::GridGlue GridGlue\n-Definition: intersection.hh:267\n-Dune::GridGlue::Intersection::indexInOutside\n-int indexInOutside(unsigned int parentId=0) const\n-Local number of codim 1 entity in outside() Entity where intersection is\n-contained in.\n-Definition: intersection.hh:417\n-Dune::GridGlue::Intersection::indexInInside\n-int indexInInside(unsigned int parentId=0) const\n-Local number of codim 1 entity in the inside() Entity where intersection is\n-contained in.\n-Definition: intersection.hh:410\n Dune::GridGlue::Intersection::self\n bool self() const\n For parallel computations: Return true if inside() entity exists locally.\n Definition: intersection.hh:393\n-Dune::GridGlue::Intersection::mydim\n-static constexpr auto mydim\n-dimension of the intersection\n-Definition: intersection.hh:291\n-Dune::GridGlue::Intersection::outsidePatch\n-static constexpr int outsidePatch\n-outside patch\n-Definition: intersection.hh:297\n-Dune::GridGlue::Intersection::type\n-Dune::GeometryType type() const\n-Type of reference element for this intersection.\n-Definition: intersection.hh:377\n Dune::GridGlue::Intersection::inside\n InsideEntity inside(unsigned int parentId=0) const\n Return element on the inside of this intersection.\n Definition: intersection.hh:319\n-Dune::GridGlue::Intersection::LocalCoordinate\n-Traits::LocalCoordinate LocalCoordinate\n-Definition: intersection.hh:284\n-Dune::GridGlue::Intersection::geometryInOutside\n-const OutsideLocalGeometry & geometryInOutside(unsigned int parentId=0) const\n-Geometric information about this intersection in local coordinates of the\n-outside() element.\n-Definition: intersection.hh:349\n-Dune::GridGlue::Intersection::unitOuterNormal\n-GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const\n-Return a unit outer normal.\n-Definition: intersection.hh:455\n-Dune::GridGlue::Intersection::ctype\n-Traits::ctype ctype\n-Definition: intersection.hh:279\n-Dune::GridGlue::Intersection::Traits\n-IntersectionTraits< P0, P1, I, O > Traits\n-Definition: intersection.hh:265\n-Dune::GridGlue::Intersection::IntersectionData\n-Traits::IntersectionData IntersectionData\n-Definition: intersection.hh:268\n-Dune::GridGlue::Intersection::GlobalCoordinate\n-Traits::GlobalCoordinate GlobalCoordinate\n-Definition: intersection.hh:285\n-Dune::GridGlue::Intersection::OutsideGeometry\n-Traits::OutsideGeometry OutsideGeometry\n-Definition: intersection.hh:276\n Dune::GridGlue::Intersection::neighbor\n size_t neighbor(unsigned int g=0) const\n Return number of embeddings into local grid0 (grid1) entities.\n Definition: intersection.hh:399\n-Dune::GridGlue::Intersection::Intersection\n-Intersection(const GridGlue *glue, const IntersectionData *i)\n-Constructor for a given Dataset.\n-Definition: intersection.hh:311\n-Dune::GridGlue::Intersection::insidePatch\n-static constexpr int insidePatch\n-inside patch\n-Definition: intersection.hh:294\n-Dune::GridGlue::Intersection::Geometry\n-Traits::Geometry Geometry\n-Definition: intersection.hh:278\n-Dune::GridGlue::Intersection::OutsideEntity\n-OutsideGridView::Traits::template Codim< 0 >::Entity OutsideEntity\n-Definition: intersection.hh:282\n Dune::GridGlue::Intersection::outside\n OutsideEntity outside(unsigned int parentId=0) const\n Return element on the outside of this intersection.\n Definition: intersection.hh:328\n-Dune::GridGlue::Intersection::geometry\n-const Geometry & geometry() const\n-Geometric information about this intersection as part of the inside grid.\n-Definition: intersection.hh:360\n-Dune::GridGlue::Intersection::InsideLocalGeometry\n-Traits::InsideLocalGeometry InsideLocalGeometry\n-Definition: intersection.hh:272\n-Dune::GridGlue::Intersection::outerNormal\n-GlobalCoordinate outerNormal(const LocalCoordinate &local) const\n-Return an outer normal (length not necessarily 1)\n-Definition: intersection.hh:427\n-Dune::GridGlue::Intersection::geometryOutside\n-const OutsideGeometry & geometryOutside() const\n-Geometric information about this intersection as part of the outside grid.\n-Definition: intersection.hh:371\n-Dune::GridGlue::Intersection::OutsideGridView\n-Traits::OutsideGridView OutsideGridView\n-Definition: intersection.hh:274\n-Dune::GridGlue::Intersection::integrationOuterNormal\n-GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const\n-Return an outer normal with the length of the integration element.\n-Definition: intersection.hh:466\n-Dune::GridGlue::Intersection::InsideGridView\n-Traits::InsideGridView InsideGridView\n-Definition: intersection.hh:271\n-Dune::GridGlue::Intersection::geometryInInside\n-const InsideLocalGeometry & geometryInInside(unsigned int parentId=0) const\n-Geometric information about this intersection in local coordinates of the\n-inside() element.\n-Definition: intersection.hh:342\n-Dune::GridGlue::Intersection::centerUnitOuterNormal\n-GlobalCoordinate centerUnitOuterNormal() const\n-Unit outer normal at the center of the intersection.\n-Definition: intersection.hh:475\n-Dune::GridGlue::Intersection::coorddim\n-static constexpr auto coorddim\n-dimension of the world space of the intersection\n-Definition: intersection.hh:288\n-Dune::GridGlue::IntersectionTraits\n-Definition: intersection.hh:234\n-Dune::GridGlue::IntersectionTraits::insidePatch\n-static constexpr int insidePatch\n-Definition: intersection.hh:249\n-Dune::GridGlue::IntersectionTraits::LocalCoordinate\n-Dune::FieldVector< ctype, mydim > LocalCoordinate\n-Definition: intersection.hh:253\n-Dune::GridGlue::IntersectionTraits::coorddim\n-static constexpr auto coorddim\n-Definition: intersection.hh:247\n-Dune::GridGlue::IntersectionTraits::OutsideGeometry\n-typename IntersectionData::template GridGeometry< outside > OutsideGeometry\n-Definition: intersection.hh:245\n-Dune::GridGlue::IntersectionTraits::ctype\n-typename GridGlue::ctype ctype\n-Definition: intersection.hh:252\n-Dune::GridGlue::IntersectionTraits::OutsideGridView\n-typename GridGlue::template GridView< outside > OutsideGridView\n-Definition: intersection.hh:239\n-Dune::GridGlue::IntersectionTraits::GlobalCoordinate\n-Dune::FieldVector< ctype, coorddim > GlobalCoordinate\n-Definition: intersection.hh:254\n-Dune::GridGlue::IntersectionTraits::OutsideLocalGeometry\n-typename IntersectionData::template GridLocalGeometry< outside >\n-OutsideLocalGeometry\n-Definition: intersection.hh:242\n-Dune::GridGlue::IntersectionTraits::InsideGridView\n-typename GridGlue::template GridView< inside > InsideGridView\n-Definition: intersection.hh:238\n-Dune::GridGlue::IntersectionTraits::Geometry\n-typename IntersectionData::template GridGeometry< inside > Geometry\n-Definition: intersection.hh:244\n-Dune::GridGlue::IntersectionTraits::InsideLocalGeometry\n-typename IntersectionData::template GridLocalGeometry< inside >\n-InsideLocalGeometry\n-Definition: intersection.hh:241\n-Dune::GridGlue::IntersectionTraits::outsidePatch\n-static constexpr int outsidePatch\n-Definition: intersection.hh:250\n-Dune::GridGlue::IntersectionTraits::mydim\n-static constexpr auto mydim\n-Definition: intersection.hh:248\n+Dune::GridGlue::GlobalId\n+Definition: gridgluecommunicate.hh:26\n+Dune::GridGlue::GlobalId::GlobalId\n+GlobalId(int i)\n+Definition: gridgluecommunicate.hh:38\n+Dune::GridGlue::GlobalId::GlobalId\n+GlobalId()\n+Definition: gridgluecommunicate.hh:30\n+Dune::GridGlue::GlobalId::GlobalId\n+GlobalId(int i, int j, unsigned int n)\n+Definition: gridgluecommunicate.hh:48\n+Dune::GridGlue::CommDataHandle\n+describes the features of a data handle for communication in parallel runs\n+using the GridGlue::commun...\n+Definition: gridgluecommunicate.hh:77\n+Dune::GridGlue::CommDataHandle::size\n+size_t size(RISType &i) const\n+Definition: gridgluecommunicate.hh:92\n+Dune::GridGlue::CommDataHandle::scatter\n+void scatter(MessageBufferImp &buff, const EntityType &e, const RISType &i,\n+size_t n)\n+Definition: gridgluecommunicate.hh:118\n+Dune::GridGlue::CommDataHandle::gather\n+void gather(MessageBufferImp &buff, const EntityType &e, const RISType &i)\n+const\n+pack data from user to message buffer\n+Definition: gridgluecommunicate.hh:104\n+Dune::GridGlue::CommDataHandle::DataType\n+DataTypeImp DataType\n+data type of data to communicate\n+Definition: gridgluecommunicate.hh:80\n+Dune::GridGlue::CommDataHandle::CommDataHandle\n+CommDataHandle()\n+Definition: gridgluecommunicate.hh:84\n+Dune::GridGlue::StreamingMessageBuffer\n+Definition: gridgluecommunicate.hh:141\n+Dune::GridGlue::StreamingMessageBuffer::j\n+size_t j\n+Definition: gridgluecommunicate.hh:183\n+Dune::GridGlue::StreamingMessageBuffer::StreamingMessageBuffer\n+StreamingMessageBuffer(DT *p)\n+Definition: gridgluecommunicate.hh:146\n+Dune::GridGlue::StreamingMessageBuffer::a\n+DT * a\n+Definition: gridgluecommunicate.hh:181\n+Dune::GridGlue::StreamingMessageBuffer::counter\n+size_t counter() const\n+Definition: gridgluecommunicate.hh:169\n+Dune::GridGlue::StreamingMessageBuffer::write\n+void write(const Y &data)\n+Definition: gridgluecommunicate.hh:155\n+Dune::GridGlue::StreamingMessageBuffer::read\n+void read(Y &data) const\n+Definition: gridgluecommunicate.hh:163\n+Dune::GridGlue::StreamingMessageBuffer::value_type\n+DT value_type\n+Definition: gridgluecommunicate.hh:143\n+Dune::GridGlue::StreamingMessageBuffer::i\n+size_t i\n+Definition: gridgluecommunicate.hh:182\n+Dune::GridGlue::StreamingMessageBuffer::clear\n+void clear()\n+Definition: gridgluecommunicate.hh:171\n+Dune::GridGlue::CommunicationOperator\n+forward gather scatter to user defined CommInfo class\n+Definition: gridgluecommunicate.hh:194\n+Dune::GridGlue::CommunicationOperator::scatter\n+static void scatter(CommInfo &commInfo, const typename CommInfo::DataType &v,\n+std::size_t i, std::size_t j=0)\n+Definition: gridgluecommunicate.hh:227\n+Dune::GridGlue::CommunicationOperator::gather\n+static const CommInfo::DataType & gather(const CommInfo &commInfo, size_t i,\n+size_t j=0)\n+Definition: gridgluecommunicate.hh:197\n+Dune::GridGlue::CommInfo\n+collects all GridGlue data requried for communication\n+Definition: gridgluecommunicate.hh:272\n+Dune::GridGlue::CommInfo::dir\n+Dune::CommunicationDirection dir\n+Definition: gridgluecommunicate.hh:288\n+Dune::GridGlue::CommInfo::value_type\n+DataTypeImp value_type\n+Definition: gridgluecommunicate.hh:273\n+Dune::GridGlue::CommInfo::DataType\n+DataTypeImp DataType\n+Definition: gridgluecommunicate.hh:275\n+Dune::GridGlue::CommInfo::data\n+::Dune::GridGlue::CommDataHandle< DataHandleImp, DataTypeImp > * data\n+Definition: gridgluecommunicate.hh:282\n+Dune::GridGlue::CommInfo::GridGlue\n+GG GridGlue\n+Definition: gridgluecommunicate.hh:274\n+Dune::GridGlue::CommInfo::currentsize\n+size_t currentsize\n+Definition: gridgluecommunicate.hh:287\n+Dune::GridGlue::CommInfo::CommInfo\n+CommInfo()\n+Definition: gridgluecommunicate.hh:277\n+Dune::GridGlue::CommInfo::mbuffer\n+mutable ::Dune::GridGlue::StreamingMessageBuffer< DataType > mbuffer\n+Definition: gridgluecommunicate.hh:286\n+Dune::GridGlue::CommInfo::gridglue\n+const GridGlue * gridglue\n+Definition: gridgluecommunicate.hh:281\n+Dune::GridGlue::CommInfo::buffer\n+std::vector< DataType > buffer\n+Definition: gridgluecommunicate.hh:285\n+Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n+>::getSize\n+static size_t getSize(const Type &commInfo, size_t i)\n+Get the number of objects at an intersection.\n+Definition: gridgluecommunicate.hh:320\n+Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n+>::IndexedType\n+DataTypeImp IndexedType\n+The datatype that should be communicated.\n+Definition: gridgluecommunicate.hh:309\n+Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n+>::Type\n+::Dune::GridGlue::CommInfo< GG, DataHandleImp, DataTypeImp > Type\n+The type of the GridGlueCommInfo.\n+Definition: gridgluecommunicate.hh:304\n+Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n+>::IndexedTypeFlag\n+VariableSize IndexedTypeFlag\n+Each intersection can communicate a different number of objects.\n+Definition: gridgluecommunicate.hh:315\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00020.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00020.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: gridgluecommunicate.hh File Reference \n+dune-grid-glue: gridgluevtkwriter.hh File Reference \n \n \n \n \n \n \n \n@@ -64,74 +64,47 @@\n \n \n \n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -6,28 +6,22 @@\n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n Classes | Namespaces\n-gridglueamirawriter.hh File Reference\n-Write all remote intersections to a AmiraMesh file. More...\n-#include \n-#include \n-#include \n+intersectionindexset.hh File Reference\n+#include \n+#include \n Go_to_the_source_code_of_this_file.\n Classes\n-class \u00a0Dune::GridGlue::GridGlueAmiraWriter\n-\u00a0 Write remote intersections to a AmiraMesh file for debugging purposes.\n- More...\n+class \u00a0Dune::GridGlue::IntersectionIndexSet<_P0,_P1_>\n \u00a0\n Namespaces\n namespace \u00a0Dune\n \u00a0\n namespace \u00a0Dune::GridGlue\n \u00a0\n-***** Detailed Description *****\n-Write all remote intersections to a AmiraMesh file.\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00023_source.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00023_source.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: gridglueamirawriter.hh Source File \n+dune-grid-glue: intersectionindexset.hh Source File \n \n \n \n \n \n \n \n@@ -62,162 +62,81 @@\n \n \n \n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,28 +5,29 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces | Functions\n-rangegenerators.hh File Reference\n-#include \n+Classes | Namespaces\n+gridglueamirawriter.hh File Reference\n+Write all remote intersections to a AmiraMesh file. More...\n+#include \n+#include \n+#include \n Go_to_the_source_code_of_this_file.\n Classes\n-struct \u00a0Dune::GridGlue::Reverse<_reverse_>\n+class \u00a0Dune::GridGlue::GridGlueAmiraWriter\n+\u00a0 Write remote intersections to a AmiraMesh file for debugging purposes.\n+ More...\n \u00a0\n Namespaces\n namespace \u00a0Dune\n \u00a0\n namespace \u00a0Dune::GridGlue\n \u00a0\n- Functions\n-template<... >\n-IteratorRange<... >\u00a0Dune::GridGlue::intersections (const GridGlue<... > &glue,\n- const Reverse<... > &reverse=!reversed)\n-\u00a0 Iterate over all intersections of a GridGlue. More...\n-\u00a0\n+***** Detailed Description *****\n+Write all remote intersections to a AmiraMesh file.\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00026_source.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00026_source.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: rangegenerators.hh Source File \n+dune-grid-glue: gridglueamirawriter.hh Source File \n \n \n \n \n \n \n \n@@ -62,74 +62,162 @@\n \n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,89 +5,166 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-rangegenerators.hh\n+gridglueamirawriter.hh\n Go_to_the_documentation_of_this_file.\n 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 3#ifndef DUNE_GRIDGLUE_ADAPTER_RANGEGENERATORS_HH\n- 4#define DUNE_GRIDGLUE_ADAPTER_RANGEGENERATORS_HH\n- 5\n- 6#include \n- 7\n- 8namespace Dune {\n- 9namespace GridGlue {\n+ 8#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n+ 9#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n 10\n- 14template\n-15struct Reverse\n- 16 : std::integral_constant\n- 17{\n-18 typedef Reverse type;\n- 19\n- 20 constexpr\n-21 Reverse operator!() const\n- 22 { return {}; }\n- 23};\n- 24\n- 25#ifdef DOXYGEN\n- 26\n-31const Reverse reversed;\n+ 11#include \n+ 12#include \n+ 13#include \n+ 14\n+ 15namespace Dune {\n+ 16namespace GridGlue {\n+ 17\n+20class GridGlueAmiraWriter\n+ 21{\n+ 22\n+ 26 template \n+ 27 static void writeIntersections(const Glue& glue, const std::string&\n+filename)\n+ 28 {\n+ 29 static_assert((side==0 || side==1), \"'side' can only be 0 or 1\");\n+ 30\n+ 31 std::ofstream fgrid;\n 32\n- 66template<...>\n- 67IteratorRange<...>\n-68intersections(const GridGlue<...>& glue, const Reverse<...>& reverse =\n-!reversed);\n- 69\n- 70#else\n- 71\n- 72namespace {\n- 73const Reverse reversed = {};\n- 74} /* namespace */\n+ 33 fgrid.open(filename.c_str());\n+ 34\n+ 35 using GridView = typename Glue::template GridView;\n+ 36 const int dim = GridView::dimension;\n+ 37 const int domdimw = GridView::dimensionworld;\n+ 38\n+ 39 // coordinates have to be in R^3 in the VTK format\n+ 40 std::string coordinatePadding;\n+ 41 for (int i=domdimw; i<3; i++)\n+ 42 coordinatePadding += \" 0\";\n+ 43\n+ 44 int overlaps = glue.size();\n+ 45\n+ 46 if (dim==3) {\n+ 47\n+ 48 fgrid << \"# HyperSurface 0.1 ASCII \\n\" << std::endl;\n+ 49 fgrid<<\"\\n\";\n+ 50 fgrid<<\"Parameters {\\n\";\n+ 51 fgrid<<\" Materials {\\n\";\n+ 52 fgrid<<\" outside {\\n\";\n+ 53 fgrid<<\" Id 0\\n\";\n+ 54 fgrid<<\" }\\n\";\n+ 55 fgrid<<\" inside {\\n\";\n+ 56 fgrid<<\" Id 1\\n\";\n+ 57 fgrid<<\" }\\n\";\n+ 58 fgrid<<\" }\\n\";\n+ 59 fgrid<<\"\\n\";\n+ 60 fgrid<<\"}\\n\";\n+ 61\n+ 62 // ////////////////////////////////////////////\n+ 63 // Write vertices\n+ 64 // ////////////////////////////////////////////\n+ 65\n+ 66 //use dim and not dim+1\n+ 67 fgrid<<\"\\nVertices \"<< overlaps*(dim)<<\"\\n\";\n+ 68 auto isEnd = glue.template iend();\n+ 69 for (auto isIt = glue.template ibegin(); isIt != isEnd; ++isIt)\n+ 70 {\n+ 71 const auto& geometry = isIt->geometry();\n+ 72 for (int i = 0; i < geometry.corners(); ++i)\n+ 73 fgrid << geometry.corner(i) << coordinatePadding << std::endl;\n+ 74 }\n 75\n- 76template\n- 77IteratorRange< typename GridGlue::template\n-IntersectionIterator >\n- 78intersections(const GridGlue& glue, const Reverse& = {})\n- 79{\n- 80 const static int side = reverse ? 1 : 0;\n- 81 return {glue.template ibegin(), glue.template iend()};\n- 82}\n- 83\n- 84#endif // DOXYGEN\n- 85\n- 86} /* namespace GridGlue */\n- 87} /* namespace Dune */\n- 88\n- 89#endif\n+ 76 // ////////////////////////////////////////////\n+ 77 // Write triangles\n+ 78 // ////////////////////////////////////////////\n+ 79\n+ 80 fgrid<<\"NBranchingPoints 0\\n\";\n+ 81 fgrid<<\"NVerticesOnCurves 0\\n\";\n+ 82 fgrid<<\"BoundaryCurves 0\\n\";\n+ 83 fgrid<<\"Patches 1\\n\";\n+ 84 fgrid<<\"{\\n\";\n+ 85 fgrid<<\"InnerRegion inside\\n\";\n+ 86 fgrid<<\"OuterRegion outside\\n\";\n+ 87 fgrid<<\"BoundaryID 0\\n\";\n+ 88 fgrid<<\"BranchingPoints 0\";\n+ 89 fgrid<<\"\\n\";\n+ 90\n+ 91 fgrid<<\"Triangles \"<();\n+ 119 for (auto isIt = glue.template ibegin(); isIt != isEnd; ++isIt) {\n+ 120 const auto& geometry = isIt->geometry();\n+ 121 for (int i = 0; i <2; ++i)\n+ 122 fgrid << geometry.corner(i) <<\" \"<<0<<\"\\n\";\n+ 123 }\n+ 124 }\n+ 125\n+ 126 fgrid.close();\n+ 127 }\n+ 128\n+ 129public:\n+ 130 template\n+131 static void write(const Glue& glue, const std::string& path, int\n+appendix=1)\n+ 132 {\n+ 133 std::ostringstream name0;\n+ 134 name0 << path;\n+ 135 name0 << \"/domain.surf\" << std::setw(3) << std::setfill('0') << appendix;\n+ 136\n+ 137 // Write extracted grid and remote intersection on the grid1-side\n+ 138 writeIntersections(glue,name0.str());\n+ 139\n+ 140 std::ostringstream name1;\n+ 141 name1 << path;\n+ 142 name1 << \"/target.surf\" << std::setw(3) << std::setfill('0') << appendix;\n+ 143\n+ 144 writeIntersections(glue, name1.str());\n+ 145 }\n+ 146\n+ 147};\n+ 148\n+ 149} // namespace GridGlue\n+ 150} // namespace Dune\n+ 151\n+ 152#endif // DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n Dune\n Definition: gridglue.hh:37\n-Dune::GridGlue::intersections\n-IteratorRange<... > intersections(const GridGlue<... > &glue, const Reverse<...\n-> &reverse=!reversed)\n-Iterate over all intersections of a GridGlue.\n-Dune::GridGlue::GridGlue\n-sequential adapter to couple two grids at specified close together boundaries\n-Definition: gridglue.hh:67\n-Dune::GridGlue::GridGlue::intersections\n-IteratorRange<... > intersections(const GridGlue<... > &glue, const Reverse<...\n-> &reverse=!reversed)\n-Iterate over all intersections of a GridGlue.\n-Dune::GridGlue::Reverse\n-Definition: rangegenerators.hh:17\n-Dune::GridGlue::Reverse::reversed\n-const Reverse< true > reversed\n-Definition: rangegenerators.hh:31\n-Dune::GridGlue::Reverse::type\n-Reverse type\n-Definition: rangegenerators.hh:18\n-Dune::GridGlue::Reverse::operator!\n-constexpr Reverse operator!() const\n-Definition: rangegenerators.hh:21\n+Dune::GridGlue::GridGlueAmiraWriter\n+Write remote intersections to a AmiraMesh file for debugging purposes.\n+Definition: gridglueamirawriter.hh:21\n+Dune::GridGlue::GridGlueAmiraWriter::write\n+static void write(const Glue &glue, const std::string &path, int appendix=1)\n+Definition: gridglueamirawriter.hh:131\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00029.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00029.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: areawriter_impl.hh File Reference \n+dune-grid-glue: areawriter.hh File Reference \n \n \n \n \n \n \n \n@@ -63,47 +63,34 @@\n \n \n
Write all remote intersections to a vtk file for debugging. \n-More...
\n-#include <fstream>\n-
#include <iomanip>\n-
#include <type_traits>\n+
#include "intersection.hh"\n
#include <vector>\n-
#include <dune/common/classname.hh>\n-
#include <dune/geometry/type.hh>\n-
#include <dune/geometry/referenceelements.hh>\n-
#include <dune/grid-glue/gridglue.hh>\n-
Go to the source code of this file.
\n-\n-Classes | |
| class | Dune::GridGlue::GridGlueVtkWriter |
| Write remote intersections to a vtk file for debugging purposes. More... | |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
\n+Functions | |
| template<typename T > | |
| void | Dune::GridGlue::printVector (const std::vector< T > &v, std::string name, int rank) |
Detailed Description
\n-Write all remote intersections to a vtk file for debugging.
\n-\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,34 +5,28 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces\n-gridgluevtkwriter.hh File Reference\n-Write all remote intersections to a vtk file for debugging. More...\n-#include intersection.hh File Reference
\n-
\n-
\n+\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n\n \n
\n
\n+
\n+\n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,23 +5,47 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces\n-intersectionindexset.hh File Reference\n+Classes | Namespaces | Macros\n+intersection.hh File Reference\n+Model of the Intersection concept provided by GridGlue. More...\n+#include \n+#include \n+#include \n+#include \n+#include \n+#include \n+#include \n #include \n-#include \n Go_to_the_source_code_of_this_file.\n Classes\n-class \u00a0Dune::GridGlue::IntersectionIndexSet<_P0,_P1_>\n+ class \u00a0Dune::GridGlue::IntersectionData<_P0,_P1_>\n+\u00a0 storage class for Dune::GridGlue::Intersection related data More...\n+\u00a0\n+struct \u00a0Dune::GridGlue::IntersectionTraits<_P0,_P1,_inside,_outside_>\n+\u00a0\n+ class \u00a0Dune::GridGlue::Intersection<_P0,_P1,_I,_O_>\n+\u00a0 The intersection of two entities of the two patches of a GridGlue.\n+ More...\n \u00a0\n Namespaces\n namespace \u00a0Dune\n \u00a0\n namespace \u00a0Dune::GridGlue\n \u00a0\n+ Macros\n+#define\u00a0ONLY_SIMPLEX_INTERSECTIONS\n+\u00a0\n+***** Detailed Description *****\n+Model of the Intersection concept provided by GridGlue.\n+ Author\n+ Christian Engwer\n+***** Macro Definition Documentation *****\n+***** \u25c6\u00a0ONLY_SIMPLEX_INTERSECTIONS *****\n+#define ONLY_SIMPLEX_INTERSECTIONS\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00014_source.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00014_source.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: intersectionindexset.hh Source File \n+dune-grid-glue: intersection.hh Source File \n \n \n \n \n \n \n \n@@ -62,81 +62,525 @@\n \n \n
\n #include <dune/grid-glue/gridglue.hh>\n-
#include <dune/grid-glue/adapter/intersection.hh>\n+\n+
Model of the Intersection concept provided by GridGlue. \n+More...
\n+#include <algorithm>\n+
#include <optional>\n+
#include <tuple>\n+
#include <dune/common/deprecated.hh>\n+
#include <dune/common/version.hh>\n+
#include <dune/geometry/affinegeometry.hh>\n+
#include <dune/geometry/referenceelements.hh>\n+
#include <dune/grid-glue/gridglue.hh>\n
Go to the source code of this file.
\n\n Classes | |
| class | Dune::GridGlue::IntersectionIndexSet< P0, P1 > |
| class | Dune::GridGlue::IntersectionData< P0, P1 > |
| storage class for Dune::GridGlue::Intersection related data More... | |
| struct | Dune::GridGlue::IntersectionTraits< P0, P1, inside, outside > |
| class | Dune::GridGlue::Intersection< P0, P1, I, O > |
| The intersection of two entities of the two patches of a GridGlue. More... | |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
\n+Macros | |
| #define | ONLY_SIMPLEX_INTERSECTIONS |
Detailed Description
\n+Model of the Intersection concept provided by GridGlue.
\n+- Author
- Christian Engwer
Macro Definition Documentation
\n+\n+◆ ONLY_SIMPLEX_INTERSECTIONS
\n+\n+\n+
\n \n+ \n+
\n+
| #define ONLY_SIMPLEX_INTERSECTIONS | \n+
\n+\n+
\n+\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,23 +5,47 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces\n-intersectionindexset.hh File Reference\n+Classes | Namespaces | Macros\n+intersection.hh File Reference\n+Model of the Intersection concept provided by GridGlue. More...\n+#include \n- \n+ \n
\n intersectionindexset.hh
intersection.hh
\n-Go to the documentation of this file.\n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,106 +5,758 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-intersectionindexset.hh\n+intersection.hh\n Go_to_the_documentation_of_this_file.\n-1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n+ 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-\n+ 2// vi: set et ts=4 sw=2 sts=2:\n+ 3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n- 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n+ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 3#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n- 4#define DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n- 5\n- 6#include \n- 7#include \n- 8\n- 9#ifndef ONLY_SIMPLEX_INTERSECTIONS\n- 10// we currently support only one intersection type. If we want to support\n-more,\n- 11// we have to think about the semantics of our IndexSet\n- 12#error Not Implemented\n- 13#endif\n- 14\n- 15namespace Dune {\n- 16 namespace GridGlue {\n+ 11#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n+ 12#define DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n+ 13\n+ 14#include \n+ 15#include \n+ 16#include \n 17\n- 18 template\n-19 class IntersectionIndexSet\n- 20 {\n- 21 friend class ::Dune::GridGlue::GridGlue;\n- 22 typedef ::Dune::GridGlue::GridGlue GridGlue;\n+ 18#include \n+ 19#include \n+ 20#include \n+ 21#include \n+ 22#include \n 23\n- 24 public:\n+24#define ONLY_SIMPLEX_INTERSECTIONS\n 25\n-27 typedef typename GridGlue::IndexType IndexType;\n-29 typedef size_t SizeType;\n- 30\n- 33 template\n-34 IndexType index (const Intersection & i) const\n- 35 {\n- 36 return i.i_->index_;\n- 37 }\n- 38\n-41 SizeType size () const\n- 42 {\n- 43 return glue_->size();\n- 44 }\n- 45\n- 46 private:\n- 47\n- 49 IntersectionIndexSet(const GridGlue * g) :\n- 50 glue_(g) {}\n+ 26namespace Dune {\n+ 27 namespace GridGlue {\n+ 28\n+ 29 // forward declaration\n+ 30 template\n+ 31 class IntersectionIndexSet;\n+ 32\n+ 36 template\n+37 class IntersectionData\n+ 38 {\n+ 39 public:\n+40 typedef ::Dune::GridGlue::GridGlue GridGlue;\n+ 41\n+42 typedef typename GridGlue::IndexType IndexType;\n+ 43\n+45 static constexpr int coorddim = GridGlue::dimworld;\n+ 46\n+ 47 private:\n+ 48 // intermediate quantities\n+ 49 template\n+ 50 static constexpr int dim() { return GridGlue::template GridView::\n+Grid::dimension - GridGlue::template GridPatch::codim; }\n 51\n- 52 const GridGlue * glue_;\n- 53 };\n- 54\n- 55 } // end namespace GridGlue\n- 56} // end namespace Dune\n- 57\n- 58#endif // DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n+ 52 public:\n+54 static constexpr int mydim = dim<0>() < dim<1>() ? dim<0>() : dim<1>();\n+ 55\n+ 56 template\n+57 using GridLocalGeometry = AffineGeometry<\n+ 58 typename GridGlue::template GridView::ctype, mydim, GridGlue::\n+template GridView::dimension>;\n+ 59\n+60 using Grid0LocalGeometry [[deprecated(\"please use GridLocalGeometry<0>\n+instead\")]] = GridLocalGeometry<0>;\n+61 using Grid1LocalGeometry [[deprecated(\"please use GridLocalGeometry<1>\n+instead\")]] = GridLocalGeometry<1>;\n+ 62\n+ 63 template\n+64 using GridGeometry = AffineGeometry<\n+ 65 typename GridGlue::template GridView::ctype, mydim, GridGlue::\n+template GridView::dimensionworld>;\n+ 66\n+67 using Grid0Geometry [[deprecated(\"please use GridGeometry<0> instead\")]] =\n+GridGeometry<0>;\n+68 using Grid1Geometry [[deprecated(\"please use GridGeometry<1> instead\")]] =\n+GridGeometry<1>;\n+ 69\n+ 70 template\n+71 using GridIndexType = typename GridGlue::template GridView::IndexSet::\n+IndexType;\n+ 72\n+73 using Grid0IndexType [[deprecated(\"please use GridIndexType<0> instead\")]] =\n+GridIndexType<0>;\n+74 using Grid1IndexType [[deprecated(\"please use GridIndexType<1> instead\")]] =\n+GridIndexType<1>;\n+ 75\n+ 77 IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned\n+int offset, bool grid0local, bool grid1local);\n+ 78\n+80 IntersectionData() = default;\n+ 81\n+ 82 /* Accessor Functions */\n+ 83\n+ 84 template\n+85 const GridLocalGeometry& localGeometry(unsigned int parentId = 0)\n+const\n+ 86 { return *std::get(sideData_).gridlocalgeom[parentId]; }\n+ 87\n+ 88 template\n+89 const GridGeometry& geometry() const\n+ 90 { return *std::get(sideData_).gridgeom; }\n+ 91\n+ 92 template\n+93 bool local() const\n+ 94 { return std::get(sideData_).gridlocal; }\n+ 95\n+ 96 template\n+97 IndexType index(unsigned int parentId = 0) const\n+ 98 { return std::get(sideData_).gridindices[parentId]; }\n+ 99\n+ 100 template\n+101 IndexType parents() const\n+ 102 { return std::get(sideData_).gridindices.size(); }\n+ 103\n+ 104 private:\n+ 105 template\n+ 106 void initializeGeometry(const GridGlue& glue, unsigned mergeindex);\n+ 107\n+ 108 /* M E M B E R V A R I A B L E S */\n+ 109\n+ 110 public:\n+112 IndexType index_;\n+ 113\n+ 114 private:\n+ 115 template\n+ 116 struct SideData {\n+ 118 bool gridlocal = false;\n+ 119\n+ 121 std::vector< GridIndexType > gridindices;\n+ 122\n+ 124 std::vector< std::optional< GridLocalGeometry > > gridlocalgeom;\n+ 125\n+ 133 std::optional< GridGeometry > gridgeom;\n+ 134 };\n+ 135\n+ 136 std::tuple< SideData<0>, SideData<1> > sideData_;\n+ 137 };\n+ 138\n+ 139 template\n+ 140 template\n+ 141 void IntersectionData::initializeGeometry(const GridGlue& glue,\n+unsigned mergeindex)\n+ 142 {\n+ 143 auto& data = std::get(sideData_);\n+ 144\n+ 145 const unsigned n_parents = glue.merger_->template parents\n+(mergeindex);\n+ 146\n+ 147 // init containers\n+ 148 data.gridindices.resize(n_parents);\n+ 149 data.gridlocalgeom.resize(n_parents);\n+ 150\n+ 151 // default values\n+ 152 data.gridindices[0] = 0;\n+ 153\n+ 154 static constexpr int nSimplexCorners = mydim + 1;\n+ 155 using ctype = typename GridGlue::ctype;\n+ 156\n+ 157 // initialize the local and the global geometries of grid `side`\n+ 158\n+ 159 // compute the coordinates of the subface's corners in codim 0 entity\n+local coordinates\n+ 160 static constexpr int elementdim = GridGlue::template GridView::\n+template Codim<0>::Geometry::mydimension;\n+ 161\n+ 162 // coordinates within the subentity that contains the remote intersection\n+ 163 std::array() >, nSimplexCorners>\n+corners_subEntity_local;\n+ 164\n+ 165 for (unsigned int par = 0; par < n_parents; ++par) {\n+ 166 for (int i = 0; i < nSimplexCorners; ++i)\n+ 167 corners_subEntity_local[i] = glue.merger_->template parentLocal\n+(mergeindex, i, par);\n+ 168\n+ 169 // Coordinates of the remote intersection corners wrt the element\n+coordinate system\n+ 170 std::array, nSimplexCorners>\n+corners_element_local;\n+ 171\n+ 172 if (data.gridlocal) {\n+ 173 data.gridindices[par] = glue.merger_->template parent\n+(mergeindex,par);\n+ 174\n+ 175 typename GridGlue::template GridPatch::LocalGeometry\n+ 176 gridLocalGeometry = glue.template patch().geometryLocal\n+(data.gridindices[par]);\n+ 177 for (std::size_t i=0; i::Geometry\n+ 196 gridWorldGeometry = glue.template patch().geometry(data.gridindices\n+[par]);\n+ 197\n+ 198 // world coordinates of the remote intersection corners\n+ 199 std::array::\n+dimensionworld>, nSimplexCorners> corners_global;\n+ 200\n+ 201 for (std::size_t i=0; i\n+213 IntersectionData::IntersectionData(const GridGlue& glue, unsigned\n+int mergeindex, unsigned int offset,\n+ 214 bool grid0local, bool grid1local)\n+ 215 : index_(mergeindex+offset)\n+ 216 {\n+ 217 // if an invalid index is given do not proceed!\n+ 218 // (happens when the parent GridGlue initializes the \"end\"-Intersection)\n+ 219 assert (0 <= mergeindex || mergeindex < glue.index__sz);\n+ 220\n+ 221 std::get<0>(sideData_).gridlocal = grid0local;\n+ 222 std::get<1>(sideData_).gridlocal = grid1local;\n+ 223\n+ 224 initializeGeometry<0>(glue, mergeindex);\n+ 225 initializeGeometry<1>(glue, mergeindex);\n+ 226 }\n+ 227\n+ 232 template\n+233 struct IntersectionTraits\n+ 234 {\n+235 using GridGlue = ::Dune::GridGlue::GridGlue;\n+236 using IntersectionData = Dune::GridGlue::IntersectionData;\n+ 237\n+238 using InsideGridView = typename GridGlue::template GridView;\n+239 using OutsideGridView = typename GridGlue::template GridView;\n+ 240\n+241 using InsideLocalGeometry = typename IntersectionData::template\n+GridLocalGeometry;\n+242 using OutsideLocalGeometry = typename IntersectionData::template\n+GridLocalGeometry;\n+ 243\n+244 using Geometry = typename IntersectionData::template GridGeometry;\n+245 using OutsideGeometry = typename IntersectionData::template\n+GridGeometry;\n+ 246\n+247 static constexpr auto coorddim = IntersectionData::coorddim;\n+248 static constexpr auto mydim = IntersectionData::mydim;\n+249 static constexpr int insidePatch = inside;\n+250 static constexpr int outsidePatch = outside;\n+ 251\n+252 using ctype = typename GridGlue::ctype;\n+253 using LocalCoordinate = Dune::FieldVector;\n+254 using GlobalCoordinate = Dune::FieldVector;\n+ 255 };\n+ 256\n+ 259 template\n+260 class Intersection\n+ 261 {\n+ 262\n+ 263 public:\n+ 264\n+265 typedef IntersectionTraits Traits;\n+ 266\n+267 typedef typename Traits::GridGlue GridGlue;\n+268 typedef typename Traits::IntersectionData IntersectionData;\n+ 269\n+ 270\n+271 typedef typename Traits::InsideGridView InsideGridView;\n+272 typedef typename Traits::InsideLocalGeometry InsideLocalGeometry;\n+ 273\n+274 typedef typename Traits::OutsideGridView OutsideGridView;\n+275 typedef typename Traits::OutsideLocalGeometry OutsideLocalGeometry;\n+276 typedef typename Traits::OutsideGeometry OutsideGeometry;\n+ 277\n+278 typedef typename Traits::Geometry Geometry;\n+279 typedef typename Traits::ctype ctype;\n+ 280\n+281 typedef typename InsideGridView::Traits::template Codim<0>::Entity\n+InsideEntity;\n+282 typedef typename OutsideGridView::Traits::template Codim<0>::Entity\n+OutsideEntity;\n+ 283\n+284 typedef typename Traits::LocalCoordinate LocalCoordinate;\n+285 typedef typename Traits::GlobalCoordinate GlobalCoordinate;\n+ 286\n+288 static constexpr auto coorddim = Traits::coorddim;\n+ 289\n+291 static constexpr auto mydim = Traits::mydim;\n+ 292\n+294 static constexpr int insidePatch = Traits::insidePatch;\n+ 295\n+297 static constexpr int outsidePatch = Traits::outsidePatch;\n+ 298\n+ 299 // typedef unsigned int IndexType;\n+ 300\n+ 301 private:\n+ 305 static constexpr int codimensionWorld = coorddim - mydim;\n+ 306\n+ 307 public:\n+ 308 /* C O N S T R U C T O R S */\n+ 309\n+311 Intersection(const GridGlue* glue, const IntersectionData* i) :\n+ 312 glue_(glue), i_(i) {}\n+ 313\n+ 314 /* F U N C T I O N A L I T Y */\n+ 315\n+ 318 InsideEntity\n+319 inside(unsigned int parentId = 0) const\n+ 320 {\n+ 321 assert(self());\n+ 322 return glue_->template patch().element(i_->template index\n+(parentId));\n+ 323 }\n+ 324\n+ 327 OutsideEntity\n+328 outside(unsigned int parentId = 0) const\n+ 329 {\n+ 330 assert(neighbor());\n+ 331 return glue_->template patch().element(i_->template index\n+(parentId));\n+ 332 }\n+ 333\n+335 bool conforming() const\n+ 336 {\n+ 337 throw Dune::NotImplemented();\n+ 338 }\n+ 339\n+342 const InsideLocalGeometry& geometryInInside(unsigned int parentId = 0)\n+const\n+ 343 {\n+ 344 return i_->template localGeometry(parentId);\n+ 345 }\n+ 346\n+349 const OutsideLocalGeometry& geometryInOutside(unsigned int parentId = 0)\n+const\n+ 350 {\n+ 351 return i_->template localGeometry(parentId);\n+ 352 }\n+ 353\n+360 const Geometry& geometry() const\n+ 361 {\n+ 362 return i_->template geometry();\n+ 363 }\n+ 364\n+371 const OutsideGeometry& geometryOutside() const // DUNE_DEPRECATED\n+ 372 {\n+ 373 return i_->template geometry();\n+ 374 }\n+ 375\n+377 Dune::GeometryType type() const\n+ 378 {\n+ 379 #ifdef ONLY_SIMPLEX_INTERSECTIONS\n+ 380 # if DUNE_VERSION_NEWER(DUNE_GEOMETRY, 2, 6)\n+ 381 return Dune::GeometryTypes::simplex(mydim);\n+ 382 # else\n+ 383 static const Dune::GeometryType type(Dune::GeometryType::simplex, mydim);\n+ 384 return type;\n+ 385 # endif\n+ 386 #else\n+ 387 #error Not Implemented\n+ 388 #endif\n+ 389 }\n+ 390\n+ 391\n+393 bool self() const\n+ 394 {\n+ 395 return i_->template local();\n+ 396 }\n+ 397\n+399 size_t neighbor(unsigned int g = 0) const\n+ 400 {\n+ 401 if (g == 0 && i_->template local()) {\n+ 402 return i_->template parents();\n+ 403 } else if (g == 1 && i_->template local()) {\n+ 404 return i_->template parents();\n+ 405 }\n+ 406 return 0;\n+ 407 }\n+ 408\n+410 int indexInInside(unsigned int parentId = 0) const\n+ 411 {\n+ 412 assert(self());\n+ 413 return glue_->template patch().indexInInside(i_->template index\n+(parentId));\n+ 414 }\n+ 415\n+417 int indexInOutside(unsigned int parentId = 0) const\n+ 418 {\n+ 419 assert(neighbor());\n+ 420 return glue_->template patch().indexInInside(i_->template index\n+(parentId));\n+ 421 }\n+ 422\n+427 GlobalCoordinate outerNormal(const LocalCoordinate &local) const\n+ 428 {\n+ 429 GlobalCoordinate normal;\n+ 430\n+ 431 if (codimensionWorld == 0)\n+ 432 DUNE_THROW(Dune::Exception, \"There is no normal vector to a full-\n+dimensional intersection\");\n+ 433 else if (codimensionWorld == 1) {\n+ 434 /* TODO: Implement the general n-ary cross product here */\n+ 435 const auto jacobianTransposed = geometry().jacobianTransposed(local);\n+ 436 if (mydim==1) {\n+ 437 normal[0] = - jacobianTransposed[0][1];\n+ 438 normal[1] = jacobianTransposed[0][0];\n+ 439 } else if (mydim==2) {\n+ 440 normal[0] = (jacobianTransposed[0][1] * jacobianTransposed[1][2] -\n+jacobianTransposed[0][2] * jacobianTransposed[1][1]);\n+ 441 normal[1] = - (jacobianTransposed[0][0] * jacobianTransposed[1][2] -\n+jacobianTransposed[0][2] * jacobianTransposed[1][0]);\n+ 442 normal[2] = (jacobianTransposed[0][0] * jacobianTransposed[1][1] -\n+jacobianTransposed[0][1] * jacobianTransposed[1][0]);\n+ 443 } else\n+ 444 DUNE_THROW(Dune::NotImplemented, \"Remote intersections don't implement the\n+'outerNormal' method for \" << mydim << \"-dimensional intersections yet\");\n+ 445 } else\n+ 446 DUNE_THROW(Dune::NotImplemented, \"Remote intersections don't implement the\n+'outerNormal' method for intersections with codim >= 2 yet\");\n+ 447\n+ 448 return normal;\n+ 449 }\n+ 450\n+455 GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const\n+ 456 {\n+ 457 GlobalCoordinate normal = outerNormal(local);\n+ 458 normal /= normal.two_norm();\n+ 459 return normal;\n+ 460 }\n+ 461\n+466 GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const\n+ 467 {\n+ 468 return (unitOuterNormal(local) *= geometry().integrationElement(local));\n+ 469 }\n+ 470\n+475 GlobalCoordinate centerUnitOuterNormal () const\n+ 476 {\n+ 477 return unitOuterNormal(ReferenceElements::general(type\n+()).position(0,0));\n+ 478 }\n+ 479\n+483 Intersection flip() const\n+ 484 {\n+ 485 return Intersection(glue_,i_);\n+ 486 }\n+ 487\n+ 488#ifdef QUICKHACK_INDEX\n+ 489 typedef typename GridGlue::IndexType IndexType;\n+ 490\n+ 491 IndexType index() const\n+ 492 {\n+ 493 return i_->index_;\n+ 494 }\n+ 495\n+ 496#endif\n+ 497\n+ 498 private:\n+ 499\n+ 500 friend class IntersectionIndexSet;\n+ 501\n+ 502 /* M E M B E R V A R I A B L E S */\n+ 503\n+ 505 const GridGlue* glue_;\n+ 506\n+ 508 const IntersectionData* i_;\n+ 509 };\n+ 510\n+ 511\n+ 512 } // end namespace GridGlue\n+ 513} // end namespace Dune\n+ 514\n+ 515#endif // DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n gridglue.hh\n Central component of the module implementing the coupling of two grids.\n-intersection.hh\n-Model of the Intersection concept provided by GridGlue.\n Dune\n Definition: gridglue.hh:37\n+Dune::GridGlue::ProjectionImplementation::inside\n+bool inside(const Coordinate &x, const Field &epsilon)\n+Definition: projection_impl.hh:111\n Dune::GridGlue::GridGlue\n sequential adapter to couple two grids at specified close together boundaries\n Definition: gridglue.hh:67\n Dune::GridGlue::GridGlue::IndexType\n unsigned int IndexType\n Definition: gridglue.hh:147\n-Dune::GridGlue::GridGlue::size\n-size_t size() const\n-Definition: gridglue.hh:393\n+Dune::GridGlue::GridGlue::ctype\n+PromotionTraits< typenameGridView< 0 >::ctype, typenameGridView< 1 >::ctype >::\n+PromotedType ctype\n+The type used for coordinates.\n+Definition: gridglue.hh:171\n+Dune::GridGlue::GridGlue::dimworld\n+static constexpr int dimworld\n+export the world dimension This is the maximum of the extractors' world\n+dimensions.\n+Definition: gridglue.hh:166\n+Dune::GridGlue::IntersectionData\n+storage class for Dune::GridGlue::Intersection related data\n+Definition: intersection.hh:38\n+Dune::GridGlue::IntersectionData::IndexType\n+GridGlue::IndexType IndexType\n+Definition: intersection.hh:42\n+Dune::GridGlue::IntersectionData::GridGeometry\n+AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim,\n+GridGlue::template GridView< side >::dimensionworld > GridGeometry\n+Definition: intersection.hh:65\n+Dune::GridGlue::IntersectionData::localGeometry\n+const GridLocalGeometry< side > & localGeometry(unsigned int parentId=0) const\n+Definition: intersection.hh:85\n+Dune::GridGlue::IntersectionData::coorddim\n+static constexpr int coorddim\n+Dimension of the world space of the intersection.\n+Definition: intersection.hh:45\n+Dune::GridGlue::IntersectionData::Grid0Geometry\n+GridGeometry< 0 > Grid0Geometry\n+Definition: intersection.hh:67\n+Dune::GridGlue::IntersectionData::GridIndexType\n+typename GridGlue::template GridView< side >::IndexSet::IndexType GridIndexType\n+Definition: intersection.hh:71\n+Dune::GridGlue::IntersectionData::index\n+IndexType index(unsigned int parentId=0) const\n+Definition: intersection.hh:97\n+Dune::GridGlue::IntersectionData::parents\n+IndexType parents() const\n+Definition: intersection.hh:101\n+Dune::GridGlue::IntersectionData::Grid1LocalGeometry\n+GridLocalGeometry< 1 > Grid1LocalGeometry\n+Definition: intersection.hh:61\n+Dune::GridGlue::IntersectionData::geometry\n+const GridGeometry< side > & geometry() const\n+Definition: intersection.hh:89\n+Dune::GridGlue::IntersectionData::local\n+bool local() const\n+Definition: intersection.hh:93\n+Dune::GridGlue::IntersectionData::GridGlue\n+::Dune::GridGlue::GridGlue< P0, P1 > GridGlue\n+Definition: intersection.hh:40\n+Dune::GridGlue::IntersectionData::mydim\n+static constexpr int mydim\n+Dimension of the intersection.\n+Definition: intersection.hh:54\n+Dune::GridGlue::IntersectionData::GridLocalGeometry\n+AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim,\n+GridGlue::template GridView< side >::dimension > GridLocalGeometry\n+Definition: intersection.hh:58\n+Dune::GridGlue::IntersectionData::Grid0LocalGeometry\n+GridLocalGeometry< 0 > Grid0LocalGeometry\n+Definition: intersection.hh:60\n Dune::GridGlue::IntersectionData::index_\n IndexType index_\n index of this intersection after GridGlue interface\n Definition: intersection.hh:112\n+Dune::GridGlue::IntersectionData::Grid1IndexType\n+GridIndexType< 1 > Grid1IndexType\n+Definition: intersection.hh:74\n+Dune::GridGlue::IntersectionData::Grid1Geometry\n+GridGeometry< 1 > Grid1Geometry\n+Definition: intersection.hh:68\n+Dune::GridGlue::IntersectionData::Grid0IndexType\n+GridIndexType< 0 > Grid0IndexType\n+Definition: intersection.hh:73\n+Dune::GridGlue::IntersectionData::IntersectionData\n+IntersectionData()=default\n+Default Constructor.\n Dune::GridGlue::Intersection\n The intersection of two entities of the two patches of a GridGlue.\n Definition: intersection.hh:261\n-Dune::GridGlue::IntersectionIndexSet\n-Definition: intersectionindexset.hh:20\n-Dune::GridGlue::IntersectionIndexSet::size\n-SizeType size() const\n-Return total number of intersections.\n-Definition: intersectionindexset.hh:41\n-Dune::GridGlue::IntersectionIndexSet::IndexType\n-GridGlue::IndexType IndexType\n-The type used for the indices.\n-Definition: intersectionindexset.hh:27\n-Dune::GridGlue::IntersectionIndexSet::index\n-IndexType index(const Intersection< P0, P1, I, O > &i) const\n-Map Dune::GridGlue::Intersection to index.\n-Definition: intersectionindexset.hh:34\n-Dune::GridGlue::IntersectionIndexSet::SizeType\n-size_t SizeType\n-The type used for the size.\n-Definition: intersectionindexset.hh:29\n+Dune::GridGlue::Intersection::OutsideLocalGeometry\n+Traits::OutsideLocalGeometry OutsideLocalGeometry\n+Definition: intersection.hh:275\n+Dune::GridGlue::Intersection::flip\n+Intersection< P0, P1, O, I > flip() const\n+Return a copy of the intersection with inside and outside switched.\n+Definition: intersection.hh:483\n+Dune::GridGlue::Intersection::conforming\n+bool conforming() const\n+Return true if intersection is conforming.\n+Definition: intersection.hh:335\n+Dune::GridGlue::Intersection::InsideEntity\n+InsideGridView::Traits::template Codim< 0 >::Entity InsideEntity\n+Definition: intersection.hh:281\n+Dune::GridGlue::Intersection::GridGlue\n+Traits::GridGlue GridGlue\n+Definition: intersection.hh:267\n+Dune::GridGlue::Intersection::indexInOutside\n+int indexInOutside(unsigned int parentId=0) const\n+Local number of codim 1 entity in outside() Entity where intersection is\n+contained in.\n+Definition: intersection.hh:417\n+Dune::GridGlue::Intersection::indexInInside\n+int indexInInside(unsigned int parentId=0) const\n+Local number of codim 1 entity in the inside() Entity where intersection is\n+contained in.\n+Definition: intersection.hh:410\n+Dune::GridGlue::Intersection::self\n+bool self() const\n+For parallel computations: Return true if inside() entity exists locally.\n+Definition: intersection.hh:393\n+Dune::GridGlue::Intersection::mydim\n+static constexpr auto mydim\n+dimension of the intersection\n+Definition: intersection.hh:291\n+Dune::GridGlue::Intersection::outsidePatch\n+static constexpr int outsidePatch\n+outside patch\n+Definition: intersection.hh:297\n+Dune::GridGlue::Intersection::type\n+Dune::GeometryType type() const\n+Type of reference element for this intersection.\n+Definition: intersection.hh:377\n+Dune::GridGlue::Intersection::inside\n+InsideEntity inside(unsigned int parentId=0) const\n+Return element on the inside of this intersection.\n+Definition: intersection.hh:319\n+Dune::GridGlue::Intersection::LocalCoordinate\n+Traits::LocalCoordinate LocalCoordinate\n+Definition: intersection.hh:284\n+Dune::GridGlue::Intersection::geometryInOutside\n+const OutsideLocalGeometry & geometryInOutside(unsigned int parentId=0) const\n+Geometric information about this intersection in local coordinates of the\n+outside() element.\n+Definition: intersection.hh:349\n+Dune::GridGlue::Intersection::unitOuterNormal\n+GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const\n+Return a unit outer normal.\n+Definition: intersection.hh:455\n+Dune::GridGlue::Intersection::ctype\n+Traits::ctype ctype\n+Definition: intersection.hh:279\n+Dune::GridGlue::Intersection::Traits\n+IntersectionTraits< P0, P1, I, O > Traits\n+Definition: intersection.hh:265\n+Dune::GridGlue::Intersection::IntersectionData\n+Traits::IntersectionData IntersectionData\n+Definition: intersection.hh:268\n+Dune::GridGlue::Intersection::GlobalCoordinate\n+Traits::GlobalCoordinate GlobalCoordinate\n+Definition: intersection.hh:285\n+Dune::GridGlue::Intersection::OutsideGeometry\n+Traits::OutsideGeometry OutsideGeometry\n+Definition: intersection.hh:276\n+Dune::GridGlue::Intersection::neighbor\n+size_t neighbor(unsigned int g=0) const\n+Return number of embeddings into local grid0 (grid1) entities.\n+Definition: intersection.hh:399\n+Dune::GridGlue::Intersection::Intersection\n+Intersection(const GridGlue *glue, const IntersectionData *i)\n+Constructor for a given Dataset.\n+Definition: intersection.hh:311\n+Dune::GridGlue::Intersection::insidePatch\n+static constexpr int insidePatch\n+inside patch\n+Definition: intersection.hh:294\n+Dune::GridGlue::Intersection::Geometry\n+Traits::Geometry Geometry\n+Definition: intersection.hh:278\n+Dune::GridGlue::Intersection::OutsideEntity\n+OutsideGridView::Traits::template Codim< 0 >::Entity OutsideEntity\n+Definition: intersection.hh:282\n+Dune::GridGlue::Intersection::outside\n+OutsideEntity outside(unsigned int parentId=0) const\n+Return element on the outside of this intersection.\n+Definition: intersection.hh:328\n+Dune::GridGlue::Intersection::geometry\n+const Geometry & geometry() const\n+Geometric information about this intersection as part of the inside grid.\n+Definition: intersection.hh:360\n+Dune::GridGlue::Intersection::InsideLocalGeometry\n+Traits::InsideLocalGeometry InsideLocalGeometry\n+Definition: intersection.hh:272\n+Dune::GridGlue::Intersection::outerNormal\n+GlobalCoordinate outerNormal(const LocalCoordinate &local) const\n+Return an outer normal (length not necessarily 1)\n+Definition: intersection.hh:427\n+Dune::GridGlue::Intersection::geometryOutside\n+const OutsideGeometry & geometryOutside() const\n+Geometric information about this intersection as part of the outside grid.\n+Definition: intersection.hh:371\n+Dune::GridGlue::Intersection::OutsideGridView\n+Traits::OutsideGridView OutsideGridView\n+Definition: intersection.hh:274\n+Dune::GridGlue::Intersection::integrationOuterNormal\n+GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const\n+Return an outer normal with the length of the integration element.\n+Definition: intersection.hh:466\n+Dune::GridGlue::Intersection::InsideGridView\n+Traits::InsideGridView InsideGridView\n+Definition: intersection.hh:271\n+Dune::GridGlue::Intersection::geometryInInside\n+const InsideLocalGeometry & geometryInInside(unsigned int parentId=0) const\n+Geometric information about this intersection in local coordinates of the\n+inside() element.\n+Definition: intersection.hh:342\n+Dune::GridGlue::Intersection::centerUnitOuterNormal\n+GlobalCoordinate centerUnitOuterNormal() const\n+Unit outer normal at the center of the intersection.\n+Definition: intersection.hh:475\n+Dune::GridGlue::Intersection::coorddim\n+static constexpr auto coorddim\n+dimension of the world space of the intersection\n+Definition: intersection.hh:288\n+Dune::GridGlue::IntersectionTraits\n+Definition: intersection.hh:234\n+Dune::GridGlue::IntersectionTraits::insidePatch\n+static constexpr int insidePatch\n+Definition: intersection.hh:249\n+Dune::GridGlue::IntersectionTraits::LocalCoordinate\n+Dune::FieldVector< ctype, mydim > LocalCoordinate\n+Definition: intersection.hh:253\n+Dune::GridGlue::IntersectionTraits::coorddim\n+static constexpr auto coorddim\n+Definition: intersection.hh:247\n+Dune::GridGlue::IntersectionTraits::OutsideGeometry\n+typename IntersectionData::template GridGeometry< outside > OutsideGeometry\n+Definition: intersection.hh:245\n+Dune::GridGlue::IntersectionTraits::ctype\n+typename GridGlue::ctype ctype\n+Definition: intersection.hh:252\n+Dune::GridGlue::IntersectionTraits::OutsideGridView\n+typename GridGlue::template GridView< outside > OutsideGridView\n+Definition: intersection.hh:239\n+Dune::GridGlue::IntersectionTraits::GlobalCoordinate\n+Dune::FieldVector< ctype, coorddim > GlobalCoordinate\n+Definition: intersection.hh:254\n+Dune::GridGlue::IntersectionTraits::OutsideLocalGeometry\n+typename IntersectionData::template GridLocalGeometry< outside >\n+OutsideLocalGeometry\n+Definition: intersection.hh:242\n+Dune::GridGlue::IntersectionTraits::InsideGridView\n+typename GridGlue::template GridView< inside > InsideGridView\n+Definition: intersection.hh:238\n+Dune::GridGlue::IntersectionTraits::Geometry\n+typename IntersectionData::template GridGeometry< inside > Geometry\n+Definition: intersection.hh:244\n+Dune::GridGlue::IntersectionTraits::InsideLocalGeometry\n+typename IntersectionData::template GridLocalGeometry< inside >\n+InsideLocalGeometry\n+Definition: intersection.hh:241\n+Dune::GridGlue::IntersectionTraits::outsidePatch\n+static constexpr int outsidePatch\n+Definition: intersection.hh:250\n+Dune::GridGlue::IntersectionTraits::mydim\n+static constexpr auto mydim\n+Definition: intersection.hh:248\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00017.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00017.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: intersection.hh File Reference \n+dune-grid-glue: gridgluecommunicate.hh File Reference \n \n \n \n \n \n \n \n@@ -65,73 +65,73 @@\n dune grid-glue adapter \n
\n \n \n+ 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n- 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+Go to the documentation of this file.\n+\n+\n-\n \n+\n \n \n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n
3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n+ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n \n-\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 50 static constexpr int dim() { return GridGlue::template GridView<side>::Grid::dimension - GridGlue::template GridPatch<side>::codim; }
\n \n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+ 58 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimension>;
\n+\n+ 60 using Grid0LocalGeometry [[deprecated("please use GridLocalGeometry<0> instead")]] = GridLocalGeometry<0>;
\n+ 61 using Grid1LocalGeometry [[deprecated("please use GridLocalGeometry<1> instead")]] = GridLocalGeometry<1>;
\n+\n+\n+\n+ 65 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimensionworld>;
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 77 IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset, bool grid0local, bool grid1local);
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 160 static constexpr int elementdim = GridGlue::template GridView<side>::template Codim<0>::Geometry::mydimension;
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 199 std::array<Dune::FieldVector<ctype, GridGlue::template GridView<side>::dimensionworld>, nSimplexCorners> corners_global;
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 213 IntersectionData<P0, P1>::IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset,
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 440 normal[0] = (jacobianTransposed[0][1] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][1]);
\n+ 441 normal[1] = - (jacobianTransposed[0][0] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][0]);
\n+ 442 normal[2] = (jacobianTransposed[0][0] * jacobianTransposed[1][1] - jacobianTransposed[0][1] * jacobianTransposed[1][0]);
\n+\n+ 444 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the 'outerNormal' method for " << mydim << "-dimensional intersections yet");
\n+\n+ 446 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the 'outerNormal' method for intersections with codim >= 2 yet");
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n Central component of the module implementing the coupling of two grids.
Model of the Intersection concept provided by GridGlue.
Definition: gridglue.hh:37
bool inside(const Coordinate &x, const Field &epsilon)
Definition: projection_impl.hh:111
sequential adapter to couple two grids at specified close together boundaries
Definition: gridglue.hh:67
PromotionTraits< typenameGridView< 0 >::ctype, typenameGridView< 1 >::ctype >::PromotedType ctype
The type used for coordinates.
Definition: gridglue.hh:171
static constexpr int dimworld
export the world dimension This is the maximum of the extractors' world dimensions.
Definition: gridglue.hh:166
storage class for Dune::GridGlue::Intersection related data
Definition: intersection.hh:38
GridGlue::IndexType IndexType
Definition: intersection.hh:42
AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim, GridGlue::template GridView< side >::dimensionworld > GridGeometry
Definition: intersection.hh:65
const GridLocalGeometry< side > & localGeometry(unsigned int parentId=0) const
Definition: intersection.hh:85
static constexpr int coorddim
Dimension of the world space of the intersection.
Definition: intersection.hh:45
GridGeometry< 0 > Grid0Geometry
Definition: intersection.hh:67
typename GridGlue::template GridView< side >::IndexSet::IndexType GridIndexType
Definition: intersection.hh:71
IndexType index(unsigned int parentId=0) const
Definition: intersection.hh:97
GridLocalGeometry< 1 > Grid1LocalGeometry
Definition: intersection.hh:61
const GridGeometry< side > & geometry() const
Definition: intersection.hh:89
::Dune::GridGlue::GridGlue< P0, P1 > GridGlue
Definition: intersection.hh:40
static constexpr int mydim
Dimension of the intersection.
Definition: intersection.hh:54
AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim, GridGlue::template GridView< side >::dimension > GridLocalGeometry
Definition: intersection.hh:58
GridLocalGeometry< 0 > Grid0LocalGeometry
Definition: intersection.hh:60
IndexType index_
index of this intersection after GridGlue interface
Definition: intersection.hh:112
GridIndexType< 1 > Grid1IndexType
Definition: intersection.hh:74
GridGeometry< 1 > Grid1Geometry
Definition: intersection.hh:68
GridIndexType< 0 > Grid0IndexType
Definition: intersection.hh:73
The intersection of two entities of the two patches of a GridGlue.
Definition: intersection.hh:261
Definition: intersectionindexset.hh:20
SizeType size() const
Return total number of intersections.
Definition: intersectionindexset.hh:41
GridGlue::IndexType IndexType
The type used for the indices.
Definition: intersectionindexset.hh:27
IndexType index(const Intersection< P0, P1, I, O > &i) const
Map Dune::GridGlue::Intersection to index.
Definition: intersectionindexset.hh:34
size_t SizeType
The type used for the size.
Definition: intersectionindexset.hh:29
Traits::OutsideLocalGeometry OutsideLocalGeometry
Definition: intersection.hh:275
Intersection< P0, P1, O, I > flip() const
Return a copy of the intersection with inside and outside switched.
Definition: intersection.hh:483
bool conforming() const
Return true if intersection is conforming.
Definition: intersection.hh:335
InsideGridView::Traits::template Codim< 0 >::Entity InsideEntity
Definition: intersection.hh:281
int indexInOutside(unsigned int parentId=0) const
Local number of codim 1 entity in outside() Entity where intersection is contained in.
Definition: intersection.hh:417
int indexInInside(unsigned int parentId=0) const
Local number of codim 1 entity in the inside() Entity where intersection is contained in.
Definition: intersection.hh:410
bool self() const
For parallel computations: Return true if inside() entity exists locally.
Definition: intersection.hh:393
static constexpr auto mydim
dimension of the intersection
Definition: intersection.hh:291
static constexpr int outsidePatch
outside patch
Definition: intersection.hh:297
Dune::GeometryType type() const
Type of reference element for this intersection.
Definition: intersection.hh:377
InsideEntity inside(unsigned int parentId=0) const
Return element on the inside of this intersection.
Definition: intersection.hh:319
Traits::LocalCoordinate LocalCoordinate
Definition: intersection.hh:284
const OutsideLocalGeometry & geometryInOutside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the outside() element.
Definition: intersection.hh:349
GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const
Return a unit outer normal.
Definition: intersection.hh:455
IntersectionTraits< P0, P1, I, O > Traits
Definition: intersection.hh:265
Traits::IntersectionData IntersectionData
Definition: intersection.hh:268
Traits::GlobalCoordinate GlobalCoordinate
Definition: intersection.hh:285
Traits::OutsideGeometry OutsideGeometry
Definition: intersection.hh:276
size_t neighbor(unsigned int g=0) const
Return number of embeddings into local grid0 (grid1) entities.
Definition: intersection.hh:399
Intersection(const GridGlue *glue, const IntersectionData *i)
Constructor for a given Dataset.
Definition: intersection.hh:311
static constexpr int insidePatch
inside patch
Definition: intersection.hh:294
OutsideGridView::Traits::template Codim< 0 >::Entity OutsideEntity
Definition: intersection.hh:282
OutsideEntity outside(unsigned int parentId=0) const
Return element on the outside of this intersection.
Definition: intersection.hh:328
const Geometry & geometry() const
Geometric information about this intersection as part of the inside grid.
Definition: intersection.hh:360
Traits::InsideLocalGeometry InsideLocalGeometry
Definition: intersection.hh:272
GlobalCoordinate outerNormal(const LocalCoordinate &local) const
Return an outer normal (length not necessarily 1)
Definition: intersection.hh:427
const OutsideGeometry & geometryOutside() const
Geometric information about this intersection as part of the outside grid.
Definition: intersection.hh:371
Traits::OutsideGridView OutsideGridView
Definition: intersection.hh:274
GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const
Return an outer normal with the length of the integration element.
Definition: intersection.hh:466
Traits::InsideGridView InsideGridView
Definition: intersection.hh:271
const InsideLocalGeometry & geometryInInside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the inside() element.
Definition: intersection.hh:342
GlobalCoordinate centerUnitOuterNormal() const
Unit outer normal at the center of the intersection.
Definition: intersection.hh:475
static constexpr auto coorddim
dimension of the world space of the intersection
Definition: intersection.hh:288
Definition: intersection.hh:234
static constexpr int insidePatch
Definition: intersection.hh:249
Dune::FieldVector< ctype, mydim > LocalCoordinate
Definition: intersection.hh:253
static constexpr auto coorddim
Definition: intersection.hh:247
typename IntersectionData::template GridGeometry< outside > OutsideGeometry
Definition: intersection.hh:245
typename GridGlue::ctype ctype
Definition: intersection.hh:252
typename GridGlue::template GridView< outside > OutsideGridView
Definition: intersection.hh:239
Dune::FieldVector< ctype, coorddim > GlobalCoordinate
Definition: intersection.hh:254
typename IntersectionData::template GridLocalGeometry< outside > OutsideLocalGeometry
Definition: intersection.hh:242
typename GridGlue::template GridView< inside > InsideGridView
Definition: intersection.hh:238
typename IntersectionData::template GridGeometry< inside > Geometry
Definition: intersection.hh:244
typename IntersectionData::template GridLocalGeometry< inside > InsideLocalGeometry
Definition: intersection.hh:241
static constexpr int outsidePatch
Definition: intersection.hh:250
\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,106 +5,758 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-intersectionindexset.hh\n+intersection.hh\n Go_to_the_documentation_of_this_file.\n-1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n+ 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-\n+ 2// vi: set et ts=4 sw=2 sts=2:\n+ 3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n- 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n+ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 3#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n- 4#define DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n- 5\n- 6#include gridgluecommunicate.hh File Reference
\n \n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n+
\n+
\n+
\n+
\n+
\n\n \n
\n-
\n+
\n \n+
\n \n Model of the Intersection concept provided by GridGlue. \n+
Describes the parallel communication interface class for Dune::GridGlue. \n More...
\n-#include <algorithm>\n-
#include <optional>\n-
#include <tuple>\n-
#include <dune/common/deprecated.hh>\n-
#include <dune/common/version.hh>\n-
#include <dune/geometry/affinegeometry.hh>\n-
#include <dune/geometry/referenceelements.hh>\n-
#include <dune/grid-glue/gridglue.hh>\n+
#include <type_traits>\n+
#include <dune/common/bartonnackmanifcheck.hh>\n+
#include <dune/common/parallel/communicator.hh>\n+
#include <dune/grid/common/datahandleif.hh>\n+
#include <dune/grid/common/gridenums.hh>\n
Go to the source code of this file.
\n\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
\n-Macros | |
| #define | ONLY_SIMPLEX_INTERSECTIONS |
\n+Typedefs | |
| typedef std::pair< int, int > | Dune::GridGlue::RankPair |
| typedef CommunicationOperator< Dune::ForwardCommunication > | Dune::GridGlue::ForwardOperator |
| typedef CommunicationOperator< Dune::BackwardCommunication > | Dune::GridGlue::BackwardOperator |
\n+Functions | |
| std::ostream & | Dune::GridGlue::operator<< (std::ostream &os, const GlobalId &id) |
Detailed Description
\n-Model of the Intersection concept provided by GridGlue.
\n+Describes the parallel communication interface class for Dune::GridGlue.
\n- Author
- Christian Engwer
Macro Definition Documentation
\n-\n-◆ ONLY_SIMPLEX_INTERSECTIONS
\n-\n-\n-
\n-\n- \n-
\n-
| #define ONLY_SIMPLEX_INTERSECTIONS | \n-
\n-\n-
\n-\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,47 +5,62 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces | Macros\n-intersection.hh File Reference\n-Model of the Intersection concept provided by GridGlue. More...\n-#include \n- \n+ \n
\n intersection.hh
gridgluecommunicate.hh
\n Go to the documentation of this file.
\n \n \n \n \n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n \n-\n \n-\n-\n-\n-\n-\n \n-\n-\n-\n \n-\n-\n-\n-\n-\n-\n-\n-\n \n-\n-\n-\n-\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n
3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 50 static constexpr int dim() { return GridGlue::template GridView<side>::Grid::dimension - GridGlue::template GridPatch<side>::codim; }
\n-\n-\n-\n-\n-\n-\n- 58 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimension>;
\n-\n- 60 using Grid0LocalGeometry [[deprecated("please use GridLocalGeometry<0> instead")]] = GridLocalGeometry<0>;
\n- 61 using Grid1LocalGeometry [[deprecated("please use GridLocalGeometry<1> instead")]] = GridLocalGeometry<1>;
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n- 65 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimensionworld>;
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 77 IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset, bool grid0local, bool grid1local);
\n-\n-\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 160 static constexpr int elementdim = GridGlue::template GridView<side>::template Codim<0>::Geometry::mydimension;
\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 199 std::array<Dune::FieldVector<ctype, GridGlue::template GridView<side>::dimensionworld>, nSimplexCorners> corners_global;
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 213 IntersectionData<P0, P1>::IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset,
\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 197 static const typename CommInfo::DataType& gather(const CommInfo& commInfo, size_t i, size_t j = 0)
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+ 227 static void scatter(CommInfo& commInfo, const typename CommInfo::DataType& v, std::size_t i, std::size_t j = 0)
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 440 normal[0] = (jacobianTransposed[0][1] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][1]);
\n- 441 normal[1] = - (jacobianTransposed[0][0] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][0]);
\n- 442 normal[2] = (jacobianTransposed[0][0] * jacobianTransposed[1][1] - jacobianTransposed[0][1] * jacobianTransposed[1][0]);
\n-\n- 444 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the 'outerNormal' method for " << mydim << "-dimensional intersections yet");
\n-\n- 446 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the 'outerNormal' method for intersections with codim >= 2 yet");
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-Central component of the module implementing the coupling of two grids.
Definition: gridglue.hh:37
bool inside(const Coordinate &x, const Field &epsilon)
Definition: projection_impl.hh:111
sequential adapter to couple two grids at specified close together boundaries
Definition: gridglue.hh:67
PromotionTraits< typenameGridView< 0 >::ctype, typenameGridView< 1 >::ctype >::PromotedType ctype
The type used for coordinates.
Definition: gridglue.hh:171
static constexpr int dimworld
export the world dimension This is the maximum of the extractors' world dimensions.
Definition: gridglue.hh:166
storage class for Dune::GridGlue::Intersection related data
Definition: intersection.hh:38
GridGlue::IndexType IndexType
Definition: intersection.hh:42
AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim, GridGlue::template GridView< side >::dimensionworld > GridGeometry
Definition: intersection.hh:65
const GridLocalGeometry< side > & localGeometry(unsigned int parentId=0) const
Definition: intersection.hh:85
static constexpr int coorddim
Dimension of the world space of the intersection.
Definition: intersection.hh:45
GridGeometry< 0 > Grid0Geometry
Definition: intersection.hh:67
typename GridGlue::template GridView< side >::IndexSet::IndexType GridIndexType
Definition: intersection.hh:71
IndexType index(unsigned int parentId=0) const
Definition: intersection.hh:97
GridLocalGeometry< 1 > Grid1LocalGeometry
Definition: intersection.hh:61
const GridGeometry< side > & geometry() const
Definition: intersection.hh:89
::Dune::GridGlue::GridGlue< P0, P1 > GridGlue
Definition: intersection.hh:40
static constexpr int mydim
Dimension of the intersection.
Definition: intersection.hh:54
AffineGeometry< typename GridGlue::template GridView< side >::ctype, mydim, GridGlue::template GridView< side >::dimension > GridLocalGeometry
Definition: intersection.hh:58
GridLocalGeometry< 0 > Grid0LocalGeometry
Definition: intersection.hh:60
IndexType index_
index of this intersection after GridGlue interface
Definition: intersection.hh:112
GridIndexType< 1 > Grid1IndexType
Definition: intersection.hh:74
GridGeometry< 1 > Grid1Geometry
Definition: intersection.hh:68
GridIndexType< 0 > Grid0IndexType
Definition: intersection.hh:73
CommunicationOperator< Dune::BackwardCommunication > BackwardOperator
Definition: gridgluecommunicate.hh:264
CommunicationOperator< Dune::ForwardCommunication > ForwardOperator
Definition: gridgluecommunicate.hh:263
std::ostream & operator<<(std::ostream &os, const GlobalId &id)
Definition: gridgluecommunicate.hh:55
The intersection of two entities of the two patches of a GridGlue.
Definition: intersection.hh:261
Traits::OutsideLocalGeometry OutsideLocalGeometry
Definition: intersection.hh:275
Intersection< P0, P1, O, I > flip() const
Return a copy of the intersection with inside and outside switched.
Definition: intersection.hh:483
bool conforming() const
Return true if intersection is conforming.
Definition: intersection.hh:335
InsideGridView::Traits::template Codim< 0 >::Entity InsideEntity
Definition: intersection.hh:281
int indexInOutside(unsigned int parentId=0) const
Local number of codim 1 entity in outside() Entity where intersection is contained in.
Definition: intersection.hh:417
int indexInInside(unsigned int parentId=0) const
Local number of codim 1 entity in the inside() Entity where intersection is contained in.
Definition: intersection.hh:410
bool self() const
For parallel computations: Return true if inside() entity exists locally.
Definition: intersection.hh:393
static constexpr auto mydim
dimension of the intersection
Definition: intersection.hh:291
static constexpr int outsidePatch
outside patch
Definition: intersection.hh:297
Dune::GeometryType type() const
Type of reference element for this intersection.
Definition: intersection.hh:377
InsideEntity inside(unsigned int parentId=0) const
Return element on the inside of this intersection.
Definition: intersection.hh:319
Traits::LocalCoordinate LocalCoordinate
Definition: intersection.hh:284
const OutsideLocalGeometry & geometryInOutside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the outside() element.
Definition: intersection.hh:349
GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const
Return a unit outer normal.
Definition: intersection.hh:455
IntersectionTraits< P0, P1, I, O > Traits
Definition: intersection.hh:265
Traits::IntersectionData IntersectionData
Definition: intersection.hh:268
Traits::GlobalCoordinate GlobalCoordinate
Definition: intersection.hh:285
Traits::OutsideGeometry OutsideGeometry
Definition: intersection.hh:276
size_t neighbor(unsigned int g=0) const
Return number of embeddings into local grid0 (grid1) entities.
Definition: intersection.hh:399
Intersection(const GridGlue *glue, const IntersectionData *i)
Constructor for a given Dataset.
Definition: intersection.hh:311
static constexpr int insidePatch
inside patch
Definition: intersection.hh:294
OutsideGridView::Traits::template Codim< 0 >::Entity OutsideEntity
Definition: intersection.hh:282
OutsideEntity outside(unsigned int parentId=0) const
Return element on the outside of this intersection.
Definition: intersection.hh:328
const Geometry & geometry() const
Geometric information about this intersection as part of the inside grid.
Definition: intersection.hh:360
Traits::InsideLocalGeometry InsideLocalGeometry
Definition: intersection.hh:272
GlobalCoordinate outerNormal(const LocalCoordinate &local) const
Return an outer normal (length not necessarily 1)
Definition: intersection.hh:427
const OutsideGeometry & geometryOutside() const
Geometric information about this intersection as part of the outside grid.
Definition: intersection.hh:371
Traits::OutsideGridView OutsideGridView
Definition: intersection.hh:274
GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const
Return an outer normal with the length of the integration element.
Definition: intersection.hh:466
Traits::InsideGridView InsideGridView
Definition: intersection.hh:271
const InsideLocalGeometry & geometryInInside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the inside() element.
Definition: intersection.hh:342
GlobalCoordinate centerUnitOuterNormal() const
Unit outer normal at the center of the intersection.
Definition: intersection.hh:475
static constexpr auto coorddim
dimension of the world space of the intersection
Definition: intersection.hh:288
Definition: intersection.hh:234
static constexpr int insidePatch
Definition: intersection.hh:249
Dune::FieldVector< ctype, mydim > LocalCoordinate
Definition: intersection.hh:253
static constexpr auto coorddim
Definition: intersection.hh:247
typename IntersectionData::template GridGeometry< outside > OutsideGeometry
Definition: intersection.hh:245
typename GridGlue::ctype ctype
Definition: intersection.hh:252
typename GridGlue::template GridView< outside > OutsideGridView
Definition: intersection.hh:239
Dune::FieldVector< ctype, coorddim > GlobalCoordinate
Definition: intersection.hh:254
typename IntersectionData::template GridLocalGeometry< outside > OutsideLocalGeometry
Definition: intersection.hh:242
typename GridGlue::template GridView< inside > InsideGridView
Definition: intersection.hh:238
typename IntersectionData::template GridGeometry< inside > Geometry
Definition: intersection.hh:244
typename IntersectionData::template GridLocalGeometry< inside > InsideLocalGeometry
Definition: intersection.hh:241
static constexpr int outsidePatch
Definition: intersection.hh:250
Definition: gridgluecommunicate.hh:26
GlobalId(int i, int j, unsigned int n)
Definition: gridgluecommunicate.hh:48
describes the features of a data handle for communication in parallel runs using the GridGlue::commun...
Definition: gridgluecommunicate.hh:77
size_t size(RISType &i) const
Definition: gridgluecommunicate.hh:92
void scatter(MessageBufferImp &buff, const EntityType &e, const RISType &i, size_t n)
Definition: gridgluecommunicate.hh:118
void gather(MessageBufferImp &buff, const EntityType &e, const RISType &i) const
pack data from user to message buffer
Definition: gridgluecommunicate.hh:104
DataTypeImp DataType
data type of data to communicate
Definition: gridgluecommunicate.hh:80
Definition: gridgluecommunicate.hh:141
StreamingMessageBuffer(DT *p)
Definition: gridgluecommunicate.hh:146
size_t counter() const
Definition: gridgluecommunicate.hh:169
void write(const Y &data)
Definition: gridgluecommunicate.hh:155
void read(Y &data) const
Definition: gridgluecommunicate.hh:163
DT value_type
Definition: gridgluecommunicate.hh:143
forward gather scatter to user defined CommInfo class
Definition: gridgluecommunicate.hh:194
static void scatter(CommInfo &commInfo, const typename CommInfo::DataType &v, std::size_t i, std::size_t j=0)
Definition: gridgluecommunicate.hh:227
static const CommInfo::DataType & gather(const CommInfo &commInfo, size_t i, size_t j=0)
Definition: gridgluecommunicate.hh:197
collects all GridGlue data requried for communication
Definition: gridgluecommunicate.hh:272
::Dune::GridGlue::CommDataHandle< DataHandleImp, DataTypeImp > * data
Definition: gridgluecommunicate.hh:282
mutable ::Dune::GridGlue::StreamingMessageBuffer< DataType > mbuffer
Definition: gridgluecommunicate.hh:286
std::vector< DataType > buffer
Definition: gridgluecommunicate.hh:285
static size_t getSize(const Type &commInfo, size_t i)
Get the number of objects at an intersection.
Definition: gridgluecommunicate.hh:320
DataTypeImp IndexedType
The datatype that should be communicated.
Definition: gridgluecommunicate.hh:309
::Dune::GridGlue::CommInfo< GG, DataHandleImp, DataTypeImp > Type
The type of the GridGlueCommInfo.
Definition: gridgluecommunicate.hh:304
VariableSize IndexedTypeFlag
Each intersection can communicate a different number of objects.
Definition: gridgluecommunicate.hh:315
\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,758 +5,443 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-intersection.hh\n+gridgluecommunicate.hh\n Go_to_the_documentation_of_this_file.\n 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-\n 2// vi: set et ts=4 sw=2 sts=2:\n 3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 11#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n- 12#define DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH\n- 13\n- 14#include \n \n- \n+Namespaces
\n+ gridgluecommunicate.hh File Reference
gridgluevtkwriter.hh File Reference
\n \n-
\n-
\n-
\n-
\n-
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n\n \n
\n
\n-
\n-
\n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,62 +5,34 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces | Typedefs | Functions\n-gridgluecommunicate.hh File Reference\n-Describes the parallel communication interface class for Dune::GridGlue.\n-More...\n+Classes | Namespaces\n+gridgluevtkwriter.hh File Reference\n+Write all remote intersections to a vtk file for debugging. More...\n+#include \n+#include \n #include \n-#include \n-#include \n-#include \n-#include \n+#include \n+#include \n+#include \n+#include \n+#include \n Go_to_the_source_code_of_this_file.\n Classes\n-struct \u00a0Dune::GridGlue::GlobalId\n-\u00a0\n- class \u00a0Dune::GridGlue::CommDataHandle<_DataHandleImp,_DataTypeImp_>\n-\u00a0 describes the features of a data handle for communication in parallel\n- runs using the GridGlue::communicate methods. More...\n-\u00a0\n- class \u00a0Dune::GridGlue::StreamingMessageBuffer<_DT_>\n-\u00a0\n- class \u00a0Dune::GridGlue::CommunicationOperator<_dir_>\n-\u00a0 forward gather scatter to user defined CommInfo class More...\n-\u00a0\n-struct \u00a0Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n-\u00a0 collects all GridGlue data requried for communication More...\n-\u00a0\n-struct \u00a0Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,\n- DataTypeImp_>_>\n-\u00a0 specialization of the CommPolicy struct, required for the\n- ParallelIndexsets More...\n+class \u00a0Dune::GridGlue::GridGlueVtkWriter\n+\u00a0 Write remote intersections to a vtk file for debugging purposes.\n+ More...\n \u00a0\n Namespaces\n namespace \u00a0Dune\n \u00a0\n namespace \u00a0Dune::GridGlue\n \u00a0\n- Typedefs\n- typedef std::pair< int, int >\u00a0Dune::GridGlue::RankPair\n-\u00a0\n-typedef CommunicationOperator< Dune::\n- ForwardCommunication >\u00a0Dune::GridGlue::ForwardOperator\n-\u00a0\n-typedef CommunicationOperator< Dune::\n- BackwardCommunication >\u00a0Dune::GridGlue::BackwardOperator\n-\u00a0\n- Functions\n-std::ostream &\u00a0Dune::GridGlue::operator<< (std::ostream &os, const GlobalId\n- &id)\n-\u00a0\n ***** Detailed Description *****\n-Describes the parallel communication interface class for Dune::GridGlue.\n- Author\n- Christian Engwer\n+Write all remote intersections to a vtk file for debugging.\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00020_source.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00020_source.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: gridgluecommunicate.hh Source File \n+dune-grid-glue: gridgluevtkwriter.hh Source File \n \n \n \n \n \n \n \n@@ -62,323 +62,314 @@\n \n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,443 +5,350 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-gridgluecommunicate.hh\n+gridgluevtkwriter.hh\n Go_to_the_documentation_of_this_file.\n 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-\n 2// vi: set et ts=4 sw=2 sts=2:\n 3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 5#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUECOMMUNICATE_HH\n- 6#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUECOMMUNICATE_HH\n- 7\n- 13#include \n- 14\n- 15#include \n- 16#include \n- 17#include \n- 18#include \n- 19\n- 20\n- 21namespace Dune {\n- 22 namespace GridGlue {\n+ 5/*\n+ 6 * Filename: GridGlueVtkWriter.hh\n+ 7 * Version: 1.0\n+ 8 * Created on: Mar 5, 2009\n+ 9 * Author: Gerrit Buse\n+ 10 * ---------------------------------\n+ 11 * Project: dune-grid-glue\n+ 12 * Description: Class thought to make graphical debugging of couplings\n+easier.\n+ 13 *\n+ 14 */\n+ 20#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUEVTKWRITER_HH\n+ 21#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUEVTKWRITER_HH\n+ 22\n 23\n-24 typedef std::pair RankPair;\n-25 struct GlobalId : public std::pair\n- 26 {\n-30 GlobalId() {\n- 31 this->first.first = 0;\n- 32 this->first.second = 0;\n- 33 this->second = 0;\n- 34 }\n-38 GlobalId(int i) {\n- 39 this->first.first = i;\n- 40 this->first.second = i;\n- 41 this->second = 0;\n- 42 }\n-48 GlobalId(int i, int j, unsigned int n) {\n- 49 this->first.first = std::min(i,j);\n- 50 this->first.second = std::max(i,j);\n- 51 this->second = n;\n- 52 }\n- 53 };\n+ 24#include \n+ 25#include \n+ 26#include \n+ 27#include \n+ 28\n+ 29#include \n+ 30#include \n+ 31#include \n+ 32\n+ 33#include \n+ 34\n+ 35namespace Dune {\n+ 36namespace GridGlue {\n+ 37\n+40class GridGlueVtkWriter\n+ 41{\n+ 42\n+ 46 template \n+ 47 static void writeExtractedPart(const Glue& glue, const std::string&\n+filename)\n+ 48 {\n+ 49 static_assert((side==0 || side==1), \"'side' can only be 0 or 1\");\n+ 50\n+ 51 std::ofstream fgrid;\n+ 52\n+ 53 fgrid.open(filename.c_str());\n 54\n-55 inline std::ostream& operator<<(std::ostream& os, const GlobalId & id)\n- 56 {\n- 57 os << \"(\"\n- 58 << id.first.first << \",\" << id.first.second << \",\"\n- 59 << id.second << \")\";\n- 60 return os;\n- 61 }\n+ 55 using GridView = typename Glue::template GridView;\n+ 56 using Extractor = typename Glue::template GridPatch;\n+ 57\n+ 58 typedef typename GridView::ctype ctype;\n+ 59\n+ 60 const int domdimw = GridView::dimensionworld;\n+ 61 const int patchDim = Extractor::dim - Extractor::codim;\n 62\n- 75 template \n-76 class CommDataHandle\n- 77 {\n- 78 public:\n-80 typedef DataTypeImp DataType;\n- 81\n- 82 protected:\n- 83 // one should not create an explicit instance of this inteface object\n-84 CommDataHandle() {}\n+ 63 // coordinates have to be in R^3 in the VTK format\n+ 64 std::string coordinatePadding;\n+ 65 for (int i=domdimw; i<3; i++)\n+ 66 coordinatePadding += \" 0\";\n+ 67\n+ 68 fgrid << \"# vtk DataFile Version 2.0\\nFilename: \" << filename << \"\\nASCII\"\n+<< std::endl;\n+ 69\n+ 70 // WRITE POINTS\n+ 71 // ----------------\n+ 72 std::vector coords;\n+ 73 glue.template patch().getCoords(coords);\n+ 74\n+ 75 fgrid << ((patchDim==3) ? \"DATASET UNSTRUCTURED_GRID\" : \"DATASET POLYDATA\")\n+<< std::endl;\n+ 76 fgrid << \"POINTS \" << coords.size() << \" \" << Dune::className() <<\n+std::endl;\n+ 77\n+ 78 for (size_t i=0; i\n-92 size_t size (RISType& i) const\n- 93 {\n- 94 CHECK_INTERFACE_IMPLEMENTATION((asImp().size(i)));\n- 95 return asImp().size(i);\n- 96 }\n+ 86 std::vector faces;\n+ 87 std::vector geometryTypes;\n+ 88 glue.template patch().getFaces(faces);\n+ 89 glue.template patch().getGeometryTypes(geometryTypes);\n+ 90\n+ 91 unsigned int faceCornerCount = 0;\n+ 92 for (size_t i=0; i\n-104 void gather (MessageBufferImp& buff, const EntityType& e, const RISType &\n-i) const\n- 105 {\n- 106 MessageBufferIF buffIF(buff);\n- 107 CHECK_AND_CALL_INTERFACE_IMPLEMENTATION((asImp().gather(buffIF,e,i)));\n- 108 }\n- 109\n- 117 template\n-118 void scatter (MessageBufferImp& buff, const EntityType& e, const RISType &\n-i, size_t n)\n- 119 {\n- 120 MessageBufferIF buffIF(buff);\n- 121 CHECK_AND_CALL_INTERFACE_IMPLEMENTATION((asImp().scatter(buffIF,e,i,n)));\n- 122 }\n- 123\n- 124 private:\n- 126 DataHandleImp& asImp () {\n- 127 return static_cast (*this);\n+ 98 for (size_t i=0; i(*this);\n- 133 }\n- 134 }; // end class CommDataHandleIF\n- 135\n- 140 template\n-141 class StreamingMessageBuffer {\n- 142 public:\n-143 typedef DT value_type;\n- 144\n- 145 // Constructor\n-146 StreamingMessageBuffer (DT *p)\n- 147 {\n- 148 a=p;\n- 149 i=0;\n- 150 j=0;\n- 151 }\n- 152\n- 153 // write data to message buffer, acts like a stream !\n- 154 template\n-155 void write (const Y& data)\n- 156 {\n- 157 static_assert(std::is_same::value, \"DataType mismatch\");\n- 158 a[i++] = data;\n- 159 }\n+ 129\n+ 130 fgrid << std::endl;\n+ 131 }\n+ 132\n+ 133 fgrid << std::endl;\n+ 134\n+ 135 // 3d VTK files need an extra section specifying the CELL_TYPES aka\n+GeometryTypes\n+ 136 if (patchDim==3) {\n+ 137\n+ 138 fgrid << \"CELL_TYPES \" << geometryTypes.size() << std::endl;\n+ 139\n+ 140 for (size_t i=0; i\n-163 void read (Y& data) const\n- 164 {\n- 165 static_assert(std::is_same::value, \"DataType mismatch\");\n- 166 data = a[j++];\n- 167 }\n- 168\n-169 size_t counter() const { return i; }\n- 170\n-171 void clear()\n- 172 {\n- 173 i = 0;\n- 174 j = 0;\n- 175 }\n+ 161 fgrid << \"CELL_DATA \" << gridSubEntityData.size() << std::endl;\n+ 162 fgrid << \"SCALARS property_coding \" << Dune::className() << \" 1\" <<\n+std::endl;\n+ 163 fgrid << \"LOOKUP_TABLE default\" << std::endl;\n+ 164\n+ 165 for (typename GridSubEntityData::const_iterator sEIt =\n+gridSubEntityData.begin();\n+ 166 sEIt != gridSubEntityData.end();\n+ 167 ++sEIt, accum += delta)\n+ 168 {\n+ 169 // \"encode\" the parent with one color...\n+ 170 fgrid << accum << std::endl;\n+ 171 }\n+ 172#endif\n+ 173 fgrid.close();\n+ 174 }\n+ 175\n 176\n- 177 // we need access to these variables in an assertion\n- 178#ifdef NDEBUG\n- 179 private:\n- 180#endif\n-181 DT *a;\n-182 size_t i;\n-183 mutable size_t j;\n- 184 };\n- 185\n- 192 template\n-193 class CommunicationOperator\n- 194 {\n- 195 public:\n- 196 template\n-197 static const typename CommInfo::DataType& gather(const CommInfo& commInfo,\n-size_t i, size_t j = 0)\n- 198 {\n- 199 // get Intersection\n- 200 typedef typename CommInfo::GridGlue::Intersection Intersection;\n- 201 Intersection ris(commInfo.gridglue->getIntersection(i));\n- 202\n- 203 // fill buffer if we have a new intersection\n- 204 if (j == 0)\n- 205 {\n- 206 commInfo.mbuffer.clear();\n- 207 if (dir == Dune::ForwardCommunication)\n- 208 {\n- 209 // read from grid0\n- 210 if(ris.self())\n- 211 commInfo.data->gather(commInfo.mbuffer, ris.inside(), ris);\n- 212 }\n- 213 else // (dir == Dune::BackwardCommunication)\n+ 180 template \n+ 181 static void writeIntersections(const Glue& glue, const std::string&\n+filename)\n+ 182 {\n+ 183 static_assert((side==0 || side==1), \"'side' can only be 0 or 1\");\n+ 184\n+ 185 std::ofstream fmerged;\n+ 186\n+ 187 fmerged.open(filename.c_str());\n+ 188\n+ 189 using GridView = typename Glue::template GridView;\n+ 190 typedef typename GridView::ctype ctype;\n+ 191\n+ 192 const int domdimw = GridView::dimensionworld;\n+ 193 const int intersectionDim = Glue::Intersection::mydim;\n+ 194\n+ 195 // coordinates have to be in R^3 in the VTK format\n+ 196 std::string coordinatePadding;\n+ 197 for (int i=domdimw; i<3; i++)\n+ 198 coordinatePadding += \" 0\";\n+ 199\n+ 200 int overlaps = glue.size();\n+ 201\n+ 202 // WRITE POINTS\n+ 203 // ----------------\n+ 204 using Extractor = typename Glue::template GridPatch<0>;\n+ 205 std::vector coords;\n+ 206 glue.template patch().getCoords(coords);\n+ 207\n+ 208 // the merged grid (i.e. the set of remote intersections\n+ 209 fmerged << \"# vtk DataFile Version 2.0\\nFilename: \" << filename <<\n+\"\\nASCII\" << std::endl;\n+ 210 fmerged << ((intersectionDim==3) ? \"DATASET UNSTRUCTURED_GRID\" : \"DATASET\n+POLYDATA\") << std::endl;\n+ 211 fmerged << \"POINTS \" << overlaps*(intersectionDim+1) << \" \" << Dune::\n+className() << std::endl;\n+ 212\n+ 213 for (const auto& intersection : intersections(glue, Reverse{}))\n 214 {\n- 215 // read from grid1\n- 216 if(ris.neighbor())\n- 217 commInfo.data->gather(commInfo.mbuffer, ris.outside(), ris.flip());\n+ 215 const auto& geometry = intersection.geometry();\n+ 216 for (int i = 0; i < geometry.corners(); ++i)\n+ 217 fmerged << geometry.corner(i) << coordinatePadding << std::endl;\n 218 }\n- 219 }\n- 220\n- 221 // return the j'th value in the buffer\n- 222 assert(j < commInfo.mbuffer.i);\n- 223 return commInfo.buffer[j];\n- 224 }\n- 225\n- 226 template\n-227 static void scatter(CommInfo& commInfo, const typename CommInfo::DataType&\n-v, std::size_t i, std::size_t j = 0)\n- 228 {\n- 229 // extract GridGlue objects...\n- 230 typedef typename CommInfo::GridGlue::Intersection Intersection;\n- 231 Intersection ris(commInfo.gridglue->getIntersection(i));\n- 232\n- 233 // get size if we have a new intersection\n- 234 if (j == 0)\n- 235 {\n- 236 commInfo.mbuffer.clear();\n- 237 commInfo.currentsize = commInfo.data->size(ris);\n- 238 }\n- 239\n- 240 // write entry to buffer\n- 241 commInfo.buffer[j] = v;\n+ 219\n+ 220 // WRITE POLYGONS\n+ 221 // ----------------\n+ 222\n+ 223 std::vector faces;\n+ 224 std::vector geometryTypes;\n+ 225 glue.template patch().getFaces(faces);\n+ 226 glue.template patch().getGeometryTypes(geometryTypes);\n+ 227\n+ 228 unsigned int faceCornerCount = 0;\n+ 229 for (size_t i=0; iscatter(commInfo.mbuffer, ris.outside(), ris.flip(),\n-commInfo.currentsize);\n+ 243 // 3d VTK files need an extra section specifying the CELL_TYPES aka\n+GeometryTypes\n+ 244 if (intersectionDim==3) {\n+ 245\n+ 246 fmerged << \"CELL_TYPES \" << overlaps << std::endl;\n+ 247\n+ 248 for (int i = 0; i < overlaps; i++)\n+ 249 fmerged << \"10\" << std::endl;\n+ 250\n 251 }\n- 252 else // (dir == Dune::BackwardCommunication)\n- 253 {\n- 254 // write to grid0\n- 255 if(ris.self())\n- 256 commInfo.data->scatter(commInfo.mbuffer, ris.inside(), ris,\n-commInfo.currentsize);\n- 257 }\n- 258 assert(commInfo.mbuffer.j <= commInfo.currentsize);\n- 259 }\n- 260 }\n- 261 };\n- 262\n-263 typedef CommunicationOperator ForwardOperator;\n-264 typedef CommunicationOperator\n-BackwardOperator;\n- 265\n- 270 template \n-271 struct CommInfo\n- 272 {\n-273 typedef DataTypeImp value_type;\n-274 typedef GG GridGlue;\n-275 typedef DataTypeImp DataType;\n- 276\n-277 CommInfo() : buffer(100), mbuffer(&buffer[0])\n- 278 {}\n- 279\n- 280 // tunnel information to the policy and the operators\n-281 const GridGlue * gridglue;\n-282 ::Dune::GridGlue::CommDataHandle * data;\n- 283\n- 284 // state variables\n-285 std::vector buffer;\n-286 mutable ::Dune::GridGlue::StreamingMessageBuffer mbuffer;\n-287 size_t currentsize;\n-288 Dune::CommunicationDirection dir;\n- 289 };\n- 290\n- 291 } // end namespace GridGlue\n+ 252\n+ 253#if 0\n+ 254 // WRITE CELL DATA\n+ 255 // ---------------\n+ 256 ctype accum = 0.0, delta = 1.0 / (ctype) (gridSubEntityData.size()-1);\n+ 257\n+ 258 fmerged << \"CELL_DATA \" << overlaps << std::endl;\n+ 259 fmerged << \"SCALARS property_coding \" << Dune::className() << \" 1\"\n+<< std::endl;\n+ 260 fmerged << \"LOOKUP_TABLE default\" << std::endl;\n+ 261\n+ 262 for (typename GridSubEntityData::const_iterator sEIt =\n+gridSubEntityData.begin();\n+ 263 sEIt != gridSubEntityData.end();\n+ 264 ++sEIt, accum += delta)\n+ 265 {\n+ 266 // ...and mark all of its merged grid parts with the same color\n+ 267 for (int j = 0; j < sEIt->first.second; ++j)\n+ 268 fmerged << accum << std::endl;\n+ 269 }\n+ 270#endif\n+ 271 fmerged.close();\n+ 272 }\n+ 273\n+ 274public:\n+ 275 template\n+276 static void write(const Glue& glue, const std::string& filenameTrunk)\n+ 277 {\n+ 278\n+ 279 // Write extracted grid and remote intersection on the grid0-side\n+ 280 writeExtractedPart(glue,\n+ 281 filenameTrunk + \"-grid0.vtk\");\n+ 282\n+ 283 writeIntersections(glue,\n+ 284 filenameTrunk + \"-intersections-grid0.vtk\");\n+ 285\n+ 286 // Write extracted grid and remote intersection on the grid1-side\n+ 287 writeExtractedPart(glue,\n+ 288 filenameTrunk + \"-grid1.vtk\");\n+ 289\n+ 290 writeIntersections(glue,\n+ 291 filenameTrunk + \"-intersections-grid1.vtk\");\n 292\n- 293#if HAVE_MPI\n- 298 template\n-299 struct CommPolicy< ::Dune::GridGlue::CommInfo >\n- 300 {\n-304 typedef ::Dune::GridGlue::CommInfo Type;\n- 305\n-309 typedef DataTypeImp IndexedType;\n- 310\n- 314 // typedef SizeOne IndexedTypeFlag;\n-315 typedef VariableSize IndexedTypeFlag;\n- 316\n-320 static size_t getSize(const Type& commInfo, size_t i)\n- 321 {\n- 322 // get Intersection\n- 323 typedef typename Type::GridGlue::Intersection Intersection;\n- 324 Intersection ris(commInfo.gridglue->getIntersection(i));\n- 325\n- 326 // ask data handle for size\n- 327 return commInfo.data->size(ris);\n- 328 }\n- 329 };\n- 330#endif\n- 331\n- 332} // end namespace Dune\n- 333#endif\n+ 293 }\n+ 294\n+ 295};\n+ 296\n+ 297} /* namespace GridGlue */\n+ 298} /* namespace Dune */\n+ 299\n+ 300#endif // DUNE_GRIDGLUE_ADAPTER_GRIDGLUEVTKWRITER_HH\n+gridglue.hh\n+Central component of the module implementing the coupling of two grids.\n Dune\n Definition: gridglue.hh:37\n-Dune::GridGlue::BackwardOperator\n-CommunicationOperator< Dune::BackwardCommunication > BackwardOperator\n-Definition: gridgluecommunicate.hh:264\n-Dune::GridGlue::ForwardOperator\n-CommunicationOperator< Dune::ForwardCommunication > ForwardOperator\n-Definition: gridgluecommunicate.hh:263\n-Dune::GridGlue::RankPair\n-std::pair< int, int > RankPair\n-Definition: gridgluecommunicate.hh:24\n-Dune::GridGlue::operator<<\n-std::ostream & operator<<(std::ostream &os, const GlobalId &id)\n-Definition: gridgluecommunicate.hh:55\n-Dune::GridGlue::Intersection\n-The intersection of two entities of the two patches of a GridGlue.\n-Definition: intersection.hh:261\n-Dune::GridGlue::Intersection::flip\n-Intersection< P0, P1, O, I > flip() const\n-Return a copy of the intersection with inside and outside switched.\n-Definition: intersection.hh:483\n-Dune::GridGlue::Intersection::self\n-bool self() const\n-For parallel computations: Return true if inside() entity exists locally.\n-Definition: intersection.hh:393\n-Dune::GridGlue::Intersection::inside\n-InsideEntity inside(unsigned int parentId=0) const\n-Return element on the inside of this intersection.\n-Definition: intersection.hh:319\n-Dune::GridGlue::Intersection::neighbor\n-size_t neighbor(unsigned int g=0) const\n-Return number of embeddings into local grid0 (grid1) entities.\n-Definition: intersection.hh:399\n-Dune::GridGlue::Intersection::outside\n-OutsideEntity outside(unsigned int parentId=0) const\n-Return element on the outside of this intersection.\n-Definition: intersection.hh:328\n-Dune::GridGlue::GlobalId\n-Definition: gridgluecommunicate.hh:26\n-Dune::GridGlue::GlobalId::GlobalId\n-GlobalId(int i)\n-Definition: gridgluecommunicate.hh:38\n-Dune::GridGlue::GlobalId::GlobalId\n-GlobalId()\n-Definition: gridgluecommunicate.hh:30\n-Dune::GridGlue::GlobalId::GlobalId\n-GlobalId(int i, int j, unsigned int n)\n-Definition: gridgluecommunicate.hh:48\n-Dune::GridGlue::CommDataHandle\n-describes the features of a data handle for communication in parallel runs\n-using the GridGlue::commun...\n-Definition: gridgluecommunicate.hh:77\n-Dune::GridGlue::CommDataHandle::size\n-size_t size(RISType &i) const\n-Definition: gridgluecommunicate.hh:92\n-Dune::GridGlue::CommDataHandle::scatter\n-void scatter(MessageBufferImp &buff, const EntityType &e, const RISType &i,\n-size_t n)\n-Definition: gridgluecommunicate.hh:118\n-Dune::GridGlue::CommDataHandle::gather\n-void gather(MessageBufferImp &buff, const EntityType &e, const RISType &i)\n+Dune::GridGlue::intersections\n+IteratorRange<... > intersections(const GridGlue<... > &glue, const Reverse<...\n+> &reverse=!reversed)\n+Iterate over all intersections of a GridGlue.\n+Dune::GridGlue::GridGlueVtkWriter\n+Write remote intersections to a vtk file for debugging purposes.\n+Definition: gridgluevtkwriter.hh:41\n+Dune::GridGlue::GridGlueVtkWriter::write\n+static void write(const Glue &glue, const std::string &filenameTrunk)\n+Definition: gridgluevtkwriter.hh:276\n+Dune::GridGlue::Reverse\n+Definition: rangegenerators.hh:17\n+Dune::GridGlue::Extractor\n+Provides codimension-independent methods for grid extraction.\n+Definition: extractor.hh:46\n+Dune::GridGlue::Extractor::codim\n+static constexpr auto codim\n+Definition: extractor.hh:52\n+Dune::GridGlue::Extractor::getCoords\n+void getCoords(std::vector< Dune::FieldVector< ctype, dimworld > > &coords)\n const\n-pack data from user to message buffer\n-Definition: gridgluecommunicate.hh:104\n-Dune::GridGlue::CommDataHandle::DataType\n-DataTypeImp DataType\n-data type of data to communicate\n-Definition: gridgluecommunicate.hh:80\n-Dune::GridGlue::CommDataHandle::CommDataHandle\n-CommDataHandle()\n-Definition: gridgluecommunicate.hh:84\n-Dune::GridGlue::StreamingMessageBuffer\n-Definition: gridgluecommunicate.hh:141\n-Dune::GridGlue::StreamingMessageBuffer::j\n-size_t j\n-Definition: gridgluecommunicate.hh:183\n-Dune::GridGlue::StreamingMessageBuffer::StreamingMessageBuffer\n-StreamingMessageBuffer(DT *p)\n-Definition: gridgluecommunicate.hh:146\n-Dune::GridGlue::StreamingMessageBuffer::a\n-DT * a\n-Definition: gridgluecommunicate.hh:181\n-Dune::GridGlue::StreamingMessageBuffer::counter\n-size_t counter() const\n-Definition: gridgluecommunicate.hh:169\n-Dune::GridGlue::StreamingMessageBuffer::write\n-void write(const Y &data)\n-Definition: gridgluecommunicate.hh:155\n-Dune::GridGlue::StreamingMessageBuffer::read\n-void read(Y &data) const\n-Definition: gridgluecommunicate.hh:163\n-Dune::GridGlue::StreamingMessageBuffer::value_type\n-DT value_type\n-Definition: gridgluecommunicate.hh:143\n-Dune::GridGlue::StreamingMessageBuffer::i\n-size_t i\n-Definition: gridgluecommunicate.hh:182\n-Dune::GridGlue::StreamingMessageBuffer::clear\n-void clear()\n-Definition: gridgluecommunicate.hh:171\n-Dune::GridGlue::CommunicationOperator\n-forward gather scatter to user defined CommInfo class\n-Definition: gridgluecommunicate.hh:194\n-Dune::GridGlue::CommunicationOperator::scatter\n-static void scatter(CommInfo &commInfo, const typename CommInfo::DataType &v,\n-std::size_t i, std::size_t j=0)\n-Definition: gridgluecommunicate.hh:227\n-Dune::GridGlue::CommunicationOperator::gather\n-static const CommInfo::DataType & gather(const CommInfo &commInfo, size_t i,\n-size_t j=0)\n-Definition: gridgluecommunicate.hh:197\n-Dune::GridGlue::CommInfo\n-collects all GridGlue data requried for communication\n-Definition: gridgluecommunicate.hh:272\n-Dune::GridGlue::CommInfo::dir\n-Dune::CommunicationDirection dir\n-Definition: gridgluecommunicate.hh:288\n-Dune::GridGlue::CommInfo::value_type\n-DataTypeImp value_type\n-Definition: gridgluecommunicate.hh:273\n-Dune::GridGlue::CommInfo::DataType\n-DataTypeImp DataType\n-Definition: gridgluecommunicate.hh:275\n-Dune::GridGlue::CommInfo::data\n-::Dune::GridGlue::CommDataHandle< DataHandleImp, DataTypeImp > * data\n-Definition: gridgluecommunicate.hh:282\n-Dune::GridGlue::CommInfo::GridGlue\n-GG GridGlue\n-Definition: gridgluecommunicate.hh:274\n-Dune::GridGlue::CommInfo::currentsize\n-size_t currentsize\n-Definition: gridgluecommunicate.hh:287\n-Dune::GridGlue::CommInfo::CommInfo\n-CommInfo()\n-Definition: gridgluecommunicate.hh:277\n-Dune::GridGlue::CommInfo::mbuffer\n-mutable ::Dune::GridGlue::StreamingMessageBuffer< DataType > mbuffer\n-Definition: gridgluecommunicate.hh:286\n-Dune::GridGlue::CommInfo::gridglue\n-const GridGlue * gridglue\n-Definition: gridgluecommunicate.hh:281\n-Dune::GridGlue::CommInfo::buffer\n-std::vector< DataType > buffer\n-Definition: gridgluecommunicate.hh:285\n-Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n->::getSize\n-static size_t getSize(const Type &commInfo, size_t i)\n-Get the number of objects at an intersection.\n-Definition: gridgluecommunicate.hh:320\n-Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n->::IndexedType\n-DataTypeImp IndexedType\n-The datatype that should be communicated.\n-Definition: gridgluecommunicate.hh:309\n-Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n->::Type\n-::Dune::GridGlue::CommInfo< GG, DataHandleImp, DataTypeImp > Type\n-The type of the GridGlueCommInfo.\n-Definition: gridgluecommunicate.hh:304\n-Dune::CommPolicy<_::Dune::GridGlue::CommInfo<_GG,_DataHandleImp,_DataTypeImp_>\n->::IndexedTypeFlag\n-VariableSize IndexedTypeFlag\n-Each intersection can communicate a different number of objects.\n-Definition: gridgluecommunicate.hh:315\n+getter for the coordinates array\n+Definition: extractor.hh:275\n+Dune::GridGlue::Extractor::dim\n+static constexpr auto dim\n+Definition: extractor.hh:51\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00023.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00023.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: gridglueamirawriter.hh File Reference \n+dune-grid-glue: intersectionindexset.hh File Reference \n \n \n \n \n \n \n \n@@ -65,41 +65,34 @@\n dune grid-glue adapter \n
\n \n Describes the parallel communication interface class for Dune::GridGlue. \n+
Write all remote intersections to a vtk file for debugging. \n More...
\n-#include <type_traits>\n-
#include <dune/common/bartonnackmanifcheck.hh>\n-
#include <dune/common/parallel/communicator.hh>\n-
#include <dune/grid/common/datahandleif.hh>\n-
#include <dune/grid/common/gridenums.hh>\n+
#include <fstream>\n+
#include <iomanip>\n+
#include <type_traits>\n+
#include <vector>\n+
#include <dune/common/classname.hh>\n+
#include <dune/geometry/type.hh>\n+
#include <dune/geometry/referenceelements.hh>\n+
#include <dune/grid-glue/gridglue.hh>\n
Go to the source code of this file.
\n\n Classes | |
| struct | Dune::GridGlue::GlobalId |
| class | Dune::GridGlue::CommDataHandle< DataHandleImp, DataTypeImp > |
| describes the features of a data handle for communication in parallel runs using the GridGlue::communicate methods. More... | |
| class | Dune::GridGlue::StreamingMessageBuffer< DT > |
| class | Dune::GridGlue::CommunicationOperator< dir > |
| forward gather scatter to user defined CommInfo class More... | |
| struct | Dune::GridGlue::CommInfo< GG, DataHandleImp, DataTypeImp > |
| collects all GridGlue data requried for communication More... | |
| struct | Dune::CommPolicy< ::Dune::GridGlue::CommInfo< GG, DataHandleImp, DataTypeImp > > |
| specialization of the CommPolicy struct, required for the ParallelIndexsets More... | |
| class | Dune::GridGlue::GridGlueVtkWriter |
| Write remote intersections to a vtk file for debugging purposes. More... | |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
\n-Typedefs | |
| typedef std::pair< int, int > | Dune::GridGlue::RankPair |
| typedef CommunicationOperator< Dune::ForwardCommunication > | Dune::GridGlue::ForwardOperator |
| typedef CommunicationOperator< Dune::BackwardCommunication > | Dune::GridGlue::BackwardOperator |
\n-Functions | |
| std::ostream & | Dune::GridGlue::operator<< (std::ostream &os, const GlobalId &id) |
Detailed Description
\n-Describes the parallel communication interface class for Dune::GridGlue.
\n-- Author
- Christian Engwer
Write all remote intersections to a vtk file for debugging.
\n\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,62 +5,34 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces | Typedefs | Functions\n-gridgluecommunicate.hh File Reference\n-Describes the parallel communication interface class for Dune::GridGlue.\n-More...\n+Classes | Namespaces\n+gridgluevtkwriter.hh File Reference\n+Write all remote intersections to a vtk file for debugging. More...\n+#include \n- \n+ \n
\n gridgluecommunicate.hh
gridgluevtkwriter.hh
\n Go to the documentation of this file.
\n \n \n \n \n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n
3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 197 static const typename CommInfo::DataType& gather(const CommInfo& commInfo, size_t i, size_t j = 0)
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 227 static void scatter(CommInfo& commInfo, const typename CommInfo::DataType& v, std::size_t i, std::size_t j = 0)
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 210 fmerged << ((intersectionDim==3) ? "DATASET UNSTRUCTURED_GRID" : "DATASET POLYDATA") << std::endl;
\n+ 211 fmerged << "POINTS " << overlaps*(intersectionDim+1) << " " << Dune::className<ctype>() << std::endl;
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+Central component of the module implementing the coupling of two grids.
Definition: gridglue.hh:37
CommunicationOperator< Dune::BackwardCommunication > BackwardOperator
Definition: gridgluecommunicate.hh:264
CommunicationOperator< Dune::ForwardCommunication > ForwardOperator
Definition: gridgluecommunicate.hh:263
std::ostream & operator<<(std::ostream &os, const GlobalId &id)
Definition: gridgluecommunicate.hh:55
The intersection of two entities of the two patches of a GridGlue.
Definition: intersection.hh:261
Intersection< P0, P1, O, I > flip() const
Return a copy of the intersection with inside and outside switched.
Definition: intersection.hh:483
bool self() const
For parallel computations: Return true if inside() entity exists locally.
Definition: intersection.hh:393
InsideEntity inside(unsigned int parentId=0) const
Return element on the inside of this intersection.
Definition: intersection.hh:319
size_t neighbor(unsigned int g=0) const
Return number of embeddings into local grid0 (grid1) entities.
Definition: intersection.hh:399
OutsideEntity outside(unsigned int parentId=0) const
Return element on the outside of this intersection.
Definition: intersection.hh:328
Definition: gridgluecommunicate.hh:26
GlobalId(int i, int j, unsigned int n)
Definition: gridgluecommunicate.hh:48
describes the features of a data handle for communication in parallel runs using the GridGlue::commun...
Definition: gridgluecommunicate.hh:77
size_t size(RISType &i) const
Definition: gridgluecommunicate.hh:92
void scatter(MessageBufferImp &buff, const EntityType &e, const RISType &i, size_t n)
Definition: gridgluecommunicate.hh:118
void gather(MessageBufferImp &buff, const EntityType &e, const RISType &i) const
pack data from user to message buffer
Definition: gridgluecommunicate.hh:104
DataTypeImp DataType
data type of data to communicate
Definition: gridgluecommunicate.hh:80
Definition: gridgluecommunicate.hh:141
StreamingMessageBuffer(DT *p)
Definition: gridgluecommunicate.hh:146
size_t counter() const
Definition: gridgluecommunicate.hh:169
void write(const Y &data)
Definition: gridgluecommunicate.hh:155
void read(Y &data) const
Definition: gridgluecommunicate.hh:163
DT value_type
Definition: gridgluecommunicate.hh:143
forward gather scatter to user defined CommInfo class
Definition: gridgluecommunicate.hh:194
static void scatter(CommInfo &commInfo, const typename CommInfo::DataType &v, std::size_t i, std::size_t j=0)
Definition: gridgluecommunicate.hh:227
static const CommInfo::DataType & gather(const CommInfo &commInfo, size_t i, size_t j=0)
Definition: gridgluecommunicate.hh:197
collects all GridGlue data requried for communication
Definition: gridgluecommunicate.hh:272
::Dune::GridGlue::CommDataHandle< DataHandleImp, DataTypeImp > * data
Definition: gridgluecommunicate.hh:282
mutable ::Dune::GridGlue::StreamingMessageBuffer< DataType > mbuffer
Definition: gridgluecommunicate.hh:286
std::vector< DataType > buffer
Definition: gridgluecommunicate.hh:285
static size_t getSize(const Type &commInfo, size_t i)
Get the number of objects at an intersection.
Definition: gridgluecommunicate.hh:320
DataTypeImp IndexedType
The datatype that should be communicated.
Definition: gridgluecommunicate.hh:309
::Dune::GridGlue::CommInfo< GG, DataHandleImp, DataTypeImp > Type
The type of the GridGlueCommInfo.
Definition: gridgluecommunicate.hh:304
VariableSize IndexedTypeFlag
Each intersection can communicate a different number of objects.
Definition: gridgluecommunicate.hh:315
IteratorRange<... > intersections(const GridGlue<... > &glue, const Reverse<... > &reverse=!reversed)
Iterate over all intersections of a GridGlue.
Write remote intersections to a vtk file for debugging purposes.
Definition: gridgluevtkwriter.hh:41
static void write(const Glue &glue, const std::string &filenameTrunk)
Definition: gridgluevtkwriter.hh:276
Definition: rangegenerators.hh:17
Provides codimension-independent methods for grid extraction.
Definition: extractor.hh:46
void getCoords(std::vector< Dune::FieldVector< ctype, dimworld > > &coords) const
getter for the coordinates array
Definition: extractor.hh:275
\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,443 +5,350 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-gridgluecommunicate.hh\n+gridgluevtkwriter.hh\n Go_to_the_documentation_of_this_file.\n 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-\n 2// vi: set et ts=4 sw=2 sts=2:\n 3// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 5#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUECOMMUNICATE_HH\n- 6#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUECOMMUNICATE_HH\n- 7\n- 13#include \n \n+ \n
\n \n Classes |\n Namespaces
\n- gridglueamirawriter.hh File Reference
intersectionindexset.hh File Reference
\n-\n-
\n-
\n-
\n+\n\n
\n
\n-\n+
\n \n Write all remote intersections to a AmiraMesh file. \n-More...
\n-#include <fstream>\n-
#include <sstream>\n-
#include <type_traits>\n+\n
Go to the source code of this file.
\n\n Classes | |
| class | Dune::GridGlue::GridGlueAmiraWriter |
| Write remote intersections to a AmiraMesh file for debugging purposes. More... | |
| class | Dune::GridGlue::IntersectionIndexSet< P0, P1 > |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
Detailed Description
\n-Write all remote intersections to a AmiraMesh file.
\n-\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -6,28 +6,22 @@\n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n Classes | Namespaces\n-gridglueamirawriter.hh File Reference\n-Write all remote intersections to a AmiraMesh file. More...\n-#include \n- \n+ \n
\n gridglueamirawriter.hh
intersectionindexset.hh
\n-Go to the documentation of this file.\n+\n \n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n \n
\n Generated by 1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,166 +5,106 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-gridglueamirawriter.hh\n+intersectionindexset.hh\n Go_to_the_documentation_of_this_file.\n- 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n+1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 8#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n- 9#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n- 10\n- 11#include \n- 12#include \n- 13#include \n+ 3#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n+ 4#define DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n+ 5\n+ 6#include \n+ 7#include \n+ 8\n+ 9#ifndef ONLY_SIMPLEX_INTERSECTIONS\n+ 10// we currently support only one intersection type. If we want to support\n+more,\n+ 11// we have to think about the semantics of our IndexSet\n+ 12#error Not Implemented\n+ 13#endif\n 14\n 15namespace Dune {\n- 16namespace GridGlue {\n+ 16 namespace GridGlue {\n 17\n-20class GridGlueAmiraWriter\n- 21{\n- 22\n- 26 template \n- 27 static void writeIntersections(const Glue& glue, const std::string&\n-filename)\n- 28 {\n- 29 static_assert((side==0 || side==1), \"'side' can only be 0 or 1\");\n+ 18 template\n+19 class IntersectionIndexSet\n+ 20 {\n+ 21 friend class ::Dune::GridGlue::GridGlue;\n+ 22 typedef ::Dune::GridGlue::GridGlue GridGlue;\n+ 23\n+ 24 public:\n+ 25\n+27 typedef typename GridGlue::IndexType IndexType;\n+29 typedef size_t SizeType;\n 30\n- 31 std::ofstream fgrid;\n- 32\n- 33 fgrid.open(filename.c_str());\n- 34\n- 35 using GridView = typename Glue::template GridView;\n- 36 const int dim = GridView::dimension;\n- 37 const int domdimw = GridView::dimensionworld;\n+ 33 template\n+34 IndexType index (const Intersection & i) const\n+ 35 {\n+ 36 return i.i_->index_;\n+ 37 }\n 38\n- 39 // coordinates have to be in R^3 in the VTK format\n- 40 std::string coordinatePadding;\n- 41 for (int i=domdimw; i<3; i++)\n- 42 coordinatePadding += \" 0\";\n- 43\n- 44 int overlaps = glue.size();\n+41 SizeType size () const\n+ 42 {\n+ 43 return glue_->size();\n+ 44 }\n 45\n- 46 if (dim==3) {\n+ 46 private:\n 47\n- 48 fgrid << \"# HyperSurface 0.1 ASCII \\n\" << std::endl;\n- 49 fgrid<<\"\\n\";\n- 50 fgrid<<\"Parameters {\\n\";\n- 51 fgrid<<\" Materials {\\n\";\n- 52 fgrid<<\" outside {\\n\";\n- 53 fgrid<<\" Id 0\\n\";\n- 54 fgrid<<\" }\\n\";\n- 55 fgrid<<\" inside {\\n\";\n- 56 fgrid<<\" Id 1\\n\";\n- 57 fgrid<<\" }\\n\";\n- 58 fgrid<<\" }\\n\";\n- 59 fgrid<<\"\\n\";\n- 60 fgrid<<\"}\\n\";\n- 61\n- 62 // ////////////////////////////////////////////\n- 63 // Write vertices\n- 64 // ////////////////////////////////////////////\n- 65\n- 66 //use dim and not dim+1\n- 67 fgrid<<\"\\nVertices \"<< overlaps*(dim)<<\"\\n\";\n- 68 auto isEnd = glue.template iend();\n- 69 for (auto isIt = glue.template ibegin(); isIt != isEnd; ++isIt)\n- 70 {\n- 71 const auto& geometry = isIt->geometry();\n- 72 for (int i = 0; i < geometry.corners(); ++i)\n- 73 fgrid << geometry.corner(i) << coordinatePadding << std::endl;\n- 74 }\n- 75\n- 76 // ////////////////////////////////////////////\n- 77 // Write triangles\n- 78 // ////////////////////////////////////////////\n- 79\n- 80 fgrid<<\"NBranchingPoints 0\\n\";\n- 81 fgrid<<\"NVerticesOnCurves 0\\n\";\n- 82 fgrid<<\"BoundaryCurves 0\\n\";\n- 83 fgrid<<\"Patches 1\\n\";\n- 84 fgrid<<\"{\\n\";\n- 85 fgrid<<\"InnerRegion inside\\n\";\n- 86 fgrid<<\"OuterRegion outside\\n\";\n- 87 fgrid<<\"BoundaryID 0\\n\";\n- 88 fgrid<<\"BranchingPoints 0\";\n- 89 fgrid<<\"\\n\";\n- 90\n- 91 fgrid<<\"Triangles \"<();\n- 119 for (auto isIt = glue.template ibegin(); isIt != isEnd; ++isIt) {\n- 120 const auto& geometry = isIt->geometry();\n- 121 for (int i = 0; i <2; ++i)\n- 122 fgrid << geometry.corner(i) <<\" \"<<0<<\"\\n\";\n- 123 }\n- 124 }\n- 125\n- 126 fgrid.close();\n- 127 }\n- 128\n- 129public:\n- 130 template\n-131 static void write(const Glue& glue, const std::string& path, int\n-appendix=1)\n- 132 {\n- 133 std::ostringstream name0;\n- 134 name0 << path;\n- 135 name0 << \"/domain.surf\" << std::setw(3) << std::setfill('0') << appendix;\n- 136\n- 137 // Write extracted grid and remote intersection on the grid1-side\n- 138 writeIntersections(glue,name0.str());\n- 139\n- 140 std::ostringstream name1;\n- 141 name1 << path;\n- 142 name1 << \"/target.surf\" << std::setw(3) << std::setfill('0') << appendix;\n- 143\n- 144 writeIntersections(glue, name1.str());\n- 145 }\n- 146\n- 147};\n- 148\n- 149} // namespace GridGlue\n- 150} // namespace Dune\n- 151\n- 152#endif // DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n+ 49 IntersectionIndexSet(const GridGlue * g) :\n+ 50 glue_(g) {}\n+ 51\n+ 52 const GridGlue * glue_;\n+ 53 };\n+ 54\n+ 55 } // end namespace GridGlue\n+ 56} // end namespace Dune\n+ 57\n+ 58#endif // DUNE_GRIDGLUE_ADAPTER_INTERSECTIONINDEXSET_HH\n+gridglue.hh\n+Central component of the module implementing the coupling of two grids.\n+intersection.hh\n+Model of the Intersection concept provided by GridGlue.\n Dune\n Definition: gridglue.hh:37\n-Dune::GridGlue::GridGlueAmiraWriter\n-Write remote intersections to a AmiraMesh file for debugging purposes.\n-Definition: gridglueamirawriter.hh:21\n-Dune::GridGlue::GridGlueAmiraWriter::write\n-static void write(const Glue &glue, const std::string &path, int appendix=1)\n-Definition: gridglueamirawriter.hh:131\n+Dune::GridGlue::GridGlue\n+sequential adapter to couple two grids at specified close together boundaries\n+Definition: gridglue.hh:67\n+Dune::GridGlue::GridGlue::IndexType\n+unsigned int IndexType\n+Definition: gridglue.hh:147\n+Dune::GridGlue::GridGlue::size\n+size_t size() const\n+Definition: gridglue.hh:393\n+Dune::GridGlue::IntersectionData::index_\n+IndexType index_\n+index of this intersection after GridGlue interface\n+Definition: intersection.hh:112\n+Dune::GridGlue::Intersection\n+The intersection of two entities of the two patches of a GridGlue.\n+Definition: intersection.hh:261\n+Dune::GridGlue::IntersectionIndexSet\n+Definition: intersectionindexset.hh:20\n+Dune::GridGlue::IntersectionIndexSet::size\n+SizeType size() const\n+Return total number of intersections.\n+Definition: intersectionindexset.hh:41\n+Dune::GridGlue::IntersectionIndexSet::IndexType\n+GridGlue::IndexType IndexType\n+The type used for the indices.\n+Definition: intersectionindexset.hh:27\n+Dune::GridGlue::IntersectionIndexSet::index\n+IndexType index(const Intersection< P0, P1, I, O > &i) const\n+Map Dune::GridGlue::Intersection to index.\n+Definition: intersectionindexset.hh:34\n+Dune::GridGlue::IntersectionIndexSet::SizeType\n+size_t SizeType\n+The type used for the size.\n+Definition: intersectionindexset.hh:29\n \n ===============================================================================\n Generated by\u00a0[doxygen] 1.9.4\n"}]}, {"source1": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00026.html", "source2": "./usr/share/doc/libdune-grid-glue-doc/doxygen/a00026.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-dune-grid-glue: rangegenerators.hh File Reference \n+dune-grid-glue: gridglueamirawriter.hh File Reference \n \n \n \n \n \n \n \n@@ -64,42 +64,42 @@\n \n
\n \n \n+ 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n+Go to the documentation of this file. 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n \n-\n+\n \n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n \n-\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+Central component of the module implementing the coupling of two grids.
Model of the Intersection concept provided by GridGlue.
Definition: gridglue.hh:37
Write remote intersections to a AmiraMesh file for debugging purposes.
Definition: gridglueamirawriter.hh:21
static void write(const Glue &glue, const std::string &path, int appendix=1)
Definition: gridglueamirawriter.hh:131
sequential adapter to couple two grids at specified close together boundaries
Definition: gridglue.hh:67
IndexType index_
index of this intersection after GridGlue interface
Definition: intersection.hh:112
The intersection of two entities of the two patches of a GridGlue.
Definition: intersection.hh:261
Definition: intersectionindexset.hh:20
SizeType size() const
Return total number of intersections.
Definition: intersectionindexset.hh:41
GridGlue::IndexType IndexType
The type used for the indices.
Definition: intersectionindexset.hh:27
IndexType index(const Intersection< P0, P1, I, O > &i) const
Map Dune::GridGlue::Intersection to index.
Definition: intersectionindexset.hh:34
size_t SizeType
The type used for the size.
Definition: intersectionindexset.hh:29
\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,166 +5,106 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-gridglueamirawriter.hh\n+intersectionindexset.hh\n Go_to_the_documentation_of_this_file.\n- 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n+1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 8#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n- 9#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUEAMIRAWRITER_HH\n- 10\n- 11#include gridglueamirawriter.hh File Reference
\n-
\n+\n+
\n+
\n+
\n\n \n
\n
\n-
\n-\n+
\n \n #include <dune/common/iteratorrange.hh>\n+\n+
Write all remote intersections to a AmiraMesh file. \n+More...
\n+#include <fstream>\n+
#include <sstream>\n+
#include <type_traits>\n
Go to the source code of this file.
\n\n Classes | |
| struct | Dune::GridGlue::Reverse< reverse > |
| class | Dune::GridGlue::GridGlueAmiraWriter |
| Write remote intersections to a AmiraMesh file for debugging purposes. More... | |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
\n-Functions | |
| template<... > | |
| IteratorRange<... > | Dune::GridGlue::intersections (const GridGlue<... > &glue, const Reverse<... > &reverse=!reversed) |
| Iterate over all intersections of a GridGlue. More... | |
Detailed Description
\n+Write all remote intersections to a AmiraMesh file.
\n+\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,28 +5,29 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-Classes | Namespaces | Functions\n-rangegenerators.hh File Reference\n-#include \n- \n+ \n
\n rangegenerators.hh
gridglueamirawriter.hh
\n Go to the documentation of this file.\n-\n-\n-\n-\n-\n-\n-\n+\n+\n
\n \n 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
\n-\n-\n-\n-\n-\n-\n-\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n Definition: gridglue.hh:37
IteratorRange<... > intersections(const GridGlue<... > &glue, const Reverse<... > &reverse=!reversed)
Iterate over all intersections of a GridGlue.
sequential adapter to couple two grids at specified close together boundaries
Definition: gridglue.hh:67
IteratorRange<... > intersections(const GridGlue<... > &glue, const Reverse<... > &reverse=!reversed)
Iterate over all intersections of a GridGlue.
Definition: rangegenerators.hh:17
constexpr Reverse<!reverse > operator!() const
Definition: rangegenerators.hh:21
Write remote intersections to a AmiraMesh file for debugging purposes.
Definition: gridglueamirawriter.hh:21
static void write(const Glue &glue, const std::string &path, int appendix=1)
Definition: gridglueamirawriter.hh:131
\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,89 +5,166 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * adapter\n-rangegenerators.hh\n+gridglueamirawriter.hh\n Go_to_the_documentation_of_this_file.\n 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 3#ifndef DUNE_GRIDGLUE_ADAPTER_RANGEGENERATORS_HH\n- 4#define DUNE_GRIDGLUE_ADAPTER_RANGEGENERATORS_HH\n- 5\n- 6#include \n \n- \n+ \n
\n areawriter_impl.hh File Reference
areawriter.hh File Reference
\n-
\n-
\n-
\n-
\n-
\n+\n\n-
\n
\n
#include <fstream>\n-
#include <vector>\n-
#include <dune/common/fvector.hh>\n-
#include <dune/geometry/type.hh>\n-
#include <dune/grid/common/mcmgmapper.hh>\n+\n
Go to the source code of this file.
\n\n-Classes | |
| struct | Dune::GridGlue::AreaWriterImplementation::FacetLayout< dimgrid > |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
| namespace | Dune::GridGlue::AreaWriterImplementation |
\n Functions | |
| template<typename GridView > | |
| void | Dune::GridGlue::AreaWriterImplementation::write_facet_geometry (const GridView &gv, std::ostream &out) |
| template<int side, typename Glue > | |
| void | Dune::GridGlue::write_glue_area_vtk (const Glue &glue, std::ostream &out) |
| template<int side, typename Glue > | |
| void | Dune::GridGlue::write_glue_area_vtk (const Glue &glue, const std::string &filename) |
| template<typename Glue > |
\n+Classes | |
| class | Dune::GridGlue::Projection< Coordinate > |
| Projection of a line (triangle) on another line (triangle). More... | |
| struct | Dune::GridGlue::Projection< Coordinate >::EdgeIntersection |
| Intersection between two edges of a triangle. More... | |
\n Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
| namespace | Dune::GridGlue::ProjectionWriterImplementation |
\n-Functions | |
| template<unsigned side, typename Coordinate , typename Corners > | |
| void | Dune::GridGlue::ProjectionWriterImplementation::write_points (const Projection< Coordinate > &projection, const Corners &corners, std::ostream &out) |
| template<unsigned side, typename Coordinate , typename Normals > | |
| void | Dune::GridGlue::ProjectionWriterImplementation::write_normals (const Projection< Coordinate > &projection, const Normals &normals, std::ostream &out) |
| template<typename Coordinate , typename Corners > | |
| void | Dune::GridGlue::ProjectionWriterImplementation::write_edge_intersection_points (const Projection< Coordinate > &projection, const Corners &corners, std::ostream &out) |
| template<typename Coordinate , typename Normals > | |
| void | Dune::GridGlue::ProjectionWriterImplementation::write_edge_intersection_normals (const Projection< Coordinate > &projection, const Normals &normals, std::ostream &out) |
| template<unsigned side, typename Coordinate > | |
| void | Dune::GridGlue::ProjectionWriterImplementation::write_success (const Projection< Coordinate > &projection, std::ostream &out) |
| template<typename Coordinate , typename Corners , typename Normals > | |
| void | Dune::GridGlue::write (const Projection< Coordinate > &projection, const Corners &corners, const Normals &normals, std::ostream &out) |
| write projection in VTK format More... | |
| template<typename Coordinate , typename Corners , typename Normals > | |
| void | Dune::GridGlue::write (const Projection< Coordinate > &projection, const Corners &corners, const Normals &normals, const std::string &filename) |
| write projection in VTK format More... | |
| template<typename Coordinate , typename Corners , typename Normals > | |
| void | Dune::GridGlue::print (const Projection< Coordinate > &projection, const Corners &corners, const Normals &normals) |
| Print information about the projection to std::cout stream. More... | |
\n Generated by
1.9.4\n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,61 +5,29 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * common\n-Namespaces | Functions\n-projectionwriter_impl.hh File Reference\n-#include \n- \n+ \n
\n projectionwriter_impl.hh
projection.hh
\n Go to the documentation of this file.\n-\n-\n-\n-\n-\n-\n-\n-\n-\n \n+\n+\n \n+\n+\n+\n+\n+\n+\n \n \n \n+\n+\n+\n
\n \n 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file LICENSE.md in module root
\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n \n-\n- 11void write_points(const Projection<Coordinate>& projection, const Corners& corners, std::ostream& out)
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 27void write_normals(const Projection<Coordinate>& projection, const Normals& normals, std::ostream& out)
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 43void write_edge_intersection_points(const Projection<Coordinate>& projection, const Corners& corners, std::ostream& out)
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n- 56void write_edge_intersection_normals(const Projection<Coordinate>& projection, const Normals& normals, std::ostream& out)
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n \n-\n-\n-\n-\n-\n-\n-\n-\n- 104 out << "CELLS " << (8 + numberOfEdgeIntersections) << " " << (26 + 3 * numberOfEdgeIntersections) << "\\n";
\n-\n-\n-\n-\n- 109 out << "CELL_TYPES " << (8 + numberOfEdgeIntersections) << "\\n" "5\\n3\\n3\\n3\\n" "5\\n3\\n3\\n3\\n";
\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n-\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+ 134 void doProjection(const std::tuple<Corners&, Corners&>& corners, const std::tuple<Normals&, Normals&>& normals);
\n+\n+\n+ 145 void doInverseProjection(const std::tuple<Corners&, Corners&>& corners, const std::tuple<Normals&, Normals&>& normals);
\n+\n+\n+ 156 void doEdgeIntersection(const std::tuple<Corners&, Corners&>& corners, const std::tuple<Normals&, Normals&>& normals);
\n+\n+\n+ 184 inline bool projectionFeasible(const Coordinate& x, const Coordinate& nx, const Coordinate& px, const Corners& corners, const Normals& normals) const;
\n+\n+\n+\n+\n+\n+\n+\n+ 213 void project(const std::tuple<Corners&, Corners&>& corners, const std::tuple<Normals&, Normals&>& normals);
\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n+\n Definition: gridglue.hh:37
void write(const Projection< Coordinate > &projection, const Corners &corners, const Normals &normals, std::ostream &out)
write projection in VTK format
Definition: projectionwriter_impl.hh:86
void print(const Projection< Coordinate > &projection, const Corners &corners, const Normals &normals)
Print information about the projection to std::cout stream.
Definition: projectionwriter_impl.hh:140
Corners::value_type interpolate(const Coordinate &x, const Corners &corners)
Definition: projection_impl.hh:70
Normals::value_type interpolate_unit_normals(const Coordinate &x, const Normals &normals)
Definition: projection_impl.hh:91
void write_normals(const Projection< Coordinate > &projection, const Normals &normals, std::ostream &out)
Definition: projectionwriter_impl.hh:27
void write_points(const Projection< Coordinate > &projection, const Corners &corners, std::ostream &out)
Definition: projectionwriter_impl.hh:11
void write_success(const Projection< Coordinate > &projection, std::ostream &out)
Definition: projectionwriter_impl.hh:72
void write_edge_intersection_points(const Projection< Coordinate > &projection, const Corners &corners, std::ostream &out)
Definition: projectionwriter_impl.hh:43
void write_edge_intersection_normals(const Projection< Coordinate > &projection, const Normals &normals, std::ostream &out)
Definition: projectionwriter_impl.hh:56
Projection of a line (triangle) on another line (triangle).
Definition: projection.hh:21
Coordinate::field_type Field
Scalar type.
Definition: projection.hh:61
Projection(const Field overlap=Field(0), const Field max_normal_product=Field(-0.1))
Definition: projection_impl.hh:130
const std::tuple< std::bitset< dim >, std::bitset< dim > > & success() const
Indicate whether projection (inverse projection) is valid for each corner or not.
Definition: projection.hh:252
static constexpr unsigned maxEdgeIntersections
maximum number of edge-edge intersections
Definition: projection.hh:54
std::array< Coordinate, dim > Images
List of corner images.
Definition: projection.hh:70
void epsilon(const Field epsilon)
Set epsilon used for floating-point comparisons.
Definition: projection_impl.hh:140
void project(const std::tuple< Corners &, Corners & > &corners, const std::tuple< Normals &, Normals & > &normals)
Do the actual projection.
Definition: projection_impl.hh:469
static constexpr unsigned dim
dimension of coordinates
Definition: projection.hh:47
unsigned numberOfEdgeIntersections() const
Number of edge intersections.
Definition: projection.hh:262
const std::tuple< Images, Preimages > & images() const
Images and preimages of corners.
Definition: projection.hh:235
const std::array< EdgeIntersection, maxEdgeIntersections > & edgeIntersections() const
Edge-edge intersections.
Definition: projection.hh:273
Intersection between two edges of a triangle.
Definition: projection.hh:29
std::array< Coordinate, 2 > local
Local coordinates of intersection and distance along normals.
Definition: projection.hh:41
std::array< unsigned, 2 > edge
Edge numbers in image and preimage triangle.
Definition: projection.hh:33
\n Generated by
1.9.4\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -5,246 +5,180 @@\n \n dune-grid-glue\u00a02.9\n \n \n * dune\n * grid-glue\n * common\n-projectionwriter_impl.hh\n+projection.hh\n Go_to_the_documentation_of_this_file.\n 1// SPDX-FileCopyrightInfo: Copyright \u00a9 DUNE Project contributors, see file\n LICENSE.md in module root\n 2// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-\n with-dune-grid-glue-exception\n- 3#include