Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
lambda_c.h

The combustion problem – the critical parameter value as a function.

The combustion problem – the critical parameter value as a function

struct alpha_c_fun {
typedef Float value_type;
typedef Float float_type;
alpha_c_fun() : _f1(0) {}
Float residue (const Float& a) const { return tanh(a) - 1/a; }
void update_derivative (const Float& a) const {
_f1 = 1/sqr(cosh(a)) + 1/sqr(a); }
Float derivative_solve (const Float& r) const { return r/_f1; }
Float dual_space_norm (const Float& r) const { return abs(r); }
mutable Float _f1;
};
Float tol = numeric_limits<Float>::epsilon();
size_t max_iter = 100;
Float ac = 1;
newton (alpha_c_fun(), ac, tol, max_iter);
return ac;
}
Float ac = alpha_c();
return 8*sqr(ac/cosh(ac));
}
see the Float page for the full documentation
Float alpha_c()
Definition lambda_c.h:36
Float lambda_c()
Definition lambda_c.h:44
Float float_type
Definition lambda_c.h:27
Float _f1
Definition lambda_c.h:34
Float derivative_solve(const Float &r) const
Definition lambda_c.h:32
Float value_type
Definition lambda_c.h:26
Float residue(const Float &a) const
Definition lambda_c.h:29
Float dual_space_norm(const Float &r) const
Definition lambda_c.h:33
void update_derivative(const Float &a) const
Definition lambda_c.h:30
Definition cavity_dg.h:29