OGS
BaseLib::GradualSubdivision Class Reference

Detailed Description

Gradual subdivision operator with a constant multiplier_

Definition at line 60 of file Subdivision.h.

#include <Subdivision.h>

Inheritance diagram for BaseLib::GradualSubdivision:
[legend]
Collaboration diagram for BaseLib::GradualSubdivision:
[legend]

Public Member Functions

 GradualSubdivision (const double L, const double dL0, const double max_dL, const double multiplier)
std::vector< double > operator() () const override
 Returns a vector of subdivided points.
Public Member Functions inherited from BaseLib::ISubdivision
virtual ~ISubdivision ()=default

Private Attributes

const double length_
const double dL0_
const double max_dL_
const double multiplier_

Constructor & Destructor Documentation

◆ GradualSubdivision()

BaseLib::GradualSubdivision::GradualSubdivision ( const double L,
const double dL0,
const double max_dL,
const double multiplier )

Constructor

Parameters
Ltotal length to be subdivided
dL0initial cell length
max_dLmaximum cell length
multipliermultiplier to cell length

Definition at line 13 of file Subdivision.cpp.

17 : length_(L), dL0_(dL0), max_dL_(max_dL), multiplier_(multiplier)
18{
19 // Check if accumulated subdivisions can ever sum up to length.
20 // Cf. geometric series formula.
21 if (multiplier < 1.0 && dL0 / (1.0 - multiplier) < L)
22 {
24 "Using dL0={:g} and multiplier={:g} the generated subdivisions can "
25 "not sum up to a total length of {:g}.",
26 dL0,
27 multiplier,
28 L);
29 }
30}
#define OGS_FATAL(...)
Definition Error.h:19

References dL0_, length_, max_dL_, multiplier_, and OGS_FATAL.

Member Function Documentation

◆ operator()()

std::vector< double > BaseLib::GradualSubdivision::operator() ( ) const
overridevirtual

Returns a vector of subdivided points.

Implements BaseLib::ISubdivision.

Definition at line 32 of file Subdivision.cpp.

33{
34 std::vector<double> vec_x;
35
36 double x = 0;
37 unsigned i = 0;
38 do
39 {
40 vec_x.push_back(x);
41 x += std::min(max_dL_,
42 dL0_ * std::pow(multiplier_, static_cast<double>(i)));
43 i++;
44 } while (x < length_);
45
46 if (vec_x.back() < length_)
47 {
48 double last_dx = vec_x[vec_x.size() - 1] - vec_x[vec_x.size() - 2];
49 if (length_ - vec_x.back() < last_dx)
50 {
51 vec_x[vec_x.size() - 1] = length_;
52 }
53 else
54 {
55 vec_x.push_back(length_);
56 }
57 }
58 return vec_x;
59}

References dL0_, length_, max_dL_, and multiplier_.

Member Data Documentation

◆ dL0_

const double BaseLib::GradualSubdivision::dL0_
private

Definition at line 80 of file Subdivision.h.

Referenced by GradualSubdivision(), and operator()().

◆ length_

const double BaseLib::GradualSubdivision::length_
private

Definition at line 79 of file Subdivision.h.

Referenced by GradualSubdivision(), and operator()().

◆ max_dL_

const double BaseLib::GradualSubdivision::max_dL_
private

Definition at line 81 of file Subdivision.h.

Referenced by GradualSubdivision(), and operator()().

◆ multiplier_

const double BaseLib::GradualSubdivision::multiplier_
private

Definition at line 82 of file Subdivision.h.

Referenced by GradualSubdivision(), and operator()().


The documentation for this class was generated from the following files: