Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 542)
by rhaas＠tapir.caltech.edu 06 Jul '13

06 Jul '13

User: rhaas Date: 2013/07/06 01:10 PM Modified: /trunk/src/ GRHydro_Con2PrimHot.F90 Log: GRHydro: Critical bugfix: Con2Prim was probably never executed in post_recover_variables or for new points in post_regrid since excision mask does typically not get initialized before MoL_PostStep! Also: Instead of aborting Con2Prim, set points to atmopshere for a point where hydro_excision_mask > 0. From: Christian Reisswig <reisswig(a)scriwalker.(none)> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2PrimHot.F90 Delta lines: +7 -4 =================================================================== --- trunk/src/GRHydro_Con2PrimHot.F90 2013-07-06 18:10:12 UTC (rev 541) +++ trunk/src/GRHydro_Con2PrimHot.F90 2013-07-06 18:10:16 UTC (rev 542) @@ -35,7 +35,7 @@ CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup ! begin EOS Omni vars - integer :: n,keytemp,anyerr,keyerr(1) + integer :: n,keytemp,anyerr,keyerr(1), reset_to_atmo real*8 :: xpress(1),xeps(1),xtemp(1),xye(1),xrho(1) n = 1;keytemp = 0;anyerr = 0;keyerr(1) = 0 xpress = 0.0d0;xeps = 0.0d0;xtemp = 0.0d0;xye = 0.0d0 @@ -89,9 +89,8 @@ do i = 1, nx - !do not compute if in atmosphere or in excised region - if ((atmosphere_mask(i,j,k) .ne. 0) .or. & - (hydro_excision_mask(i,j,k) .ne. 0)) cycle + !do not compute if in atmosphere + if (atmosphere_mask(i,j,k) .gt. 0) cycle epsnegative = .false. @@ -121,8 +120,12 @@ GRHydro_Y_e_min) endif + reset_to_atmo = 0 IF_BELOW_ATMO(dens(i,j,k), sqrt(det)*GRHydro_rho_min, GRHydro_atmo_tolerance, r(i,j,k)) then + reset_to_atmo = 1 + endif + if (reset_to_atmo .gt. 0 .or. hydro_excision_mask(i,j,k) .gt. 0) then SET_ATMO_MIN(dens(i,j,k), sqrt(det)*GRHydro_rho_min, r(i,j,k)) SET_ATMO_MIN(rho(i,j,k), GRHydro_rho_min, r(i,j,k)) scon(i,j,k,:) = 0.d0

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 541)
by rhaas＠tapir.caltech.edu 06 Jul '13

06 Jul '13

User: rhaas Date: 2013/07/06 01:10 PM Modified: /trunk/src/ GRHydro_Prim2Con.F90 Log: GRHydro: limit temperature to valid range in prim2con add some more error checking From: Christian David Ott <cott(a)zwicky-b.(none)> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Prim2Con.F90 Delta lines: +45 -3 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2013-07-06 18:10:08 UTC (rev 540) +++ trunk/src/GRHydro_Prim2Con.F90 2013-07-06 18:10:12 UTC (rev 541) @@ -322,14 +322,56 @@ *dvelz*dvelz + 2*gxy*dvelx*dvely + 2*gxz*dvelx *dvelz + 2*gyz& *dvely*dvelz)) + ! make sure Y_e and temperature within allowed range ye = max(min(ye,GRHydro_Y_e_max),GRHydro_Y_e_min) + temp = max(GRHydro_hot_atmo_temp,min(temp,GRHydro_max_temp)) call EOS_Omni_press(handle,lkeytemp,GRHydro_eos_rf_prec,n,& drho,deps,temp,ye,dpress,keyerr,anyerr) ! error handling if(anyerr.ne.0) then if(reconstruct_temper.ne.0) then - + if(keytemp.eq.1) then + !$OMP CRITICAL + call CCTK_WARN(1,"EOS error in prim2con_hot:") + write(warnline,"(i8,4i5,1P10E15.6)") cctk_iteration,GRHydro_Reflevel,ii,jj,kk,x,y,z,r + call CCTK_WARN(1,warnline) + write(warnline,"(1P10E15.6)") drho,deps,temp,ye + call CCTK_WARN(1,warnline) + write(warnline,"(A7,i8,A10,i8)") "keyerr: ",keyerr(1)," keytemp:",lkeytemp + call CCTK_WARN(0,warnline) + !$OMP END CRITICAL + else + if(keyerr(1).eq.668 .and. temp(1).lt.10.0d0) then + !$OMP CRITICAL + write(warnline,"(A18,i7,4i4,1P10E15.6)") "p2c resetting T: ",cctk_iteration, & + GRHydro_reflevel,ii,jj,kk,drho,deps,temp,ye + call CCTK_WARN(1,warnline) + !$OMP END CRITICAL + temp = GRHydro_hot_atmo_temp + lkeytemp = 1 + call EOS_Omni_press(handle,lkeytemp,GRHydro_eos_rf_prec,n,& + drho,deps,temp,ye,dpress,keyerr,anyerr) + lkeytemp = 0 + if(anyerr.ne.0) then + !$OMP CRITICAL + write(warnline,"(A7,i8,A10,i8)") "keyerr: ",keyerr(1)," keytemp:",lkeytemp + call CCTK_WARN(1,warnline) + call CCTK_WARN(0,"Fatal EOS error in p2c!") + !$OMP END CRITICAL + endif + else + !$OMP CRITICAL + call CCTK_WARN(1,"EOS error in prim2con_hot:") + write(warnline,"(i8,4i5,1P10E15.6)") cctk_iteration,GRHydro_Reflevel,ii,jj,kk,x,y,z,r + call CCTK_WARN(1,warnline) + write(warnline,"(1P10E15.6)") drho,deps,temp,ye + call CCTK_WARN(1,warnline) + write(warnline,"(A7,i8,A10,i8)") "keyerr: ",keyerr(1)," keytemp:",lkeytemp + call CCTK_WARN(0,warnline) + !$OMP END CRITICAL + endif + endif else if(GRHydro_reflevel.lt.GRHydro_c2p_warn_from_reflevel) then ! in this case (coarse grid error that is hopefully restricted @@ -341,7 +383,7 @@ lkeytemp=0 if(anyerr.ne.0) then !$OMP CRITICAL - call CCTK_WARN(1,"EOS error in prim2con_hot: lev 2") + call CCTK_WARN(1,"EOS error in prim2con_hot:") write(warnline,"(i8,3i5,1P10E15.6)") cctk_iteration,ii,jj,kk,x,y,z,r call CCTK_WARN(1,warnline) write(warnline,"(1P10E15.6)") drho,deps,temp,ye @@ -391,7 +433,7 @@ endif endif endif - + vlowx = gxx*dvelx + gxy*dvely + gxz*dvelz vlowy = gxy*dvelx + gyy*dvely + gyz*dvelz vlowz = gxz*dvelx + gyz*dvely + gzz*dvelz

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 540)
by rhaas＠tapir.caltech.edu 06 Jul '13

06 Jul '13

User: rhaas Date: 2013/07/06 01:10 PM Modified: /trunk/src/ Con2Prim_fortran_interfaces.F90 Log: GRHydro: change C2P interface From: Christian David Ott <cott(a)zwicky-b.(none)> File Changes: Directory: /trunk/src/ ====================== File [modified]: Con2Prim_fortran_interfaces.F90 Delta lines: +3 -1 =================================================================== --- trunk/src/Con2Prim_fortran_interfaces.F90 2013-07-06 18:10:06 UTC (rev 539) +++ trunk/src/Con2Prim_fortran_interfaces.F90 2013-07-06 18:10:08 UTC (rev 540) @@ -6,7 +6,8 @@ interface subroutine Con2Prim_pt(& - handle, & + cctk_iteration,ii,jj,kk,& + handle, & dens, & sx, sy, sz, & tau, & @@ -23,6 +24,7 @@ GRHydro_reflevel, GRHydro_C2P_failed) implicit none + CCTK_INT cctk_iteration, ii, jj, kk CCTK_INT handle CCTK_REAL dens CCTK_REAL sx, sy, sz

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 539)
by rhaas＠tapir.caltech.edu 06 Jul '13

06 Jul '13

User: rhaas Date: 2013/07/06 01:10 PM Added: /trunk/src/ GRHydro_Con2PrimHot.F90 Modified: /trunk/src/ GRHydro_Con2Prim.F90, make.code.defn Log: GRHydro: call dedicated HOT C2P routine in Con2Prim From: Christian David Ott <cott(a)zwicky-b.(none)> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: <Unavailable> File [added]: GRHydro_Con2PrimHot.F90 Delta lines: <Unavailable> File [modified]: make.code.defn Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 538)
by rhaas＠tapir.caltech.edu 06 Jul '13

06 Jul '13

User: rhaas Date: 2013/07/06 01:10 PM Modified: /trunk/ param.ccl Log: GRHydro: add GRHydro_max_temp parameter for max allowed T From: Christian David Ott <cott(a)zwicky-b.(none)> File Changes: Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +7 -2 =================================================================== --- trunk/param.ccl 2013-07-06 18:09:53 UTC (rev 537) +++ trunk/param.ccl 2013-07-06 18:10:02 UTC (rev 538) @@ -370,9 +370,14 @@ REAL GRHydro_hot_atmo_temp "Temperature of the hot atmosphere in MeV" STEERABLE=ALWAYS { - 0.0:* :: "Larger than 0 MeV" -} 0.5e0 + (0.0:* :: "Larger than 0 MeV" +} 0.1e0 +REAL GRHydro_max_temp "maximum temperature we allow" STEERABLE=ALWAYS +{ + (0.0:* :: "Larger than 0 MeV" +} 90.0e0 + REAL GRHydro_hot_atmo_Y_e "Y_e of the hot atmosphere" STEERABLE=ALWAYS { 0.0:* :: "Larger than 0"

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 537)
by rhaas＠tapir.caltech.edu 06 Jul '13

06 Jul '13

User: rhaas Date: 2013/07/06 01:09 PM Modified: /trunk/ interface.ccl Log: GRHydro: add EOS_Omni multi-var call for dpde and dprho From: Christian David Ott <cott(a)zwicky-b.(none)> File Changes: Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: +14 -0 =================================================================== --- trunk/interface.ccl 2013-06-13 19:00:01 UTC (rev 536) +++ trunk/interface.ccl 2013-07-06 18:09:53 UTC (rev 537) @@ -342,7 +342,21 @@ USES FUNCTION EOS_Omni_DPressByDRho +void FUNCTION EOS_Omni_dpderho_dpdrhoe(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_REAL IN rf_precision, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY dpderho, \ + CCTK_REAL OUT ARRAY dpdrhoe, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) +USES FUNCTION EOS_Omni_dpderho_dpdrhoe + void FUNCTION EOS_Omni_cs2(CCTK_INT IN eoskey, \ CCTK_INT IN havetemp, \ CCTK_REAL IN rf_precision, \

1 0

[svn:einsteintoolkit] www/ (Rev. 1038)
by knarf＠cct.lsu.edu 05 Jul '13

05 Jul '13

User: knarf Date: 2013/07/05 08:39 PM Added: /about/members/ latex.pl Log: small latex script to take members list and generate a small latex list Directory Changes: Directory: /svn:executable/ =========================== + * File Changes: Directory: /about/members/ ========================== File [added]: latex.pl Delta lines: +17 -0 =================================================================== --- about/members/latex.pl (rev 0) +++ about/members/latex.pl 2013-07-06 01:39:42 UTC (rev 1038) @@ -0,0 +1,17 @@ +#!/usr/bin/perl -w +use strict; + +my $itemopen = 0; +foreach my $line (<>) { + if ($line =~ /^([^ ]+) (.*)/) { + if ($itemopen) { + print " \\end{cititemize2}\n"; + } + print " \\item $2\n \\begin{cititemize2}\n"; + $itemopen = 1; + } elsif ($line =~ /^ (.*)/) { + $line =~ s/ http[^ ]*//g; + print " \\item$line"; + } +} +print " \\end{cititemize2}\n"; Property changes on: about/members/latex.pl ___________________________________________________________________

1 0

[svn:einsteintoolkit] www/download/ (Rev. 1037)
by rhaas＠tapir.caltech.edu 05 Jul '13

05 Jul '13

User: rhaas Date: 2013/07/05 12:47 PM Modified: /download/ index.php Log: tell curl to ignore github's ssl certificate File Changes: Directory: /download/ ===================== File [modified]: index.php Delta lines: +2 -2 =================================================================== --- download/index.php 2013-06-17 00:20:15 UTC (rev 1036) +++ download/index.php 2013-07-05 17:47:11 UTC (rev 1037) @@ -22,7 +22,7 @@ notes</a> for more information. <pre> -curl -O https://raw.github.com/gridaphobe/CRL/master/GetComponents +curl -kO https://raw.github.com/gridaphobe/CRL/master/GetComponents chmod a+x GetComponents @@ -40,7 +40,7 @@ as the released version. <pre> -curl -O https://raw.github.com/gridaphobe/CRL/master/GetComponents +curl -kO https://raw.github.com/gridaphobe/CRL/master/GetComponents chmod a+x GetComponents

1 0

[svn:einsteintoolkit] incoming/Symmetry/ParitySymmetry/test/ (Rev. 93)
by yosef＠astro.rit.edu 01 Jul '13

01 Jul '13

User: zlochower Date: 2013/07/01 01:04 PM Modified: /Symmetry/ParitySymmetry/test/ test.ccl Log: Set tolerances to 1e-9 File Changes: Directory: /Symmetry/ParitySymmetry/test/ ========================================= File [modified]: test.ccl Delta lines: +2 -0 =================================================================== --- Symmetry/ParitySymmetry/test/test.ccl 2013-07-01 17:31:01 UTC (rev 92) +++ Symmetry/ParitySymmetry/test/test.ccl 2013-07-01 18:04:48 UTC (rev 93) @@ -2,4 +2,6 @@ TEST KerrSchild-parity { NPROCS 2 + ABSTOL 1.0e-9 + RELTOL 1.0e-9 }

1 0

[svn:einsteintoolkit] incoming/Symmetry/ParitySymmetry/ (Rev. 92)
by yosef＠astro.rit.edu 01 Jul '13

01 Jul '13

User: zlochower Date: 2013/07/01 12:31 PM Modified: /Symmetry/ParitySymmetry/ README /Symmetry/ParitySymmetry/test/ KerrSchild-parity.par Log: Describe the thorn and give correct title to test File Changes: Directory: /Symmetry/ParitySymmetry/ ==================================== File [modified]: README Delta lines: +12 -10 =================================================================== --- Symmetry/ParitySymmetry/README 2013-06-23 02:29:19 UTC (rev 91) +++ Symmetry/ParitySymmetry/README 2013-07-01 17:31:01 UTC (rev 92) @@ -1,4 +1,4 @@ -Cactus Code Thorn RotatingSymmetry180 +Cactus Code Thorn ParitySymmetry Author(s) : Erik Schnetter <schnetter(a)cct.lsu.edu> Maintainer(s): Erik Schnetter <schnetter(a)cct.lsu.edu> Licence : GPL @@ -6,17 +6,19 @@ 1. Purpose -Provide a 180 degree rotational symmetry boundary condition. This is -the kind of symmetry that exists e.g. in a binary black hole evolution -when the two black holes are identical. The same setup also admits an -additional bitant symmetry. +Provides parity symmetry boundary condition. That is, there is a +symmetry (only) under simulatneous reflections about all coordinate +directions: (x,y,z) -> (-x, -y, -z). This is a the kind of +symmetry that exists e.g. in certain binary black hole evolution when +the momentum of the two black holes are equal in magnitude and +opposite in direction, the spins are equal in magnitude and direction, +and the masses are identical. -Currently the symmetry is fixed to be located at the lower x boundary, -so that there is a rotation about the z axis. +Currently the symmetry is fixed to be located at the lower z boundary. The coordinate thorn does not know about this symmetry. You will have -to set up the domain such that you have an upper bitant in the x +to set up the domain such that you have an upper bitant in the z direction. -This symmetry boundary condition can be mixed with arbitrary symmetry -conditions in the z direction, e.g. bitant or periodicity. +This symmetry boundary condition should not be mixed with other symmetry +conditions. Directory: /Symmetry/ParitySymmetry/test/ ========================================= File [modified]: KerrSchild-parity.par Delta lines: +1 -1 =================================================================== --- Symmetry/ParitySymmetry/test/KerrSchild-parity.par 2013-06-23 02:29:19 UTC (rev 91) +++ Symmetry/ParitySymmetry/test/KerrSchild-parity.par 2013-07-01 17:31:01 UTC (rev 92) @@ -1,4 +1,4 @@ -Cactus::cctk_run_title = "Kerr-Schild with a 180 degree rotational symmetry" +Cactus::cctk_run_title = "Kerr-Schild with parity symmetry" Cactus::cctk_full_warnings = yes Cactus::highlight_warning_messages = no

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits