3.3.0
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
boxgeometryhelper.hh
Go to the documentation of this file.
1// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2// vi: set et ts=4 sw=4 sts=4:
3/*****************************************************************************
4 * See the file COPYING for full copying permissions. *
5 * *
6 * This program is free software: you can redistribute it and/or modify *
7 * it under the terms of the GNU General Public License as published by *
8 * the Free Software Foundation, either version 3 of the License, or *
9 * (at your option) any later version. *
10 * *
11 * This program is distributed in the hope that it will be useful, *
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
14 * GNU General Public License for more details. *
15 * *
16 * You should have received a copy of the GNU General Public License *
17 * along with this program. If not, see <http://www.gnu.org/licenses/>. *
18 *****************************************************************************/
25#ifndef DUMUX_DISCRETIZATION_BOX_GEOMETRY_HELPER_HH
26#define DUMUX_DISCRETIZATION_BOX_GEOMETRY_HELPER_HH
27
28#include <array>
29
30#include <dumux/common/math.hh>
31
32namespace Dumux {
33
35template<class GridView, int dim, class ScvType, class ScvfType>
36class BoxGeometryHelper;
37
39template <class GridView, class ScvType, class ScvfType>
40class BoxGeometryHelper<GridView, 1, ScvType, ScvfType>
41{
42private:
43 using Scalar = typename GridView::ctype;
44 using GlobalPosition = typename Dune::FieldVector<Scalar, GridView::dimensionworld>;
45 using ScvCornerStorage = typename ScvType::Traits::CornerStorage;
46 using ScvfCornerStorage = typename ScvfType::Traits::CornerStorage;
47
48 using Element = typename GridView::template Codim<0>::Entity;
49 using Intersection = typename GridView::Intersection;
50
53 static constexpr int maxPoints = 3;
54
55public:
56
57 BoxGeometryHelper(const typename Element::Geometry& geometry)
58 : elementGeometry_(geometry), corners_(geometry.corners()),
59 p_({geometry.center(), geometry.corner(0), geometry.corner(1)}) {}
60
62 ScvCornerStorage getScvCorners(unsigned int localScvIdx) const
63 {
65 static const std::uint8_t map[2][2] =
66 {
67 {1, 0},
68 {2, 0}
69 };
70
71 return ScvCornerStorage{ {p_[map[localScvIdx][0]],
72 p_[map[localScvIdx][1]]} };
73 }
74
76 ScvfCornerStorage getScvfCorners(unsigned int localScvfIdx) const
77 {
78 return ScvfCornerStorage{{p_[0]}};
79 }
80
82 ScvfCornerStorage getBoundaryScvfCorners(const Intersection& is,
83 const typename Intersection::Geometry& geometry,
84 unsigned int indexInIntersection) const
85 {
86 return ScvfCornerStorage{{geometry.corner(0)}};
87 }
88
90 GlobalPosition normal(const ScvfCornerStorage& scvfCorners,
91 const std::vector<unsigned int>& scvIndices) const
92 {
93 auto normal = p_[2] - p_[1];
94 normal /= normal.two_norm();
95 return normal;
96 }
97
99 Scalar scvVolume(const ScvCornerStorage& scvCorners) const
100 {
101 return (scvCorners[1] - scvCorners[0]).two_norm();
102 }
103
105 Scalar scvfArea(const ScvfCornerStorage& scvfCorners) const
106 {
107 return 1.0;
108 }
109
110protected:
111 const typename Element::Geometry& elementGeometry_;
112 std::size_t corners_; // number of element corners
113 std::array<GlobalPosition, maxPoints> p_; // the points needed for construction of the geometries
114};
115
117template <class GridView, class ScvType, class ScvfType>
118class BoxGeometryHelper<GridView, 2, ScvType, ScvfType>
119{
120 using Scalar = typename GridView::ctype;
121 using GlobalPosition = typename Dune::FieldVector<Scalar, GridView::dimensionworld>;
122 using ScvCornerStorage = typename ScvType::Traits::CornerStorage;
123 using ScvfCornerStorage = typename ScvfType::Traits::CornerStorage;
124
125 using Element = typename GridView::template Codim<0>::Entity;
126 using Intersection = typename GridView::Intersection;
127
128 static constexpr auto dim = GridView::dimension;
129 static constexpr auto dimWorld = GridView::dimensionworld;
130
133 static constexpr int maxPoints = 9;
134
135public:
136
137 BoxGeometryHelper(const typename Element::Geometry& geometry)
138 : elementGeometry_(geometry)
139 , corners_(geometry.corners())
140 {
141 const auto refElement = referenceElement(geometry);
142
143 // the element center
144 p_[0] = geometry.center();
145
146 // vertices
147 for (int i = 0; i < corners_; ++i)
148 p_[i+1] = geometry.corner(i);
149
150 // face midpoints
151 for (int i = 0; i < refElement.size(1); ++i)
152 p_[i+corners_+1] = geometry.global(refElement.position(i, 1));
153 }
154
156 ScvCornerStorage getScvCorners(unsigned int localScvIdx) const
157 {
158 // proceed according to number of corners of the element
159 switch (corners_)
160 {
161 case 3: // triangle
162 {
164 static const std::uint8_t vo = 1;
165 static const std::uint8_t fo = 4;
166 static const std::uint8_t map[3][4] =
167 {
168 {vo+0, fo+0, fo+1, 0},
169 {vo+1, fo+2, fo+0, 0},
170 {vo+2, fo+1, fo+2, 0}
171 };
172
173 return ScvCornerStorage{ {p_[map[localScvIdx][0]],
174 p_[map[localScvIdx][1]],
175 p_[map[localScvIdx][2]],
176 p_[map[localScvIdx][3]]} };
177 }
178 case 4: // quadrilateral
179 {
181 static const std::uint8_t vo = 1;
182 static const std::uint8_t fo = 5;
183 static const std::uint8_t map[4][4] =
184 {
185 {vo+0, fo+2, fo+0, 0},
186 {vo+1, fo+1, fo+2, 0},
187 {vo+2, fo+0, fo+3, 0},
188 {vo+3, fo+3, fo+1, 0}
189 };
190
191 return ScvCornerStorage{ {p_[map[localScvIdx][0]],
192 p_[map[localScvIdx][1]],
193 p_[map[localScvIdx][2]],
194 p_[map[localScvIdx][3]]} };
195 }
196 default:
197 DUNE_THROW(Dune::NotImplemented, "Box scv geometries for dim=" << dim
198 << " dimWorld=" << dimWorld
199 << " corners=" << corners_);
200 }
201 }
202
203
205 ScvfCornerStorage getScvfCorners(unsigned int localScvfIdx) const
206 {
207 // proceed according to number of corners
208 switch (corners_)
209 {
210 case 3: // triangle
211 {
213 static const std::uint8_t fo = 4;
214 static const std::uint8_t map[3][2] =
215 {
216 {0, fo+0},
217 {fo+1, 0},
218 {0, fo+2}
219 };
220
221 return ScvfCornerStorage{ {p_[map[localScvfIdx][0]],
222 p_[map[localScvfIdx][1]]} };
223 }
224 case 4: // quadrilateral
225 {
227 static const std::uint8_t fo = 5;
228 static const std::uint8_t map[4][2] =
229 {
230 {fo+0, 0},
231 {0, fo+1},
232 {0, fo+2},
233 {fo+3, 0}
234 };
235
236 return ScvfCornerStorage{ {p_[map[localScvfIdx][0]],
237 p_[map[localScvfIdx][1]]} };
238 }
239 default:
240 DUNE_THROW(Dune::NotImplemented, "Box scvf geometries for dim=" << dim
241 << " dimWorld=" << dimWorld
242 << " corners=" << corners_);
243 }
244 }
245
247 ScvfCornerStorage getBoundaryScvfCorners(const Intersection& is,
248 const typename Intersection::Geometry& isGeom,
249 unsigned int indexInIntersection) const
250 {
251 const auto refElement = referenceElement(elementGeometry_);
252
253 const auto vIdxLocal = refElement.subEntity(is.indexInInside(), 1, indexInIntersection, dim);
254 if (indexInIntersection == 0)
255 return ScvfCornerStorage({p_[vIdxLocal+1], isGeom.center()});
256 else if (indexInIntersection == 1)
257 return ScvfCornerStorage({isGeom.center(), p_[vIdxLocal+1]});
258 else
259 DUNE_THROW(Dune::InvalidStateException, "local index exceeds the number of corners of 2d intersections");
260 }
261
263 template <int w = dimWorld>
264 typename std::enable_if<w == 3, GlobalPosition>::type
265 normal(const ScvfCornerStorage& scvfCorners,
266 const std::vector<unsigned int>& scvIndices) const
267 {
268 const auto v1 = elementGeometry_.corner(1) - elementGeometry_.corner(0);
269 const auto v2 = elementGeometry_.corner(2) - elementGeometry_.corner(0);
270 const auto v3 = Dumux::crossProduct(v1, v2);
271 const auto t = scvfCorners[1] - scvfCorners[0];
272 GlobalPosition normal = Dumux::crossProduct(v3, t);
273 normal /= normal.two_norm();
274
276 const auto v = elementGeometry_.corner(scvIndices[1]) - elementGeometry_.corner(scvIndices[0]);
277 const auto s = v*normal;
278 if (std::signbit(s))
279 normal *= -1;
280
281 return normal;
282 }
283
285 template <int w = dimWorld>
286 typename std::enable_if<w == 2, GlobalPosition>::type
287 normal(const ScvfCornerStorage& scvfCorners,
288 const std::vector<unsigned int>& scvIndices) const
289 {
291 const auto t = scvfCorners[1] - scvfCorners[0];
292 GlobalPosition normal({-t[1], t[0]});
293 normal /= normal.two_norm();
294
296 const auto v = elementGeometry_.corner(scvIndices[1]) - elementGeometry_.corner(scvIndices[0]);
297 const auto s = v*normal;
298 if (std::signbit(s))
299 normal *= -1;
300
301 return normal;
302 }
303
305 template <int w = dimWorld>
306 typename std::enable_if<w == 3, Scalar>::type
307 scvVolume(const ScvCornerStorage& p) const
308 {
309 return 0.5*Dumux::crossProduct(p[3]-p[0], p[2]-p[1]).two_norm();
310 }
311
313 template <int w = dimWorld>
314 typename std::enable_if<w == 2, Scalar>::type
315 scvVolume(const ScvCornerStorage& p) const
316 {
319 using std::abs;
320 return 0.5*abs(Dumux::crossProduct(p[3]-p[0], p[2]-p[1]));
321 }
322
324 Scalar scvfArea(const ScvfCornerStorage& p) const
325 {
326 return (p[1]-p[0]).two_norm();
327 }
328
329protected:
330 const typename Element::Geometry& elementGeometry_;
331 std::size_t corners_; // number of element corners
332 std::array<GlobalPosition, maxPoints> p_; // the points needed for construction of the geometries
333};
334
336template <class GridView, class ScvType, class ScvfType>
337class BoxGeometryHelper<GridView, 3, ScvType, ScvfType>
338{
339 using Scalar = typename GridView::ctype;
340 using GlobalPosition = typename Dune::FieldVector<Scalar, GridView::dimensionworld>;
341 using ScvCornerStorage = typename ScvType::Traits::CornerStorage;
342 using ScvfCornerStorage = typename ScvfType::Traits::CornerStorage;
343
344 using Element = typename GridView::template Codim<0>::Entity;
345 using Intersection = typename GridView::Intersection;
346
347 static constexpr auto dim = GridView::dimension;
348 static constexpr auto dimWorld = GridView::dimensionworld;
349
352 static constexpr int maxPoints = 27;
353public:
354 BoxGeometryHelper(const typename Element::Geometry& geometry)
355 : elementGeometry_(geometry)
356 , corners_(geometry.corners())
357 {
358 const auto refElement = referenceElement(geometry);
359
360 // the element center
361 p_[0] = geometry.center();
362
363 // vertices
364 for (int i = 0; i < corners_; ++i)
365 p_[i+1] = geometry.corner(i);
366
367 // edge midpoints
368 for (int i = 0; i < refElement.size(dim-1); ++i)
369 p_[i+corners_+1] = geometry.global(refElement.position(i, dim-1));
370
371 // face midpoints
372 for (int i = 0; i < refElement.size(1); ++i)
373 p_[i+corners_+1+refElement.size(dim-1)] = geometry.global(refElement.position(i, 1));
374 }
375
377 ScvCornerStorage getScvCorners(unsigned int localScvIdx) const
378 {
379 // proceed according to number of corners of the element
380 switch (corners_)
381 {
382 case 4: // tetrahedron
383 {
385 static const std::uint8_t vo = 1;
386 static const std::uint8_t eo = 5;
387 static const std::uint8_t fo = 11;
388 static const std::uint8_t map[4][8] =
389 {
390 {vo+0, eo+0, eo+1, fo+0, eo+3, fo+1, fo+2, 0},
391 {vo+1, eo+2, eo+0, fo+0, eo+4, fo+3, fo+1, 0},
392 {vo+2, eo+1, eo+2, fo+0, eo+5, fo+2, fo+3, 0},
393 {vo+3, eo+3, eo+5, fo+2, eo+4, fo+1, fo+3, 0}
394 };
395
396 return ScvCornerStorage{ {p_[map[localScvIdx][0]],
397 p_[map[localScvIdx][1]],
398 p_[map[localScvIdx][2]],
399 p_[map[localScvIdx][3]],
400 p_[map[localScvIdx][4]],
401 p_[map[localScvIdx][5]],
402 p_[map[localScvIdx][6]],
403 p_[map[localScvIdx][7]]} };
404 }
405 case 6: // prism
406 {
408 static const std::uint8_t vo = 1;
409 static const std::uint8_t eo = 7;
410 static const std::uint8_t fo = 16;
411 static const std::uint8_t map[6][8] =
412 {
413 {vo+0, eo+3, eo+4, fo+3, eo+0, fo+0, fo+1, 0},
414 {vo+1, eo+5, eo+3, fo+3, eo+1, fo+2, fo+0, 0},
415 {vo+2, eo+4, eo+5, fo+3, eo+2, fo+1, fo+2, 0},
416 {vo+3, eo+7, eo+6, fo+4, eo+0, fo+1, fo+0, 0},
417 {vo+4, eo+6, eo+8, fo+4, eo+1, fo+0, fo+2, 0},
418 {vo+5, eo+8, eo+7, fo+4, eo+2, fo+2, fo+1, 0}
419 };
420
421 return ScvCornerStorage{ {p_[map[localScvIdx][0]],
422 p_[map[localScvIdx][1]],
423 p_[map[localScvIdx][2]],
424 p_[map[localScvIdx][3]],
425 p_[map[localScvIdx][4]],
426 p_[map[localScvIdx][5]],
427 p_[map[localScvIdx][6]],
428 p_[map[localScvIdx][7]]} };
429 }
430 case 8: // hexahedron
431 {
433 static const std::uint8_t vo = 1;
434 static const std::uint8_t eo = 9;
435 static const std::uint8_t fo = 21;
436 static const std::uint8_t map[8][8] =
437 {
438 {vo+0, eo+6, eo+4, fo+4, eo+0, fo+2, fo+0, 0},
439 {vo+1, eo+5, eo+6, fo+4, eo+1, fo+1, fo+2, 0},
440 {vo+2, eo+4, eo+7, fo+4, eo+2, fo+0, fo+3, 0},
441 {vo+3, eo+7, eo+5, fo+4, eo+3, fo+3, fo+1, 0},
442 {vo+4, eo+8, eo+10, fo+5, eo+0, fo+0, fo+2, 0},
443 {vo+5, eo+10, eo+9, fo+5, eo+1, fo+2, fo+1, 0},
444 {vo+6, eo+11, eo+8, fo+5, eo+2, fo+3, fo+0, 0},
445 {vo+7, eo+9, eo+11, fo+5, eo+3, fo+1, fo+3, 0},
446 };
447
448 return ScvCornerStorage{ {p_[map[localScvIdx][0]],
449 p_[map[localScvIdx][1]],
450 p_[map[localScvIdx][2]],
451 p_[map[localScvIdx][3]],
452 p_[map[localScvIdx][4]],
453 p_[map[localScvIdx][5]],
454 p_[map[localScvIdx][6]],
455 p_[map[localScvIdx][7]]} };
456 }
457 default:
458 DUNE_THROW(Dune::NotImplemented, "Box scv geometries for dim=" << dim
459 << " dimWorld=" << dimWorld
460 << " corners=" << corners_);
461 }
462 }
463
465 ScvfCornerStorage getScvfCorners(unsigned int localScvfIdx) const
466 {
467 // proceed according to number of corners of the element
468 switch (corners_)
469 {
470 case 4: // tetrahedron
471 {
473 static const std::uint8_t eo = 5;
474 static const std::uint8_t fo = 11;
475 static const std::uint8_t map[6][4] =
476 {
477 {eo+0, fo+0, fo+1, 0},
478 {fo+0, eo+1, 0, fo+2},
479 {eo+2, fo+0, fo+3, 0},
480 {fo+2, eo+3, 0, fo+1},
481 {fo+3, 0, eo+4, fo+1},
482 {eo+5, fo+2, fo+3, 0}
483 };
484
485 return ScvfCornerStorage{ {p_[map[localScvfIdx][0]],
486 p_[map[localScvfIdx][1]],
487 p_[map[localScvfIdx][2]],
488 p_[map[localScvfIdx][3]]} };
489 }
490 case 6: // prism
491 {
493 static const std::uint8_t eo = 7;
494 static const std::uint8_t fo = 16;
495 static const std::uint8_t map[9][4] =
496 {
497 {eo+0, fo+0, fo+1, 0},
498 {eo+1, fo+2, fo+0, 0},
499 {eo+2, fo+1, fo+2, 0},
500 {eo+3, fo+0, fo+3, 0},
501 {eo+4, fo+3, fo+1, 0},
502 {eo+5, fo+2, fo+3, 0},
503 {eo+6, fo+4, fo+0, 0},
504 {eo+7, fo+1, fo+4, 0},
505 {eo+8, fo+4, fo+2, 0}
506 };
507
508 return ScvfCornerStorage{ {p_[map[localScvfIdx][0]],
509 p_[map[localScvfIdx][1]],
510 p_[map[localScvfIdx][2]],
511 p_[map[localScvfIdx][3]]} };
512 }
513 case 8: // hexahedron
514 {
516 static const std::uint8_t eo = 9;
517 static const std::uint8_t fo = 21;
518 static const std::uint8_t map[12][4] =
519 {
520 {fo+0, eo+0, 0, fo+2},
521 {fo+1, 0, eo+1, fo+2},
522 {fo+3, eo+2, 0, fo+0},
523 {eo+3, fo+3, fo+1, 0},
524 {fo+4, eo+4, 0, fo+0},
525 {eo+5, fo+4, fo+1, 0},
526 {eo+6, fo+4, fo+2, 0},
527 {fo+4, eo+7, 0, fo+3},
528 { 0, fo+0, fo+5, eo+8},
529 {eo+9, fo+1, fo+5, 0},
530 {eo+10, fo+2, fo+5, 0},
531 {eo+11, fo+5, fo+3, 0}
532 };
533
534 return ScvfCornerStorage{ {p_[map[localScvfIdx][0]],
535 p_[map[localScvfIdx][1]],
536 p_[map[localScvfIdx][2]],
537 p_[map[localScvfIdx][3]]} };
538 }
539 default:
540 DUNE_THROW(Dune::NotImplemented, "Box scv geometries for dim=" << dim
541 << " dimWorld=" << dimWorld
542 << " corners=" << corners_);
543 }
544 }
545
547 ScvfCornerStorage getBoundaryScvfCorners(const Intersection& is,
548 const typename Intersection::Geometry& geometry,
549 unsigned int indexInIntersection) const
550 {
551 const auto refElement = referenceElement(elementGeometry_);
552 const auto faceRefElem = referenceElement(geometry);
553
554 GlobalPosition pi[9];
555 auto corners = geometry.corners();
556
557 // the facet center
558 pi[0] = geometry.center();
559
560 // corners
561 const auto idxInInside = is.indexInInside();
562 for (int i = 0; i < corners; ++i)
563 {
564 const auto vIdxLocal = refElement.subEntity(idxInInside, 1, i, dim);
565 pi[i+1] = elementGeometry_.corner(vIdxLocal);
566 }
567
568 // edge midpoints
569 for (int i = 0; i < faceRefElem.size(1); ++i)
570 {
571 const auto edgeIdxLocal = refElement.subEntity(idxInInside, 1, i, dim-1);
572 pi[i+corners+1] = p_[edgeIdxLocal+corners_+1];
573 }
574
575 // procees according to number of corners
576 switch (corners)
577 {
578 case 3: // triangle
579 {
581 static const std::uint8_t vo = 1;
582 static const std::uint8_t eo = 4;
583 static const std::uint8_t map[3][4] =
584 {
585 {vo+0, eo+0, eo+1, 0},
586 {vo+1, eo+2, eo+0, 0},
587 {vo+2, eo+1, eo+2, 0}
588 };
589
590 return ScvfCornerStorage{ {pi[map[indexInIntersection][0]],
591 pi[map[indexInIntersection][1]],
592 pi[map[indexInIntersection][2]],
593 pi[map[indexInIntersection][3]]} };
594 }
595 case 4: // quadrilateral
596 {
598 static const std::uint8_t vo = 1;
599 static const std::uint8_t eo = 5;
600 static const std::uint8_t map[4][4] =
601 {
602 {vo+0, eo+2, eo+0, 0},
603 {vo+1, eo+1, eo+2, 0},
604 {vo+2, eo+0, eo+3, 0},
605 {vo+3, eo+3, eo+1, 0}
606 };
607
608 return ScvfCornerStorage{ {pi[map[indexInIntersection][0]],
609 pi[map[indexInIntersection][1]],
610 pi[map[indexInIntersection][2]],
611 pi[map[indexInIntersection][3]]} };
612 }
613 default:
614 DUNE_THROW(Dune::NotImplemented, "Box scvf boundary geometries for dim=" << dim
615 << " dimWorld=" << dimWorld
616 << " corners=" << corners);
617 }
618 }
619
621 GlobalPosition normal(const ScvfCornerStorage& p,
622 const std::vector<unsigned int>& scvIndices) const
623 {
624 auto normal = Dumux::crossProduct(p[1]-p[0], p[2]-p[0]);
625 normal /= normal.two_norm();
626
627 const auto v = elementGeometry_.corner(scvIndices[1]) - elementGeometry_.corner(scvIndices[0]);
628 const auto s = v*normal;
629 if (std::signbit(s))
630 normal *= -1;
631
632 return normal;
633 }
634
636 Scalar scvVolume(const ScvCornerStorage& p) const
637 {
638 // after Grandy 1997, Efficient computation of volume of hexahedron
639 const auto v = p[7]-p[0];
640 return 1.0/6.0 * ( Dumux::tripleProduct(v, p[1]-p[0], p[3]-p[5])
641 + Dumux::tripleProduct(v, p[4]-p[0], p[5]-p[6])
642 + Dumux::tripleProduct(v, p[2]-p[0], p[6]-p[3]));
643 }
644
646 Scalar scvfArea(const ScvfCornerStorage& p) const
647 {
648 // after Wolfram alpha quadrilateral area
649 return 0.5*Dumux::crossProduct(p[3]-p[0], p[2]-p[1]).two_norm();
650 }
651
652protected:
653 const typename Element::Geometry& elementGeometry_;
654 std::size_t corners_; // number of element corners
655 std::array<GlobalPosition, maxPoints> p_; // the points needed for construction of the scv/scvf geometries
656};
657
658} // end namespace Dumux
659
660#endif
math.hh
Define some often used mathematical functions.
Dumux::normal
Vector normal(const Vector &v)
Create a vector normal to the given one (v is expected to be non-zero)
Definition: geometry/normal.hh:36
Dumux::crossProduct
Dune::FieldVector< Scalar, 3 > crossProduct(const Dune::FieldVector< Scalar, 3 > &vec1, const Dune::FieldVector< Scalar, 3 > &vec2)
Cross product of two vectors in three-dimensional Euclidean space.
Definition: math.hh:640
Dumux::tripleProduct
Scalar tripleProduct(const Dune::FieldVector< Scalar, 3 > &vec1, const Dune::FieldVector< Scalar, 3 > &vec2, const Dune::FieldVector< Scalar, 3 > &vec3)
Triple product of three vectors in three-dimensional Euclidean space retuning scalar.
Definition: math.hh:669
Dumux
Definition: adapt.hh:29
Dumux::BoxGeometryHelper
Create sub control volumes and sub control volume face geometries.
Definition: boxgeometryhelper.hh:36
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::scvfArea
Scalar scvfArea(const ScvfCornerStorage &scvfCorners) const
get scvf area
Definition: boxgeometryhelper.hh:105
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::BoxGeometryHelper
BoxGeometryHelper(const typename Element::Geometry &geometry)
Definition: boxgeometryhelper.hh:57
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::getScvCorners
ScvCornerStorage getScvCorners(unsigned int localScvIdx) const
Create a vector with the scv corners.
Definition: boxgeometryhelper.hh:62
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::p_
std::array< GlobalPosition, maxPoints > p_
Definition: boxgeometryhelper.hh:113
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::normal
GlobalPosition normal(const ScvfCornerStorage &scvfCorners, const std::vector< unsigned int > &scvIndices) const
get scvf normal vector
Definition: boxgeometryhelper.hh:90
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::getScvfCorners
ScvfCornerStorage getScvfCorners(unsigned int localScvfIdx) const
Create a vector with the corners of sub control volume faces.
Definition: boxgeometryhelper.hh:76
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::getBoundaryScvfCorners
ScvfCornerStorage getBoundaryScvfCorners(const Intersection &is, const typename Intersection::Geometry &geometry, unsigned int indexInIntersection) const
Create the sub control volume face geometries on the boundary.
Definition: boxgeometryhelper.hh:82
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::elementGeometry_
const Element::Geometry & elementGeometry_
Reference to the element geometry.
Definition: boxgeometryhelper.hh:111
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::scvVolume
Scalar scvVolume(const ScvCornerStorage &scvCorners) const
get scv volume
Definition: boxgeometryhelper.hh:99
Dumux::BoxGeometryHelper< GridView, 1, ScvType, ScvfType >::corners_
std::size_t corners_
Definition: boxgeometryhelper.hh:112
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::normal
std::enable_if< w==2, GlobalPosition >::type normal(const ScvfCornerStorage &scvfCorners, const std::vector< unsigned int > &scvIndices) const
get scvf normal vector for dim == 2, dimworld == 2
Definition: boxgeometryhelper.hh:287
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::scvfArea
Scalar scvfArea(const ScvfCornerStorage &p) const
get scvf area
Definition: boxgeometryhelper.hh:324
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::p_
std::array< GlobalPosition, maxPoints > p_
Definition: boxgeometryhelper.hh:332
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::getScvfCorners
ScvfCornerStorage getScvfCorners(unsigned int localScvfIdx) const
Create a vector with the corners of sub control volume faces.
Definition: boxgeometryhelper.hh:205
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::elementGeometry_
const Element::Geometry & elementGeometry_
Reference to the element geometry.
Definition: boxgeometryhelper.hh:330
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::getBoundaryScvfCorners
ScvfCornerStorage getBoundaryScvfCorners(const Intersection &is, const typename Intersection::Geometry &isGeom, unsigned int indexInIntersection) const
Create the sub control volume face geometries on the boundary.
Definition: boxgeometryhelper.hh:247
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::BoxGeometryHelper
BoxGeometryHelper(const typename Element::Geometry &geometry)
Definition: boxgeometryhelper.hh:137
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::scvVolume
std::enable_if< w==3, Scalar >::type scvVolume(const ScvCornerStorage &p) const
get scv volume for dim == 2, dimworld == 3
Definition: boxgeometryhelper.hh:307
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::corners_
std::size_t corners_
Definition: boxgeometryhelper.hh:331
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::getScvCorners
ScvCornerStorage getScvCorners(unsigned int localScvIdx) const
Create a vector with the scv corners.
Definition: boxgeometryhelper.hh:156
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::scvVolume
std::enable_if< w==2, Scalar >::type scvVolume(const ScvCornerStorage &p) const
get scv volume for dim == 2, dimworld == 2
Definition: boxgeometryhelper.hh:315
Dumux::BoxGeometryHelper< GridView, 2, ScvType, ScvfType >::normal
std::enable_if< w==3, GlobalPosition >::type normal(const ScvfCornerStorage &scvfCorners, const std::vector< unsigned int > &scvIndices) const
get scvf normal vector for dim == 2, dimworld == 3
Definition: boxgeometryhelper.hh:265
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::scvfArea
Scalar scvfArea(const ScvfCornerStorage &p) const
get scvf area
Definition: boxgeometryhelper.hh:646
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::elementGeometry_
const Element::Geometry & elementGeometry_
Reference to the element geometry.
Definition: boxgeometryhelper.hh:653
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::normal
GlobalPosition normal(const ScvfCornerStorage &p, const std::vector< unsigned int > &scvIndices) const
get scvf normal vector
Definition: boxgeometryhelper.hh:621
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::p_
std::array< GlobalPosition, maxPoints > p_
Definition: boxgeometryhelper.hh:655
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::BoxGeometryHelper
BoxGeometryHelper(const typename Element::Geometry &geometry)
Definition: boxgeometryhelper.hh:354
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::scvVolume
Scalar scvVolume(const ScvCornerStorage &p) const
get scv volume
Definition: boxgeometryhelper.hh:636
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::getScvCorners
ScvCornerStorage getScvCorners(unsigned int localScvIdx) const
Create a vector with the scv corners.
Definition: boxgeometryhelper.hh:377
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::getBoundaryScvfCorners
ScvfCornerStorage getBoundaryScvfCorners(const Intersection &is, const typename Intersection::Geometry &geometry, unsigned int indexInIntersection) const
Create the sub control volume face geometries on the boundary.
Definition: boxgeometryhelper.hh:547
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::corners_
std::size_t corners_
Definition: boxgeometryhelper.hh:654
Dumux::BoxGeometryHelper< GridView, 3, ScvType, ScvfType >::getScvfCorners
ScvfCornerStorage getScvfCorners(unsigned int localScvfIdx) const
Create a vector with the scvf corners.
Definition: boxgeometryhelper.hh:465
Generated by doxygen 1.9.3