23 char bhole_pseudo_viriel_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Bhole_binaire/bhole_pseudo_viriel.C,v 1.5 2014/10/13 08:52:40 j_novak Exp $" ;
86 double erreur = (2*(*devel_psi[1])(nz-1, 0, 0, 0)
87 + (*devel_n[1])(nz-1, 0, 0, 0))/fabs ((*devel_n[1])(nz-1, 0, 0, 0)) ;
95 int nz_un = hole1.mp.get_mg()->get_nzone() ;
96 int nz_deux = hole2.mp.get_mg()->get_nzone() ;
98 Valeur** devel_psi_un (hole1.psi_auto().asymptot(1)) ;
99 Valeur** devel_psi_deux (hole2.psi_auto().asymptot(1)) ;
100 Valeur** devel_n_un (hole1.n_auto().asymptot(1)) ;
101 Valeur** devel_n_deux (hole2.n_auto().asymptot(1)) ;
104 (2*(*devel_psi_un[1])(nz_un-1, 0, 0, 0)+
105 2*(*devel_psi_deux[1])(nz_deux-1, 0, 0, 0)+
106 (*devel_n_deux[1])(nz_deux-1, 0, 0, 0) +
107 (*devel_n_un[1])(nz_un-1, 0, 0, 0))
108 / fabs ((*devel_n_deux[1])(nz_deux-1, 0, 0, 0) +
109 (*devel_n_un[1])(nz_un-1, 0, 0, 0)) ;
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tenseur psi_auto
Part of generated by the hole.
Values and coefficients of a (real-value) function.
double viriel_seul() const
Computes the viriel error, that is the difference between the ADM and the Komar masses, calculated by the asymptotic behaviours of respectively and N .
Map_af & mp
Affine mapping.
int get_nzone() const
Returns the number of domains.
double viriel() const
Computes the viriel error, that is the difference between the ADM and the Komar masses, calculated by the asymptotic behaviours of respectively and N .
Tenseur n_auto
Part of N generated by the hole.