29char star_bhns_equilibrium_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star_bhns/star_bhns_equilibrium.C,v 1.4 2014/10/13 08:53:40 j_novak Exp $" ;
62#include "utilitaires.h"
68 const double&
sepa,
bool kerrschild,
205 for (
int i=0;
i<3;
i++) {
231 cout <<
"-----------------------------------------------" <<
endl ;
307 cout <<
"k_b, j_b, alpha_r: " <<
k_b <<
" " <<
j_b <<
" "
393 cout <<
"******* FROZEN MAPPING *********" <<
endl ;
438 for (
int i=1;
i<nz-4;
i++) {
492 cout <<
"Max. abs(enthalpy) at the boundary between domains nzet-1"
493 <<
" and nzet : " <<
endl ;
575 "Relative error in the resolution of the equation for lapconf_auto : "
577 for (
int l=0;
l<nz;
l++) {
634 "Relative error in the resolution of the equation for confo_auto : "
636 for (
int l=0;
l<nz;
l++) {
658 for (
int i=1;
i<=3;
i++) {
667 for (
int i=1;
i<=3;
i++) {
673 for (
int i=1;
i<=3;
i++) {
699 for (
int i=1;
i<=3;
i++) {
708 for (
int i=1;
i<=3;
i++) {
716 for (
int i=1;
i<=3;
i++) {
729 for (
int i=0;
i<3;
i++) {
735 for (
int i=0;
i<3 ;
i++) {
749 for (
int i=1;
i<=3;
i++)
754 for (
int i=0;
i<3;
i++) {
771 "Relative error in the resolution of the equation for shift_auto : "
773 cout <<
"x component : " ;
774 for (
int l=0;
l<nz;
l++) {
778 cout <<
"y component : " ;
779 for (
int l=0;
l<nz;
l++) {
783 cout <<
"z component : " ;
784 for (
int l=0;
l<nz;
l++) {
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Time evolution with partial storage (*** under development ***).
Radial mapping of rather general form.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord y
y coordinate centered on the grid
virtual void reevaluate(const Map *mp_prev, int nzet, Cmp &uu) const =0
Recomputes the values of a Cmp at the collocation points after a change in the mapping.
virtual void resize(int l, double lambda)=0
Rescales the outer boundary of one domain.
virtual void adapt(const Cmp &ent, const Param &par, int nbr=0)=0
Adaptation of the mapping to a given scalar field.
Coord x
x coordinate centered on the grid
virtual double val_r(int l, double xi, double theta, double pphi) const =0
Returns the value of the radial coordinate r for a given in a given domain.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Tensor field of valence 0 (or component of a tensorial field).
const Scalar & dsdy() const
Returns of *this , where .
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
const Scalar & dsdx() const
Returns of *this , where .
const Scalar & dsdr() const
Returns of *this .
Valeur & set_spectral_va()
Returns va (read/write version)
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star ).
Sym_tensor taij_auto
Part of the extrinsic curvature tensor generated by shift_auto , lapse_auto , and confo_auto .
Scalar lapconf_auto
Lapconf function generated by the star.
double velo_pot_bhns(const double &mass_bh, const double &sepa, bool kerrschild, int mermax, double precis, double relax)
Computes the non-translational part of the velocity scalar potential by solving the continuity equat...
void equilibrium_bhns(double ent_c, const double &mass_bh, const double &sepa, bool kerrschild, int mer, int mermax_ns, int mermax_potvit, int mermax_poisson, int filter_r, int filter_r_s, int filter_p_s, double relax_poisson, double relax_potvit, double thres_adapt, double resize_ns, const Tbl &fact_resize, Tbl &diff)
Computes an equilibrium configuration.
Vector ssjm1_wshift
Effective source at the previous step for the resolution of the vector Poisson equation for by means...
Scalar confo_tot
Total conformal factor.
Vector d_lapconf_auto
Derivative of the lapconf function generated by the star .
Vector d_confo_auto
Derivative of the conformal factor generated by the star .
double radius_p(double phi)
Radius of the star to the direction of and .
Scalar gam0
Lorentz factor between the co-orbiting observer and the Eulerian one.
void hydro_euler_bhns(bool kerrschild, const double &mass_bh, const double &sepa)
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
Scalar gam
Lorentz factor between the fluid and the co-orbiting observer.
Scalar ssjm1_lapconf
Effective source at the previous step for the resolution of the Poisson equation for lapconf_auto .
bool irrotational
true for an irrotational star, false for a corotating one
Scalar pot_centri
Centrifugal potential.
double chi_rp(double radius, double phi)
Sensitive indicator of the mass-shedding to the direction of , , .
Scalar lapconf_comp
Lapconf function generated by the companion black hole.
Scalar taij_quad_auto
Part of the scalar generated by .
double phi_min()
Azimuthal angle when the indicator of the mass-shedding takes its minimum chi_min.
Vector shift_auto
Shift vector generated by the star.
Scalar ssjm1_khi
Effective source at the previous step for the resolution of the Poisson equation for the scalar by m...
Scalar confo_auto
Conformal factor generated by the star.
Scalar ssjm1_confo
Effective source at the previous step for the resolution of the Poisson equation for confo_auto .
Scalar loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.
Scalar ener_euler
Total energy density in the Eulerian frame.
void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
Scalar s_euler
Trace of the stress scalar in the Eulerian frame.
Scalar press
Fluid pressure.
double ray_eq_pi() const
Coordinate radius at , [r_unit].
Vector u_euler
Fluid 3-velocity with respect to the Eulerian observer.
Map & mp
Mapping associated with the star.
int nzet
Number of domains of *mp occupied by the star.
double ray_pole() const
Coordinate radius at [r_unit].
Tensor handling *** DEPRECATED : use class Tensor instead ***.
void smooth(int nzet, Valeur &uuva) const
Changes the function *this as a smooth one when there exists a discontinuity between the nucleus and ...
Tensor field of valence 1.
const Scalar & divergence(const Metric &) const
The divergence of this with respect to a Metric .
Scalar & set(int)
Read/write access to a component.
Cmp sqrt(const Cmp &)
Square root.
Cmp exp(const Cmp &)
Exponential.
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
Tbl max(const Cmp &)
Maximum values of a Cmp in each domain.
Cmp pow(const Cmp &, int)
Power .
Cmp abs(const Cmp &)
Absolute value.
Cmp log(const Cmp &)
Neperian logarithm.
const Tensor & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
Standard units of space, time and mass.