Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

New changesets in Carpet Mercurial
by Carpet mercurial repositories 07 Jun '12

07 Jun '12

changeset: 3573:0fc51a62a32b tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Thu Jun 07 17:20:59 2012 -0500 summary: Carpet: Update setting CPU affinity

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 07 Jun '12

07 Jun '12

changeset: 3572:48f8d1c9a36b tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Thu Jun 07 15:19:38 2012 -0500 summary: CarpetLib: Output error message for illegal bbox setups

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 07 Jun '12

07 Jun '12

changeset: 3571:ac912892433c tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Thu Jun 07 14:57:49 2012 -0500 summary: Carpet: Correct error in "balanced" domain decomposition

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 07 Jun '12

07 Jun '12

changeset: 3570:bed7a2e1b4a9 tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Thu Jun 07 14:53:59 2012 -0500 summary: Carpet: Add short instructions to backtraces

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 07 Jun '12

07 Jun '12

changeset: 3569:18d226c24a35 tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Thu Jun 07 14:50:29 2012 -0500 summary: Carpet: Iterate over timelevels in post_recover_variables

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 07 Jun '12

07 Jun '12

changeset: 3568:eefe46d6982a tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Thu Jun 07 13:42:03 2012 -0500 summary: Carpet: Separate creating and initialising Timers

1 0

[svn:einsteintoolkit] GRHydro/trunk/test/ (Rev. 353)
by rhaas＠tapir.caltech.edu 06 Jun '12

06 Jun '12

User: rhaas Date: 2012/06/06 09:37 AM Added: /trunk/test/ test_tov_carpet_refined_nosync.par /trunk/test/test_tov_carpet_refined_nosync/ rho.x.asc, rho.y.asc, rho.z.asc Log: GRHydro: add files missing from r340 Thanks to Ian Hinder for pointing these (and the last one) out. File Changes: Directory: /trunk/test/ ======================= File [added]: test_tov_carpet_refined_nosync.par Delta lines: <Unavailable> Directory: /trunk/test/test_tov_carpet_refined_nosync/ ====================================================== File [added]: rho.x.asc Delta lines: <Unavailable> File [added]: rho.y.asc Delta lines: <Unavailable> File [added]: rho.z.asc Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/ (Rev. 352)
by rhaas＠tapir.caltech.edu 06 Jun '12

06 Jun '12

User: rhaas Date: 2012/06/06 09:35 AM Added: /trunk/src/ GRHydro_LastMoLPostStep.c Log: GRHydro: add file missing from r343 Directory Changes: Directory: /svn:mime-type/ ========================== + text/x-csrc File Changes: Directory: /trunk/src/ ====================== File [added]: GRHydro_LastMoLPostStep.c Delta lines: +41 -0 =================================================================== --- trunk/src/GRHydro_LastMoLPostStep.c (rev 0) +++ trunk/src/GRHydro_LastMoLPostStep.c 2012-06-06 14:35:52 UTC (rev 352) @@ -0,0 +1,41 @@ +// GRHydro_LastMoLPostStep.c +// +// Compute is this is the last MoL PostStep call. Code taken from Christian +// Reisswig's rejected MoL changes. +// +// Roland Haas +// Sun Jun 3 17:35:53 PDT 2012 + +#include <stdio.h> +#include <string.h> +#include <assert.h> + +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" + +void GRHydro_SetLastMoLPostStep(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + + const CCTK_INT *MoL_Intermediate_Step = + CCTK_VarDataPtr(cctkGH,0,"MoL::MoL_Intermediate_Step"); + if(NULL == MoL_Intermediate_Step) + { + CCTK_WARN(CCTK_WARN_ABORT, + "Could not get data pointer for MoL::MoL_Intermediate_Step"); + } + + // If counter becomes zero, the only thing left to do is to call PostStep! + if (*MoL_Intermediate_Step == 0) + *InLastMoLPostStep = 1; + else + *InLastMoLPostStep = 0; +} + +void GRHydro_ClearLastMoLPostStep(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + + *InLastMoLPostStep = 0; +} Property changes on: trunk/src/GRHydro_LastMoLPostStep.c ___________________________________________________________________

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 351)
by rhaas＠tapir.caltech.edu 05 Jun '12

05 Jun '12

User: rhaas Date: 2012/06/05 03:51 PM Modified: /trunk/ param.ccl /trunk/src/ GRHydro_PreLoop.F90, GRHydro_Scalars.F90, GRHydro_TVDReconstruct.F90, make.code.defn Log: GRHydro: Fix slope limiters and get rid of broken cases; there is a related patch in TRAC at the moment Patch by Josh Faber. Some additional checks by Roland Haas. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_PreLoop.F90 Delta lines: +0 -9 =================================================================== --- trunk/src/GRHydro_PreLoop.F90 2012-06-05 20:50:25 UTC (rev 350) +++ trunk/src/GRHydro_PreLoop.F90 2012-06-05 20:51:49 UTC (rev 351) @@ -39,9 +39,6 @@ DECLARE_CCTK_FUNCTIONS MINMOD = .false. - MINMOD2 = .false. - MINMOD3 = .false. - MC1 = .false. MC2 = .false. SUPERBEE = .false. HLLE = .false. @@ -49,12 +46,6 @@ if (CCTK_EQUALS(tvd_limiter,"minmod")) then MINMOD = .true. - else if (CCTK_EQUALS(tvd_limiter,"minmod2")) then - MINMOD2 = .true. - else if (CCTK_EQUALS(tvd_limiter,"minmod3")) then - MINMOD3 = .true. - else if (CCTK_EQUALS(tvd_limiter,"vanleerMC")) then - MC1 = .true. else if (CCTK_EQUALS(tvd_limiter,"vanleerMC2")) then MC2 = .true. else if (CCTK_EQUALS(tvd_limiter,"Superbee")) then File [modified]: GRHydro_Scalars.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Scalars.F90 2012-06-05 20:50:25 UTC (rev 350) +++ trunk/src/GRHydro_Scalars.F90 2012-06-05 20:51:49 UTC (rev 351) @@ -16,7 +16,7 @@ implicit none - LOGICAL :: MINMOD, MINMOD2, MINMOD3, MC1, MC2, SUPERBEE, PPM3, PPM4 + LOGICAL :: MINMOD, MC2, SUPERBEE, PPM3, PPM4 LOGICAL :: ANALYTICAL LOGICAL :: FAST LOGICAL :: HLLE, LLF File [modified]: GRHydro_TVDReconstruct.F90 Delta lines: +18 -14 =================================================================== --- trunk/src/GRHydro_TVDReconstruct.F90 2012-06-05 20:50:25 UTC (rev 350) +++ trunk/src/GRHydro_TVDReconstruct.F90 2012-06-05 20:51:49 UTC (rev 351) @@ -79,22 +79,26 @@ dloc = orig(i+xoffset, j+yoffset, k+zoffset) - orig(i, j, k) if (MINMOD) then - delta = minmod_func(dupw,dloc) + delta = minmod_func(dupw,dloc) else if (MC2) then -!!$ This is an alternative equivalent implementation -!!$ of vanLeeer MC slopelimiter - if (dupw*dloc < 0.d0) then - delta=0.d0 - else - delta=sign(min(2.d0*abs(dupw),2.d0*abs(dloc),& - 0.5d0*(abs(dupw)+abs(dloc))),dupw+dloc) - end if +!!$ This is the van Leer MC slope limiter + if (dupw*dloc < 0.d0) then + delta=0.d0 + else + delta=sign(min(2.d0*abs(dupw),2.d0*abs(dloc),& + 0.5d0*(abs(dupw)+abs(dloc))),dupw+dloc) + end if + else if (SUPERBEE) then + if (dupw*dloc < 0.d0) then + delta=0.d0 + else + delta=sign(max(min(2.d0*abs(dupw),abs(dloc)),& + min(2.d0*abs(dloc),abs(dupw))),dupw+dloc) + end if else - delta = 0.5d0*(dupw + dloc) - if (abs(dupw) < myfloor ) dupw = myfloor*sign(1.d0, dupw) - if (abs(dloc) < myfloor ) dloc = myfloor*sign(1.d0, dloc) - ratio = dupw / dloc - call slopelimiter(ratio, delta) + call CCTK_WARN(0, "Type of limiter not recognized") + ! NOTREACHED + delta = 0d0 end if hdelta = 0.5d0 * delta bextm(i, j, k) = orig(i, j, k) - hdelta File [modified]: make.code.defn Delta lines: +0 -1 =================================================================== --- trunk/src/make.code.defn 2012-06-05 20:50:25 UTC (rev 350) +++ trunk/src/make.code.defn 2012-06-05 20:51:49 UTC (rev 351) @@ -35,7 +35,6 @@ GRHydro_RoeSolver.F90 \ GRHydro_Scalars.F90 \ GRHydro_Shift.F90 \ - GRHydro_SlopeLimiter.F90 \ GRHydro_Source.F90 \ GRHydro_Startup.F90 \ GRHydro_TVDReconstruct.F90 \ Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +0 -3 =================================================================== --- trunk/param.ccl 2012-06-05 20:50:25 UTC (rev 350) +++ trunk/param.ccl 2012-06-05 20:51:49 UTC (rev 351) @@ -116,10 +116,7 @@ keyword tvd_limiter "Which slope limiter to use" { "minmod" :: "Minmod" - "vanleerMC" :: "Van Leer MC - Ian" "vanleerMC2" :: "Van Leer MC - Luca" - "minmod2" :: "Minmod - Ian" - "minmod3" :: "Minmod - Luca (or vice versa). Expect I've made a typo" "Superbee" :: "Superbee" } "minmod"

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 350)
by rhaas＠tapir.caltech.edu 05 Jun '12

05 Jun '12

User: rhaas Date: 2012/06/05 03:50 PM Modified: /trunk/src/ GRHydro_RegisterVars.cc Log: GRHydro: register TmunuBase variables as MoL constraint. Patch by Christian Reisswig. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_RegisterVars.cc Delta lines: +3 -0 =================================================================== --- trunk/src/GRHydro_RegisterVars.cc 2012-06-05 20:49:41 UTC (rev 349) +++ trunk/src/GRHydro_RegisterVars.cc 2012-06-05 20:50:25 UTC (rev 350) @@ -59,6 +59,9 @@ register_constrained("HydroBase::eps"); register_constrained("HydroBase::w_lorentz"); register_constrained("HydroBase::vel"); + register_constrained("TmunuBase::stress_energy_scalar"); + register_constrained("TmunuBase::stress_energy_vector"); + register_constrained("TmunuBase::stress_energy_tensor"); if (general_coordinates) { register_constrained("GRHydro::lvel");

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits