29 Float shift = -2*floor((x0+1)/2);
30 Float xs = x0 + shift;
31 check_macro (xs >= -1 && xs <= 1,
"invalid xs="<<xs);
35 h0(b1*t1),
t(t1),
a(a1),
b(b1),
c(c1) {}
see the Float page for the full documentation
see the point page for the full documentation
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
The Burgers problem: the Harten exact solution at t=0.
harten(Float t1=0, Float a1=1, Float b1=0.5, Float c1=0)
Float operator()(const point &x) const