30 char bhole_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Bhole/bhole.C,v 1.15 2014/10/13 08:52:39 j_novak Exp $" ;
171 #include "utilitaires.h"
172 #include "et_bin_nsbh.h"
173 #include "graphique.h"
178 rayon ((mp.get_alpha())[0]),
179 omega(0), omega_local(0), rot_state(COROT), boost (new double[3]), regul(0),
180 n_auto(mpi), n_comp(mpi), n_tot(mpi),
181 psi_auto(mpi), psi_comp(mpi), psi_tot(mpi),
182 grad_n_tot (mpi, 1, COV, mpi.get_bvect_cart()),
183 grad_psi_tot (mpi, 1, COV, mpi.get_bvect_cart()),
184 shift_auto(mpi, 1, CON, mpi.get_bvect_cart()),
185 taij_auto(mpi, 2, CON, mpi.get_bvect_cart()),
186 taij_comp(mpi, 2, CON, mpi.get_bvect_cart()),
187 taij_tot (mpi, 2, CON, mpi.get_bvect_cart()),
188 tkij_auto (mpi, 2, CON, mpi.get_bvect_cart()),
189 tkij_tot (mpi, 2, CON, mpi.get_bvect_cart()),
191 for (
int i=0 ; i<3 ; i++)
198 mp(source.mp), rayon(source.rayon), omega (source.omega), omega_local(source.omega_local), rot_state(source.rot_state),
199 boost (new double [3]), regul(source.regul), n_auto(source.n_auto),
200 n_comp(source.n_comp), n_tot(source.n_tot), psi_auto(source.psi_auto),
201 psi_comp(source.psi_comp), psi_tot(source.psi_tot),
202 grad_n_tot (source.grad_n_tot), grad_psi_tot (source.grad_psi_tot),
203 shift_auto(source.shift_auto),
204 taij_auto (source.taij_auto),
205 taij_comp(source.taij_comp), taij_tot(source.taij_tot),
206 tkij_auto(source.tkij_auto), tkij_tot(source.tkij_tot), decouple(source.decouple) {
208 for (
int i=0 ; i<3 ; i++)
220 assert (&
mp == &source.
mp) ;
226 for (
int i=0 ; i<3 ; i++)
386 trace = grad(0, 0)+grad(1, 1)+grad(2, 2) ;
390 for (
int i=0 ; i<3 ; i++) {
391 for (
int j=i+1 ; j<3 ; j++)
396 for (
int i=0 ; i<3 ; i++)
397 for (
int j=0 ; j<3 ; j++)
497 mp(mpi), rayon ((mp.get_alpha())[0]),
498 boost (new double[3]), n_auto(mpi), n_comp(mpi), n_tot(mpi),
499 psi_auto(mpi), psi_comp(mpi), psi_tot(mpi),
500 grad_n_tot (mpi, 1, COV, mpi.get_bvect_cart()),
501 grad_psi_tot (mpi, 1, COV, mpi.get_bvect_cart()),
502 shift_auto(mpi, 1, CON, mpi.get_bvect_cart()),
503 taij_auto(mpi, 2, CON, mpi.get_bvect_cart()),
504 taij_comp(mpi, 2, CON, mpi.get_bvect_cart()),
505 taij_tot (mpi, 2, CON, mpi.get_bvect_cart()),
506 tkij_auto (mpi, 2, CON, mpi.get_bvect_cart()) ,
507 tkij_tot (mpi, 2, CON, mpi.get_bvect_cart()) ,
560 double Bhole::local_momentum()
const {
571 vecteur.set_etat_qcq() ;
572 for (
int i=0 ; i<3 ; i++)
574 vecteur.set_std_base() ;
579 integrant.std_base_scal() ;