Users August 2010

users@lists.einsteintoolkit.org

8 participants
19 discussions

Improve finite differencing in Exact
by Erik Schnetter 16 Nov '10

16 Nov '10

The enclosed patch improves the finite differencing that is used in Exact to calculate the extrinsic curvature. These finite differences are independent of the Cactus grid spacing. This is prompted by a remark by Peter Diener that the current finite differencing accuracy for the extrinsic curvature from thorn Exact is only about nine digits. It introduces new parameters "exact_order" and "exact_eps" that choose the order (2nd or 4th) and stencil size (default 1e-6). It also introduces new functions Exact__metric_deriv and Exact__metric_deriv2 that evaluate these finite differences. Some finite differencing terms for the spatial derivative of the shift have to be rewritten, which may lead to differences on the order of the finite differencing accuracy. The code could be rewritten to evaluate each stencil point only once in all calls to the new finite differencing operators. This would reduce the number of evaluations. Since thorn Exact is often very fast, this does not seem necessary at this point. Unfortunately, this patch causes five of the AHFinderDirect test cases to fail. -erik -- Erik Schnetter <schnetter(a)cct.lsu.edu> http://www.cct.lsu.edu/~eschnett/

2 2

fix WeylScal4::teukolskyID testsuite
by Roland Haas 09 Sep '10

09 Sep '10

Hello all, Ian Hinder found that the testsuite teukolskyID for WeylScal4 was poisoned. It turns out I had forgotten to add MoL to the list of ActiveThorns in the testsuite. Attached please find a patch which rectifies that. Ok to apply (both development and release branch)? Expect some more patches while Tanja and I try and add some extra statements to WeylScal4 to make sure MoL is present (ie. something like inherit from MoL maybe) and adapt to a possible future change in how MoL_PseudoEvolution is scheduled. Ian also pointed out some issues with how (outer hopefully) boundaries are handled which will have to be looked at eventually. Yours, Roland

4 9

Meeting reminder
by Frank Loeffler 06 Sep '10

06 Sep '10

Hi, This is a short reminder about the regular EinsteinToolkit phone call later today (at 10am CDT). Instructions how to join the call can be found here: http://einsteintoolkit.org/info/support/ Frank Loeffler

3 4

Meeting minutes
by Roland Haas 30 Aug '10

30 Aug '10

eeting minutes: Present were: Tanja, Roland, Frank, Peter, Scott, Josh, Bruno, Yosef * Frank reported of the Carpet workshop (see http://ccrg.rit.edu/~carpet/wiki/Main_Page) * GRHydro support for thorn staticconformal has been removed * new features to Carpet will be only added to the Mercurial version, users who use the current git checkout from the Einstein Toolkit will not benefit from them * loose ends from the Carpet workshop should be tied up within one or two weeks (3d output without buffer zones) * Scott and Josh reported on progress adding MHD to GRHydro which is summarized in the Carpet workshop wiki Yours, Roland -- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://keys.gnupg.net.

1 0

Removing dependence of GRHydro on StaticConformal
by Christian D. Ott 29 Aug '10

29 Aug '10

Hi, I suggest to remove the dependence of GRHydro on StaticConformal and to remove all code in GRHydro that assumes the presence of a static conformal factor. The code that would allow the use of a static conformal factor was never completed and has been disfunctional for years. Objections? - Christian

2 2

[cott@tapir.caltech.edu: [Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 152)]
by Christian D. Ott 27 Aug '10

27 Aug '10

----- Forwarded message from cott(a)tapir.caltech.edu ----- To: commits(a)einsteintoolkit.org From: cott(a)tapir.caltech.edu Date: Fri, 27 Aug 2010 20:30:51 GMT Subject: [Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 152) User: cott Date: 2010/08/27 03:30 PM Modified: /trunk/ interface.ccl, interface.ccl.omni /trunk/src/ GRHydro_Boundaries.F90, GRHydro_Con2Prim.F90, GRHydro_EoSChangeGamma.F90, GRHydro_FluxSplit.F90, GRHydro_HLLE.F90, GRHydro_HLLEPoly.F90, GRHydro_Marquina.F90, GRHydro_PPM.F90, GRHydro_ParamCheck.F90, GRHydro_Prim2Con.F90, GRHydro_Reconstruct.F90, GRHydro_ReconstructPoly.F90, GRHydro_RiemannSolve.F90, GRHydro_RoeSolver.F90, GRHydro_Source.F90, GRHydro_Tmunu.F90, GRHydro_UpdateMask.F90, Utils.F90 Log: * remove dependence on StaticConformal * change calculation of the determinant of the 3-metric from a subroutine call to a macro. File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Boundaries.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2Prim.F90 Delta lines: <Unavailable> File [modified]: GRHydro_EoSChangeGamma.F90 Delta lines: <Unavailable> File [modified]: GRHydro_FluxSplit.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLE.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLEPoly.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Marquina.F90 Delta lines: <Unavailable> File [modified]: GRHydro_PPM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ParamCheck.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2Con.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Reconstruct.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ReconstructPoly.F90 Delta lines: <Unavailable> File [modified]: GRHydro_RiemannSolve.F90 Delta lines: <Unavailable> File [modified]: GRHydro_RoeSolver.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Source.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Tmunu.F90 Delta lines: <Unavailable> File [modified]: GRHydro_UpdateMask.F90 Delta lines: <Unavailable> File [modified]: Utils.F90 Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: interface.ccl Delta lines: <Unavailable> File [modified]: interface.ccl.omni Delta lines: <Unavailable> _______________________________________________ Commits mailing list Commits(a)einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/commits ----- End forwarded message -----

1 0

Carpet 3d output
by Michael Koppitz 25 Aug '10

25 Aug '10

Hi all, I made the observation that Carpet's output is quite large compared to the grid structure: on 100^3 gridpoints per refinement level dumping admbase::metric (6 gridfunctions) in 3d hdf5 I would expect for 1 ref-level 10^6*6*8b = 48Mb. I got 80Mb. I assume that is because Carpet also includes x, y, z, and t in the dump, which brings it to 10 GF's and 80Mb expected. But for 2 levels I got 250Mb. For 3 levels I got 450Mb and for 4 levels I got 620Mb. This seems much even if I take ghost- and buffer-zones into account. Am I missing something, or is Carpet doing something unexpected? thank you Michael

2 1

Meeting minutes
by Roland Haas 17 Aug '10

17 Aug '10

Meeting minutes: Next meeting: There will be no phone call next week (Monday August 23rd). Present were: Frank, Josh, Tanja, Roland, Yosef, Peter, Druno * CactusTest: no objection to include it once it passes on a set of common machines. Will include it without further discussion at that point. * MHD: no progress * EOS_Omni: ** Frank has added some code to GRHydro to make progress towards using EOS_Onmi only Yours, Roland -- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://keys.gnupg.net.

3 4

Trying EOS_Omni with simfactory
by Frank Loeffler 16 Aug '10

16 Aug '10

Hi, I keep forgetting this myself, and it might also be interesting for others: if you try EOS_Omni currently and want to use simfactory, this is how you can do this with one command (point to your thornlist, including EOS_Omni and _not_ including the other EOS thorns): ./simfactory/sim build --thornlist=../thornlists/Omni.th --attach=FPPFLAGS=-DUSE_EOS_OMNI=1 --attach=CPPFLAGS=-DUSE_EOS_OMNI=1 Frank

1 0

Re: [Users] [Commits] [svn:einsteintoolkit] manifest/trunk/ (Rev. 35)
by Frank Loeffler 12 Aug '10

12 Aug '10

On Wed, Aug 11, 2010 at 01:37:59AM +0000, schnetter(a)cct.lsu.edu wrote: > Add CactusTest arrangement, which contains many good test cases. Do all the testcases pass? Frank

2 4

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Users August 2010