│ │ │ │ Classes
Model of the Intersection concept provided by GridGlue.
Implement iterators over GridGlue intersections.
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
│ │ │ │
│ │ │ │ -
│ │ │ │ +Macros
│ │ │ │ + gridglue.cc File Reference
intersection.hh File Reference
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ + │ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │ Generated by
1.9.8
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -1,28 +1,46 @@
│ │ │ │ │ dune-grid-glue 2.10
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _g_r_i_d_-_g_l_u_e
│ │ │ │ │ * _a_d_a_p_t_e_r
│ │ │ │ │ -_N_a_m_e_s_p_a_c_e_s | _F_u_n_c_t_i_o_n_s
│ │ │ │ │ -gridglue.cc File Reference
│ │ │ │ │ -#include "_i_n_t_e_r_s_e_c_t_i_o_n_._h_h"
│ │ │ │ │ -#include
│ │ │ │ │ -#include
│ │ │ │ │ -#include "_._._/_g_r_i_d_g_l_u_e_._h_h"
│ │ │ │ │ -#include "_._._/_c_o_m_m_o_n_/_r_i_n_g_c_o_m_m_._h_h"
│ │ │ │ │ -#include
│ │ │ │ │ +_C_l_a_s_s_e_s | _N_a_m_e_s_p_a_c_e_s | _M_a_c_r_o_s
│ │ │ │ │ +intersection.hh File Reference
│ │ │ │ │ +Model of the Intersection concept provided by GridGlue. _M_o_r_e_._._.
│ │ │ │ │ +#include
│ │ │ │ │ +#include
│ │ │ │ │ +#include
│ │ │ │ │ +#include
│ │ │ │ │ +#include
│ │ │ │ │ +#include
│ │ │ │ │ +#include
│ │ │ │ │ +#include <_d_u_n_e_/_g_r_i_d_-_g_l_u_e_/_g_r_i_d_g_l_u_e_._h_h>
│ │ │ │ │ _G_o_ _t_o_ _t_h_e_ _s_o_u_r_c_e_ _c_o_d_e_ _o_f_ _t_h_i_s_ _f_i_l_e_.
│ │ │ │ │ +CCllaasssseess
│ │ │ │ │ + class _D_u_n_e_:_:_G_r_i_d_G_l_u_e_:_:_I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a_<_ _P_0_,_ _P_1_ _>
│ │ │ │ │ + storage class for _D_u_n_e_:_:_G_r_i_d_G_l_u_e_:_:_I_n_t_e_r_s_e_c_t_i_o_n related data _M_o_r_e_._._.
│ │ │ │ │ +
│ │ │ │ │ +struct _D_u_n_e_:_:_G_r_i_d_G_l_u_e_:_:_I_n_t_e_r_s_e_c_t_i_o_n_T_r_a_i_t_s_<_ _P_0_,_ _P_1_,_ _i_n_s_i_d_e_,_ _o_u_t_s_i_d_e_ _>
│ │ │ │ │ +
│ │ │ │ │ + class _D_u_n_e_:_:_G_r_i_d_G_l_u_e_:_:_I_n_t_e_r_s_e_c_t_i_o_n_<_ _P_0_,_ _P_1_,_ _I_,_ _O_ _>
│ │ │ │ │ + The intersection of two entities of the two patches of a _G_r_i_d_G_l_u_e.
│ │ │ │ │ + _M_o_r_e_._._.
│ │ │ │ │ +
│ │ │ │ │ NNaammeessppaacceess
│ │ │ │ │ namespace _D_u_n_e
│ │ │ │ │
│ │ │ │ │ namespace _D_u_n_e_:_:_G_r_i_d_G_l_u_e
│ │ │ │ │
│ │ │ │ │ -FFuunnccttiioonnss
│ │ │ │ │ -template
│ │ │ │ │ -void _D_u_n_e_:_:_G_r_i_d_G_l_u_e_:_:_p_r_i_n_t_V_e_c_t_o_r (const std::vector< T > &v, std::string name,
│ │ │ │ │ - int rank)
│ │ │ │ │ +MMaaccrrooss
│ │ │ │ │ +#define _O_N_L_Y___S_I_M_P_L_E_X___I_N_T_E_R_S_E_C_T_I_O_N_S
│ │ │ │ │
│ │ │ │ │ +********** DDeettaaiilleedd DDeessccrriippttiioonn **********
│ │ │ │ │ +Model of the Intersection concept provided by GridGlue.
│ │ │ │ │ + Author
│ │ │ │ │ + Christian Engwer
│ │ │ │ │ +********** MMaaccrroo DDeeffiinniittiioonn DDooccuummeennttaattiioonn **********
│ │ │ │ │ +********** _?◆_? OONNLLYY__SSIIMMPPLLEEXX__IINNTTEERRSSEECCTTIIOONNSS **********
│ │ │ │ │ +#define ONLY_SIMPLEX_INTERSECTIONS
│ │ │ │ │ ===============================================================================
│ │ │ │ │ Generated by _[_d_o_x_y_g_e_n_] 1.9.8
│ │ │ ├── ./usr/share/doc/libdune-grid-glue-doc/doxygen/a00008_source.html
│ │ │ │ @@ -1,15 +1,15 @@
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ -dune-grid-glue: gridglue.cc Source File
│ │ │ │ +dune-grid-glue: intersection.hh Source File
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ @@ -74,502 +74,569 @@
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ #include "intersection.hh"│ │ │ │ -
#include <vector>│ │ │ │ -
#include <iterator>│ │ │ │ -
#include "../gridglue.hh"│ │ │ │ -
#include "../common/ringcomm.hh"│ │ │ │ -
#include <dune/common/unused.hh>│ │ │ │ + │ │ │ │ +
Model of the Intersection concept provided by GridGlue. │ │ │ │ +More...
│ │ │ │ +#include <algorithm>│ │ │ │ +
#include <optional>│ │ │ │ +
#include <tuple>│ │ │ │ +
#include <dune/common/deprecated.hh>│ │ │ │ +
#include <dune/common/version.hh>│ │ │ │ +
#include <dune/geometry/affinegeometry.hh>│ │ │ │ +
#include <dune/geometry/referenceelements.hh>│ │ │ │ +
#include <dune/grid-glue/gridglue.hh>│ │ │ │
Go to the source code of this file.
│ │ │ ││ │ │ │ +Classes | |
| 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... | |
│ │ │ │ Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
│ │ │ │ -Functions | |
| template<typename T > | |
| void | Dune::GridGlue::printVector (const std::vector< T > &v, std::string name, int rank) |
│ │ │ │ +Macros | |
| #define | ONLY_SIMPLEX_INTERSECTIONS |
Detailed Description
│ │ │ │ +Model of the Intersection concept provided by GridGlue.
│ │ │ │ + │ │ │ │ +Macro Definition Documentation
│ │ │ │ + │ │ │ │ +◆ ONLY_SIMPLEX_INTERSECTIONS
│ │ │ │ + │ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
| #define ONLY_SIMPLEX_INTERSECTIONS | │ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +│ │ │ │ Generated by
1.9.8
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -1,28 +1,46 @@
│ │ │ │ │ dune-grid-glue 2.10
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _g_r_i_d_-_g_l_u_e
│ │ │ │ │ * _a_d_a_p_t_e_r
│ │ │ │ │ -_N_a_m_e_s_p_a_c_e_s | _F_u_n_c_t_i_o_n_s
│ │ │ │ │ -gridglue.cc File Reference
│ │ │ │ │ -#include "_i_n_t_e_r_s_e_c_t_i_o_n_._h_h"
│ │ │ │ │ -#include
│ │ │ │ -
│ │ │ │ +
│ │ │ │
│ │ │ │ gridglue.cc
intersection.hh
│ │ │ │ Go to the documentation of this file.
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
3// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file LICENSE.md in module root
│ │ │ │ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 21GridGlue<P0, P1>::GridGlue(const std::shared_ptr< const GridPatch<0> >& gp0, const std::shared_ptr< const GridPatch<1> >& gp1, const std::shared_ptr<Merger>& merger)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 47 std::vector<IntersectionData> dummy(1); // we need size 1, as we always store data for the end-intersection
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ - 50 static constexpr int dim() { return GridGlue::template GridView<side>::Grid::dimension - GridGlue::template GridPatch<side>::codim; }
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ + 58 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimension>;
│ │ │ │ +
│ │ │ │ + 60 using Grid0LocalGeometry [[deprecated("please use GridLocalGeometry<0> instead")]] = GridLocalGeometry<0>;
│ │ │ │ + 61 using Grid1LocalGeometry [[deprecated("please use GridLocalGeometry<1> instead")]] = GridLocalGeometry<1>;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ + 65 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimensionworld>;
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ + 77 IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset, bool grid0local, bool grid1local);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ + 195 std::array<Dune::FieldVector<ctype, GridGlue::template GridView<side>::dimensionworld>, nSimplexCorners> corners_global;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ + 209 IntersectionData<P0, P1>::IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset,
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ + 318 typedef Dune::EnumRange <Dune::PartitionType, Dune::InteriorEntity, Dune::GhostEntity> AllFlags;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 345 DUNE_THROW(Dune::NotImplemented, "GridGlue::communicate for interface " << iftype << " not implemented");
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -Model of the Intersection concept provided by GridGlue.
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ + 431 normal[0] = (jacobianTransposed[0][1] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][1]);
│ │ │ │ + 432 normal[1] = - (jacobianTransposed[0][0] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][0]);
│ │ │ │ + 433 normal[2] = (jacobianTransposed[0][0] * jacobianTransposed[1][1] - jacobianTransposed[0][1] * jacobianTransposed[1][0]);
│ │ │ │ +
│ │ │ │ + 435 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the 'outerNormal' method for " << mydim << "-dimensional intersections yet");
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +Central component of the module implementing the coupling of two grids.
Definition gridglue.hh:37
IteratorRange<... > intersections(const GridGlue<... > &glue, const Reverse<... > &reverse=!reversed)
Iterate over all intersections of a GridGlue.
void printVector(const std::vector< T > &v, std::string name, int rank)
Definition gridglue.cc:168
void MPI_AllApply(MPI_Comm mpicomm, OP &&op, const Args &... data)
apply an operator locally to a difstributed data set
Definition ringcomm.hh:297
sequential adapter to couple two grids at specified close together boundaries
Definition gridglue.hh:67
void mergePatches(const std::vector< Dune::FieldVector< ctype, dimworld > > &patch0coords, const std::vector< unsigned int > &patch0entities, const std::vector< Dune::GeometryType > &patch0types, const int patch0rank, const std::vector< Dune::FieldVector< ctype, dimworld > > &patch1coords, const std::vector< unsigned int > &patch1entities, const std::vector< Dune::GeometryType > &patch1types, const int patch1rank)
after building the merged grid the intersection can be updated through this method (for internal use)
Definition gridglue.cc:179
void communicate(Dune::GridGlue::CommDataHandle< DataHandleImp, DataTypeImp > &data, Dune::InterfaceType iftype, Dune::CommunicationDirection dir) const
Communicate information on the MergedGrid of a GridGlue.
Definition gridglue.cc:300
void extractGrid(const Extractor &extractor, std::vector< Dune::FieldVector< ctype, dimworld > > &coords, std::vector< unsigned int > &faces, std::vector< Dune::GeometryType > &geometryTypes) const
Definition gridglue.cc:265
std::conditional_t< side==0, P0, std::conditional_t< side==1, P1, void > > GridPatch
Definition gridglue.hh:96
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
Definition gridgluecommunicate.hh:26
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
Definition gridgluecommunicate.hh:141
forward gather scatter to user defined CommInfo class
Definition gridgluecommunicate.hh:194
collects all GridGlue data required for communication
Definition gridgluecommunicate.hh:272
::Dune::GridGlue::CommDataHandle< DataHandleImp, DataTypeImp > * data
Definition gridgluecommunicate.hh:282
Provides codimension-independent methods for grid extraction.
Definition extractor.hh:46
void getFaces(std::vector< VertexVector > &faces) const
Get the corners of the extracted subentities.
Definition extractor.hh:285
void getGeometryTypes(std::vector< Dune::GeometryType > &geometryTypes) const
Get the list of geometry types.
Definition extractor.hh:274
void getCoords(std::vector< Dune::FieldVector< ctype, dimworld > > &coords) const
getter for the coordinates array
Definition extractor.hh:256
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:257
Traits::OutsideLocalGeometry OutsideLocalGeometry
Definition intersection.hh:271
Intersection< P0, P1, O, I > flip() const
Return a copy of the intersection with inside and outside switched.
Definition intersection.hh:474
bool conforming() const
Return true if intersection is conforming.
Definition intersection.hh:331
InsideGridView::Traits::template Codim< 0 >::Entity InsideEntity
Definition intersection.hh:277
int indexInOutside(unsigned int parentId=0) const
Local number of codim 1 entity in outside() Entity where intersection is contained in.
Definition intersection.hh:408
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:401
bool self() const
For parallel computations: Return true if inside() entity exists locally.
Definition intersection.hh:384
static constexpr auto mydim
dimension of the intersection
Definition intersection.hh:287
static constexpr int outsidePatch
outside patch
Definition intersection.hh:293
Dune::GeometryType type() const
Type of reference element for this intersection.
Definition intersection.hh:373
InsideEntity inside(unsigned int parentId=0) const
Return element on the inside of this intersection.
Definition intersection.hh:315
Traits::LocalCoordinate LocalCoordinate
Definition intersection.hh:280
const OutsideLocalGeometry & geometryInOutside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the outside() element.
Definition intersection.hh:345
GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const
Return a unit outer normal.
Definition intersection.hh:446
IntersectionTraits< P0, P1, I, O > Traits
Definition intersection.hh:261
Traits::IntersectionData IntersectionData
Definition intersection.hh:264
Traits::GlobalCoordinate GlobalCoordinate
Definition intersection.hh:281
Traits::OutsideGeometry OutsideGeometry
Definition intersection.hh:272
size_t neighbor(unsigned int g=0) const
Return number of embeddings into local grid0 (grid1) entities.
Definition intersection.hh:390
Intersection(const GridGlue *glue, const IntersectionData *i)
Constructor for a given Dataset.
Definition intersection.hh:307
static constexpr int insidePatch
inside patch
Definition intersection.hh:290
OutsideGridView::Traits::template Codim< 0 >::Entity OutsideEntity
Definition intersection.hh:278
OutsideEntity outside(unsigned int parentId=0) const
Return element on the outside of this intersection.
Definition intersection.hh:324
const Geometry & geometry() const
Geometric information about this intersection as part of the inside grid.
Definition intersection.hh:356
Traits::InsideLocalGeometry InsideLocalGeometry
Definition intersection.hh:268
GlobalCoordinate outerNormal(const LocalCoordinate &local) const
Return an outer normal (length not necessarily 1)
Definition intersection.hh:418
const OutsideGeometry & geometryOutside() const
Geometric information about this intersection as part of the outside grid.
Definition intersection.hh:367
Traits::OutsideGridView OutsideGridView
Definition intersection.hh:270
GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const
Return an outer normal with the length of the integration element.
Definition intersection.hh:457
Traits::InsideGridView InsideGridView
Definition intersection.hh:267
const InsideLocalGeometry & geometryInInside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the inside() element.
Definition intersection.hh:338
GlobalCoordinate centerUnitOuterNormal() const
Unit outer normal at the center of the intersection.
Definition intersection.hh:466
static constexpr auto coorddim
dimension of the world space of the intersection
Definition intersection.hh:284
Definition intersection.hh:230
static constexpr int insidePatch
Definition intersection.hh:245
Dune::FieldVector< ctype, mydim > LocalCoordinate
Definition intersection.hh:249
static constexpr auto coorddim
Definition intersection.hh:243
typename IntersectionData::template GridGeometry< outside > OutsideGeometry
Definition intersection.hh:241
typename GridGlue::ctype ctype
Definition intersection.hh:248
typename GridGlue::template GridView< outside > OutsideGridView
Definition intersection.hh:235
Dune::FieldVector< ctype, coorddim > GlobalCoordinate
Definition intersection.hh:250
typename IntersectionData::template GridLocalGeometry< outside > OutsideLocalGeometry
Definition intersection.hh:238
typename GridGlue::template GridView< inside > InsideGridView
Definition intersection.hh:234
typename IntersectionData::template GridGeometry< inside > Geometry
Definition intersection.hh:240
typename IntersectionData::template GridLocalGeometry< inside > InsideLocalGeometry
Definition intersection.hh:237
static constexpr int outsidePatch
Definition intersection.hh:246
│ │ │ │ Generated by
1.9.8
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -1,584 +1,743 @@
│ │ │ │ │ dune-grid-glue 2.10
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _g_r_i_d_-_g_l_u_e
│ │ │ │ │ * _a_d_a_p_t_e_r
│ │ │ │ │ -gridglue.cc
│ │ │ │ │ +intersection.hh
│ │ │ │ │ _G_o_ _t_o_ _t_h_e_ _d_o_c_u_m_e_n_t_a_t_i_o_n_ _o_f_ _t_h_i_s_ _f_i_l_e_.
│ │ │ │ │ 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
│ │ │ │ │ 2// vi: set et ts=4 sw=2 sts=2:
│ │ │ │ │ 3// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file
│ │ │ │ │ LICENSE.md in module root
│ │ │ │ │ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-
│ │ │ │ │ dune-grid-glue-exception
│ │ │ │ │ -5/* IMPLEMENTATION OF CLASS G R I D G L U E */
│ │ │ │ │ -6
│ │ │ │ │ -7#include "_i_n_t_e_r_s_e_c_t_i_o_n_._h_h"
│ │ │ │ │ -8#include -
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │ +Namespaces
│ │ │ │ + intersection.hh File Reference
gridgluevtkwriter.hh File Reference
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │
│ │ │ │ +
│ │ │ │
│ │ │ │ Model of the Intersection concept provided by GridGlue. │ │ │ │ +
Write all remote intersections to a vtk file for debugging. │ │ │ │ More...
│ │ │ │ -#include <algorithm>│ │ │ │ -
#include <optional>│ │ │ │ -
#include <tuple>│ │ │ │ -
#include <dune/common/deprecated.hh>│ │ │ │ -
#include <dune/common/version.hh>│ │ │ │ -
#include <dune/geometry/affinegeometry.hh>│ │ │ │ +
#include <fstream>│ │ │ │ +
#include <iomanip>│ │ │ │ +
#include <type_traits>│ │ │ │ +
#include <vector>│ │ │ │ +
#include <dune/common/classname.hh>│ │ │ │ +
#include <dune/geometry/type.hh>│ │ │ │
#include <dune/geometry/referenceelements.hh>│ │ │ │
#include <dune/grid-glue/gridglue.hh>│ │ │ │
Go to the source code of this file.
│ │ │ ││ │ │ │ Classes | |
| 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... | |
| class | Dune::GridGlue::GridGlueVtkWriter |
| Write remote intersections to a vtk file for debugging purposes. More... | |
│ │ │ │ Namespaces | |
| namespace | Dune |
| namespace | Dune::GridGlue |
│ │ │ │ -Macros | |
| #define | ONLY_SIMPLEX_INTERSECTIONS |
Detailed Description
│ │ │ │ -Model of the Intersection concept provided by GridGlue.
│ │ │ │ - │ │ │ │ -Macro Definition Documentation
│ │ │ │ - │ │ │ │ -◆ ONLY_SIMPLEX_INTERSECTIONS
│ │ │ │ - │ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
| #define ONLY_SIMPLEX_INTERSECTIONS | │ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -Write all remote intersections to a vtk file for debugging.
│ │ │ │ +│ │ │ │ Generated by
1.9.8
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -1,46 +1,33 @@
│ │ │ │ │ dune-grid-glue 2.10
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _g_r_i_d_-_g_l_u_e
│ │ │ │ │ * _a_d_a_p_t_e_r
│ │ │ │ │ -_C_l_a_s_s_e_s | _N_a_m_e_s_p_a_c_e_s | _M_a_c_r_o_s
│ │ │ │ │ -intersection.hh File Reference
│ │ │ │ │ -Model of the Intersection concept provided by GridGlue. _M_o_r_e_._._.
│ │ │ │ │ -#include
│ │ │ │ -
│ │ │ │ +
│ │ │ │
│ │ │ │ intersection.hh
gridgluevtkwriter.hh
│ │ │ │ Go to the documentation of this file.
│ │ │ │
│ │ │ │
│ │ │ │ Generated by 1.9.8
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -1,743 +1,349 @@
│ │ │ │ │ dune-grid-glue 2.10
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _g_r_i_d_-_g_l_u_e
│ │ │ │ │ * _a_d_a_p_t_e_r
│ │ │ │ │ -intersection.hh
│ │ │ │ │ +gridgluevtkwriter.hh
│ │ │ │ │ _G_o_ _t_o_ _t_h_e_ _d_o_c_u_m_e_n_t_a_t_i_o_n_ _o_f_ _t_h_i_s_ _f_i_l_e_.
│ │ │ │ │ 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
│ │ │ │ │ 2// vi: set et ts=4 sw=2 sts=2:
│ │ │ │ │ 3// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file
│ │ │ │ │ LICENSE.md in module root
│ │ │ │ │ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-
│ │ │ │ │ dune-grid-glue-exception
│ │ │ │ │ -11#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH
│ │ │ │ │ -12#define DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH
│ │ │ │ │ -13
│ │ │ │ │ -14#include
│ │ │ │ │ -15#include
│ │ │ │ │ -16#include
│ │ │ │ │ -17
│ │ │ │ │ -18#include
│ │ │ │ │ -19#include
│ │ │ │ │ -20#include
│ │ │ │ │ -21#include
│ │ │ │ │ -22#include <_d_u_n_e_/_g_r_i_d_-_g_l_u_e_/_g_r_i_d_g_l_u_e_._h_h>
│ │ │ │ │ +5/*
│ │ │ │ │ +6 * Filename: GridGlueVtkWriter.hh
│ │ │ │ │ +7 * Version: 1.0
│ │ │ │ │ +8 * Created on: Mar 5, 2009
│ │ │ │ │ +9 * Author: Gerrit Buse
│ │ │ │ │ +10 * ---------------------------------
│ │ │ │ │ +11 * Project: dune-grid-glue
│ │ │ │ │ +12 * Description: Class thought to make graphical debugging of couplings
│ │ │ │ │ +easier.
│ │ │ │ │ +13 *
│ │ │ │ │ +14 */
│ │ │ │ │ +20#ifndef DUNE_GRIDGLUE_ADAPTER_GRIDGLUEVTKWRITER_HH
│ │ │ │ │ +21#define DUNE_GRIDGLUE_ADAPTER_GRIDGLUEVTKWRITER_HH
│ │ │ │ │ +22
│ │ │ │ │ 23
│ │ │ │ │ -_2_4#define ONLY_SIMPLEX_INTERSECTIONS
│ │ │ │ │ -25
│ │ │ │ │ -26namespace _D_u_n_e {
│ │ │ │ │ -27 namespace GridGlue {
│ │ │ │ │ +24#include
│ │ │ │ │ +25#include
│ │ │ │ │ +26#include
│ │ │ │ │ +27#include
│ │ │ │ │ 28
│ │ │ │ │ -29 // forward declaration
│ │ │ │ │ -30 template
│ │ │ │ │ -31 class IntersectionIndexSet;
│ │ │ │ │ +29#include
│ │ │ │ │ +30#include
│ │ │ │ │ +31#include
│ │ │ │ │ 32
│ │ │ │ │ -36 template
│ │ │ │ │ -_3_7 class _I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a
│ │ │ │ │ -38 {
│ │ │ │ │ -39 public:
│ │ │ │ │ -_4_0 typedef ::Dune::GridGlue::GridGlue _G_r_i_d_G_l_u_e;
│ │ │ │ │ -41
│ │ │ │ │ -_4_2 typedef typename _G_r_i_d_G_l_u_e_:_:_I_n_d_e_x_T_y_p_e _I_n_d_e_x_T_y_p_e;
│ │ │ │ │ -43
│ │ │ │ │ -_4_5 static constexpr int _c_o_o_r_d_d_i_m = _G_r_i_d_G_l_u_e_:_:_d_i_m_w_o_r_l_d;
│ │ │ │ │ -46
│ │ │ │ │ -47 private:
│ │ │ │ │ -48 // intermediate quantities
│ │ │ │ │ -49 template
│ │ │ │ │ -50 static constexpr int dim() { return GridGlue::template GridView::
│ │ │ │ │ -Grid::dimension - GridGlue::template GridPatch::codim; }
│ │ │ │ │ -51
│ │ │ │ │ -52 public:
│ │ │ │ │ -_5_4 static constexpr int _m_y_d_i_m = dim<0>() < dim<1>() ? dim<0>() : dim<1>();
│ │ │ │ │ -55
│ │ │ │ │ -56 template
│ │ │ │ │ -_5_7 using _G_r_i_d_L_o_c_a_l_G_e_o_m_e_t_r_y = AffineGeometry<
│ │ │ │ │ -58 typename GridGlue::template GridView::ctype, _m_y_d_i_m, GridGlue::template
│ │ │ │ │ -GridView::dimension>;
│ │ │ │ │ +33#include <_d_u_n_e_/_g_r_i_d_-_g_l_u_e_/_g_r_i_d_g_l_u_e_._h_h>
│ │ │ │ │ +34
│ │ │ │ │ +35namespace _D_u_n_e {
│ │ │ │ │ +36namespace GridGlue {
│ │ │ │ │ +37
│ │ │ │ │ +_4_0class _G_r_i_d_G_l_u_e_V_t_k_W_r_i_t_e_r
│ │ │ │ │ +41{
│ │ │ │ │ +42
│ │ │ │ │ +46 template
│ │ │ │ │ +47 static void writeExtractedPart(const Glue& glue, const std::string&
│ │ │ │ │ +filename)
│ │ │ │ │ +48 {
│ │ │ │ │ +49 static_assert((side==0 || side==1), "'side' can only be 0 or 1");
│ │ │ │ │ +50
│ │ │ │ │ +51 std::ofstream fgrid;
│ │ │ │ │ +52
│ │ │ │ │ +53 fgrid.open(filename.c_str());
│ │ │ │ │ +54
│ │ │ │ │ +55 using GridView = typename Glue::template GridView;
│ │ │ │ │ +56 using _E_x_t_r_a_c_t_o_r = typename Glue::template GridPatch;
│ │ │ │ │ +57
│ │ │ │ │ +58 typedef typename GridView::ctype ctype;
│ │ │ │ │ 59
│ │ │ │ │ -_6_0 using _G_r_i_d_0_L_o_c_a_l_G_e_o_m_e_t_r_y [[deprecated("please use GridLocalGeometry<0>
│ │ │ │ │ -instead")]] = _G_r_i_d_L_o_c_a_l_G_e_o_m_e_t_r_y_<_0_>;
│ │ │ │ │ -_6_1 using _G_r_i_d_1_L_o_c_a_l_G_e_o_m_e_t_r_y [[deprecated("please use GridLocalGeometry<1>
│ │ │ │ │ -instead")]] = _G_r_i_d_L_o_c_a_l_G_e_o_m_e_t_r_y_<_1_>;
│ │ │ │ │ +60 const int domdimw = GridView::dimensionworld;
│ │ │ │ │ +61 const int patchDim = _E_x_t_r_a_c_t_o_r_:_:_d_i_m - _E_x_t_r_a_c_t_o_r_:_:_c_o_d_i_m;
│ │ │ │ │ 62
│ │ │ │ │ -63 template
│ │ │ │ │ -_6_4 using _G_r_i_d_G_e_o_m_e_t_r_y = AffineGeometry<
│ │ │ │ │ -65 typename GridGlue::template GridView::ctype, _m_y_d_i_m, GridGlue::template
│ │ │ │ │ -GridView::dimensionworld>;
│ │ │ │ │ -66
│ │ │ │ │ -_6_7 using _G_r_i_d_0_G_e_o_m_e_t_r_y [[deprecated("please use GridGeometry<0> instead")]] =
│ │ │ │ │ -_G_r_i_d_G_e_o_m_e_t_r_y_<_0_>;
│ │ │ │ │ -_6_8 using _G_r_i_d_1_G_e_o_m_e_t_r_y [[deprecated("please use GridGeometry<1> instead")]] =
│ │ │ │ │ -_G_r_i_d_G_e_o_m_e_t_r_y_<_1_>;
│ │ │ │ │ +63 // coordinates have to be in R^3 in the VTK format
│ │ │ │ │ +64 std::string coordinatePadding;
│ │ │ │ │ +65 for (int i=domdimw; i<3; i++)
│ │ │ │ │ +66 coordinatePadding += " 0";
│ │ │ │ │ +67
│ │ │ │ │ +68 fgrid << "# vtk DataFile Version 2.0\nFilename: " << filename << "\nASCII"
│ │ │ │ │ +<< std::endl;
│ │ │ │ │ 69
│ │ │ │ │ -70 template
│ │ │ │ │ -_7_1 using _G_r_i_d_I_n_d_e_x_T_y_p_e = typename GridGlue::template GridView::IndexSet::
│ │ │ │ │ -IndexType;
│ │ │ │ │ -72
│ │ │ │ │ -_7_3 using _G_r_i_d_0_I_n_d_e_x_T_y_p_e [[deprecated("please use GridIndexType<0> instead")]] =
│ │ │ │ │ -_G_r_i_d_I_n_d_e_x_T_y_p_e_<_0_>;
│ │ │ │ │ -_7_4 using _G_r_i_d_1_I_n_d_e_x_T_y_p_e [[deprecated("please use GridIndexType<1> instead")]] =
│ │ │ │ │ -_G_r_i_d_I_n_d_e_x_T_y_p_e_<_1_>;
│ │ │ │ │ -75
│ │ │ │ │ -77 _I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a(const _G_r_i_d_G_l_u_e& glue, unsigned int mergeindex, unsigned int
│ │ │ │ │ -offset, bool grid0local, bool grid1local);
│ │ │ │ │ -78
│ │ │ │ │ -_8_0 _I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a() = default;
│ │ │ │ │ -81
│ │ │ │ │ -82 /* Accessor Functions */
│ │ │ │ │ -83
│ │ │ │ │ -84 template
│ │ │ │ │ -_8_5 const _G_r_i_d_L_o_c_a_l_G_e_o_m_e_t_r_y_<_s_i_d_e_>& _l_o_c_a_l_G_e_o_m_e_t_r_y(unsigned int parentId = 0)
│ │ │ │ │ -const
│ │ │ │ │ -86 { return *std::get(sideData_).gridlocalgeom[parentId]; }
│ │ │ │ │ -87
│ │ │ │ │ -88 template
│ │ │ │ │ -_8_9 const _G_r_i_d_G_e_o_m_e_t_r_y_<_s_i_d_e_>& _g_e_o_m_e_t_r_y() const
│ │ │ │ │ -90 { return *std::get(sideData_).gridgeom; }
│ │ │ │ │ -91
│ │ │ │ │ -92 template
│ │ │ │ │ -_9_3 bool _l_o_c_a_l() const
│ │ │ │ │ -94 { return std::get(sideData_).gridlocal; }
│ │ │ │ │ -95
│ │ │ │ │ -96 template
│ │ │ │ │ -_9_7 _I_n_d_e_x_T_y_p_e _i_n_d_e_x(unsigned int parentId = 0) const
│ │ │ │ │ -98 { return std::get(sideData_).gridindices[parentId]; }
│ │ │ │ │ +70 // WRITE POINTS
│ │ │ │ │ +71 // ----------------
│ │ │ │ │ +72 std::vector coords;
│ │ │ │ │ +73 glue.template patch()._g_e_t_C_o_o_r_d_s(coords);
│ │ │ │ │ +74
│ │ │ │ │ +75 fgrid << ((patchDim==3) ? "DATASET UNSTRUCTURED_GRID" : "DATASET POLYDATA")
│ │ │ │ │ +<< std::endl;
│ │ │ │ │ +76 fgrid << "POINTS " << coords.size() << " " << Dune::className() <<
│ │ │ │ │ +std::endl;
│ │ │ │ │ +77
│ │ │ │ │ +78 for (size_t i=0; i faces;
│ │ │ │ │ +87 std::vector geometryTypes;
│ │ │ │ │ +88 glue.template patch().getFaces(faces);
│ │ │ │ │ +89 glue.template patch().getGeometryTypes(geometryTypes);
│ │ │ │ │ +90
│ │ │ │ │ +91 unsigned int faceCornerCount = 0;
│ │ │ │ │ +92 for (size_t i=0; i
│ │ │ │ │ -_1_0_1 _I_n_d_e_x_T_y_p_e _p_a_r_e_n_t_s() const
│ │ │ │ │ -102 { return std::get(sideData_).gridindices.size(); }
│ │ │ │ │ -103
│ │ │ │ │ -104 private:
│ │ │ │ │ -105 template
│ │ │ │ │ -106 void initializeGeometry(const _G_r_i_d_G_l_u_e& glue, unsigned mergeindex);
│ │ │ │ │ +100 fgrid << faces[i].size();
│ │ │ │ │ +101
│ │ │ │ │ +102 // vtk expects the vertices to by cyclically ordered
│ │ │ │ │ +103 // therefore unfortunately we have to deal with several element types on a
│ │ │ │ │ +case-by-case basis
│ │ │ │ │ +104 if (geometryTypes[i].isSimplex()) {
│ │ │ │ │ +105 for (int j=0; j
│ │ │ │ │ -116 struct SideData {
│ │ │ │ │ -118 bool gridlocal = false;
│ │ │ │ │ +108 } else if (geometryTypes[i].isQuadrilateral()) {
│ │ │ │ │ +109 fgrid << " " << faces[i][0] << " " << faces[i][1]
│ │ │ │ │ +110 << " " << faces[i][3] << " " << faces[i][2];
│ │ │ │ │ +111
│ │ │ │ │ +112 } else if (geometryTypes[i].isPyramid()) {
│ │ │ │ │ +113 fgrid << " " << faces[i][0] << " " << faces[i][1]
│ │ │ │ │ +114 << " " << faces[i][3] << " " << faces[i][2] << " " << faces[i][4];
│ │ │ │ │ +115
│ │ │ │ │ +116 } else if (geometryTypes[i].isPrism()) {
│ │ │ │ │ +117 fgrid << " " << faces[i][0] << " " << faces[i][2] << " " << faces[i][1]
│ │ │ │ │ +118 << " " << faces[i][3] << " " << faces[i][5] << " " << faces[i][4];
│ │ │ │ │ 119
│ │ │ │ │ -121 std::vector< GridIndexType > gridindices;
│ │ │ │ │ -122
│ │ │ │ │ -124 std::vector< std::optional< GridLocalGeometry > > gridlocalgeom;
│ │ │ │ │ +120 } else if (geometryTypes[i].isHexahedron()) {
│ │ │ │ │ +121 fgrid << " " << faces[i][0] << " " << faces[i][1]
│ │ │ │ │ +122 << " " << faces[i][3] << " " << faces[i][2]
│ │ │ │ │ +123 << " " << faces[i][4] << " " << faces[i][5]
│ │ │ │ │ +124 << " " << faces[i][7] << " " << faces[i][6];
│ │ │ │ │ 125
│ │ │ │ │ -133 std::optional< GridGeometry > gridgeom;
│ │ │ │ │ -134 };
│ │ │ │ │ -135
│ │ │ │ │ -136 std::tuple< SideData<0>, SideData<1> > sideData_;
│ │ │ │ │ -137 };
│ │ │ │ │ -138
│ │ │ │ │ -139 template
│ │ │ │ │ -140 template
│ │ │ │ │ -141 void IntersectionData::initializeGeometry(const GridGlue& glue,
│ │ │ │ │ -unsigned mergeindex)
│ │ │ │ │ -142 {
│ │ │ │ │ -143 auto& data = std::get(sideData_);
│ │ │ │ │ -144
│ │ │ │ │ -145 const unsigned n_parents = glue.merger_->template parents
│ │ │ │ │ -(mergeindex);
│ │ │ │ │ -146
│ │ │ │ │ -147 // init containers
│ │ │ │ │ -148 data.gridindices.resize(n_parents);
│ │ │ │ │ -149 data.gridlocalgeom.resize(n_parents);
│ │ │ │ │ -150
│ │ │ │ │ -151 // default values
│ │ │ │ │ -152 data.gridindices[0] = 0;
│ │ │ │ │ +126 } else {
│ │ │ │ │ +127 DUNE_THROW(Dune::NotImplemented, "Geometry type " << geometryTypes[i] << "
│ │ │ │ │ +not supported yet");
│ │ │ │ │ +128 }
│ │ │ │ │ +129
│ │ │ │ │ +130 fgrid << std::endl;
│ │ │ │ │ +131 }
│ │ │ │ │ +132
│ │ │ │ │ +133 fgrid << std::endl;
│ │ │ │ │ +134
│ │ │ │ │ +135 // 3d VTK files need an extra section specifying the CELL_TYPES aka
│ │ │ │ │ +GeometryTypes
│ │ │ │ │ +136 if (patchDim==3) {
│ │ │ │ │ +137
│ │ │ │ │ +138 fgrid << "CELL_TYPES " << geometryTypes.size() << std::endl;
│ │ │ │ │ +139
│ │ │ │ │ +140 for (size_t i=0; i::
│ │ │ │ │ -template Codim<0>::Geometry::mydimension;
│ │ │ │ │ -161
│ │ │ │ │ -162 // coordinates within the subentity that contains the remote intersection
│ │ │ │ │ -163 std::array() >, nSimplexCorners>
│ │ │ │ │ -corners_subEntity_local;
│ │ │ │ │ +154 }
│ │ │ │ │ +155
│ │ │ │ │ +156#if 0
│ │ │ │ │ +157 // WRITE CELL DATA
│ │ │ │ │ +158 // ---------------
│ │ │ │ │ +159 ctype accum = 0.0, delta = 1.0 / (ctype) (gridSubEntityData.size()-1);
│ │ │ │ │ +160
│ │ │ │ │ +161 fgrid << "CELL_DATA " << gridSubEntityData.size() << std::endl;
│ │ │ │ │ +162 fgrid << "SCALARS property_coding " << Dune::className() << " 1" <<
│ │ │ │ │ +std::endl;
│ │ │ │ │ +163 fgrid << "LOOKUP_TABLE default" << std::endl;
│ │ │ │ │ 164
│ │ │ │ │ -165 for (unsigned int par = 0; par < n_parents; ++par) {
│ │ │ │ │ -166 for (int i = 0; i < nSimplexCorners; ++i)
│ │ │ │ │ -167 corners_subEntity_local[i] = glue.merger_->template parentLocal
│ │ │ │ │ -(mergeindex, i, par);
│ │ │ │ │ -168
│ │ │ │ │ -169 // Coordinates of the remote intersection corners wrt the element
│ │ │ │ │ -coordinate system
│ │ │ │ │ -170 std::array, nSimplexCorners>
│ │ │ │ │ -corners_element_local;
│ │ │ │ │ -171
│ │ │ │ │ -172 if (data.gridlocal) {
│ │ │ │ │ -173 data.gridindices[par] = glue.merger_->template parent
│ │ │ │ │ -(mergeindex,par);
│ │ │ │ │ -174
│ │ │ │ │ -175 typename GridGlue::template GridPatch::LocalGeometry
│ │ │ │ │ -176 gridLocalGeometry = glue.template patch().geometryLocal
│ │ │ │ │ -(data.gridindices[par]);
│ │ │ │ │ -177 for (std::size_t i=0; i
│ │ │ │ │ +181 static void writeIntersections(const Glue& glue, const std::string&
│ │ │ │ │ +filename)
│ │ │ │ │ +182 {
│ │ │ │ │ +183 static_assert((side==0 || side==1), "'side' can only be 0 or 1");
│ │ │ │ │ +184
│ │ │ │ │ +185 std::ofstream fmerged;
│ │ │ │ │ +186
│ │ │ │ │ +187 fmerged.open(filename.c_str());
│ │ │ │ │ 188
│ │ │ │ │ -189 // Add world geometry only for 0th parent
│ │ │ │ │ -190 if (par == 0) {
│ │ │ │ │ -191 typename GridGlue::template GridPatch::Geometry
│ │ │ │ │ -192 gridWorldGeometry = glue.template patch().geometry(data.gridindices
│ │ │ │ │ -[par]);
│ │ │ │ │ -193
│ │ │ │ │ -194 // world coordinates of the remote intersection corners
│ │ │ │ │ -195 std::array::
│ │ │ │ │ -dimensionworld>, nSimplexCorners> corners_global;
│ │ │ │ │ -196
│ │ │ │ │ -197 for (std::size_t i=0; i
│ │ │ │ │ -_2_0_9 _I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a_<_P_0_,_ _P_1_>_:_:_I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a(const _G_r_i_d_G_l_u_e& glue, unsigned
│ │ │ │ │ -int mergeindex, unsigned int offset,
│ │ │ │ │ -210 bool grid0local, bool grid1local)
│ │ │ │ │ -211 : index_(mergeindex+offset)
│ │ │ │ │ -212 {
│ │ │ │ │ -213 // if an invalid index is given do not proceed!
│ │ │ │ │ -214 // (happens when the parent GridGlue initializes the "end"-Intersection)
│ │ │ │ │ -215 assert (0 <= mergeindex || mergeindex < glue.index__sz);
│ │ │ │ │ -216
│ │ │ │ │ -217 std::get<0>(sideData_).gridlocal = grid0local;
│ │ │ │ │ -218 std::get<1>(sideData_).gridlocal = grid1local;
│ │ │ │ │ +189 using GridView = typename Glue::template GridView;
│ │ │ │ │ +190 typedef typename GridView::ctype ctype;
│ │ │ │ │ +191
│ │ │ │ │ +192 const int domdimw = GridView::dimensionworld;
│ │ │ │ │ +193 const int intersectionDim = Glue::Intersection::mydim;
│ │ │ │ │ +194
│ │ │ │ │ +195 // coordinates have to be in R^3 in the VTK format
│ │ │ │ │ +196 std::string coordinatePadding;
│ │ │ │ │ +197 for (int i=domdimw; i<3; i++)
│ │ │ │ │ +198 coordinatePadding += " 0";
│ │ │ │ │ +199
│ │ │ │ │ +200 int overlaps = glue.size();
│ │ │ │ │ +201
│ │ │ │ │ +202 // WRITE POINTS
│ │ │ │ │ +203 // ----------------
│ │ │ │ │ +204 using _E_x_t_r_a_c_t_o_r = typename Glue::template GridPatch<0>;
│ │ │ │ │ +205 std::vector coords;
│ │ │ │ │ +206 glue.template patch()._g_e_t_C_o_o_r_d_s(coords);
│ │ │ │ │ +207
│ │ │ │ │ +208 // the merged grid (i.e. the set of remote intersections
│ │ │ │ │ +209 fmerged << "# vtk DataFile Version 2.0\nFilename: " << filename <<
│ │ │ │ │ +"\nASCII" << std::endl;
│ │ │ │ │ +210 fmerged << ((intersectionDim==3) ? "DATASET UNSTRUCTURED_GRID" : "DATASET
│ │ │ │ │ +POLYDATA") << std::endl;
│ │ │ │ │ +211 fmerged << "POINTS " << overlaps*(intersectionDim+1) << " " << Dune::
│ │ │ │ │ +className() << std::endl;
│ │ │ │ │ +212
│ │ │ │ │ +213 for (const auto& intersection : _i_n_t_e_r_s_e_c_t_i_o_n_s(glue, _R_e_v_e_r_s_e_<_s_i_d_e_ _=_=_ _1_>{}))
│ │ │ │ │ +214 {
│ │ │ │ │ +215 const auto& geometry = intersection.geometry();
│ │ │ │ │ +216 for (int i = 0; i < geometry.corners(); ++i)
│ │ │ │ │ +217 fmerged << geometry.corner(i) << coordinatePadding << std::endl;
│ │ │ │ │ +218 }
│ │ │ │ │ 219
│ │ │ │ │ -220 initializeGeometry<0>(glue, mergeindex);
│ │ │ │ │ -221 initializeGeometry<1>(glue, mergeindex);
│ │ │ │ │ -222 }
│ │ │ │ │ -223
│ │ │ │ │ -228 template
│ │ │ │ │ -_2_2_9 struct _I_n_t_e_r_s_e_c_t_i_o_n_T_r_a_i_t_s
│ │ │ │ │ -230 {
│ │ │ │ │ -_2_3_1 using _G_r_i_d_G_l_u_e = _:_:_D_u_n_e_:_:_G_r_i_d_G_l_u_e_:_:_G_r_i_d_G_l_u_e_<_P_0_,_ _P_1_>;
│ │ │ │ │ -_2_3_2 using _I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a = _D_u_n_e_:_:_G_r_i_d_G_l_u_e_:_:_I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a_<_P_0_,_ _P_1_>;
│ │ │ │ │ -233
│ │ │ │ │ -_2_3_4 using _I_n_s_i_d_e_G_r_i_d_V_i_e_w = typename GridGlue::template GridView;
│ │ │ │ │ -_2_3_5 using _O_u_t_s_i_d_e_G_r_i_d_V_i_e_w = typename GridGlue::template GridView;
│ │ │ │ │ -236
│ │ │ │ │ -_2_3_7 using _I_n_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y = typename IntersectionData::template
│ │ │ │ │ -GridLocalGeometry;
│ │ │ │ │ -_2_3_8 using _O_u_t_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y = typename IntersectionData::template
│ │ │ │ │ -GridLocalGeometry;
│ │ │ │ │ -239
│ │ │ │ │ -_2_4_0 using _G_e_o_m_e_t_r_y = typename IntersectionData::template GridGeometry;
│ │ │ │ │ -_2_4_1 using _O_u_t_s_i_d_e_G_e_o_m_e_t_r_y = typename IntersectionData::template
│ │ │ │ │ -GridGeometry;
│ │ │ │ │ +220 // WRITE POLYGONS
│ │ │ │ │ +221 // ----------------
│ │ │ │ │ +222
│ │ │ │ │ +223 std::vector faces;
│ │ │ │ │ +224 std::vector geometryTypes;
│ │ │ │ │ +225 glue.template patch().getFaces(faces);
│ │ │ │ │ +226 glue.template patch().getGeometryTypes(geometryTypes);
│ │ │ │ │ +227
│ │ │ │ │ +228 unsigned int faceCornerCount = 0;
│ │ │ │ │ +229 for (size_t i=0; i;
│ │ │ │ │ -_2_5_0 using _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e = Dune::FieldVector;
│ │ │ │ │ -251 };
│ │ │ │ │ +248 for (int i = 0; i < overlaps; i++)
│ │ │ │ │ +249 fmerged << "10" << std::endl;
│ │ │ │ │ +250
│ │ │ │ │ +251 }
│ │ │ │ │ 252
│ │ │ │ │ -255 template
│ │ │ │ │ -_2_5_6 class _I_n_t_e_r_s_e_c_t_i_o_n
│ │ │ │ │ -257 {
│ │ │ │ │ -258
│ │ │ │ │ -259 public:
│ │ │ │ │ -260
│ │ │ │ │ -_2_6_1 typedef _I_n_t_e_r_s_e_c_t_i_o_n_T_r_a_i_t_s_<_P_0_,_P_1_,_I_,_O_> _T_r_a_i_t_s;
│ │ │ │ │ -262
│ │ │ │ │ -_2_6_3 typedef typename _T_r_a_i_t_s_:_:_G_r_i_d_G_l_u_e _G_r_i_d_G_l_u_e;
│ │ │ │ │ -_2_6_4 typedef typename _T_r_a_i_t_s_:_:_I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a _I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a;
│ │ │ │ │ -265
│ │ │ │ │ -266
│ │ │ │ │ -_2_6_7 typedef typename _T_r_a_i_t_s_:_:_I_n_s_i_d_e_G_r_i_d_V_i_e_w _I_n_s_i_d_e_G_r_i_d_V_i_e_w;
│ │ │ │ │ -_2_6_8 typedef typename _T_r_a_i_t_s_:_:_I_n_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y _I_n_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y;
│ │ │ │ │ -269
│ │ │ │ │ -_2_7_0 typedef typename _T_r_a_i_t_s_:_:_O_u_t_s_i_d_e_G_r_i_d_V_i_e_w _O_u_t_s_i_d_e_G_r_i_d_V_i_e_w;
│ │ │ │ │ -_2_7_1 typedef typename _T_r_a_i_t_s_:_:_O_u_t_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y _O_u_t_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y;
│ │ │ │ │ -_2_7_2 typedef typename _T_r_a_i_t_s_:_:_O_u_t_s_i_d_e_G_e_o_m_e_t_r_y _O_u_t_s_i_d_e_G_e_o_m_e_t_r_y;
│ │ │ │ │ +253#if 0
│ │ │ │ │ +254 // WRITE CELL DATA
│ │ │ │ │ +255 // ---------------
│ │ │ │ │ +256 ctype accum = 0.0, delta = 1.0 / (ctype) (gridSubEntityData.size()-1);
│ │ │ │ │ +257
│ │ │ │ │ +258 fmerged << "CELL_DATA " << overlaps << std::endl;
│ │ │ │ │ +259 fmerged << "SCALARS property_coding " << Dune::className() << " 1"
│ │ │ │ │ +<< std::endl;
│ │ │ │ │ +260 fmerged << "LOOKUP_TABLE default" << std::endl;
│ │ │ │ │ +261
│ │ │ │ │ +262 for (typename GridSubEntityData::const_iterator sEIt =
│ │ │ │ │ +gridSubEntityData.begin();
│ │ │ │ │ +263 sEIt != gridSubEntityData.end();
│ │ │ │ │ +264 ++sEIt, accum += delta)
│ │ │ │ │ +265 {
│ │ │ │ │ +266 // ...and mark all of its merged grid parts with the same color
│ │ │ │ │ +267 for (int j = 0; j < sEIt->first.second; ++j)
│ │ │ │ │ +268 fmerged << accum << std::endl;
│ │ │ │ │ +269 }
│ │ │ │ │ +270#endif
│ │ │ │ │ +271 fmerged.close();
│ │ │ │ │ +272 }
│ │ │ │ │ 273
│ │ │ │ │ -_2_7_4 typedef typename _T_r_a_i_t_s_:_:_G_e_o_m_e_t_r_y _G_e_o_m_e_t_r_y;
│ │ │ │ │ -_2_7_5 typedef typename _T_r_a_i_t_s_:_:_c_t_y_p_e _c_t_y_p_e;
│ │ │ │ │ -276
│ │ │ │ │ -_2_7_7 typedef typename InsideGridView::Traits::template Codim<0>::Entity
│ │ │ │ │ -_I_n_s_i_d_e_E_n_t_i_t_y;
│ │ │ │ │ -_2_7_8 typedef typename OutsideGridView::Traits::template Codim<0>::Entity
│ │ │ │ │ -_O_u_t_s_i_d_e_E_n_t_i_t_y;
│ │ │ │ │ -279
│ │ │ │ │ -_2_8_0 typedef typename _T_r_a_i_t_s_:_:_L_o_c_a_l_C_o_o_r_d_i_n_a_t_e _L_o_c_a_l_C_o_o_r_d_i_n_a_t_e;
│ │ │ │ │ -_2_8_1 typedef typename _T_r_a_i_t_s_:_:_G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e;
│ │ │ │ │ +274public:
│ │ │ │ │ +275 template
│ │ │ │ │ +_2_7_6 static void _w_r_i_t_e(const Glue& glue, const std::string& filenameTrunk)
│ │ │ │ │ +277 {
│ │ │ │ │ +278
│ │ │ │ │ +279 // Write extracted grid and remote intersection on the grid0-side
│ │ │ │ │ +280 writeExtractedPart(glue,
│ │ │ │ │ +281 filenameTrunk + "-grid0.vtk");
│ │ │ │ │ 282
│ │ │ │ │ -_2_8_4 static constexpr auto _c_o_o_r_d_d_i_m = _T_r_a_i_t_s_:_:_c_o_o_r_d_d_i_m;
│ │ │ │ │ +283 writeIntersections(glue,
│ │ │ │ │ +284 filenameTrunk + "-intersections-grid0.vtk");
│ │ │ │ │ 285
│ │ │ │ │ -_2_8_7 static constexpr auto _m_y_d_i_m = _T_r_a_i_t_s_:_:_m_y_d_i_m;
│ │ │ │ │ -288
│ │ │ │ │ -_2_9_0 static constexpr int _i_n_s_i_d_e_P_a_t_c_h = _T_r_a_i_t_s_:_:_i_n_s_i_d_e_P_a_t_c_h;
│ │ │ │ │ -291
│ │ │ │ │ -_2_9_3 static constexpr int _o_u_t_s_i_d_e_P_a_t_c_h = _T_r_a_i_t_s_:_:_o_u_t_s_i_d_e_P_a_t_c_h;
│ │ │ │ │ +286 // Write extracted grid and remote intersection on the grid1-side
│ │ │ │ │ +287 writeExtractedPart(glue,
│ │ │ │ │ +288 filenameTrunk + "-grid1.vtk");
│ │ │ │ │ +289
│ │ │ │ │ +290 writeIntersections(glue,
│ │ │ │ │ +291 filenameTrunk + "-intersections-grid1.vtk");
│ │ │ │ │ +292
│ │ │ │ │ +293 }
│ │ │ │ │ 294
│ │ │ │ │ -295 // typedef unsigned int IndexType;
│ │ │ │ │ +295};
│ │ │ │ │ 296
│ │ │ │ │ -297 private:
│ │ │ │ │ -301 static constexpr int codimensionWorld = _c_o_o_r_d_d_i_m - _m_y_d_i_m;
│ │ │ │ │ -302
│ │ │ │ │ -303 public:
│ │ │ │ │ -304 /* C O N S T R U C T O R S */
│ │ │ │ │ -305
│ │ │ │ │ -_3_0_7 _I_n_t_e_r_s_e_c_t_i_o_n(const _G_r_i_d_G_l_u_e* glue, const _I_n_t_e_r_s_e_c_t_i_o_n_D_a_t_a* i) :
│ │ │ │ │ -308 glue_(glue), i_(i) {}
│ │ │ │ │ -309
│ │ │ │ │ -310 /* F U N C T I O N A L I T Y */
│ │ │ │ │ -311
│ │ │ │ │ -314 _I_n_s_i_d_e_E_n_t_i_t_y
│ │ │ │ │ -_3_1_5 _i_n_s_i_d_e(unsigned int parentId = 0) const
│ │ │ │ │ -316 {
│ │ │ │ │ -317 assert(_s_e_l_f());
│ │ │ │ │ -318 return glue_->template patch().element(i_->template index(parentId));
│ │ │ │ │ -319 }
│ │ │ │ │ -320
│ │ │ │ │ -323 _O_u_t_s_i_d_e_E_n_t_i_t_y
│ │ │ │ │ -_3_2_4 _o_u_t_s_i_d_e(unsigned int parentId = 0) const
│ │ │ │ │ -325 {
│ │ │ │ │ -326 assert(_n_e_i_g_h_b_o_r());
│ │ │ │ │ -327 return glue_->template patch().element(i_->template index(parentId));
│ │ │ │ │ -328 }
│ │ │ │ │ -329
│ │ │ │ │ -_3_3_1 bool _c_o_n_f_o_r_m_i_n_g() const
│ │ │ │ │ -332 {
│ │ │ │ │ -333 throw Dune::NotImplemented();
│ │ │ │ │ -334 }
│ │ │ │ │ -335
│ │ │ │ │ -_3_3_8 const _I_n_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y& _g_e_o_m_e_t_r_y_I_n_I_n_s_i_d_e(unsigned int parentId = 0)
│ │ │ │ │ -const
│ │ │ │ │ -339 {
│ │ │ │ │ -340 return i_->template localGeometry(parentId);
│ │ │ │ │ -341 }
│ │ │ │ │ -342
│ │ │ │ │ -_3_4_5 const _O_u_t_s_i_d_e_L_o_c_a_l_G_e_o_m_e_t_r_y& _g_e_o_m_e_t_r_y_I_n_O_u_t_s_i_d_e(unsigned int parentId = 0)
│ │ │ │ │ -const
│ │ │ │ │ -346 {
│ │ │ │ │ -347 return i_->template localGeometry(parentId);
│ │ │ │ │ -348 }
│ │ │ │ │ -349
│ │ │ │ │ -_3_5_6 const _G_e_o_m_e_t_r_y& _g_e_o_m_e_t_r_y() const
│ │ │ │ │ -357 {
│ │ │ │ │ -358 return i_->template geometry();
│ │ │ │ │ -359 }
│ │ │ │ │ -360
│ │ │ │ │ -_3_6_7 const _O_u_t_s_i_d_e_G_e_o_m_e_t_r_y& _g_e_o_m_e_t_r_y_O_u_t_s_i_d_e() const // DUNE_DEPRECATED
│ │ │ │ │ -368 {
│ │ │ │ │ -369 return i_->template geometry();
│ │ │ │ │ -370 }
│ │ │ │ │ -371
│ │ │ │ │ -_3_7_3 Dune::GeometryType _t_y_p_e() const
│ │ │ │ │ -374 {
│ │ │ │ │ -375 #ifdef ONLY_SIMPLEX_INTERSECTIONS
│ │ │ │ │ -376 return Dune::GeometryTypes::simplex(_m_y_d_i_m);
│ │ │ │ │ -377 #else
│ │ │ │ │ -378 #error Not Implemented
│ │ │ │ │ -379 #endif
│ │ │ │ │ -380 }
│ │ │ │ │ -381
│ │ │ │ │ -382
│ │ │ │ │ -_3_8_4 bool _s_e_l_f() const
│ │ │ │ │ -385 {
│ │ │ │ │ -386 return i_->template local();
│ │ │ │ │ -387 }
│ │ │ │ │ -388
│ │ │ │ │ -_3_9_0 size_t _n_e_i_g_h_b_o_r(unsigned int g = 0) const
│ │ │ │ │ -391 {
│ │ │ │ │ -392 if (g == 0 && i_->template local()) {
│ │ │ │ │ -393 return i_->template parents();
│ │ │ │ │ -394 } else if (g == 1 && i_->template local()) {
│ │ │ │ │ -395 return i_->template parents();
│ │ │ │ │ -396 }
│ │ │ │ │ -397 return 0;
│ │ │ │ │ -398 }
│ │ │ │ │ -399
│ │ │ │ │ -_4_0_1 int _i_n_d_e_x_I_n_I_n_s_i_d_e(unsigned int parentId = 0) const
│ │ │ │ │ -402 {
│ │ │ │ │ -403 assert(_s_e_l_f());
│ │ │ │ │ -404 return glue_->template patch().indexInInside(i_->template index
│ │ │ │ │ -(parentId));
│ │ │ │ │ -405 }
│ │ │ │ │ -406
│ │ │ │ │ -_4_0_8 int _i_n_d_e_x_I_n_O_u_t_s_i_d_e(unsigned int parentId = 0) const
│ │ │ │ │ -409 {
│ │ │ │ │ -410 assert(_n_e_i_g_h_b_o_r());
│ │ │ │ │ -411 return glue_->template patch().indexInInside(i_->template index
│ │ │ │ │ -(parentId));
│ │ │ │ │ -412 }
│ │ │ │ │ -413
│ │ │ │ │ -_4_1_8 _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e _o_u_t_e_r_N_o_r_m_a_l(const _L_o_c_a_l_C_o_o_r_d_i_n_a_t_e &local) const
│ │ │ │ │ -419 {
│ │ │ │ │ -420 _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e normal;
│ │ │ │ │ -421
│ │ │ │ │ -422 if (codimensionWorld == 0)
│ │ │ │ │ -423 DUNE_THROW(Dune::Exception, "There is no normal vector to a full-
│ │ │ │ │ -dimensional intersection");
│ │ │ │ │ -424 else if (codimensionWorld == 1) {
│ │ │ │ │ -425 /* TODO: Implement the general n-ary cross product here */
│ │ │ │ │ -426 const auto jacobianTransposed = _g_e_o_m_e_t_r_y().jacobianTransposed(local);
│ │ │ │ │ -427 if (_m_y_d_i_m==1) {
│ │ │ │ │ -428 normal[0] = - jacobianTransposed[0][1];
│ │ │ │ │ -429 normal[1] = jacobianTransposed[0][0];
│ │ │ │ │ -430 } else if (_m_y_d_i_m==2) {
│ │ │ │ │ -431 normal[0] = (jacobianTransposed[0][1] * jacobianTransposed[1][2] -
│ │ │ │ │ -jacobianTransposed[0][2] * jacobianTransposed[1][1]);
│ │ │ │ │ -432 normal[1] = - (jacobianTransposed[0][0] * jacobianTransposed[1][2] -
│ │ │ │ │ -jacobianTransposed[0][2] * jacobianTransposed[1][0]);
│ │ │ │ │ -433 normal[2] = (jacobianTransposed[0][0] * jacobianTransposed[1][1] -
│ │ │ │ │ -jacobianTransposed[0][1] * jacobianTransposed[1][0]);
│ │ │ │ │ -434 } else
│ │ │ │ │ -435 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the
│ │ │ │ │ -'outerNormal' method for " << _m_y_d_i_m << "-dimensional intersections yet");
│ │ │ │ │ -436 } else
│ │ │ │ │ -437 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the
│ │ │ │ │ -'outerNormal' method for intersections with codim >= 2 yet");
│ │ │ │ │ -438
│ │ │ │ │ -439 return normal;
│ │ │ │ │ -440 }
│ │ │ │ │ -441
│ │ │ │ │ -_4_4_6 _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e _u_n_i_t_O_u_t_e_r_N_o_r_m_a_l(const _L_o_c_a_l_C_o_o_r_d_i_n_a_t_e &local) const
│ │ │ │ │ -447 {
│ │ │ │ │ -448 _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e normal = _o_u_t_e_r_N_o_r_m_a_l(local);
│ │ │ │ │ -449 normal /= normal.two_norm();
│ │ │ │ │ -450 return normal;
│ │ │ │ │ -451 }
│ │ │ │ │ -452
│ │ │ │ │ -_4_5_7 _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e _i_n_t_e_g_r_a_t_i_o_n_O_u_t_e_r_N_o_r_m_a_l(const _L_o_c_a_l_C_o_o_r_d_i_n_a_t_e &local) const
│ │ │ │ │ -458 {
│ │ │ │ │ -459 return (_u_n_i_t_O_u_t_e_r_N_o_r_m_a_l(local) *= _g_e_o_m_e_t_r_y().integrationElement(local));
│ │ │ │ │ -460 }
│ │ │ │ │ -461
│ │ │ │ │ -_4_6_6 _G_l_o_b_a_l_C_o_o_r_d_i_n_a_t_e _c_e_n_t_e_r_U_n_i_t_O_u_t_e_r_N_o_r_m_a_l () const
│ │ │ │ │ -467 {
│ │ │ │ │ -468 return _u_n_i_t_O_u_t_e_r_N_o_r_m_a_l(ReferenceElements::general(_t_y_p_e
│ │ │ │ │ -()).position(0,0));
│ │ │ │ │ -469 }
│ │ │ │ │ -470
│ │ │ │ │ -_4_7_4 _I_n_t_e_r_s_e_c_t_i_o_n_<_P_0_,_P_1_,_O_,_I_> _f_l_i_p() const
│ │ │ │ │ -475 {
│ │ │ │ │ -476 return _I_n_t_e_r_s_e_c_t_i_o_n_<_P_0_,_P_1_,_O_,_I_>(glue_,i_);
│ │ │ │ │ -477 }
│ │ │ │ │ -478
│ │ │ │ │ -479#ifdef QUICKHACK_INDEX
│ │ │ │ │ -480 typedef typename _G_r_i_d_G_l_u_e_:_:_I_n_d_e_x_T_y_p_e IndexType;
│ │ │ │ │ -481
│ │ │ │ │ -482 IndexType index() const
│ │ │ │ │ -483 {
│ │ │ │ │ -484 return i_->_i_n_d_e_x__;
│ │ │ │ │ -485 }
│ │ │ │ │ -486
│ │ │ │ │ -487#endif
│ │ │ │ │ -488
│ │ │ │ │ -489 private:
│ │ │ │ │ -490
│ │ │ │ │ -491 friend class IntersectionIndexSet
3// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file LICENSE.md in module root
│ │ │ │ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │
│ │ │ │ 50 static constexpr int dim() { return GridGlue::template GridView<side>::Grid::dimension - GridGlue::template GridPatch<side>::codim; }
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 58 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimension>;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 60 using Grid0LocalGeometry [[deprecated("please use GridLocalGeometry<0> instead")]] = GridLocalGeometry<0>;
│ │ │ │ - 61 using Grid1LocalGeometry [[deprecated("please use GridLocalGeometry<1> instead")]] = GridLocalGeometry<1>;
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 65 typename GridGlue::template GridView<side>::ctype, mydim, GridGlue::template GridView<side>::dimensionworld>;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 77 IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset, bool grid0local, bool grid1local);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 160 static constexpr int elementdim = GridGlue::template GridView<side>::template Codim<0>::Geometry::mydimension;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 195 std::array<Dune::FieldVector<ctype, GridGlue::template GridView<side>::dimensionworld>, nSimplexCorners> corners_global;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 209 IntersectionData<P0, P1>::IntersectionData(const GridGlue& glue, unsigned int mergeindex, unsigned int offset,
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ + 210 fmerged << ((intersectionDim==3) ? "DATASET UNSTRUCTURED_GRID" : "DATASET POLYDATA") << std::endl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ - 431 normal[0] = (jacobianTransposed[0][1] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][1]);
│ │ │ │ - 432 normal[1] = - (jacobianTransposed[0][0] * jacobianTransposed[1][2] - jacobianTransposed[0][2] * jacobianTransposed[1][0]);
│ │ │ │ - 433 normal[2] = (jacobianTransposed[0][0] * jacobianTransposed[1][1] - jacobianTransposed[0][1] * jacobianTransposed[1][0]);
│ │ │ │ -
│ │ │ │ - 435 DUNE_THROW(Dune::NotImplemented, "Remote intersections don't implement the 'outerNormal' method for " << mydim << "-dimensional intersections yet");
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ Central component of the module implementing the coupling of two grids.
Definition gridglue.hh:37
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:257
Traits::OutsideLocalGeometry OutsideLocalGeometry
Definition intersection.hh:271
Intersection< P0, P1, O, I > flip() const
Return a copy of the intersection with inside and outside switched.
Definition intersection.hh:474
bool conforming() const
Return true if intersection is conforming.
Definition intersection.hh:331
InsideGridView::Traits::template Codim< 0 >::Entity InsideEntity
Definition intersection.hh:277
int indexInOutside(unsigned int parentId=0) const
Local number of codim 1 entity in outside() Entity where intersection is contained in.
Definition intersection.hh:408
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:401
bool self() const
For parallel computations: Return true if inside() entity exists locally.
Definition intersection.hh:384
static constexpr auto mydim
dimension of the intersection
Definition intersection.hh:287
static constexpr int outsidePatch
outside patch
Definition intersection.hh:293
Dune::GeometryType type() const
Type of reference element for this intersection.
Definition intersection.hh:373
InsideEntity inside(unsigned int parentId=0) const
Return element on the inside of this intersection.
Definition intersection.hh:315
Traits::LocalCoordinate LocalCoordinate
Definition intersection.hh:280
const OutsideLocalGeometry & geometryInOutside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the outside() element.
Definition intersection.hh:345
GlobalCoordinate unitOuterNormal(const LocalCoordinate &local) const
Return a unit outer normal.
Definition intersection.hh:446
IntersectionTraits< P0, P1, I, O > Traits
Definition intersection.hh:261
Traits::IntersectionData IntersectionData
Definition intersection.hh:264
Traits::GlobalCoordinate GlobalCoordinate
Definition intersection.hh:281
Traits::OutsideGeometry OutsideGeometry
Definition intersection.hh:272
size_t neighbor(unsigned int g=0) const
Return number of embeddings into local grid0 (grid1) entities.
Definition intersection.hh:390
Intersection(const GridGlue *glue, const IntersectionData *i)
Constructor for a given Dataset.
Definition intersection.hh:307
static constexpr int insidePatch
inside patch
Definition intersection.hh:290
OutsideGridView::Traits::template Codim< 0 >::Entity OutsideEntity
Definition intersection.hh:278
OutsideEntity outside(unsigned int parentId=0) const
Return element on the outside of this intersection.
Definition intersection.hh:324
const Geometry & geometry() const
Geometric information about this intersection as part of the inside grid.
Definition intersection.hh:356
Traits::InsideLocalGeometry InsideLocalGeometry
Definition intersection.hh:268
GlobalCoordinate outerNormal(const LocalCoordinate &local) const
Return an outer normal (length not necessarily 1)
Definition intersection.hh:418
const OutsideGeometry & geometryOutside() const
Geometric information about this intersection as part of the outside grid.
Definition intersection.hh:367
Traits::OutsideGridView OutsideGridView
Definition intersection.hh:270
GlobalCoordinate integrationOuterNormal(const LocalCoordinate &local) const
Return an outer normal with the length of the integration element.
Definition intersection.hh:457
Traits::InsideGridView InsideGridView
Definition intersection.hh:267
const InsideLocalGeometry & geometryInInside(unsigned int parentId=0) const
Geometric information about this intersection in local coordinates of the inside() element.
Definition intersection.hh:338
GlobalCoordinate centerUnitOuterNormal() const
Unit outer normal at the center of the intersection.
Definition intersection.hh:466
static constexpr auto coorddim
dimension of the world space of the intersection
Definition intersection.hh:284
Definition intersection.hh:230
static constexpr int insidePatch
Definition intersection.hh:245
Dune::FieldVector< ctype, mydim > LocalCoordinate
Definition intersection.hh:249
static constexpr auto coorddim
Definition intersection.hh:243
typename IntersectionData::template GridGeometry< outside > OutsideGeometry
Definition intersection.hh:241
typename GridGlue::ctype ctype
Definition intersection.hh:248
typename GridGlue::template GridView< outside > OutsideGridView
Definition intersection.hh:235
Dune::FieldVector< ctype, coorddim > GlobalCoordinate
Definition intersection.hh:250
typename IntersectionData::template GridLocalGeometry< outside > OutsideLocalGeometry
Definition intersection.hh:238
typename GridGlue::template GridView< inside > InsideGridView
Definition intersection.hh:234
typename IntersectionData::template GridGeometry< inside > Geometry
Definition intersection.hh:240
typename IntersectionData::template GridLocalGeometry< inside > InsideLocalGeometry
Definition intersection.hh:237
static constexpr int outsidePatch
Definition intersection.hh:246
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:256
│ │ │ │ Generated by
1.9.8
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -1,743 +1,349 @@
│ │ │ │ │ dune-grid-glue 2.10
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _g_r_i_d_-_g_l_u_e
│ │ │ │ │ * _a_d_a_p_t_e_r
│ │ │ │ │ -intersection.hh
│ │ │ │ │ +gridgluevtkwriter.hh
│ │ │ │ │ _G_o_ _t_o_ _t_h_e_ _d_o_c_u_m_e_n_t_a_t_i_o_n_ _o_f_ _t_h_i_s_ _f_i_l_e_.
│ │ │ │ │ 1// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
│ │ │ │ │ 2// vi: set et ts=4 sw=2 sts=2:
│ │ │ │ │ 3// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file
│ │ │ │ │ LICENSE.md in module root
│ │ │ │ │ 4// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-
│ │ │ │ │ dune-grid-glue-exception
│ │ │ │ │ -11#ifndef DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH
│ │ │ │ │ -12#define DUNE_GRIDGLUE_ADAPTER_INTERSECTION_HH
│ │ │ │ │ -13
│ │ │ │ │ -14#include