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
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))) {}
};
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))) {}
};
see the Float page for the full documentation
see the point page for the full documentation
point operator()(const point &x) const
point operator()(const point &x) const