30 char et_bin_bhns_extr_excurv_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_bhns_extr_excurv.C,v 1.4 2014/10/13 08:52:55 j_novak Exp $" ;
56 #include "et_bin_bhns_extr.h"
104 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
109 xx_con.
set(0) = xx + sepa ;
115 xsr_con = xx_con / r_bh ;
119 msr = ggrav * mass / r_bh ;
123 lapse_bh2 = 1. / (1.+2.*msr) ;
138 for (
int i=0; i<3; i++) {
139 for (
int j=0; j<3; j++) {
140 tmp1.
set(i, j) = -2.*lapse_bh2()%msr()%xsr_con(i)%xsr_con(j) ;
142 tmp2.
set(i, j) = -3.*lapse_bh2()%xsr_con(i)%xsr_con(j)
143 -4.*lapse_bh2()*msr()%xsr_con(i)%xsr_con(j) ;
145 tmp3.
set(i, j) = xsr_con(i)%shift_auto_local(j) ;
158 for (
int i=0; i<3; i++)
159 tmp4.
set() += xsr_con(i) % shift_auto_local(i) ;
172 for (
int i=0; i<3; i++) {
173 for (
int j=i; j<3; j++) {
175 + tmp1dn(i, j) + tmp1dn(j, i)
176 + 2.*lapse_bh2()%msr()/r_bh()%( tmp3(i, j) + tmp3(j, i)
177 + tmp4() % tmp2(i, j) )
178 -
double(2)/double(3) * tmp1(i, j)
179 * (divn() - lapse_bh2() % msr() / r_bh() % tmp4()) ;
182 * (divn() - lapse_bh2() % msr() / r_bh() % tmp4()) ;
194 xx_cov.
set(0) = xx + sepa ;
200 xsr_cov = xx_cov / r_bh ;
208 for (
int i=0; i<3; i++) {
209 for (
int j=0; j<3; j++) {
215 for (
int i=0; i<3; i++) {
216 for (
int j=0; j<3; j++) {
217 tmp5.
set(i, j) = 2.*msr()%xsr_cov(i)%xsr_cov(j) ;
221 for (
int i=0; i<3; i++) {
222 for (
int j=0; j<3; j++) {
223 for (
int k=0; k<3; k++) {
237 for (
int i=0; i<3; i++) {
238 for (
int j=0; j<3; j++) {
239 tmp7.
set() += tmp5(i,j) % tmp6(i,j) ;
249 for (
int i=0; i<3; i++) {
250 for (
int j=0; j<3; j++) {
260 for (
int i=0; i<3; i++) {
261 for (
int j=0; j<3; j++) {
304 for (
int i=0; i<3; i++) {
305 for (
int j=i; j<3; j++) {
322 for (
int i=0; i<3; i++) {
323 for (
int j=0; j<3; j++) {