class Bin_BH

Binary black hole configuration on a Cartesian grid.

Documentation

Binary black hole configuration on a Cartesian grid.

A binary black hole system is constructed on a Cartesian grid from data stored in a file resulting from a computation by Grandclement, Gourgoulhon and Bonazzola, Phys. Rev. D 65, 044021 (2002).

All the quantities are in units derived from the length scale defined by the coordinate radius of black hole 1 apparent horizon (throat).

Importation of Lorene data is performed by means of the constructor Bin_BH::Bin_BH(int, const double*, const double*, const double*, const char*). This constructor takes general arrays for the location of the Cartesian coordinates , i.e. it does not assume that the grid is a uniform one. Note also that these arrays are 1-D, as well as all the metric fields, in order to be use with any ordering of the 3-D storage.

This class is very simple, with all data members being public. A typical example of use is the following one

 *	    // Define the Cartesian grid by means of the arrays xg, yg, zg:
 *	    for (int i=0; i<nb_points; i++) {
 *           xg[i] = ...
 *           yg[i] = ...
 *           zg[i] = ...
 *	    }
 *
 *	    // Read the file containing the spectral data and evaluate
 *	    //  all the fields on the Cartesian grid :
 *
 *	    Bin_BH binary_system(nb_points, xg, yg, zg, fill, datafile) ;
 *
 *	    // Extract what you need :
 *
 *	    double* gamma_xx = binary_system.g_xx ; // metric coefficient g_xx
 *
 *	    double* shift_x = binary_system.beta_x ; // x comp. of shift vector
 *
 *	    ...
 *
 *	    // Save everything in an ASCII file :
 *
 *	    ofstream file_ini("ini.d") ;
 *	    binary_system.save_form(file_ini) ;
 *	    file_ini.close() ;
 *
 *  


Champs publics

[more]double omega
Orbital angular velocity [unit: ]
[more]double dist
Distance between the coordinate centers of two black holes [unit: ]
[more]double radius2
Coordinate radius of the apparent horizon (throat) of black hole 2 [unit: ].
[more]int np
Total number of grid points
[more]double* xx
1-D array storing the values of coordinate x of the np grid points [unit: ]
[more]double* yy
1-D array storing the values of coordinate y of the np grid points [unit: ]
[more]double* zz
1-D array storing the values of coordinate z of the np grid points [unit: ]
[more]double* nnn
Lapse function at the np grid points (1-D array)
[more]double* beta_x
Component of the shift vector of corotating coordinates [unit: ]
[more]double* beta_y
Component of the shift vector of corotating coordinates [unit: ]
[more]double* beta_z
Component of the shift vector of corotating coordinates [unit: ]
[more]double* g_xx
Metric coefficient at the grid points (1-D array)
[more]double* g_xy
Metric coefficient at the grid points (1-D array)
[more]double* g_xz
Metric coefficient at the grid points (1-D array)
[more]double* g_yy
Metric coefficient at the grid points (1-D array)
[more]double* g_yz
Metric coefficient at the grid points (1-D array)
[more]double* g_zz
Metric coefficient at the grid points (1-D array)
[more]double* k_xx
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]
[more]double* k_xy
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]
[more]double* k_xz
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]
[more]double* k_yy
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]
[more]double* k_yz
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]
[more]double* k_zz
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]
[more]double* dpsi_x
First derivative of the conformal factor [unit: ]
[more]double* dpsi_y
First derivative of the conformal factor [unit: ]
[more]double* dpsi_z
First derivative of the conformal factor [unit: ]
[more]double* d2psi_xx
Second derivative of the conformal factor [unit: ]
[more]double* d2psi_xy
Second derivative of the conformal factor [unit: ]
[more]double* d2psi_xz
Second derivative of the conformal factor [unit: ]
[more]double* d2psi_yy
Second derivative of the conformal factor [unit: ]
[more]double* d2psi_yz
Second derivative of the conformal factor [unit: ]
[more]double* d2psi_zz
Second derivative of the conformal factor [unit: ]

Méthodes publiques

[more] Bin_BH(int nbpoints, const double* xi, const double* yi, const double* zi, int fill, const char* filename, bool mdiff=false)
Constructor from Lorene spectral data.
[more] Bin_BH(FILE* )
Constructor from a binary file (previously created by save_bin)
[more] Bin_BH(ifstream& )
Constructor from a formatted file (previously created by save_form)
[more] ~Bin_BH()
Destructor
[more]void save_bin(FILE* ) const
Save in a binary file.
[more]void save_form(ofstream& ) const
Save in a formatted file.

Méthodes privées

[more]void alloc_memory()
Allocate the memory for the arrays g_ij, k_ij, etc

odouble omega
Orbital angular velocity [unit: ]

odouble dist
Distance between the coordinate centers of two black holes [unit: ]

odouble radius2
Coordinate radius of the apparent horizon (throat) of black hole 2 [unit: ]. NB: The coordinate radius of black hole 1 is 1 by definition of the length unit.

oint np
Total number of grid points

odouble* xx
1-D array storing the values of coordinate x of the np grid points [unit: ]

odouble* yy
1-D array storing the values of coordinate y of the np grid points [unit: ]

odouble* zz
1-D array storing the values of coordinate z of the np grid points [unit: ]

odouble* nnn
Lapse function at the np grid points (1-D array)

odouble* beta_x
Component of the shift vector of corotating coordinates [unit: ]

odouble* beta_y
Component of the shift vector of corotating coordinates [unit: ]

odouble* beta_z
Component of the shift vector of corotating coordinates [unit: ]

odouble* g_xx
Metric coefficient at the grid points (1-D array)

odouble* g_xy
Metric coefficient at the grid points (1-D array)

odouble* g_xz
Metric coefficient at the grid points (1-D array)

odouble* g_yy
Metric coefficient at the grid points (1-D array)

odouble* g_yz
Metric coefficient at the grid points (1-D array)

odouble* g_zz
Metric coefficient at the grid points (1-D array)

odouble* k_xx
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]

odouble* k_xy
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]

odouble* k_xz
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]

odouble* k_yy
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]

odouble* k_yz
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]

odouble* k_zz
Component of the extrinsic curvature at the grid points (1-D array) [unit: ]

odouble* dpsi_x
First derivative of the conformal factor [unit: ]

odouble* dpsi_y
First derivative of the conformal factor [unit: ]

odouble* dpsi_z
First derivative of the conformal factor [unit: ]

odouble* d2psi_xx
Second derivative of the conformal factor [unit: ]

odouble* d2psi_xy
Second derivative of the conformal factor [unit: ]

odouble* d2psi_xz
Second derivative of the conformal factor [unit: ]

odouble* d2psi_yy
Second derivative of the conformal factor [unit: ]

odouble* d2psi_yz
Second derivative of the conformal factor [unit: ]

odouble* d2psi_zz
Second derivative of the conformal factor [unit: ]

o Bin_BH(int nbpoints, const double* xi, const double* yi, const double* zi, int fill, const char* filename, bool mdiff=false)
Constructor from Lorene spectral data.

This constructor takes general arrays xi, yi, zi for the location of the Cartesian coordinates , i.e. it does not assume that the grid is a uniform one. These arrays are 1-D to deal with any ordering of a 3-D storage.

Paramètres :
nbpoints - [input] Total number of grid points
xi - [input] 1-D array (size nbpoints) storing the values of coordinate x of the grid points [unit: ]
yi - [input] 1-D array (size nbpoints) storing the values of coordinate y of the grid points [unit: ]
zi - [input] 1-D array (size nbpoints) storing the values of coordinate z of the grid points [unit: ]
fill - [input] sets how the hole "interiors" must be filled:
fill = 0 : all the fields are set to zero
fill = 1 : the fields are extrapolated from their values "outside" the holes, by means of parabolas along radial directions
filename - [input] Name of the (binary) file containing the result of a computation by means of the multi-domain spectral method.

o Bin_BH(FILE* )
Constructor from a binary file (previously created by save_bin)

o Bin_BH(ifstream& )
Constructor from a formatted file (previously created by save_form)

o ~Bin_BH()
Destructor

ovoid alloc_memory()
Allocate the memory for the arrays g_ij, k_ij, etc

ovoid save_bin(FILE* ) const
Save in a binary file. This file can be subsenquently read by the evolution code, or by the constructor Bin_BH::Bin_BH(FILE* ).

ovoid save_form(ofstream& ) const
Save in a formatted file. This file can be subsenquently read by the evolution code, or by the constructor Bin_BH::Bin_BH(ifstream& ).


Cette classe n'a pas de classe descendante.
Friends:
ostream& operator<<(ostream& , const Bin_BH& )
Version :
$Id: Bin_BH.html,v 1.2 2010/07/14 16:56:37 e_gourgoulhon Exp $

Index alphabétique Hiérarchie des classes HTML ou Java



This page was generated with the help of DOC++.