Commits November 2013

commits@lists.einsteintoolkit.org

10 participants
131 discussions

[svn:einsteintoolkit] NullSHRExtract/branches/tapir/src/ (Rev. 36)
by bela＠caltech.edu 18 Nov '13

18 Nov '13

User: szilagyi Date: 2013/11/18 04:58 PM Modified: /branches/tapir/src/ NullSHRE_InitialWT.F90 Log: remove unused symbol File Changes: Directory: /branches/tapir/src/ =============================== File [modified]: NullSHRE_InitialWT.F90 Delta lines: +0 -4 =================================================================== --- branches/tapir/src/NullSHRE_InitialWT.F90 2013-11-18 22:57:06 UTC (rev 35) +++ branches/tapir/src/NullSHRE_InitialWT.F90 2013-11-18 22:58:04 UTC (rev 36) @@ -34,10 +34,6 @@ DECLARE_CCTK_FUNCTIONS CCTK_REAL :: rl_min, elld_min -#define dt cctk_delta_time -#define dd null_delta(1) - - !--------------------------------------------------------------------------- ! world-tube, and analytic part of the coordinate transformation !---------------------------------------------------------------------------

1 0

[svn:einsteintoolkit] NullSHRExtract/branches/tapir/src/ (Rev. 35)
by bela＠caltech.edu 18 Nov '13

18 Nov '13

User: szilagyi Date: 2013/11/18 04:57 PM Modified: /branches/tapir/src/ NullSHRE_modEta.F90 Log: Fix a bug in d(elld)/dt. This commit is also present in the master branch File Changes: Directory: /branches/tapir/src/ =============================== File [modified]: NullSHRE_modEta.F90 Delta lines: +1 -1 =================================================================== --- branches/tapir/src/NullSHRE_modEta.F90 2013-11-16 01:10:11 UTC (rev 34) +++ branches/tapir/src/NullSHRE_modEta.F90 2013-11-18 22:57:06 UTC (rev 35) @@ -92,7 +92,7 @@ delld(4)%d = delld(4)%d & - g(i,j)%d * dbeta(i,4)%d * sa(j)%d& - g(i,j)%d * beta(i)%d * dsa(j,4)%d& - - beta(i)%d * sa(j)%d !*dg(i,j,4)%d + - beta(i)%d * sa(j)%d *dg(i,j,4)%d end do end do

1 0

[svn:einsteintoolkit] NullGrid/branches/tapir/ (Rev. 13)
by bela＠caltech.edu 18 Nov '13

18 Nov '13

User: szilagyi Date: 2013/11/18 04:29 PM Modified: /branches/tapir/ interface.ccl, schedule.ccl /branches/tapir/src/ NullGrid_InitCoord.F90 Log: add mechanism of saving the last few values of cctk_time File Changes: Directory: /branches/tapir/src/ =============================== File [modified]: NullGrid_InitCoord.F90 Delta lines: +24 -0 =================================================================== --- branches/tapir/src/NullGrid_InitCoord.F90 2013-10-22 18:26:54 UTC (rev 12) +++ branches/tapir/src/NullGrid_InitCoord.F90 2013-11-18 22:29:43 UTC (rev 13) @@ -33,7 +33,31 @@ end function NullGrid_RegisterCoords +subroutine NullGrid_InitTime(CCTK_ARGUMENTS) + use cctk + implicit none + + DECLARE_CCTK_ARGUMENTS + + null_time = cctk_time + null_time_p = cctk_time - cctk_delta_time + null_time_p_p = cctk_time - 2.0*cctk_delta_time + +end subroutine NullGrid_InitTime + + +subroutine NullGrid_CopyTime(CCTK_ARGUMENTS) + + use cctk + implicit none + + DECLARE_CCTK_ARGUMENTS + + null_time = cctk_time + +end subroutine NullGrid_CopyTime + subroutine NullGrid_Coord(CCTK_ARGUMENTS) use cctk Directory: /branches/tapir/ =========================== File [modified]: interface.ccl Delta lines: +5 -0 =================================================================== --- branches/tapir/interface.ccl 2013-10-22 18:26:54 UTC (rev 12) +++ branches/tapir/interface.ccl 2013-11-18 22:29:43 UTC (rev 13) @@ -23,6 +23,11 @@ null_dx, qsize } "scalars" +CCTK_REAL NullTime TYPE=SCALAR TIMELEVELS=3 +{ + null_time +} "histor of last 3 time-levels" + ################################################################## # INTEGER CONSTANTS DESCRIBING THE GRID ################################################################## File [modified]: schedule.ccl Delta lines: +12 -1 =================================================================== --- branches/tapir/schedule.ccl 2013-10-22 18:26:54 UTC (rev 12) +++ branches/tapir/schedule.ccl 2013-11-18 22:29:43 UTC (rev 13) @@ -2,7 +2,7 @@ # $Header$ storage: RadCrd NullSclr StCrd StCmCrd NullGrParArrI NullGrParArrR -storage: EG_mask, guard_mask, EQ_mask, EV_mask +storage: EG_mask, guard_mask, EQ_mask, EV_mask, NullTime[3] schedule NullGrid_RegisterCoords at STARTUP { @@ -28,3 +28,14 @@ OPTIONS: global } "setup variables" +schedule NullGrid_InitTime at CCTK_INITIAL +{ + LANG: Fortran + OPTIONS: global +} "initialize null time variable" + +schedule NullGrid_CopyTime at CCTK_PRESTEP +{ + LANG: Fortran + OPTIONS: global +} "copy value of time to local var"

1 0

[svn:einsteintoolkit] www/about/releases/ (Rev. 1074)
by rhaas＠tapir.caltech.edu 18 Nov '13

18 Nov '13

User: rhaas Date: 2013/11/18 11:46 AM Modified: /about/releases/ ET_2013_11_announcement.txt Log: add bugfix in NullSHRE_modEta.F90 (rev 34 of NullSHRExtract) File Changes: Directory: /about/releases/ =========================== File [modified]: ET_2013_11_announcement.txt Delta lines: +1 -0 =================================================================== --- about/releases/ET_2013_11_announcement.txt 2013-11-18 17:11:26 UTC (rev 1073) +++ about/releases/ET_2013_11_announcement.txt 2013-11-18 17:46:48 UTC (rev 1074) @@ -86,6 +86,7 @@ * TOVSolver: - remove dependency on GRHydro, in particular remove atmosphere handling - remove option to use different K and Gamma for different stars +* PITTNull: fix bug in computing L which affects CCE runs * Kranc: Too many improvements/rewrites to list here * Multipole: various improvements, especially convergence testing * LocalInterp: add 0th order interpolation

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 1073)
by rhaas＠tapir.caltech.edu 18 Nov '13

18 Nov '13

User: rhaas Date: 2013/11/18 11:11 AM Modified: /release-info/ lonestar__1_6.log, lonestar__2_6.log Log: update lonestar results File Changes: Directory: /release-info/ ========================= File [modified]: lonestar__1_6.log Delta lines: <Unavailable> File [modified]: lonestar__2_6.log Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/test/ (Rev. 588)
by rhaas＠tapir.caltech.edu 18 Nov '13

18 Nov '13

User: rhaas Date: 2013/11/18 11:03 AM Modified: /trunk/test/ tov_slowsector.par Log: GRHydro: remove misleading comment File Changes: Directory: /trunk/test/ ======================= File [modified]: tov_slowsector.par Delta lines: +0 -1 =================================================================== --- trunk/test/tov_slowsector.par 2013-11-15 17:09:36 UTC (rev 587) +++ trunk/test/tov_slowsector.par 2013-11-18 17:03:43 UTC (rev 588) @@ -76,7 +76,6 @@ GRHydro::tvd_limiter = "VanLeerMC2" GRHydro::rho_abs_min = 1.e-7 GRHydro::initial_atmosphere_factor = 0.999 -# these do not match TOVSolver to make the star actually move EOS_Omni::poly_gamma = 2.0 EOS_Omni::gl_gamma = 2.0 EOS_Omni::poly_k = 100

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 1072)
by bcmsma＠astro.rit.edu 17 Nov '13

17 Nov '13

User: bmundim Date: 2013/11/17 07:58 AM Modified: /release-info/ loewe__1_24.log, loewe__2_12.log Log: loewe results. File Changes: Directory: /release-info/ ========================= File [modified]: loewe__1_24.log Delta lines: +22 -58 =================================================================== --- release-info/loewe__1_24.log 2013-11-15 21:30:26 UTC (rev 1071) +++ release-info/loewe__1_24.log 2013-11-17 13:58:01 UTC (rev 1072) @@ -242,12 +242,12 @@ (LSUThorns/PeriodicCarpet/test/ml-gw1d-small-amr.par) Requires 2 processors - KerrSchild-rotating-180 in RotatingSymmetry180 - (CactusNumerical/RotatingSymmetry180/test/KerrSchild-rotating-180.par) + KerrSchild-rotating-180-EE in RotatingSymmetry180 + (CactusNumerical/RotatingSymmetry180/test/KerrSchild-rotating-180-EE.par) Requires 2 processors - KerrSchild-rotating-90 in RotatingSymmetry90 - (CactusNumerical/RotatingSymmetry90/test/KerrSchild-rotating-90.par) + KerrSchild-rotating-90-EE in RotatingSymmetry90 + (CactusNumerical/RotatingSymmetry90/test/KerrSchild-rotating-90-EE.par) Requires 2 processors @@ -265,12 +265,6 @@ Test created 4 extra files: Multipole (test_44) Test created 4 extra files: - RotatingSymmetry180 (Kerr-rotating-180) - Test created 4 extra files: - RotatingSymmetry180 (Kerr-rotating-180-staggered) - Test created 4 extra files: - RotatingSymmetry180 (Kerr-staggered) - Test created 4 extra files: TestLoop (testloop) Test created 8 extra files: @@ -314,9 +308,9 @@ EOS_Omni [1] - EinsteinExact_Test [6] + EinsteinExact_Test [7] - Exact [5] + Exact [4] Extract [2] @@ -350,9 +344,9 @@ RotatingDBHIVP [2] - RotatingSymmetry180 [8] + RotatingSymmetry180 [4] - RotatingSymmetry90 [4] + RotatingSymmetry90 [2] SlabTest [1] @@ -563,6 +557,8 @@ GaugeWave + KS-tilted-EE + KerrSchild Minkowski @@ -575,8 +571,6 @@ Exact: - KS-tilted - Schwarzschild_EF bowl-evolution @@ -717,30 +711,18 @@ RotatingSymmetry180: - Kerr - Kerr-EE - Kerr-rotating-180 - Kerr-rotating-180-EE - Kerr-rotating-180-staggered - Kerr-rotating-180-staggered-EE - Kerr-staggered - Kerr-staggered-EE RotatingSymmetry90: - Kerr-rotating-90 - Kerr-rotating-90-EE - Kerr-rotating-90-staggered - Kerr-rotating-90-staggered-EE SlabTest: @@ -1175,6 +1157,8 @@ Success: 1 files identical EinsteinExact_Test: GaugeWave Success: 6 files compared, 4 differ in the last digits + EinsteinExact_Test: KS-tilted-EE + Success: 235 files compared, 34 differ in the last digits EinsteinExact_Test: KerrSchild Success: 6 files compared, 4 differ in the last digits EinsteinExact_Test: Minkowski @@ -1185,8 +1169,6 @@ Success: 6 files compared, 6 differ in the last digits EinsteinExact_Test: Vaidya2 Success: 6 files compared, 5 differ in the last digits - Exact: KS-tilted - Failure: 235 files compared, 62 differ, 5 differ significantly Exact: Schwarzschild_EF Success: 112 files compared, 4 differ in the last digits Exact: bowl-evolution @@ -1266,7 +1248,7 @@ Multipole: test_rads Success: 39 files compared, 27 differ in the last digits Multipole: test_simpson - Failure: 15 files compared, 11 differ, 1 differ significantly + Success: 15 files compared, 11 differ in the last digits Multipole: test_trapezoidal Success: 15 files compared, 10 differ in the last digits Multipole: test_vars @@ -1293,28 +1275,16 @@ Failure: 96 files compared, 92 differ, 6 differ significantly RotatingDBHIVP: test_rdbh Failure: 120 files missing, 0 files compared, 0 differ - RotatingSymmetry180: Kerr - Success: 144 files compared, 13 differ in the last digits RotatingSymmetry180: Kerr-EE Success: 204 files compared, 15 differ in the last digits - RotatingSymmetry180: Kerr-rotating-180 - Success: 140 files compared, 11 differ in the last digits RotatingSymmetry180: Kerr-rotating-180-EE Success: 204 files compared, 14 differ in the last digits - RotatingSymmetry180: Kerr-rotating-180-staggered - Failure: 140 files compared, 29 differ, 1 differ significantly RotatingSymmetry180: Kerr-rotating-180-staggered-EE Success: 204 files compared, 34 differ in the last digits - RotatingSymmetry180: Kerr-staggered - Failure: 140 files compared, 27 differ, 1 differ significantly RotatingSymmetry180: Kerr-staggered-EE Success: 204 files compared, 36 differ in the last digits - RotatingSymmetry90: Kerr-rotating-90 - Success: 140 files compared, 5 differ in the last digits RotatingSymmetry90: Kerr-rotating-90-EE Success: 200 files compared, 14 differ in the last digits - RotatingSymmetry90: Kerr-rotating-90-staggered - Failure: 140 files compared, 29 differ, 1 differ significantly RotatingSymmetry90: Kerr-rotating-90-staggered-EE Success: 200 files compared, 31 differ in the last digits SlabTest: slabtest @@ -1431,20 +1401,20 @@ Summary for configuration sim - Time -> Tue Nov 12 15:15:28 CET 2013 - Host -> node6-077 + Time -> Sun Nov 17 14:11:34 CET 2013 + Host -> node3-023 Processes -> 1 User -> mundim - Total available tests -> 250 + Total available tests -> 244 Unrunnable tests -> 60 - Runnable tests -> 190 + Runnable tests -> 184 Total number of thorns -> 194 Number of tested thorns -> 56 - Number of tests passed -> 181 + Number of tests passed -> 180 Number passed only to - set tolerance -> 118 - Number failed -> 9 + set tolerance -> 117 + Number failed -> 4 Tests passed: @@ -1513,6 +1483,7 @@ E2xeon_test_dbh (from DistortedBHIVP) eos_omni (from EOS_Omni) GaugeWave (from EinsteinExact_Test) + KS-tilted-EE (from EinsteinExact_Test) KerrSchild (from EinsteinExact_Test) Minkowski (from EinsteinExact_Test) ModifiedSchwarzschildBL (from EinsteinExact_Test) @@ -1555,6 +1526,7 @@ test_44 (from Multipole) test_driscollhealy (from Multipole) test_rads (from Multipole) + test_simpson (from Multipole) test_trapezoidal (from Multipole) test_vars (from Multipole) nancount (from NaNChecker) @@ -1566,13 +1538,10 @@ qlm-ks-shifted (from QuasiLocalMeasures) qlm-ks-tilted (from QuasiLocalMeasures) qlm-minkowski (from QuasiLocalMeasures) - Kerr (from RotatingSymmetry180) Kerr-EE (from RotatingSymmetry180) - Kerr-rotating-180 (from RotatingSymmetry180) Kerr-rotating-180-EE (from RotatingSymmetry180) Kerr-rotating-180-staggered-EE (from RotatingSymmetry180) Kerr-staggered-EE (from RotatingSymmetry180) - Kerr-rotating-90 (from RotatingSymmetry90) Kerr-rotating-90-EE (from RotatingSymmetry90) Kerr-rotating-90-staggered-EE (from RotatingSymmetry90) slabtest (from SlabTest) @@ -1632,15 +1601,10 @@ Tests failed: - KS-tilted (from Exact) E2xeon_test_axioddbh (from IDAxiOddBrillBH) test_axioddbh (from IDAxiOddBrillBH) - test_simpson (from Multipole) E2xeon_test_rdbh (from RotatingDBHIVP) test_rdbh (from RotatingDBHIVP) - Kerr-rotating-180-staggered (from RotatingSymmetry180) - Kerr-staggered (from RotatingSymmetry180) - Kerr-rotating-90-staggered (from RotatingSymmetry90) ======================================================================== File [modified]: loewe__2_12.log Delta lines: +24 -61 =================================================================== --- release-info/loewe__2_12.log 2013-11-15 21:30:26 UTC (rev 1071) +++ release-info/loewe__2_12.log 2013-11-17 13:58:01 UTC (rev 1072) @@ -150,6 +150,10 @@ (EinsteinExact/EinsteinExact_Test/test/GaugeWave.par) Requires 1 processors + KS-tilted-EE in EinsteinExact_Test + (EinsteinExact/EinsteinExact_Test/test/KS-tilted-EE.par) + Requires 1 processors + KerrSchild in EinsteinExact_Test (EinsteinExact/EinsteinExact_Test/test/KerrSchild.par) Requires 1 processors @@ -211,12 +215,6 @@ Test created 4 extra files: Multipole (test_44) Test created 4 extra files: - RotatingSymmetry180 (Kerr-rotating-180) - Test created 4 extra files: - RotatingSymmetry180 (Kerr-rotating-180-staggered) - Test created 4 extra files: - RotatingSymmetry180 (Kerr-staggered) - Test created 4 extra files: ------------------------------------------------------------------------ @@ -256,7 +254,7 @@ EOS_Omni [1] - Exact [5] + Exact [4] Extract [2] @@ -300,9 +298,9 @@ RotatingDBHIVP [2] - RotatingSymmetry180 [9] + RotatingSymmetry180 [5] - RotatingSymmetry90 [5] + RotatingSymmetry90 [3] SlabTest [1] @@ -489,8 +487,6 @@ Exact: - KS-tilted - Schwarzschild_EF bowl-evolution @@ -703,35 +699,23 @@ RotatingSymmetry180: - Kerr - Kerr-EE - Kerr-rotating-180 - Kerr-rotating-180-EE - Kerr-rotating-180-staggered - Kerr-rotating-180-staggered-EE - Kerr-staggered - Kerr-staggered-EE - KerrSchild-rotating-180 + KerrSchild-rotating-180-EE RotatingSymmetry90: - Kerr-rotating-90 - Kerr-rotating-90-EE - Kerr-rotating-90-staggered - Kerr-rotating-90-staggered-EE - KerrSchild-rotating-90 + KerrSchild-rotating-90-EE SlabTest: @@ -1129,8 +1113,6 @@ Success: 100 files compared, 44 differ in the last digits EOS_Omni: eos_omni Success: 1 files identical - Exact: KS-tilted - Failure: 235 files compared, 70 differ, 5 differ significantly Exact: Schwarzschild_EF Success: 112 files compared, 4 differ in the last digits Exact: bowl-evolution @@ -1268,7 +1250,7 @@ Multipole: test_rads Success: 39 files compared, 27 differ in the last digits Multipole: test_simpson - Failure: 15 files compared, 11 differ, 1 differ significantly + Success: 15 files compared, 11 differ in the last digits Multipole: test_trapezoidal Success: 15 files compared, 10 differ in the last digits Multipole: test_vars @@ -1299,34 +1281,22 @@ Failure: 96 files compared, 92 differ, 6 differ significantly RotatingDBHIVP: test_rdbh Failure: 120 files missing, 0 files compared, 0 differ - RotatingSymmetry180: Kerr - Success: 144 files compared, 17 differ in the last digits RotatingSymmetry180: Kerr-EE Success: 204 files compared, 24 differ in the last digits - RotatingSymmetry180: Kerr-rotating-180 - Success: 140 files compared, 13 differ in the last digits RotatingSymmetry180: Kerr-rotating-180-EE Success: 204 files compared, 18 differ in the last digits - RotatingSymmetry180: Kerr-rotating-180-staggered - Failure: 140 files compared, 31 differ, 1 differ significantly RotatingSymmetry180: Kerr-rotating-180-staggered-EE Success: 204 files compared, 35 differ in the last digits - RotatingSymmetry180: Kerr-staggered - Failure: 140 files compared, 28 differ, 1 differ significantly RotatingSymmetry180: Kerr-staggered-EE Success: 204 files compared, 39 differ in the last digits - RotatingSymmetry180: KerrSchild-rotating-180 - Failure: 120 files compared, 104 differ, 13 differ significantly - RotatingSymmetry90: Kerr-rotating-90 - Success: 140 files compared, 6 differ in the last digits + RotatingSymmetry180: KerrSchild-rotating-180-EE + Success: 120 files compared, 83 differ in the last digits RotatingSymmetry90: Kerr-rotating-90-EE Success: 200 files compared, 16 differ in the last digits - RotatingSymmetry90: Kerr-rotating-90-staggered - Failure: 140 files compared, 29 differ, 1 differ significantly RotatingSymmetry90: Kerr-rotating-90-staggered-EE Success: 200 files compared, 33 differ in the last digits - RotatingSymmetry90: KerrSchild-rotating-90 - Failure: 120 files compared, 104 differ, 8 differ significantly + RotatingSymmetry90: KerrSchild-rotating-90-EE + Success: 120 files compared, 75 differ in the last digits SlabTest: slabtest Success: 1 files identical SphericalHarmonicRecon: regression_test @@ -1336,7 +1306,7 @@ TOVSolver: test_one_static_max Success: 232 files compared, 5 differ in the last digits TOVSolver: test_tov_carpet - Success: 116 files compared, 16 differ in the last digits + Success: 116 files compared, 17 differ in the last digits TOVSolver: test_two_av Success: 232 files compared, 7 differ in the last digits TOVSolver: test_two_max @@ -1427,20 +1397,20 @@ Summary for configuration sim - Time -> Tue Nov 12 15:23:30 CET 2013 - Host -> node6-081 + Time -> Sun Nov 17 14:10:30 CET 2013 + Host -> node4-067 Processes -> 2 User -> mundim - Total available tests -> 250 - Unrunnable tests -> 47 - Runnable tests -> 203 + Total available tests -> 244 + Unrunnable tests -> 48 + Runnable tests -> 196 Total number of thorns -> 194 Number of tested thorns -> 56 Number of tests passed -> 190 Number passed only to set tolerance -> 127 - Number failed -> 13 + Number failed -> 6 Tests passed: @@ -1565,6 +1535,7 @@ test_44 (from Multipole) test_driscollhealy (from Multipole) test_rads (from Multipole) + test_simpson (from Multipole) test_trapezoidal (from Multipole) test_vars (from Multipole) nancount (from NaNChecker) @@ -1578,15 +1549,14 @@ qlm-ks-shifted (from QuasiLocalMeasures) qlm-ks-tilted (from QuasiLocalMeasures) qlm-minkowski (from QuasiLocalMeasures) - Kerr (from RotatingSymmetry180) Kerr-EE (from RotatingSymmetry180) - Kerr-rotating-180 (from RotatingSymmetry180) Kerr-rotating-180-EE (from RotatingSymmetry180) Kerr-rotating-180-staggered-EE (from RotatingSymmetry180) Kerr-staggered-EE (from RotatingSymmetry180) - Kerr-rotating-90 (from RotatingSymmetry90) + KerrSchild-rotating-180-EE (from RotatingSymmetry180) Kerr-rotating-90-EE (from RotatingSymmetry90) Kerr-rotating-90-staggered-EE (from RotatingSymmetry90) + KerrSchild-rotating-90-EE (from RotatingSymmetry90) slabtest (from SlabTest) regression_test (from SphericalHarmonicRecon) test_one_boost_max (from TOVSolver) @@ -1637,19 +1607,12 @@ Tests failed: - KS-tilted (from Exact) balsara4_1d (from GRHydro) balsara4_2d (from GRHydro) E2xeon_test_axioddbh (from IDAxiOddBrillBH) test_axioddbh (from IDAxiOddBrillBH) - test_simpson (from Multipole) E2xeon_test_rdbh (from RotatingDBHIVP) test_rdbh (from RotatingDBHIVP) - Kerr-rotating-180-staggered (from RotatingSymmetry180) - Kerr-staggered (from RotatingSymmetry180) - KerrSchild-rotating-180 (from RotatingSymmetry180) - Kerr-rotating-90-staggered (from RotatingSymmetry90) - KerrSchild-rotating-90 (from RotatingSymmetry90) ========================================================================

1 0

[svn:einsteintoolkit] EinsteinExamples/trunk/par/ (Rev. 37)
by rhaas＠tapir.caltech.edu 16 Nov '13

16 Nov '13

User: rhaas Date: 2013/11/16 07:38 PM Modified: /trunk/par/ static_tov.par Log: update static_tov.par file to TOVSolver ref 137 File Changes: Directory: /trunk/par/ ====================== File [modified]: static_tov.par Delta lines: +2 -2 =================================================================== --- trunk/par/static_tov.par 2013-05-27 22:30:09 UTC (rev 36) +++ trunk/par/static_tov.par 2013-11-17 01:38:06 UTC (rev 37) @@ -181,8 +181,8 @@ ADMBase::initial_dtshift = "zero" TOVSolver::TOV_Rho_Central[0] = 1.28e-3 -TOVSolver::TOV_Gamma[0] = 2.0 -TOVSolver::TOV_K[0] = 100.0 +TOVSolver::TOV_Gamma = 2.0 +TOVSolver::TOV_K = 100.0 IOBasic::outInfo_every = 1

1 0

[svn:einsteintoolkit] NullSHRExtract/trunk/src/ (Rev. 34)
by reisswig＠tapir.caltech.edu 15 Nov '13

15 Nov '13

User: reisswig Date: 2013/11/15 07:10 PM Modified: /trunk/src/ NullSHRE_modEta.F90 Log: Fixed a bug in modEta.F90. There was a missing time derivative of g in time derivative of $L = \alpha - g_{ij} \beta^{i} s^{j}$. File Changes: Directory: /trunk/src/ ====================== File [modified]: NullSHRE_modEta.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/NullSHRE_modEta.F90 2013-11-15 22:46:10 UTC (rev 33) +++ trunk/src/NullSHRE_modEta.F90 2013-11-16 01:10:11 UTC (rev 34) @@ -92,7 +92,7 @@ delld(4)%d = delld(4)%d & - g(i,j)%d * dbeta(i,4)%d * sa(j)%d& - g(i,j)%d * beta(i)%d * dsa(j,4)%d& - - beta(i)%d * sa(j)%d !*dg(i,j,4)%d + - beta(i)%d * sa(j)%d *dg(i,j,4)%d end do end do

1 0

[svn:einsteintoolkit] NullSHRExtract/branches/tapir/src/ (Rev. 33)
by bela＠caltech.edu 15 Nov '13

15 Nov '13

User: szilagyi Date: 2013/11/15 04:46 PM Modified: /branches/tapir/src/ NullSHRE_WTSphHarm.F90 Log: add output for time-deriv of lapse,shift File Changes: Directory: /branches/tapir/src/ =============================== File [modified]: NullSHRE_WTSphHarm.F90 Delta lines: +66 -46 =================================================================== --- branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-11-15 20:09:38 UTC (rev 32) +++ branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-11-15 22:46:10 UTC (rev 33) @@ -85,12 +85,13 @@ type (gf2d), dimension(2) :: angular_tensor_1_n, angular_tensor_1_s logical, parameter :: output_gij = .false. - logical, parameter :: output_lapse = .false. - logical, parameter :: output_shift = .false. + logical, parameter :: output_lapse = .true. + logical, parameter :: output_shift = .true. logical, parameter :: output_dgij = .false. - logical, parameter :: output_ell = .true. + logical, parameter :: output_ell = .false. logical, parameter :: output_w0parts = .false. logical, parameter :: output_u1parts = .false. + logical, parameter :: output_bondi = .false. CCTK_INT :: a, b angular_tensor_1_n(1)%d => tmp1(:,:,1) @@ -140,6 +141,16 @@ call NullDecomp_WriteCoefsForRealFunc('shiftU_z_wt', cctkGH, truncate, & null_lsh, zeta, beta_n(3)%d, beta_s(3)%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('dshiftU_x_dt_wt', cctkGH, truncate, & + null_lsh, zeta, dbeta_n(1,4)%d, dbeta_s(1,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dshiftU_y_dt_wt', cctkGH, truncate, & + null_lsh, zeta, dbeta_n(2,4)%d, dbeta_s(2,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dshiftU_z_dt_wt', cctkGH, truncate, & + null_lsh, zeta, dbeta_n(3,4)%d, dbeta_s(3,4)%d, cctk_time) + + end if if(output_lapse) then @@ -149,6 +160,9 @@ call NullDecomp_WriteCoefsForRealFunc('alpha_wt', cctkGH, truncate, & null_lsh, zeta, alpha_n%d, alpha_s%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('dalpha_dt_wt', cctkGH, truncate, & + null_lsh, zeta, dalpha_n(4)%d, dalpha_s(4)%d, cctk_time) + end if if(output_dgij) then @@ -490,68 +504,74 @@ end if - ! xb, rb, r_l, r_t - call NullDecomp_WriteCoefsForRealFunc('X_wt', cctkGH, truncate, & - null_lsh, zeta, x_wt_n%d, x_wt_s%d, cctk_time) + if(output_bondi) then - call NullDecomp_WriteCoefsForRealFunc('r0_wt', cctkGH, truncate, & - null_lsh, zeta, r0_n%d, r0_s%d, cctk_time) + ! xb, rb, r_l, r_t + call NullDecomp_WriteCoefsForRealFunc('X_wt', cctkGH, truncate, & + null_lsh, zeta, x_wt_n%d, x_wt_s%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('r_l_wt', cctkGH, truncate, & - null_lsh, zeta, dr0_n(1)%d, dr0_s(1)%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('r0_wt', cctkGH, truncate, & + null_lsh, zeta, r0_n%d, r0_s%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('r_t_wt', cctkGH, truncate, & - null_lsh, zeta, dr0_n(4)%d, dr0_s(4)%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('r_l_wt', cctkGH, truncate, & + null_lsh, zeta, dr0_n(1)%d, dr0_s(1)%d, cctk_time) - call NullSHRE_ConvertAngOneForm(ip_n, dr0_n(2:3), angular_tensor_1_n) - call NullSHRE_ConvertAngOneForm(ip_s, dr0_s(2:3), angular_tensor_1_s) + call NullDecomp_WriteCoefsForRealFunc('r_t_wt', cctkGH, truncate, & + null_lsh, zeta, dr0_n(4)%d, dr0_s(4)%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('r_theta_wt', cctkGH, truncate, & - null_lsh, zeta, angular_tensor_1_n(1)%d, angular_tensor_1_s(1)%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('r_phi_wt', cctkGH, truncate, & - null_lsh, zeta, angular_tensor_1_n(2)%d, angular_tensor_1_s(2)%d, cctk_time) + call NullSHRE_ConvertAngOneForm(ip_n, dr0_n(2:3), angular_tensor_1_n) + call NullSHRE_ConvertAngOneForm(ip_s, dr0_s(2:3), angular_tensor_1_s) - call NullDecomp_WriteCoefsForRealFunc('AffineInvMetric_lambda_lambda', cctkGH, truncate, & - null_lsh, zeta, etaup0_n(1,1)%d, etaup0_s(1,1)%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('r_theta_wt', cctkGH, truncate, & + null_lsh, zeta, angular_tensor_1_n(1)%d, angular_tensor_1_s(1)%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('r_phi_wt', cctkGH, truncate, & + null_lsh, zeta, angular_tensor_1_n(2)%d, angular_tensor_1_s(2)%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('AffineInvMetric_lambda_u', cctkGH, truncate, & - null_lsh, zeta, etaup0_n(1,4)%d, etaup0_s(1,4)%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc(& + 'AffineInvMetric_lambda_lambda', cctkGH, truncate, & + null_lsh, zeta, etaup0_n(1,1)%d, etaup0_s(1,1)%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) + call NullDecomp_WriteCoefsForRealFunc(& + 'AffineInvMetric_lambda_u', cctkGH, truncate, & + null_lsh, zeta, etaup0_n(1,4)%d, etaup0_s(1,4)%d, cctk_time) - ! J, beta, U, Q, W + 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) - call NullDecomp_WriteCoefsForComplexFunc('J_wt', cctkGH, truncate, & - null_lsh, 2, zeta, j_wt_n%d, j_wt_s%d, cctk_time) + ! J, beta, U, Q, W - call NullDecomp_WriteCoefsForComplexFunc('Jl_wt', cctkGH, truncate, & - null_lsh, 2, zeta, j_l_n%d, j_l_s%d, cctk_time) + call NullDecomp_WriteCoefsForComplexFunc('J_wt', cctkGH, truncate, & + null_lsh, 2, zeta, j_wt_n%d, j_wt_s%d, cctk_time) - call NullDecomp_WriteCoefsForComplexFunc('Q_wt', cctkGH, truncate, & - null_lsh, 1, zeta, q_wt_n%d, q_wt_s%d, cctk_time) + call NullDecomp_WriteCoefsForComplexFunc('Jl_wt', cctkGH, truncate, & + null_lsh, 2, zeta, j_l_n%d, j_l_s%d, cctk_time) - call NullDecomp_WriteCoefsForComplexFunc('U_wt', cctkGH, truncate, & - null_lsh, 1, zeta, u_wt_n%d, u_wt_s%d, cctk_time) + call NullDecomp_WriteCoefsForComplexFunc('Q_wt', cctkGH, truncate, & + null_lsh, 1, zeta, q_wt_n%d, q_wt_s%d, cctk_time) - call NullDecomp_WriteCoefsForComplexFunc('Ul_wt', cctkGH, truncate, & - null_lsh, 1, zeta, u_l_n%d, u_l_s%d, cctk_time) + call NullDecomp_WriteCoefsForComplexFunc('U_wt', cctkGH, truncate, & + null_lsh, 1, zeta, u_wt_n%d, u_wt_s%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('B_wt', cctkGH, truncate, & - null_lsh, zeta, beta_wt_n%d, beta_wt_s%d, cctk_time) + call NullDecomp_WriteCoefsForComplexFunc('Ul_wt', cctkGH, truncate, & + null_lsh, 1, zeta, u_l_n%d, u_l_s%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('Bl_wt', cctkGH, truncate, & - null_lsh, zeta, beta_l_n%d, beta_l_s%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('B_wt', cctkGH, truncate, & + null_lsh, zeta, beta_wt_n%d, beta_wt_s%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('W_wt', cctkGH, truncate, & - null_lsh, zeta, w_wt_n%d, w_wt_s%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('Bl_wt', cctkGH, truncate, & + null_lsh, zeta, beta_l_n%d, beta_l_s%d, cctk_time) - call NullDecomp_WriteCoefsForRealFunc('Wl_wt', cctkGH, truncate, & - null_lsh, zeta, w_l_n%d, w_l_s%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('W_wt', cctkGH, truncate, & + null_lsh, zeta, w_wt_n%d, w_wt_s%d, cctk_time) + call NullDecomp_WriteCoefsForRealFunc('Wl_wt', cctkGH, truncate, & + null_lsh, zeta, w_l_n%d, w_l_s%d, cctk_time) + + end if + ! bits of W if(output_w0parts) then

1 0

← Newer
1
...
5
6
7
8
9
10
11
...
14
Older →

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits November 2013