Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
lambda_c.h
Go to the documentation of this file.
1
28 alpha_c_fun() : _f1(0) {}
29 Float residue (const Float& a) const { return tanh(a) - 1/a; }
30 void update_derivative (const Float& a) const {
31 _f1 = 1/sqr(cosh(a)) + 1/sqr(a); }
32 Float derivative_solve (const Float& r) const { return r/_f1; }
33 Float dual_space_norm (const Float& r) const { return abs(r); }
34 mutable Float _f1;
35};
37 Float tol = numeric_limits<Float>::epsilon();
38 size_t max_iter = 100;
40 Float ac = 1;
41 newton (alpha_c_fun(), ac, tol, max_iter);
42 return ac;
43}
45 Float ac = alpha_c();
46 return 8*sqr(ac/cosh(ac));
47}
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