[Commits] [svn:einsteintoolkit] NullSHRExtract/branches/tapir/src/ (Rev. 32)
bela at caltech.edu
bela at caltech.edu
Fri Nov 15 14:09:39 CST 2013
User: szilagyi
Date: 2013/11/15 02:09 PM
Modified:
/branches/tapir/src/
NullSHRE_WTSphHarm.F90
Log:
add more diagnostic output (time-deriv of ell)
File Changes:
Directory: /branches/tapir/src/
===============================
File [modified]: NullSHRE_WTSphHarm.F90
Delta lines: +136 -103
===================================================================
--- branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-11-15 20:03:03 UTC (rev 31)
+++ branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-11-15 20:09:38 UTC (rev 32)
@@ -72,12 +72,12 @@
CCTK_REAL, dimension(null_lsh(1), null_lsh(2)) :: fx_s,fy_s,fz_s
CCTK_REAL, dimension(null_lsh(1), null_lsh(2)) :: &
- w_1_n, w_2_n, w_3_n, w_4_n, w_5_n, &
- w_1_s, w_2_s, w_3_s, w_4_s, w_5_s
+ w_1_n, w_2_n, w_3_n, w_4_n, w_5_n, &
+ w_1_s, w_2_s, w_3_s, w_4_s, w_5_s
CCTK_COMPLEX, dimension(null_lsh(1), null_lsh(2)) :: &
- ul_1_n, ul_2_n, ul_3_n, ul_4_n, ul_5_n, &
- ul_1_s, ul_2_s, ul_3_s, ul_4_s, ul_5_s
+ ul_1_n, ul_2_n, ul_3_n, ul_4_n, ul_5_n, &
+ ul_1_s, ul_2_s, ul_3_s, ul_4_s, ul_5_s
CCTK_REAL, dimension(null_lsh(1),null_lsh(2),4),target :: tmp1
CCTK_REAL, dimension(null_lsh(1),null_lsh(2),6),target :: &
@@ -88,7 +88,9 @@
logical, parameter :: output_lapse = .false.
logical, parameter :: output_shift = .false.
logical, parameter :: output_dgij = .false.
- logical, parameter :: output_ell = .false.
+ logical, parameter :: output_ell = .true.
+ logical, parameter :: output_w0parts = .false.
+ logical, parameter :: output_u1parts = .false.
CCTK_INT :: a, b
angular_tensor_1_n(1)%d => tmp1(:,:,1)
@@ -420,19 +422,7 @@
call NullDecomp_WriteCoefsForRealFunc('ell_t_wt', cctkGH, truncate, &
null_lsh, zeta, ell_n(4)%d, ell_s(4)%d, cctk_time)
- call NullDecomp_WriteCoefsForRealFunc('dell_x_dt_wt', cctkGH, truncate, &
- null_lsh, zeta, j1_n(1,4)%d, j1_s(1,4)%d, cctk_time)
- call NullDecomp_WriteCoefsForRealFunc('dell_y_dt_wt', cctkGH, truncate, &
- null_lsh, zeta, j1_n(2,4)%d, j1_s(2,4)%d, cctk_time)
-
- call NullDecomp_WriteCoefsForRealFunc('dell_z_dt_wt', cctkGH, truncate, &
- null_lsh, zeta, j1_n(3,4)%d, j1_s(3,4)%d, cctk_time)
-
- call NullDecomp_WriteCoefsForRealFunc('dell_t_dt_wt', cctkGH, truncate, &
- null_lsh, zeta, j1_n(4,4)%d, j1_s(4,4)%d, cctk_time)
-
-
call NullDecomp_WriteCoefsForRealFunc('s_x_wt', cctkGH, truncate, &
null_lsh, zeta, sa_n(1)%d, sa_s(1)%d, cctk_time)
@@ -458,6 +448,46 @@
call NullDecomp_WriteCoefsForRealFunc('elldenom_wt', cctkGH, truncate, &
null_lsh, zeta, elld_n%d, elld_s%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dell_x_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, j1_n(1,4)%d, j1_s(1,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dell_y_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, j1_n(2,4)%d, j1_s(2,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dell_z_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, j1_n(3,4)%d, j1_s(3,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dell_t_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, j1_n(4,4)%d, j1_s(4,4)%d, cctk_time)
+
+
+ call NullDecomp_WriteCoefsForRealFunc('ds_x_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, dsa_n(1,4)%d, dsa_s(1,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('ds_y_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, dsa_n(2,4)%d, dsa_s(2,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('ds_z_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, dsa_n(3,4)%d, dsa_s(3,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dn_x_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, dna_n(1,4)%d, dna_s(1,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dn_y_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, dna_n(2,4)%d, dna_s(2,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dn_z_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, dna_n(3,4)%d, dna_s(3,4)%d, cctk_time)
+
+ call NullDecomp_WriteCoefsForRealFunc('dn_t_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, dna_n(4,4)%d, dna_s(4,4)%d, cctk_time)
+
+
+ call NullDecomp_WriteCoefsForRealFunc('delldenom_dt_wt', cctkGH, truncate, &
+ null_lsh, zeta, delld_n(4)%d, delld_s(4)%d, cctk_time)
+
+
end if
! xb, rb, r_l, r_t
@@ -488,10 +518,10 @@
null_lsh, zeta, etaup0_n(1,4)%d, etaup0_s(1,4)%d, cctk_time)
call NullDecomp_WriteCoefsForComplexFunc(&
- 'AffineInvMetric_rA_qA', cctkGH, truncate, null_lsh, 1, zeta, &
- etaup0_n(1,2)%d*qa(2)%d + etaup0_n(1,3)%d*qa(3)%d,&
- etaup0_s(1,2)%d*qa(2)%d + etaup0_s(1,3)%d*qa(3)%d,&
- cctk_time)
+ 'AffineInvMetric_rA_qA', cctkGH, truncate, null_lsh, 1, zeta, &
+ etaup0_n(1,2)%d*qa(2)%d + etaup0_n(1,3)%d*qa(3)%d,&
+ etaup0_s(1,2)%d*qa(2)%d + etaup0_s(1,3)%d*qa(3)%d,&
+ cctk_time)
! J, beta, U, Q, W
@@ -524,113 +554,116 @@
! bits of W
+ if(output_w0parts) then
+ w_1_n = dr0_n(1)%d * etaup0_n(1,1)%d / r0_n%d
+ w_2_n = - 2.d0 * dr0_n(4)%d / r0_n%d
+ w_3_n = &
+ + 2.d0 * ( dr0_n(2)%d * etaup0_n(1,2)%d &
+ + dr0_n(3)%d * etaup0_n(1,3)%d ) / r0_n%d
+ w_4_n = &
+ + ( dr0_n(2)%d ** 2 * etaup0_n(2,2)%d &
+ + 2.d0 * dr0_n(2)%d * dr0_n(3)%d * etaup0_n(2,3)%d &
+ + dr0_n(3)%d ** 2 * etaup0_n(3,3)%d &
+ ) / (r0_n%d*dr0_n(1)%d)
+ w_5_n = - 1.d0 / r0_n%d
- w_1_n = dr0_n(1)%d * etaup0_n(1,1)%d / r0_n%d
- w_2_n = - 2.d0 * dr0_n(4)%d / r0_n%d
- w_3_n = &
- + 2.d0 * ( dr0_n(2)%d * etaup0_n(1,2)%d &
- + dr0_n(3)%d * etaup0_n(1,3)%d ) / r0_n%d
- w_4_n = &
- + ( dr0_n(2)%d ** 2 * etaup0_n(2,2)%d &
- + 2.d0 * dr0_n(2)%d * dr0_n(3)%d * etaup0_n(2,3)%d &
- + dr0_n(3)%d ** 2 * etaup0_n(3,3)%d &
- ) / (r0_n%d*dr0_n(1)%d)
- w_5_n = - 1.d0 / r0_n%d
+ w_1_s = dr0_s(1)%d * etaup0_s(1,1)%d / r0_s%d
+ w_2_s = - 2.d0 * dr0_s(4)%d / r0_s%d
+ w_3_s = &
+ + 2.d0 * ( dr0_s(2)%d * etaup0_s(1,2)%d &
+ + dr0_s(3)%d * etaup0_s(1,3)%d ) / r0_s%d
+ w_4_s = &
+ + ( dr0_s(2)%d ** 2 * etaup0_s(2,2)%d &
+ + 2.d0 * dr0_s(2)%d * dr0_s(3)%d * etaup0_s(2,3)%d &
+ + dr0_s(3)%d ** 2 * etaup0_s(3,3)%d &
+ ) / (r0_s%d*dr0_s(1)%d)
+ w_5_s = - 1.d0 / r0_s%d
- w_1_s = dr0_s(1)%d * etaup0_s(1,1)%d / r0_s%d
- w_2_s = - 2.d0 * dr0_s(4)%d / r0_s%d
- w_3_s = &
- + 2.d0 * ( dr0_s(2)%d * etaup0_s(1,2)%d &
- + dr0_s(3)%d * etaup0_s(1,3)%d ) / r0_s%d
- w_4_s = &
- + ( dr0_s(2)%d ** 2 * etaup0_s(2,2)%d &
- + 2.d0 * dr0_s(2)%d * dr0_s(3)%d * etaup0_s(2,3)%d &
- + dr0_s(3)%d ** 2 * etaup0_s(3,3)%d &
- ) / (r0_s%d*dr0_s(1)%d)
- w_5_s = - 1.d0 / r0_s%d
+ call NullDecomp_WriteCoefsForRealFunc(&
+ 'W_part_rl_etaup_ll_over_r', cctkGH, truncate, &
+ null_lsh, zeta, w_1_n, w_1_s, cctk_time)
- call NullDecomp_WriteCoefsForRealFunc(&
- 'W_part_rl_etaup_ll_over_r', cctkGH, truncate, &
- null_lsh, zeta, w_1_n, w_1_s, cctk_time)
+ call NullDecomp_WriteCoefsForRealFunc(&
+ 'W_part_negative_2ru_over_r', cctkGH, truncate, &
+ null_lsh, zeta, w_2_n, w_2_s, cctk_time)
- call NullDecomp_WriteCoefsForRealFunc(&
- 'W_part_negative_2ru_over_r', cctkGH, truncate, &
- null_lsh, zeta, w_2_n, w_2_s, cctk_time)
+ call NullDecomp_WriteCoefsForRealFunc(&
+ 'W_part_2rA_etaup_lA_over_r', cctkGH, truncate, &
+ null_lsh, zeta, w_3_n, w_3_s, cctk_time)
- call NullDecomp_WriteCoefsForRealFunc(&
- 'W_part_2rA_etaup_lA_over_r', cctkGH, truncate, &
- null_lsh, zeta, w_3_n, w_3_s, cctk_time)
+ call NullDecomp_WriteCoefsForRealFunc(&
+ 'W_part_rA_rB_etaup_AB_over_rl_over_r', cctkGH, truncate, &
+ null_lsh, zeta, w_4_n, w_4_s, cctk_time)
- call NullDecomp_WriteCoefsForRealFunc(&
- 'W_part_rA_rB_etaup_AB_over_rl_over_r', cctkGH, truncate, &
- null_lsh, zeta, w_4_n, w_4_s, cctk_time)
+ call NullDecomp_WriteCoefsForRealFunc(&
+ 'W_part_negative_one_over_r', cctkGH, truncate, &
+ null_lsh, zeta, w_5_n, w_5_s, cctk_time)
- call NullDecomp_WriteCoefsForRealFunc(&
- 'W_part_negative_one_over_r', cctkGH, truncate, &
- null_lsh, zeta, w_5_n, w_5_s, cctk_time)
+ end if
+ ! write (*,*) 'diff_n = ' , &
+ ! maxval(abs( -w_wt_n%d + w_1_n + w_2_n + w_3_n + w_4_n + w_5_n))
-! write (*,*) 'diff_n = ' , &
-! maxval(abs( -w_wt_n%d + w_1_n + w_2_n + w_3_n + w_4_n + w_5_n))
-
-! write (*,*) 'diff_s = ' , &
-! maxval(abs( -w_wt_s%d + w_1_s + w_2_s + w_3_s + w_4_s + w_5_s))
-
+ ! write (*,*) 'diff_s = ' , &
+ ! maxval(abs( -w_wt_s%d + w_1_s + w_2_s + w_3_s + w_4_s + w_5_s))
+
! bits of U_{,\lambda}
+ if(output_w0parts) then
- ul_1_s = - etaup1_s(1,2)%d * qa(2)%d - etaup1_s(1,3)%d * qa(3)%d
- ul_1_n = - etaup1_n(1,2)%d * qa(2)%d - etaup1_n(1,3)%d * qa(3)%d
+ ul_1_s = - etaup1_s(1,2)%d * qa(2)%d - etaup1_s(1,3)%d * qa(3)%d
+ ul_1_n = - etaup1_n(1,2)%d * qa(2)%d - etaup1_n(1,3)%d * qa(3)%d
- call NullDecomp_WriteCoefsForComplexFunc(&
- 'Ul_part_negative_detaup_dlambda_lA_qA', cctkGH, truncate, &
- null_lsh, 1, zeta, ul_1_n, ul_1_s, cctk_time)
+ call NullDecomp_WriteCoefsForComplexFunc(&
+ 'Ul_part_negative_detaup_dlambda_lA_qA', cctkGH, truncate, &
+ null_lsh, 1, zeta, ul_1_n, ul_1_s, cctk_time)
- ul_2_s = 0
- ul_2_n = 0
- ul_3_s = 0
- ul_3_n = 0
- do a = 2, 3
- do b = 2, 3
- ul_2_s = ul_2_s - dr1_s(b)%d * etaup0_s(a,b)%d * qa(a)%d / dr0_s(1)%d
- ul_2_n = ul_2_n - dr1_n(b)%d * etaup0_n(a,b)%d * qa(a)%d / dr0_n(1)%d
- ul_3_s = ul_3_s - dr0_s(b)%d * etaup1_s(a,b)%d * qa(a)%d / dr0_s(1)%d
- ul_3_n = ul_3_n - dr0_n(b)%d * etaup1_n(a,b)%d * qa(a)%d / dr0_n(1)%d
- end do
- end do
+ ul_2_s = 0
+ ul_2_n = 0
+ ul_3_s = 0
+ ul_3_n = 0
+ do a = 2, 3
+ do b = 2, 3
+ ul_2_s = ul_2_s - dr1_s(b)%d * etaup0_s(a,b)%d * qa(a)%d / dr0_s(1)%d
+ ul_2_n = ul_2_n - dr1_n(b)%d * etaup0_n(a,b)%d * qa(a)%d / dr0_n(1)%d
+ ul_3_s = ul_3_s - dr0_s(b)%d * etaup1_s(a,b)%d * qa(a)%d / dr0_s(1)%d
+ ul_3_n = ul_3_n - dr0_n(b)%d * etaup1_n(a,b)%d * qa(a)%d / dr0_n(1)%d
+ end do
+ end do
- call NullDecomp_WriteCoefsForComplexFunc(&
- 'Ul_part_negative_etaup_AB_r1_B_qA_over_r1', cctkGH, truncate, &
- null_lsh, 1, zeta, ul_2_n, ul_2_s, cctk_time)
+ call NullDecomp_WriteCoefsForComplexFunc(&
+ 'Ul_part_negative_etaup_AB_r1_B_qA_over_r1', cctkGH, truncate, &
+ null_lsh, 1, zeta, ul_2_n, ul_2_s, cctk_time)
- call NullDecomp_WriteCoefsForComplexFunc(&
- 'Ul_part_negative_detaup_dlambda_AB_r0_B_qA_over_r1', cctkGH, truncate, &
- null_lsh, 1, zeta, ul_3_n, ul_3_s, cctk_time)
+ call NullDecomp_WriteCoefsForComplexFunc(&
+ 'Ul_part_negative_detaup_dlambda_AB_r0_B_qA_over_r1', cctkGH, truncate, &
+ null_lsh, 1, zeta, ul_3_n, ul_3_s, cctk_time)
- ul_4_s = 2.d0 * beta_l_s%d * u_wt_s%d
- ul_4_n = 2.d0 * beta_l_n%d * u_wt_n%d
+ ul_4_s = 2.d0 * beta_l_s%d * u_wt_s%d
+ ul_4_n = 2.d0 * beta_l_n%d * u_wt_n%d
- call NullDecomp_WriteCoefsForComplexFunc(&
- 'Ul_part_twice_beta_l_U', cctkGH, truncate, &
- null_lsh, 1, zeta, ul_4_n, ul_4_s, cctk_time)
+ call NullDecomp_WriteCoefsForComplexFunc(&
+ 'Ul_part_twice_beta_l_U', cctkGH, truncate, &
+ null_lsh, 1, zeta, ul_4_n, ul_4_s, cctk_time)
- ul_5_s = 2.d0 * beta_l_s%d * (etaup0_s(1,2)%d * qa(2)%d +&
- & etaup0_s(1,3)%d * qa(3)%d)
- ul_5_n = 2.d0 * beta_l_n%d * (etaup0_n(1,2)%d * qa(2)%d +&
- & etaup0_n(1,3)%d * qa(3)%d)
+ ul_5_s = 2.d0 * beta_l_s%d * (etaup0_s(1,2)%d * qa(2)%d +&
+ & etaup0_s(1,3)%d * qa(3)%d)
+ ul_5_n = 2.d0 * beta_l_n%d * (etaup0_n(1,2)%d * qa(2)%d +&
+ & etaup0_n(1,3)%d * qa(3)%d)
- call NullDecomp_WriteCoefsForComplexFunc(&
- 'Ul_part_twice_beta_l_etaup_lA_qA', cctkGH, truncate, &
- null_lsh, 1, zeta, ul_5_n, ul_5_s, cctk_time)
+ call NullDecomp_WriteCoefsForComplexFunc(&
+ 'Ul_part_twice_beta_l_etaup_lA_qA', cctkGH, truncate, &
+ null_lsh, 1, zeta, ul_5_n, ul_5_s, cctk_time)
+ end if
-! write (*,*) 'diff_ul = ' , &
-! maxval(abs( -u_l_n%d + ul_1_n + ul_2_n + ul_3_n + ul_4_n + ul_5_n)), &
-! maxval(abs( -u_l_s%d + ul_1_s + ul_2_s + ul_3_s + ul_4_s + ul_5_s))
-
+ ! write (*,*) 'diff_ul = ' , &
+ ! maxval(abs( -u_l_n%d + ul_1_n + ul_2_n + ul_3_n + ul_4_n + ul_5_n)), &
+ ! maxval(abs( -u_l_s%d + ul_1_s + ul_2_s + ul_3_s + ul_4_s + ul_5_s))
+
end subroutine NullSHRE_WTSphHarm
More information about the Commits
mailing list