The tensor transport-diffusion benchmark – right-hand-side and exact solution.
The tensor transport-diffusion benchmark – right-hand-side and exact solution
}
}
s(0,0) = c0*(1+cos(2*
t));
s(1,1) = c0*(1-cos(2*
t));
s(0,1) =
return s;
}
s(0,0) = c1*(1+cos(2*
t)) - 2*c0*sin(2*
t);
s(1,1) = c1*(1-cos(2*
t)) + 2*c0*sin(2*
t);
s(0,1) =
s(1,0) = c1*sin(2*
t) + 2*c0*cos(2*
t);
return s;
}
for (size_t i = 0; i < 2; i++)
for (size_t j = 0; j < 2; j++)
for (size_t k = 0; k < 2; k++)
A(i,j,k) = s(i,j)*vec[k];
}
}
{
}
};
}
};
}
};
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
rheolef::details::is_vec dot
tensor_basic< T > exp(const tensor_basic< T > &a, size_t d)
T norm2(const vec< T, M > &x)
norm2(x): see the expression page for the full documentation
tensor operator()(const point &x) const
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 df_dt(const point &x) const
tensor time_derivative(const point &x) const