Rheolef  7.1
an efficient C++ finite element environment
sinusrad_laplace.h
Go to the documentation of this file.
1 struct f {
26  Float operator() (const point& x) const {
27  Float r = sqrt(sqr(x[0])+sqr(x[1])+sqr(x[2]));
28  return (n == 2) ? 4 : n*n*pow(r,n-2);
29  }
30  f(size_t=0) : n(2) {}
31  size_t n;
32 };
33 struct g {
34  Float operator() (const point& x) const {
35  Float r = sqrt(sqr(x[0])+sqr(x[1])+sqr(x[2]));
36  return 1-pow(r,n);
37  }
38  g(size_t=0) : n(2) {}
39  size_t n;
40 };
f::n
size_t n
Definition: sinusrad_laplace.h:31
rheolef::pow
space_mult_list< T, M > pow(const space_basic< T, M > &X, size_t n)
Definition: space_mult.h:120
g::g
g(size_t=0)
Definition: sinusrad_laplace.h:38
f::operator()
point operator()(const point &x) const
Definition: cavity_dg.h:30
phi::r
Float r
Definition: phi.h:54
Float
see the Float page for the full documentation
point
see the point page for the full documentation
g::operator()
point operator()(const point &x) const
Definition: cavity_dg.h:26
f::f
f(size_t=0)
Definition: sinusrad_laplace.h:30
g
Definition: cavity_dg.h:25
g::n
size_t n
Definition: sinusrad_laplace.h:39
f
Definition: cavity_dg.h:29