57 #include "et_bin_bhns_extr.h" 58 #include "utilitaires.h" 78 double ppp = M_PI/2. ;
86 while ( dp > 1.e-15 ) {
92 diff_prev = ( dff.
val_point(0,1.,M_PI/2.,ppp)
96 - x_max *
cos(ppp) - y_max *
sin(ppp) )
99 * ( - x_max *
sin(ppp) + y_max *
cos(ppp) ) ;
101 if ( diff_prev > 0. ) {
108 while ( diff > 1.e-15 ) {
111 ptmp = ppp + mm * dp ;
113 diff = ss * ( ( dff.
val_point(0,1.,M_PI/2.,ptmp)
115 * ( 1. + ff0.
val_point(0,1.,M_PI/2.,ptmp)
117 - x_max *
cos(ptmp) - y_max *
sin(ptmp) )
118 - ( 1. + ff0.
val_point(0,1.,M_PI/2.,ptmp)
120 * ( - x_max *
sin(ptmp) + y_max *
cos(ptmp) ) ) ;
123 ppp += ss * (mm - 1) * dp ;
double phi_longest_rad(double x_max, double y_max) const
Searches the position (phi) of the longest radius of NS from the position of the maximum enthalpy...
Radial mapping of rather general form.
const Valeur & dsdp() const
Returns of *this.
const Valeur & get_gg() const
Returns a (constant) reference to the function .
Values and coefficients of a (real-value) function.
Cmp cos(const Cmp &)
Cosine.
const Valeur & get_ff() const
Returns a (constant) reference to the function .
double val_point(int l, double x, double theta, double phi) const
Computes the value of the field represented by *this at an arbitrary point, by means of the spectral ...
Map & mp
Mapping associated with the star.
Cmp sin(const Cmp &)
Sine.