Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
sinusprod_helmholtz.h

The sinus product function – right-hand-side and boundary condition for the Helmholtz problem

The sinus product function – right-hand-side and boundary condition for the Helmholtz problem

struct f {
Float operator() (const point& x) const {
switch (d) {
case 1: return (1+d*pi*pi)*sin(pi*x[0]);
case 2: return (1+d*pi*pi)*sin(pi*x[0])*sin(pi*x[1]);
default: return (1+d*pi*pi)*sin(pi*x[0])*sin(pi*x[1])*sin(pi*x[2]);
}}
f(size_t d1) : d(d1), pi(acos(Float(-1.0))) {}
size_t d; const Float pi;
};
struct g {
Float operator() (const point& x) const {
switch (d) {
case 1: return -pi;
case 2: return -pi*(sin(pi*x[0]) + sin(pi*x[1]));
default: return -pi*( sin(pi*x[0])*sin(pi*x[1])
+ sin(pi*x[1])*sin(pi*x[2])
+ sin(pi*x[2])*sin(pi*x[0]));
}}
g(size_t d1) : d(d1), pi(acos(Float(-1.0))) {}
size_t d; const Float pi;
};
see the Float page for the full documentation
see the point page for the full documentation
Definition cavity_dg.h:29
point operator()(const point &x) const
Definition cavity_dg.h:30
f()
Definition taylor.h:34
const Float pi
Definition cavity_dg.h:25
size_t d
point operator()(const point &x) const
Definition cavity_dg.h:26
g()
Definition taylor.h:29
const Float pi