3.5-git
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
normalaxis.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 *****************************************************************************/
24#ifndef DUMUX_DISCRETIZATION_FACECENTERED_STAGGERED_NORMAL_AXIS_HH
25#define DUMUX_DISCRETIZATION_FACECENTERED_STAGGERED_NORMAL_AXIS_HH
26
27#include <cstddef>
28#include <algorithm>
29#include <numeric>
30
31namespace Dumux {
32
38template<class Vector>
39inline static std::size_t normalAxis(const Vector& v)
40{
41 using std::abs;
42
43 constexpr auto eps = 1e-8;
44 assert(std::any_of(v.begin(), v.end(), [=](auto x){ return abs(x) > eps; }));
45
46 const auto result = std::distance(
47 std::begin(v), std::find_if(v.begin(), v.end(), [eps=eps](auto x){ return abs(x) > eps; })
48 );
49
50 // make sure there is only one non-zero entry
51 assert(v[result] == std::accumulate(v.begin(), v.end(), 0.0));
52
53 return result;
54}
55
56} // end namespace Dumux
57
58#endif
static ctype distance(const Dune::FieldVector< ctype, dimWorld > &a, const Dune::FieldVector< ctype, dimWorld > &b)
Compute the shortest distance between two points.
Definition: distance.hh:292
Definition: adapt.hh:29
static std::size_t normalAxis(const Vector &v)
Definition: normalaxis.hh:39