31char et_bin_bhns_extr_phi_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_bhns_extr_phi.C,v 1.3 2014/10/13 08:52:55 j_novak Exp $" ;
54#include "et_bin_bhns_extr.h"
55#include "utilitaires.h"
75 double ppp = M_PI/2. ;
83 while ( dp > 1.e-15 ) {
89 diff_prev = ( dff.
val_point(0,1.,M_PI/2.,ppp)
93 - x_max *
cos(ppp) - y_max *
sin(ppp) )
96 * ( - x_max *
sin(ppp) + y_max *
cos(ppp) ) ;
98 if ( diff_prev > 0. ) {
105 while ( diff > 1.e-15 ) {
108 ptmp = ppp + mm * dp ;
110 diff = ss * ( ( dff.
val_point(0,1.,M_PI/2.,ptmp)
112 * ( 1. + ff0.
val_point(0,1.,M_PI/2.,ptmp)
114 - x_max *
cos(ptmp) - y_max *
sin(ptmp) )
115 - ( 1. + ff0.
val_point(0,1.,M_PI/2.,ptmp)
117 * ( - x_max *
sin(ptmp) + y_max *
cos(ptmp) ) ) ;
120 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.
Map & mp
Mapping associated with the star.
Radial mapping of rather general form.
const Valeur & get_ff() const
Returns a (constant) reference to the function .
const Valeur & get_gg() const
Returns a (constant) reference to the function .
Values and coefficients of a (real-value) 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 ...
const Valeur & dsdp() const
Returns of *this.
Cmp sin(const Cmp &)
Sine.
Cmp cos(const Cmp &)
Cosine.