1#ifndef _RHEOLEF_INTERPOLATE_RTK_POLYNOM_ICC
2#define _RHEOLEF_INTERPOLATE_RTK_POLYNOM_ICC
30#include "rheolef/basis.h"
40 case 0:
return point(1,0);
41 case 1:
return point(0,1);
42 case 2:
return point(x[0],x[1]);
44 case 3:
return point(x[0],-x[1]);
45 case 4:
return point(0,x[0]);
46 case 5:
return point(x[1],0);
47 case 6:
return point(sqr(x[0]),x[0]*x[1]);
48 default:
return point(x[0]*x[1],sqr(x[1]));
55 case 0:
return point(1,0);
56 case 1:
return point(0,1);
57 case 2:
return point(x[0],0);
58 case 3:
return point(0,x[1]);
60 case 4:
return point(0,x[0]);
61 case 5:
return point(x[1],0);
62 case 6:
return point(x[0]*x[1],0);
63 case 7:
return point(0,x[0]*x[1]);
64 case 8:
return point(sqr(x[0]),0);
65 case 9:
return point(0,sqr(x[1]));
66 case 10:
return point(sqr(x[0])*x[1],0);
67 default:
return point(0,x[0]*sqr(x[1]));
84 case 6:
return 3*x[0];
85 default:
return 3*x[1];
101 case 8:
return 2*x[0];
102 case 9:
return 2*x[1];
103 case 10:
return 2*x[0]*x[1];
104 default:
return 2*x[0]*x[1];
128 case 2:
return point(cos(
w*(x[0]+2*x[1])),
129 sin(
w*(x[0]-2*x[1])));
130 default:
return point(cos(
w*(x[0]+2*x[1]+x[2])),
131 sin(
w*(x[0]-2*x[1]-x[2])),
132 sin(
w*(x[0]+2*x[1]-x[2])));
137 case 2:
return -
w*( sin(
w*( x[0]+2*x[1]))
138 + 2*cos(
w*(-x[0]+2*x[1])));
139 default:
return -
w*( sin(
w*( x[0]+2*x[1]+x[2]))
140 + 2*cos(
w*(-x[0]+2*x[1]+x[2]))
141 + cos(
w*(-x[0]-2*x[1]+x[2])));
see the Float page for the full documentation
see the point page for the full documentation
div_psi(char hat_K_name, size_t i)
Float operator()(const point &x) const
Float operator()(const point &x) const
div_u(size_t d, Float w=acos(Float(-1)))
point operator()(const point &x) const
Float div(const point &x) const
static size_t n_index(char hat_K_name, size_t k)
psi(char hat_K_name1, size_t i1)
point operator()(const point &x) const
Float div(const point &x) const
u_exact(size_t d1, Float w1=acos(Float(-1)))