Rheolef  7.2
an efficient C++ finite element environment
combustion_exact.icc
Go to the documentation of this file.
1 #include "lambda2alpha.h"
26 struct u_exact {
27  Float operator() (const point& x) const {
28  return 2*log(cosh(a)/cosh(a*(1-2*x[0]))); }
29  u_exact (Float lambda, bool is_upper)
30  : a(lambda2alpha(lambda,is_upper)) {}
31  u_exact (Float a1) : a(a1) {}
33 };
34 struct grad_u {
35  point operator() (const point& x) const {
36  return point(4*a*tanh(a*(1-2*x[0]))); }
37  grad_u (Float lambda, bool is_upper)
38  : a(lambda2alpha(lambda,is_upper)) {}
39  grad_u (Float a1) : a(a1) {}
41 };
see the Float page for the full documentation
see the point page for the full documentation
The combustion problem – inversion of the parameter function.
Float lambda2alpha(Float lambda, bool up=false)
Definition: lambda2alpha.h:26
point operator()(const point &x) const
grad_u(Float a1)
grad_u(Float lambda, bool is_upper)
point operator()(const point &x) const
u_exact(Float a1)
u_exact(Float lambda, bool is_upper)