30 return 0.5*exp( -
nu*
t - norm2(x-
xt)/sqr(
r0));
33 return ( -
nu + 2*dot(
vt,x-
xt)/sqr(
r0))*
f(x);
38 s(0,0) = c0*(1+cos(2*
t));
39 s(1,1) = c0*(1-cos(2*
t));
48 s(0,0) = c1*(1+cos(2*
t)) - 2*c0*sin(2*
t);
49 s(1,1) = c1*(1-cos(2*
t)) + 2*c0*sin(2*
t);
51 s(1,0) = c1*sin(2*
t) + 2*c0*cos(2*
t);
58 for (
size_t i = 0; i < 2; i++)
59 for (
size_t j = 0; j < 2; j++)
60 for (
size_t k = 0; k < 2; k++)
61 A(i,j,k) = s(i,j)*vec[k];
67 return (norm2(x-
xt)/sqr(
r0) - 1)*(4/sqr(
r0))*s;
73 x0[0]*sin(
t) +
x0[1]*cos(
t));
75 x0[0]*cos(
t) -
x0[1]*sin(
t));
see the Float page for the full documentation
see the point page for the full documentation
see the tensor3 page for the full documentation
see the tensor page for the full documentation
tensor operator()(const point &x) const
chi(Float eps, Float nu, Float t=0)
grad_sigma_g(Float nu, Float t=0)
tensor3 operator()(const point &x) const
tensor laplacian(const point &x) const
tensor operator()(const point &x) const
tensor3 grad(const point &x) const
Float f(const point &x) const
Float df_dt(const point &x) const
tensor time_derivative(const point &x) const
sigma_exact(Float nu1, Float t1=0)