Helper class to integrate over an elliptic domain.
More...
#include <dumux/multidomain/embedded/cylinderintegration.hh>
template<class Scalar>
class Dumux::EmbeddedCoupling::EllipseIntegration< Scalar >
Helper class to integrate over an elliptic domain.
- Note
- This is mostly useful if the integral is known and the integral mass has to be locally distributed to some non-matching domain
-
The algorithm is based on evenly spaced area elements with rejection sampling. To improve the quality of the integral, increase the sample size.
-
The area integral of a constant function is exact See Koch et al. (2020) https://doi.org/10.1016/j.jcp.2020.109369
◆ EllipseIntegration()
Constructor.
- Parameters
-
relCharLength | characteristic relative integration length (real number between 0 and 1) |
- Note
- half the characteristic length means 2*2=4 times more integration points
◆ integrationElement()
obtain ith integration element
◆ integrationPoint()
obtain ith integration point
◆ setGeometry()
set geometry of an ellipse
- Parameters
-
center | the center position |
firstAxis | first ellipse axis (length corresponding to axis length) |
secondAxis | second ellipse axis (length corresponding to axis length) |
verbosity | the verbosity level (default: 0 -> no terminal output) |
◆ size()
number of integration points
The documentation for this class was generated from the following file: