Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
combustion_exact.icc
Go to the documentation of this file.
1
25#include "lambda2alpha.h"
26struct 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};
34struct 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)
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)