Commits November 2010

commits@lists.einsteintoolkit.org

9 participants
258 discussions

[svn:einsteintoolkit] manifest/trunk/ (Rev. 45)
by bcmsma＠astro.rit.edu 15 Nov '10

15 Nov '10

User: bmundim Date: 2010/11/15 11:49 AM Modified: /trunk/ einsteintoolkit.bib Log: Noexcision belongs to EinsteinInitialData instead. File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.bib Delta lines: +1 -1 =================================================================== --- trunk/einsteintoolkit.bib 2010-11-15 17:25:56 UTC (rev 44) +++ trunk/einsteintoolkit.bib 2010-11-15 17:49:27 UTC (rev 45) @@ -59,7 +59,7 @@ } @Article{Brown:2007pg, - required-for = "EinsteinAnalysis/NoExcision", + required-for = "EinsteinInitialData/NoExcision", author = "David Brown and Olivier Sarbach and Erik Schnetter and Manuel Tiglio and Peter Diener and Ian Hawke and Denis Pollney",

1 0

[svn:einsteintoolkit] tools/trunk/bibtex/ (Rev. 18)
by knarf＠cct.lsu.edu 15 Nov '10

15 Nov '10

User: knarf Date: 2010/11/15 11:41 AM Added: /trunk/bibtex/ bib2xhtml, parse_bibtex.pl Log: bibtex utils File Changes: Directory: /trunk/bibtex/ ========================= File [added]: bib2xhtml Delta lines: <Unavailable> File [added]: parse_bibtex.pl Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] manifest/trunk/ (Rev. 44)
by bcmsma＠astro.rit.edu 15 Nov '10

15 Nov '10

User: bmundim Date: 2010/11/15 11:25 AM Modified: /trunk/ einsteintoolkit.bib Log: TwoPunctures is located in the EinsteinInitialData arrangement. File Changes: Directory: /trunk/ ================== File [modified]: einsteintoolkit.bib Delta lines: +1 -1 =================================================================== --- trunk/einsteintoolkit.bib 2010-11-08 19:27:40 UTC (rev 43) +++ trunk/einsteintoolkit.bib 2010-11-15 17:25:56 UTC (rev 44) @@ -25,7 +25,7 @@ } @Article{Ansorg:2004ds, - required-for = "EinsteinAnalysis/TwoPunctures", + required-for = "EinsteinInitialData/TwoPunctures", author = "Ansorg, Marcus and Br{\"u}gmann, Bernd and Tichy, Wolfgang", title = "{A single-domain spectral method for black hole

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 403)
by knarf＠cct.lsu.edu 14 Nov '10

14 Nov '10

User: knarf Date: 2010/11/14 10:22 PM Modified: /release-info/ philip__1.log, philip__2.log Log: update File Changes: Directory: /release-info/ ========================= File [modified]: philip__1.log Delta lines: +6 -10 =================================================================== --- release-info/philip__1.log 2010-11-15 04:21:29 UTC (rev 402) +++ release-info/philip__1.log 2010-11-15 04:22:41 UTC (rev 403) @@ -1003,7 +1003,7 @@ WaveToyFreeF90: test_rad Success: 6 files compared, 1 differ in the last digits WeylScal4: teukolsky - Failure: 8 files compared, 8 differ, 6 differ significantly + Success: 8 files compared, 6 differ in the last digits WeylScal4: teukolskyID Success: 8 files compared, 6 differ in the last digits @@ -1012,8 +1012,8 @@ Summary for configuration sim - Time -> Thu Nov 4 23:16:28 CDT 2010 - Host -> philip033 + Time -> Sun Nov 14 18:09:54 CST 2010 + Host -> philip034 User -> knarf Total available tests -> 172 @@ -1021,14 +1021,10 @@ Runnable tests -> 130 Total number of thorns -> 155 Number of tested thorns -> 38 - Number of tests passed -> 129 + Number of tests passed -> 130 Number passed only to - set tolerance -> 80 - Number failed -> 1 + set tolerance -> 81 + Number failed -> 0 - Tests failed: - - teukolsky (from WeylScal4) - ======================================================================== File [modified]: philip__2.log Delta lines: +2 -2 =================================================================== --- release-info/philip__2.log 2010-11-15 04:21:29 UTC (rev 402) +++ release-info/philip__2.log 2010-11-15 04:22:41 UTC (rev 403) @@ -1003,8 +1003,8 @@ Summary for configuration sim - Time -> Thu Nov 4 23:30:40 CDT 2010 - Host -> philip033 + Time -> Sun Nov 14 18:24:20 CST 2010 + Host -> philip034 User -> knarf Total available tests -> 172

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 402)
by knarf＠cct.lsu.edu 14 Nov '10

14 Nov '10

User: knarf Date: 2010/11/14 10:21 PM Modified: /release-info/ topf__1.log, topf__2.log Log: update File Changes: Directory: /release-info/ ========================= File [modified]: topf__1.log Delta lines: +5 -6 =================================================================== --- release-info/topf__1.log 2010-11-12 21:25:14 UTC (rev 401) +++ release-info/topf__1.log 2010-11-15 04:21:29 UTC (rev 402) @@ -814,7 +814,7 @@ CarpetReduce: staggered Success: 5 files identical DistortedBHIVP: E2xeon_test_dbh - Failure: 100 files compared, 44 differ, 44 differ significantly + Success: 100 files compared, 44 differ in the last digits Exact: Schwarzschild_EF Success: 112 files compared, 3 differ in the last digits Exact: bowl-evolution @@ -1025,7 +1025,7 @@ Summary for configuration sim - Time -> Wed Nov 3 22:23:36 CDT 2010 + Time -> Sun Nov 14 21:49:44 CST 2010 Host -> topf.wg User -> knarf @@ -1034,14 +1034,13 @@ Runnable tests -> 132 Total number of thorns -> 157 Number of tested thorns -> 39 - Number of tests passed -> 128 + Number of tests passed -> 129 Number passed only to - set tolerance -> 76 - Number failed -> 4 + set tolerance -> 77 + Number failed -> 3 Tests failed: - E2xeon_test_dbh (from DistortedBHIVP) halfsphere_pugh_ppm (from Hydro_InitExcision) x_flip_pugh_ppm (from Hydro_InitExcision) x_pugh_ppm (from Hydro_InitExcision) File [modified]: topf__2.log Delta lines: +6 -8 =================================================================== --- release-info/topf__2.log 2010-11-12 21:25:14 UTC (rev 401) +++ release-info/topf__2.log 2010-11-15 04:21:29 UTC (rev 402) @@ -761,7 +761,7 @@ ADM: test_ADM_1 Success: 32 files compared, 1 differ in the last digits ADM: test_ADM_2 - Failure: 32 files compared, 32 differ, 32 differ significantly + Success: 32 files identical ADM: test_ADM_3 Success: 32 files compared, 1 differ in the last digits ADM: test_ADM_4 @@ -809,7 +809,7 @@ AHFinderDirect: recoverML Failure: 5 files compared, 5 differ, 5 differ significantly DistortedBHIVP: E2xeon_test_dbh - Failure: 100 files compared, 44 differ, 44 differ significantly + Success: 100 files compared, 44 differ in the last digits Exact: Schwarzschild_EF Success: 112 files compared, 3 differ in the last digits Exact: bowl-evolution @@ -1016,7 +1016,7 @@ Summary for configuration sim - Time -> Wed Nov 3 22:54:24 CDT 2010 + Time -> Sun Nov 14 22:20:03 CST 2010 Host -> topf.wg User -> knarf @@ -1025,14 +1025,13 @@ Runnable tests -> 126 Total number of thorns -> 157 Number of tested thorns -> 36 - Number of tests passed -> 108 + Number of tests passed -> 110 Number passed only to - set tolerance -> 62 - Number failed -> 18 + set tolerance -> 63 + Number failed -> 16 Tests failed: - test_ADM_2 (from ADM) Kerr (from AHFinderDirect) Kerr-definition-expansion (from AHFinderDirect) Kerr-definition-expansion-product (from AHFinderDirect) @@ -1046,7 +1045,6 @@ checkpointML (from AHFinderDirect) misner1.2-025 (from AHFinderDirect) recoverML (from AHFinderDirect) - E2xeon_test_dbh (from DistortedBHIVP) halfsphere_pugh_ppm (from Hydro_InitExcision) x_flip_pugh_ppm (from Hydro_InitExcision) x_pugh_ppm (from Hydro_InitExcision)

1 0

[svn:einsteintoolkit] incoming/EOS_Omni/ (Rev. 19)
by cott＠tapir.caltech.edu 14 Nov '10

14 Nov '10

User: cott Date: 2010/11/14 12:17 AM Modified: /EOS_Omni/ param.ccl /EOS_Omni/src/nuc_eos/ readtable.F90 Log: * fix a stupid issue with hdf5 File Changes: Directory: /EOS_Omni/ ===================== File [modified]: param.ccl Delta lines: +1 -1 =================================================================== --- EOS_Omni/param.ccl 2010-11-14 05:04:22 UTC (rev 18) +++ EOS_Omni/param.ccl 2010-11-14 06:17:54 UTC (rev 19) @@ -58,7 +58,7 @@ ################ Finite-Temperature Nuclear EOS -BOOLEAN nuceos_read_table "Read in EOS table?" +BOOLEAN nuceos_read_table "Read in EOS table?" STEERABLE=ALWAYS { } "No" Directory: /EOS_Omni/src/nuc_eos/ ================================= File [modified]: readtable.F90 Delta lines: +2 -1 =================================================================== --- EOS_Omni/src/nuc_eos/readtable.F90 2010-11-14 05:04:22 UTC (rev 18) +++ EOS_Omni/src/nuc_eos/readtable.F90 2010-11-14 06:17:54 UTC (rev 19) @@ -213,7 +213,8 @@ call h5fclose_f (file_id,error) - call h5close_f (error) +! not a good idea to call this function as it will screw recovery +! call h5close_f (error) ! set min-max values:

1 0

[svn:einsteintoolkit] GRHydro/branches/hot_and_MHD_temp_dev/src/ (Rev. 175)
by knarf＠cct.lsu.edu 13 Nov '10

13 Nov '10

User: knarf Date: 2010/11/13 11:28 PM Modified: /branches/hot_and_MHD_temp_dev/src/ GRHydro_Marquina.F90 Log: fix hot part of eigenproblem_marquina File Changes: Directory: /branches/hot_and_MHD_temp_dev/src/ ============================================== File [modified]: GRHydro_Marquina.F90 Delta lines: +2 -2 =================================================================== --- branches/hot_and_MHD_temp_dev/src/GRHydro_Marquina.F90 2010-11-14 05:15:59 UTC (rev 174) +++ branches/hot_and_MHD_temp_dev/src/GRHydro_Marquina.F90 2010-11-14 05:28:48 UTC (rev 175) @@ -256,7 +256,7 @@ marquinaflux(5)) else xtemp = temperature(i,j,k) - call eigenproblem_marquina(GRHydro_eos_handle,& + call eigenproblem_marquina_hot(GRHydro_eos_handle,& primm_i(1),primm_i(3), & primm_i(4),primm_i(2),primm_i(5),primp(1), & primp(3),primp(4),primp(2),primp(5), & @@ -285,7 +285,7 @@ marquinaflux(5)) else xtemp = temperature(i,j,k) - call eigenproblem_marquina(GRHydro_eos_handle,& + call eigenproblem_marquina_hot(GRHydro_eos_handle,& primm_i(1),primm_i(4), & primm_i(2),primm_i(3),primm_i(5),primp(1), & primp(4),primp(2),primp(3),primp(5), &

1 0

[svn:einsteintoolkit] GRHydro/branches/hot_and_MHD_temp_dev/src/ (Rev. 174)
by knarf＠cct.lsu.edu 13 Nov '10

13 Nov '10

User: knarf Date: 2010/11/13 11:15 PM Modified: /branches/hot_and_MHD_temp_dev/src/ GRHydro_Con2Prim.F90 Log: fix type of GRHydro_C2P_failed: this is a CCTK_REAL File Changes: Directory: /branches/hot_and_MHD_temp_dev/src/ ============================================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: +4 -2 =================================================================== --- branches/hot_and_MHD_temp_dev/src/GRHydro_Con2Prim.F90 2010-11-11 23:05:07 UTC (rev 173) +++ branches/hot_and_MHD_temp_dev/src/GRHydro_Con2Prim.F90 2010-11-14 05:15:59 UTC (rev 174) @@ -264,7 +264,8 @@ y, z, r, GRHydro_rho_min CCTK_REAL s2, c0, c1, c2, c3, c4, f, df, ftol, v2, w, vlowx, vlowy, vlowz - CCTK_INT count, i, handle, GRHydro_reflevel, GRHydro_C2P_failed + CCTK_INT count, i, handle, GRHydro_reflevel + CCTK_REAL GRHydro_C2P_failed CCTK_REAL udens, usx, usy, usz, utau, pold, pnew, epsold, epsnew, w2, & w2mhalf, temp1, drhobydpress, depsbydpress, dpressbydeps, dpressbydrho, pmin, epsmin character(len=200) warnline @@ -486,7 +487,8 @@ y, z, r, GRHydro_rho_min CCTK_REAL temp, ye CCTK_REAL s2, c0, c1, c2, c3, c4, f, df, ftol, v2, w, vlowx, vlowy, vlowz - CCTK_INT ii,jj,kk,count, i, handle, GRHydro_reflevel, GRHydro_C2P_failed + CCTK_INT ii,jj,kk,count, i, handle, GRHydro_reflevel + CCTK_REAL GRHydro_C2P_failed CCTK_REAL udens, usx, usy, usz, utau, pold, pnew, epsold, epsnew, w2, & w2mhalf, temp1, drhobydpress, depsbydpress, dpressbydeps, dpressbydrho, pmin, epsmin character(len=200) warnline

1 0

[svn:einsteintoolkit] incoming/EOS_Omni/src/nuc_eos/ (Rev. 18)
by knarf＠cct.lsu.edu 13 Nov '10

13 Nov '10

User: knarf Date: 2010/11/13 11:04 PM Modified: /EOS_Omni/src/nuc_eos/ bisection.F90 Log: correct compile error File Changes: Directory: /EOS_Omni/src/nuc_eos/ ================================= File [modified]: bisection.F90 Delta lines: +1 -1 =================================================================== --- EOS_Omni/src/nuc_eos/bisection.F90 2010-11-13 04:13:28 UTC (rev 17) +++ EOS_Omni/src/nuc_eos/bisection.F90 2010-11-14 05:04:22 UTC (rev 18) @@ -19,7 +19,7 @@ integer bcount,i,itmax,maxbcount - real*8 bivar + real*8 bivar(*) bcount = 0 maxbcount = 80

1 0

[svn:einsteintoolkit] incoming/EOS_Omni/src/ (Rev. 17)
by cott＠tapir.caltech.edu 12 Nov '10

12 Nov '10

User: cott Date: 2010/11/12 10:13 PM Modified: /EOS_Omni/src/ EOS_Omni_SingleVarCalls.F90 /EOS_Omni/src/nuc_eos/ nuc_eos.F90 Log: * small optimization: introduce call to nuc_eos that just returns the pressure (and eps) File Changes: Directory: /EOS_Omni/src/nuc_eos/ ================================= File [modified]: nuc_eos.F90 Delta lines: +101 -0 =================================================================== --- EOS_Omni/src/nuc_eos/nuc_eos.F90 2010-11-12 00:35:41 UTC (rev 16) +++ EOS_Omni/src/nuc_eos/nuc_eos.F90 2010-11-13 04:13:28 UTC (rev 17) @@ -264,6 +264,87 @@ end subroutine nuc_eos_short + +subroutine nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& + keytemp,keyerr,rfeps) + + use eosmodule + implicit none + + real*8, intent(in) :: xrho,xye + real*8, intent(inout) :: xtemp,xenr + real*8, intent(in) :: rfeps + real*8, intent(out) :: xprs + integer, intent(in) :: keytemp + integer, intent(out) :: keyerr + + ! local variables + real*8 :: xcs2,xdpderho,xdpdrhoe + real*8 :: lr,lt,y,xx,xeps,leps,xs + real*8 :: d1,d2,d3,ff(8) + integer :: keyerrt = 0 + + if(xrho.gt.eos_rhomax) then + stop "nuc_eos: rho > rhomax" + endif + + if(xrho.lt.eos_rhomin*1.2d0) then + call nuc_low_eos(xrho,xenr,xprs,xcs2,xdpderho,xdpdrhoe,keytemp) + return + endif + + if(xye.gt.eos_yemax) then + stop "nuc_eos: ye > yemax" + endif + + if(xye.lt.eos_yemin) then + stop "nuc_eos: ye < yemin" + endif + + if(keytemp.eq.1) then + if(xtemp.gt.eos_tempmax) then + stop "nuc_eos: temp > tempmax" + endif + + if(xtemp.lt.eos_tempmin) then + call nuc_low_eos(xrho,xenr,xprs,xcs2,xdpderho,xdpdrhoe,keytemp) + return + endif + endif + + lr = log10(xrho) + lt = log10(xtemp) + y = xye + xeps = xenr + energy_shift + leps = log10(max(xeps,1.0d0)) + + keyerr = 0 + + if(keytemp.eq.0) then + !need to find temperature based on xeps + call findtemp(lr,lt,y,leps,keyerrt,rfeps) + if(keyerrt.ne.0) then + keyerr = keyerrt + return + endif + xtemp = 10.0d0**lt + + elseif(keytemp.eq.2) then + stop "eos_nuc_press does not support keytemp.eq.2" + endif + + ! have temperature, proceed: + call findall_press_eps(lr,lt,y,ff) + xprs = 10.0d0**ff(1) + + xeps = 10.0d0**ff(2) - energy_shift + if((keytemp.eq.1).or.(keytemp.eq.2)) then + xenr = xeps + endif + +end subroutine nuc_eos_press_eps + + subroutine findthis(lr,lt,y,value,array,d1,d2,d3) use eosmodule @@ -321,3 +402,23 @@ ff(:) = ffx(:,1) end subroutine findall_short + + +subroutine findall_press_eps(lr,lt,y,ff) + + use eosmodule + implicit none + + real*8 ffx(8,1) + real*8 ff(8) + real*8 lr,lt,y + integer i + integer :: nvarsx = 2 + + +! Ewald's interpolator + call intp3d_many(lr,lt,y,ffx,1,alltables(:,:,:,1:2), & + nrho,ntemp,nye,nvarsx,logrho,logtemp,ye) + ff(:) = ffx(:,1) + +end subroutine findall_press_eps Directory: /EOS_Omni/src/ ========================= File [modified]: EOS_Omni_SingleVarCalls.F90 Delta lines: +1 -2 =================================================================== --- EOS_Omni/src/EOS_Omni_SingleVarCalls.F90 2010-11-12 00:35:41 UTC (rev 16) +++ EOS_Omni/src/EOS_Omni_SingleVarCalls.F90 2010-11-13 04:13:28 UTC (rev 17) @@ -96,8 +96,7 @@ xtemp = temp(i) xye = ye(i) xenr = eps(i) * inv_eps_gf - call nuc_eos_short(xrho,xtemp,xye,xenr,xprs,& - xent,xcs2,xdedt,xdpderho,xdpdrhoe,xmunu,& + call nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& keytemp,keyerr(i),rf_precision) if(keyerr(i).ne.0) then

1 0

← Newer
1
...
16
17
18
19
20
21
22
...
26
Older →

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits November 2010