31char strot_dirac_diff_faitomeg_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/strot_dirac_diff_faitomeg.C,v 1.2 2014/10/13 08:53:39 j_novak Exp $" ;
40#include "star_rot_dirac_diff.h"
41#include "utilitaires.h"
45double strot_dirac_diff_fzero(
double omeg,
const Param& par) ;
148double strot_dirac_diff_fzero(
double omeg,
const Param&
par) {
153 int l =
par.get_int(0) ;
154 int k =
par.get_int(1) ;
155 int j =
par.get_int(2) ;
156 int i =
par.get_int(3) ;
203 for (
int l=
nzet+1;
l<nz;
l++) {
Time evolution with partial storage (*** under development ***).
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
virtual const Sym_tensor & cov() const
Read-only access to the covariant representation.
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).
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.
Tbl & set_domain(int l)
Read/write of the value in a given domain.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
virtual void allocate_all()
Sets the logical state to ETATQCQ (ordinary state) and performs the memory allocation of all the elem...
const Tbl & domain(int l) const
Read-only of the value in a given domain.
Class for relativistic differentially rotating stars in Dirac gauge and maximal slicing.
void fait_omega_field(double omeg_min, double omeg_max, double precis, int nitermax)
Computes (member omega_field ).
Scalar omega_field
Field .
double funct_omega(double omeg) const
Evaluates , where F is the function defining the rotation profile.
void fait_prim_field()
Computes the member prim_field from omega_field .
double omega_min
Minimum value of .
double(* primfrot)(double, const Tbl &)
Primitive of the function , which vanishes at the stellar center.
Tbl par_frot
Parameters of the function .
double omega_max
Maximum value of .
Scalar psi4
Conformal factor .
Map & mp
Mapping associated with the star.
int nzet
Number of domains of *mp occupied by the star.
Tbl min(const Cmp &)
Minimum values of a Cmp in each domain.
Tbl max(const Cmp &)
Maximum values of a Cmp in each domain.
double zerosec(double(*f)(double, const Param &), const Param &par, double a, double b, double precis, int nitermax, int &niter, bool abort=true)
Finding the zero a function.