Rheolef  7.1
an efficient C++ finite element environment
rotating-hill-statio.h
Go to the documentation of this file.
1 #include "rotating-hill.h"
22 struct f: phi { // - d phi/dt
23  Float operator() (const point& x) const {
24  return (4*nu/t0 - 4*nu*dist2(x,x0t())/sqr(t0+4*nu*t)
25  -2*dot(d_x0t_dt(),x-x0t())/(t0+4*nu*t))*phi::operator()(x); }
26  f (size_t d1, Float nu1, Float t1=0) : phi(d1,nu1,t1) {}
27 };
see the Float page for the full documentation
see the point page for the full documentation
rheolef::details::is_vec dot
T dist2(const point_basic< T > &x, const point_basic< T > &y)
Definition: point.h:293
Convection-diffusion equation – the rotating hill benchmark.
Definition: cavity_dg.h:29
f(size_t d1, Float nu1, Float t1=0)
point operator()(const point &x) const
Definition: cavity_dg.h:30
Definition: nu.h:26
Definition: phi.h:25
Float t
Definition: rotating-hill.h:46
point x0t() const
Definition: rotating-hill.h:38
static constexpr Float t0
Definition: rotating-hill.h:47
point d_x0t_dt() const
Definition: rotating-hill.h:42
Float operator()(const Float &x) const
Definition: phi.h:27