/usr/share/doc/rheolef-doc/examples/sinusprod_helmholtz.icc is in rheolef-doc 6.7-6.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | 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;
};
|