GCC Code Coverage Report

Directory:	../../../builds/dumux-repositories/
File:	/builds/dumux-repositories/dumux/dumux/geometry/normal.hh
Date:	2024-05-04 19:09:25

	Exec	Total	Coverage
Lines:	17	17	100.0%
Functions:	2	3	66.7%
Branches:	75	132	56.8%

  
      Line
      Branch
      Exec
      Source
    
      // -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
    
      // vi: set et ts=4 sw=4 sts=4:
    
      //
    
      // SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
    
      // SPDX-License-Identifier: GPL-3.0-or-later
    
      //
    
      /*!
    
       * \file
    
       * \ingroup Geometry
    
       * \brief Helper functions for normals
    
       */
    
      #ifndef DUMUX_GEOMETRY_NORMAL_HH
    
      #define DUMUX_GEOMETRY_NORMAL_HH
    
      #include <algorithm>
    
      #include <dune/common/float_cmp.hh>
    
      namespace Dumux {
    
      /*!
    
       * \ingroup Geometry
    
       * \brief Create a vector normal to the given one (v is expected to be non-zero)
    
       * \note This returns some orthogonal vector with arbitrary length
    
       */
    
      template<class Vector>
    
      6059809
      inline Vector normal(const Vector& v)
    
      {
    
          static_assert(Vector::size() > 1, "normal expects a coordinate dimension > 1");
    
          if constexpr (Vector::size() == 2)
    
        20/40✓ Branch 1 taken 400000 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 400000 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 400000 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 400000 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 400000 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 400000 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 400000 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 400000 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 400000 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 400000 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 400000 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 400000 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 400000 times.
✗ Branch 38 not taken.
✓ Branch 40 taken 400000 times.
✗ Branch 41 not taken.
✓ Branch 43 taken 400000 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 400000 times.
✗ Branch 47 not taken.
✓ Branch 49 taken 400000 times.
✗ Branch 50 not taken.
✓ Branch 52 taken 400000 times.
✗ Branch 53 not taken.
✓ Branch 55 taken 400000 times.
✗ Branch 56 not taken.
✓ Branch 58 taken 400000 times.
✗ Branch 59 not taken.

      8000960
              return Vector({-v[1], v[0]});
    
        50/84✓ Branch 0 taken 719480 times.
✓ Branch 1 taken 2069252 times.
✓ Branch 2 taken 719480 times.
✓ Branch 3 taken 2069252 times.
✓ Branch 4 taken 416300 times.
✓ Branch 5 taken 1652952 times.
✓ Branch 6 taken 416300 times.
✓ Branch 7 taken 1652952 times.
✓ Branch 8 taken 375214 times.
✓ Branch 9 taken 1277738 times.
✓ Branch 10 taken 375214 times.
✓ Branch 11 taken 1277738 times.
✓ Branch 12 taken 334572 times.
✓ Branch 13 taken 943166 times.
✓ Branch 14 taken 334572 times.
✓ Branch 15 taken 943166 times.
✓ Branch 16 taken 38794 times.
✓ Branch 17 taken 21015 times.
✓ Branch 18 taken 38794 times.
✓ Branch 19 taken 21015 times.
✓ Branch 20 taken 111954 times.
✓ Branch 21 taken 63495 times.
✓ Branch 22 taken 111954 times.
✓ Branch 23 taken 63495 times.
✓ Branch 24 taken 63495 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 63495 times.
✗ Branch 27 not taken.
✓ Branch 28 taken 779316 times.
✓ Branch 29 taken 1220684 times.
✓ Branch 30 taken 779316 times.
✓ Branch 31 taken 1220684 times.
✓ Branch 32 taken 621936 times.
✓ Branch 33 taken 598748 times.
✓ Branch 34 taken 621936 times.
✓ Branch 35 taken 598748 times.
✓ Branch 36 taken 443230 times.
✓ Branch 37 taken 155518 times.
✓ Branch 38 taken 443230 times.
✓ Branch 39 taken 155518 times.
✓ Branch 40 taken 155518 times.
✗ Branch 41 not taken.
✓ Branch 42 taken 155518 times.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✗ Branch 51 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 69 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✓ Branch 72 taken 950558 times.
✓ Branch 73 taken 1049442 times.
✓ Branch 74 taken 950558 times.
✓ Branch 75 taken 1049442 times.
✓ Branch 76 taken 798972 times.
✓ Branch 77 taken 250470 times.
✓ Branch 78 taken 798972 times.
✓ Branch 79 taken 250470 times.
✓ Branch 80 taken 250470 times.
✗ Branch 81 not taken.
✓ Branch 82 taken 250470 times.
✗ Branch 83 not taken.

      34954221
          const auto it = std::find_if(v.begin(), v.end(), [](const auto& x) { return Dune::FloatCmp::ne(x, 0.0); });
    
        3/6✗ Branch 0 not taken.
✓ Branch 1 taken 59809 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 59809 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 59809 times.

      18179427
          const auto index = std::distance(v.begin(), it);
    
        2/2✓ Branch 0 taken 43156 times.
✓ Branch 1 taken 16653 times.

      6059809
          if (index != Vector::size()-1)
    
          {
    
      5590326
              Vector normal(0.0);
    
      11180652
              normal[index] = -v[index+1];
    
      11180652
              normal[index+1] = v[index];
    
      5590326
              return normal;
    
          }
    
          else
    
          {
    
      469483
              Vector normal(0.0);
    
      938966
              normal[index-1] = -v[index];
    
      938966
              normal[index] = v[index-1];
    
      469483
              return normal;
    
          }
    
      }
    
      /*!
    
       * \ingroup Geometry
    
       * \brief Create a vector normal to the given one (v is expected to be non-zero)
    
       * \note This returns some orthogonal vector with unit length
    
       */
    
      template<class Vector>
    
      59569
      inline Vector unitNormal(const Vector& v)
    
      {
    
      59569
          auto normal = Dumux::normal(v);
    
      119138
          normal /= normal.two_norm();
    
      59569
          return normal;
    
      }
    
      } // end namespace Dumux
    
      #endif

Line	Branch	Exec	Source
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			// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
5			// SPDX-License-Identifier: GPL-3.0-or-later
6			//
7			/*!
8			* \file
9			* \ingroup Geometry
10			* \brief Helper functions for normals
11			*/
12			#ifndef DUMUX_GEOMETRY_NORMAL_HH
13			#define DUMUX_GEOMETRY_NORMAL_HH
14
15			#include <algorithm>
16			#include <dune/common/float_cmp.hh>
17
18			namespace Dumux {
19
20			/*!
21			* \ingroup Geometry
22			* \brief Create a vector normal to the given one (v is expected to be non-zero)
23			* \note This returns some orthogonal vector with arbitrary length
24			*/
25			template<class Vector>
26		6059809	inline Vector normal(const Vector& v)
27			{
28			static_assert(Vector::size() > 1, "normal expects a coordinate dimension > 1");
29
30			if constexpr (Vector::size() == 2)
31	20/40 ✓ Branch 1 taken 400000 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 400000 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 400000 times. ✗ Branch 8 not taken. ✓ Branch 10 taken 400000 times. ✗ Branch 11 not taken. ✓ Branch 13 taken 400000 times. ✗ Branch 14 not taken. ✓ Branch 16 taken 400000 times. ✗ Branch 17 not taken. ✓ Branch 19 taken 400000 times. ✗ Branch 20 not taken. ✓ Branch 22 taken 400000 times. ✗ Branch 23 not taken. ✓ Branch 25 taken 400000 times. ✗ Branch 26 not taken. ✓ Branch 28 taken 400000 times. ✗ Branch 29 not taken. ✓ Branch 31 taken 400000 times. ✗ Branch 32 not taken. ✓ Branch 34 taken 400000 times. ✗ Branch 35 not taken. ✓ Branch 37 taken 400000 times. ✗ Branch 38 not taken. ✓ Branch 40 taken 400000 times. ✗ Branch 41 not taken. ✓ Branch 43 taken 400000 times. ✗ Branch 44 not taken. ✓ Branch 46 taken 400000 times. ✗ Branch 47 not taken. ✓ Branch 49 taken 400000 times. ✗ Branch 50 not taken. ✓ Branch 52 taken 400000 times. ✗ Branch 53 not taken. ✓ Branch 55 taken 400000 times. ✗ Branch 56 not taken. ✓ Branch 58 taken 400000 times. ✗ Branch 59 not taken.	8000960	return Vector({-v[1], v[0]});
32
33	50/84 ✓ Branch 0 taken 719480 times. ✓ Branch 1 taken 2069252 times. ✓ Branch 2 taken 719480 times. ✓ Branch 3 taken 2069252 times. ✓ Branch 4 taken 416300 times. ✓ Branch 5 taken 1652952 times. ✓ Branch 6 taken 416300 times. ✓ Branch 7 taken 1652952 times. ✓ Branch 8 taken 375214 times. ✓ Branch 9 taken 1277738 times. ✓ Branch 10 taken 375214 times. ✓ Branch 11 taken 1277738 times. ✓ Branch 12 taken 334572 times. ✓ Branch 13 taken 943166 times. ✓ Branch 14 taken 334572 times. ✓ Branch 15 taken 943166 times. ✓ Branch 16 taken 38794 times. ✓ Branch 17 taken 21015 times. ✓ Branch 18 taken 38794 times. ✓ Branch 19 taken 21015 times. ✓ Branch 20 taken 111954 times. ✓ Branch 21 taken 63495 times. ✓ Branch 22 taken 111954 times. ✓ Branch 23 taken 63495 times. ✓ Branch 24 taken 63495 times. ✗ Branch 25 not taken. ✓ Branch 26 taken 63495 times. ✗ Branch 27 not taken. ✓ Branch 28 taken 779316 times. ✓ Branch 29 taken 1220684 times. ✓ Branch 30 taken 779316 times. ✓ Branch 31 taken 1220684 times. ✓ Branch 32 taken 621936 times. ✓ Branch 33 taken 598748 times. ✓ Branch 34 taken 621936 times. ✓ Branch 35 taken 598748 times. ✓ Branch 36 taken 443230 times. ✓ Branch 37 taken 155518 times. ✓ Branch 38 taken 443230 times. ✓ Branch 39 taken 155518 times. ✓ Branch 40 taken 155518 times. ✗ Branch 41 not taken. ✓ Branch 42 taken 155518 times. ✗ Branch 43 not taken. ✗ Branch 44 not taken. ✗ Branch 45 not taken. ✗ Branch 46 not taken. ✗ Branch 47 not taken. ✗ Branch 48 not taken. ✗ Branch 49 not taken. ✗ Branch 50 not taken. ✗ Branch 51 not taken. ✗ Branch 52 not taken. ✗ Branch 53 not taken. ✗ Branch 54 not taken. ✗ Branch 55 not taken. ✗ Branch 56 not taken. ✗ Branch 57 not taken. ✗ Branch 58 not taken. ✗ Branch 59 not taken. ✗ Branch 60 not taken. ✗ Branch 61 not taken. ✗ Branch 62 not taken. ✗ Branch 63 not taken. ✗ Branch 64 not taken. ✗ Branch 65 not taken. ✗ Branch 66 not taken. ✗ Branch 67 not taken. ✗ Branch 68 not taken. ✗ Branch 69 not taken. ✗ Branch 70 not taken. ✗ Branch 71 not taken. ✓ Branch 72 taken 950558 times. ✓ Branch 73 taken 1049442 times. ✓ Branch 74 taken 950558 times. ✓ Branch 75 taken 1049442 times. ✓ Branch 76 taken 798972 times. ✓ Branch 77 taken 250470 times. ✓ Branch 78 taken 798972 times. ✓ Branch 79 taken 250470 times. ✓ Branch 80 taken 250470 times. ✗ Branch 81 not taken. ✓ Branch 82 taken 250470 times. ✗ Branch 83 not taken.	34954221	const auto it = std::find_if(v.begin(), v.end(), [](const auto& x) { return Dune::FloatCmp::ne(x, 0.0); });
34	3/6 ✗ Branch 0 not taken. ✓ Branch 1 taken 59809 times. ✗ Branch 2 not taken. ✓ Branch 3 taken 59809 times. ✗ Branch 4 not taken. ✓ Branch 5 taken 59809 times.	18179427	const auto index = std::distance(v.begin(), it);
35	2/2 ✓ Branch 0 taken 43156 times. ✓ Branch 1 taken 16653 times.	6059809	if (index != Vector::size()-1)
36			{
37		5590326	Vector normal(0.0);
38		11180652	normal[index] = -v[index+1];
39		11180652	normal[index+1] = v[index];
40		5590326	return normal;
41			}
42			else
43			{
44		469483	Vector normal(0.0);
45		938966	normal[index-1] = -v[index];
46		938966	normal[index] = v[index-1];
47		469483	return normal;
48
49			}
50			}
51
52			/*!
53			* \ingroup Geometry
54			* \brief Create a vector normal to the given one (v is expected to be non-zero)
55			* \note This returns some orthogonal vector with unit length
56			*/
57			template<class Vector>
58		59569	inline Vector unitNormal(const Vector& v)
59			{
60		59569	auto normal = Dumux::normal(v);
61		119138	normal /= normal.two_norm();
62		59569	return normal;
63			}
64
65			} // end namespace Dumux
66
67			#endif
68