55 #include "et_bin_bhns_extr.h" 77 while ( fabs(dent(0).val_point(rr, M_PI/2., pp)) > 1.e-15 ||
78 fabs(dent(1).val_point(rr, M_PI/2., pp)) > 5.e-15) {
85 while ( dy > 1.e-15 ) {
91 rr =
sqrt( xxp*xxp + yyp*yyp ) ;
97 pp =
acos( xxp / rr ) ;
100 dval_y = dent(1).val_point(rr, M_PI/2., pp) ;
109 while( diff_dent_y > 1.e-15 ) {
112 ytmp = yyp + ss * nn * dy ;
114 rr =
sqrt( xxp*xxp + ytmp*ytmp ) ;
125 pp =
acos( xxp / rr ) ;
128 diff_dent_y = ss * dent(1).val_point(rr, M_PI/2., pp) ;
131 yyp += ss * (nn - 1) * dy ;
135 while ( dx > 1.e-15 ) {
141 rr =
sqrt( xxp*xxp + yyp*yyp ) ;
147 pp =
acos( xxp / rr ) ;
150 dval_x = dent(0).val_point(rr, M_PI/2., pp) ;
159 while( diff_dent_x > 1.e-15 ) {
162 xtmp = xxp + ss * mm * dx ;
164 rr =
sqrt( xtmp*xtmp + yyp*yyp ) ;
175 pp =
acos( xtmp / rr ) ;
178 diff_dent_x = ss * dent(0).val_point(rr, M_PI/2., pp) ;
181 xxp += ss * (mm - 1) * dx ;
Cmp sqrt(const Cmp &)
Square root.
void ent_max_search(double &xx, double &yy) const
Searches the position of the maximum enthalpy.
Tenseur ent
Log-enthalpy (relativistic case) or specific enthalpy (Newtonian case)
Cmp acos(const Cmp &)
Arccosine.
Tensor handling *** DEPRECATED : use class Tensor instead ***.
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates)