Returns the type of sampling in the direction: SYM : : symmetry with respect to the equatorial plane NONSYM : : no symmetry with respect to the equatorial plane.
Returns the type of sampling in the direction: SYM : : symmetry with respect to the transformation NONSYM : : no symmetry with respect to the transformation
A multi-domain grid is a set of 3D grids for the implementation of the multi-domain spectral method described in Bonazzola, Gourgoulhon & Marck, Phys. Rev. D58, 104020 (1998). Each domain is represented by a 3D mono grid (Grille3d) and is called a zone. For each direction, the Number of Degrees of Freedom (NDF) is apriori independent of the zone. However, some methods or routines may refuse to work if the NDF of some direction is not identical in all the zones. This holds for the type of sampling (symmetry) too.
[input] Array (size: nz ) of number of degree of freedom (NDF) in r-direction
typr[]
[input] Array (size: nz ) of type of sampling in r -direction
nbt[]
[input] Array (size: nz ) of NDF in -direction
typt
[input] Type of sampling in -direction
nbp[]
[input] Array (size: nz) of NDF in -direction
typp
[input] Type of sampling in -direction
base_r
[input] Types of r bases in each domain, to define the collocation points: BASE_CHEB , BASE_LEG or BASE_JAC02. If the pointer is null, BASE_CHEB is set.
This provides a multi-grid with the same number of degrees of freedom in all the domains.
The domain of index l = 0 will be a nucleus: , rarefied sampling (type RARE) near the origin;
domains of indices will be shells: , dense sampling (type FIN) near -1 and 1;
if compact == true, the domains of index l = nz-1 will be the outermost compactified shell: , dense sampling (type UNSURR) near -1 and 1 for a 1/r discretization.
Parameters
nz
[input] Number of domains (zones).
nbr
[input] Number of degree of freedom (NDF) in r -direction in each domain
nbt
[input] Number of degree of freedom (NDF) in -direction in each domain
nbp
[input] Number of degree of freedom (NDF) in -direction in each domain
typt
[input] Type of sampling in -direction: SYM for a sampling in (symmetry with respect to the equatorial plane), NONSYM for a sampling in
typp
[input] Type of sampling in -direction: SYM for a sampling in (symmetry with respect to a translation in ) NONSYM for a sampling in
compact
[input] true for the last domain to have a 1/r sampling (UNSURR ) instead of a r sampling (FIN ).
legendre
[input] true if the nucleus and shells have Legendre-type collocation points. If false, Chebyshev ones shall be used.
Simplified constructor for a standard multi-grid only with shells and Chebyshev.
This provides a multi-grid with the same number of degrees of freedom in all the domains.
ALL DOMAINS ARE TREATED AS SHELLS , dense sampling (type FIN) near -1 and 1;
and only Chebyshev bases are considered.
Parameters
nz
[input] Number of domains (zones).
nbr
[input] Number of degree of freedom (NDF) in r -direction in each domain
nbt
[input] Number of degree of freedom (NDF) in -direction in each domain
nbp
[input] Number of degree of freedom (NDF) in -direction in each domain
typt
[input] Type of sampling in -direction: SYM for a sampling in (symmetry with respect to the equatorial plane), NONSYM for a sampling in
typp
[input] Type of sampling in -direction: SYM for a sampling in (symmetry with respect to a translation in ) NONSYM for a sampling in
If the boolean flag read_base is false (default) the spectral basis is not saved, for compatibility with older versions. If it is set to true , the basis is saved.
Returns the type of sampling in the direction: SYM : : symmetry with respect to the transformation NONSYM : : no symmetry with respect to the transformation
Returns the type of sampling in the direction: SYM : : symmetry with respect to the equatorial plane NONSYM : : no symmetry with respect to the equatorial plane.
By default, if save_base is false, the spectral decomposition basis colloc_r is not saved (for compatibility reasons). If save_base is true, then colloc_r is save, too. The same value for the flag read_base must be given in the contructor reading the file.