Rheolef  7.1
an efficient C++ finite element environment
cosinusrad_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  Float sin_over_ar = (r == 0) ? 1 : sin(a*r)/(a*r);
29  return sqr(a)*((d-1)*sin_over_ar + cos(a*r)); }
30  f(size_t d1) : d(d1), a(acos(Float(-1.0))) {}
31  size_t d; Float a;
32 };
33 struct g {
34  Float operator() (const point& x) const {
35  return cos(a*sqrt(sqr(x[0])+sqr(x[1])+sqr(x[2]))); }
36  g(size_t=0) : a(acos(Float(-1.0))) {}
38 };
f::f
f(size_t d1)
Definition: cosinusrad_laplace.h:30
f::d
size_t d
Definition: cosinusprod_dirichlet.h:30
g::g
g(size_t=0)
Definition: cosinusrad_laplace.h:36
a
Definition: diffusion_isotropic.h:25
f::operator()
point operator()(const point &x) const
Definition: cavity_dg.h:30
g::a
Float a
Definition: cosinusrad_laplace.h:37
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
g
Definition: cavity_dg.h:25
f::a
Float a
Definition: cosinusrad_laplace.h:31
f
Definition: cavity_dg.h:29