Users

users@lists.einsteintoolkit.org

3 participants
3984 discussions

Debugging patch for SetMask_SphericalSurface
by Erik Schnetter 02 Aug '10

02 Aug '10

Christian The enclosed patch may help debug SetMask_SphericalSurface. It checks a few consistency conditions. -erik -- Erik Schnetter <schnetter(a)cct.lsu.edu> http://www.cct.lsu.edu/~eschnett/ $ svn diff Index: param.ccl =================================================================== --- param.ccl (revision 107) +++ param.ccl (working copy) @@ -2,7 +2,7 @@ STRING SetMask_MaskName "Full variable name of the mask to be set" { - ".+" :: "Any full Cactus variable name, GF of type CCTK_BYTE" + ".+" :: "Any full Cactus variable name, GF of type CCTK_INT" } "HydroBase::hydro_excision_mask" #BOOLEAN SetMask_ResetAll "'yes' (default): Overwrite complete mask, 'no': Only set excision points" Index: src/SetMask.c =================================================================== --- src/SetMask.c (revision 107) +++ src/SetMask.c (working copy) @@ -10,9 +10,18 @@ DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; - CCTK_INT *mask = (CCTK_INT*) CCTK_VarDataPtr(cctkGH, 0, SetMask_MaskName); + int varindex = CCTK_VarIndex(SetMask_MaskName); + if (varindex<0) + CCTK_WARN(0, "No such variable"); + int grouptype = CCTK_GroupTypeFromVarI(varindex); + if (grouptype!=CCTK_GF) + CCTK_WARN(0, "Wrong group type, must be GF"); + int vartype = CCTK_VarTypeI(varindex); + if (vartype!=CCTK_VARIABLE_INT) + CCTK_WARN(0, "Wrong variable type, must be CCTK_INT"); + CCTK_INT *mask = (CCTK_INT*) CCTK_VarDataPtrI(cctkGH, 0, varindex); if (!mask) - CCTK_WARN(0, "No such variable, or no storage enabled"); + CCTK_WARN(0, "No storage enabled"); for (int smi = 0; smi < 10; smi++) {

1 0

Carpet developer workshop at RIT on Aug 23-25
by Bruno Coutinho Mundim 02 Aug '10

02 Aug '10

Dear Colleagues, We would like to invite you to attend a Carpet developer workshop to be hosted by the Center for Computational Relativity and Gravitation at the Rochester Institute of Technology on August 23th, 24th and 25th (2010). This workshop is aimed at experienced Carpet users that want to gain further expertise on the internal workings of the Carpet driver, and, thereby, to be able to contribute to the infrastructure development. This workshop will be led by Erik Schnetter, the main architect for the Carpet AMR infrastructure, and will count with the presence of an HPC specialist. One important goal of this workshop will be to discuss and implement ways to improve the performance for Cactus-based applications for current and future high-end petascale architectures. Due to space limitations, we encourage you to contact the organizers as soon as possible to register your presence. Sincerely, The organizers: Manuela Campanelli (manuela(a)astro.rit.edu) Bruno C. Mundim (bcmsma(a)astro.rit.edu) Erik Schnetter (schnetter(a)cct.lsu.edu)

1 1

Meeting reminder
by Frank Loeffler 02 Aug '10

02 Aug '10

Hi, this is a short reminder about the regular EinsteinToolkit phone call later today (at 10am CDT, in about 15 minutes). Instructions how to join the call can be found here: http://einsteintoolkit.org/info/support/ I don't expect this to be a long meeting, also due to the expected absence of a couple of people because of the TeraGrid conference. Frank Loeffler

2 1

prim2con in HydroBase_Prim2ConInitial
by Frank Loeffler 29 Jul '10

29 Jul '10

Hi, GRHydro currently does not convert the primitive variables the ID thorns populate to the conserved quantities. However, since poststep is now scheduled also for initial data, it converts the conserved quantities to the primivite ones before evolution. This means ID thorns have to setup the conserved quantities to make this work. Currently e.g. TOVSolver does this, but it really should not. HydroBase defines a group HydroBase_Prim2ConInitial exactly for this purpose: to have any evolution thorn scheduling prim2con routines it might need to use ID given in only primary variables defined in HydroBase. GRHydro does not schedule anything there. I suggest to schedule the corresponding GRHydro prim2con routines in this group. I tried this and all testsuites still pass - but ID thorns no longer have to specify the conserved quantities anymore. Are there objections? Frank

1 0

simfactory: adding new machine?
by Bernard Kelly 29 Jul '10

29 Jul '10

Hi all. I'm trying out einsteintoolkit + simfactory on a non-standard machine, and am failing rather quickly. I followed the simfactory tutorial advice: https://svn.cct.lsu.edu/repos/numrel/simfactory/tutorial/index.html ... which said to edit the section of udb.pm entitled "Customize your local machine here". Since this heading no longer appears, I imagine the tutorial is slightly out of date. Anyway, changing "Options for new machines" along these lines doesn't seem to work -- if I run ./sim, I get the following error message: ------------------------ Uncaught exception from user code: Wile setting user option [discover][nickname]="discover": Machine "discover" does not exist at udb.pm line 123. ------------------------ What's the proper method for adding a new machine? Do I actually have to add to the "mdb.pm" file? Thanks, Bernard -- ----------------------------------------------- Bernard J. Kelly NASA Goddard Space Flight Center, Code 663 8800 Greenbelt Road Greenbelt, MD 20771, U.S.A. phone: +1 (301) 286-7243 -----------------------------------------------

4 4

Meeting minutes
by Roland Haas 26 Jul '10

26 Jul '10

Meeting minutes: Present were: Frank, Roland, Gabrielle, Bruno, Ian, Eloisa * Frank reported on patches applied to the release version of the toolkit * we agreed to use EOS_Omni in the current set of initial data tests * a partial list of thorns/modules for which there are papers to be cited is: ** Kranc ** Carpet ** AHFinderDirect ** Cactus flesh ** TwoPunctures ** GRHydro (citing the Whisky paper on whiskycode.org) ** McLachlan * a draft of the citation policy has been send to the maintainers and will have to be discussed on the mailing list * Gabrielle will draft a website to list these publications We have one item for next the meeting/mailing list discussions: * decide on dates/locations for toolkit workshops/conferences 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

Meeting reminder
by Frank Loeffler 19 Jul '10

19 Jul '10

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

2 1

EOS_Omni
by Christian D. Ott 19 Jul '10

19 Jul '10

Dear All, you can give the new EOS for GRHydro a try by recompiling with -DUSE_EOS_OMNI in CPPFLAGS and FPPFLAGS. This makes GRHydro _independent_ of EOS_Base and all EOS thorns except EOS_Omni (dependence via aliased functions). Let's say you want to evolve a Gamma=2 polytrope with K=100 with a gamma-law. For this, you would set: ActiveThorns = "EOS_Omni" grhydro::grhydro_poly_eoskey = 1 grhydro::grhydro_eoskey = 2 eos_omni::gl_gamma = 2.0 eos_omni::gl_k = 100.0d0 eos_omni::poly_gamma = 2.0 eos_omni::poly_gamma_ini = 2.0 eos_omni::poly_k = 100.0d0 Note that gl_ controls the gamma-law settings and that poly settings are still needed (for fallback purposes in GRHydro if c2p fails). Without the preprocessor defines, the code behaves as usual and is wedded to EOS_Base and EOS_Polytrope. All test suites pass. A test suite using EOS_Omni has been added. If EOS_Omni is activated and EOS_Omni parameters are added, all standard testsuites pass. This confirms that EOS_Omni produces the same results as EOS_Base + EOS_Polytrope/EOS_Ideal_Fluid/EOS_Hybrid. Currently, finite-temperature EOS support is lacking. This requires messing with various routines in GRHydro to include temperature in calls etc. This will be the next step in the EOS_Omni development. Since EOS_Omni already provides the entire functionality of the EOS_Base interface, I would suggest to remove this old interface from the code soon (once people have tried out and are happy with EOS_Omni). Please report any problems you may find with EOS_Omni. Best, - Christian

5 10

CarpetIOHDF5: checkpoint initial data
by Baiotti Luca 19 Jul '10

19 Jul '10

Hello, I found a minor unexpected behaviour (or at least inaccurate error message) of CarpetIOHDF5. When setting among other parameters: Cactus::terminate = time Cactus::cctk_final_time = 0.0 IO::checkpoint_on_terminate= yes IO::checkpoint_ID = no the below message is printed: "INFO (CarpetIOHDF5): Termination checkpoint already dumped as last evolution checkpoint at iteration 0" but no checkpoint is dumped. I would think that the checkpoint should be dumped with this choice of parameters, but this is just my point of view. Anyway, the info message is wrong. Luca

2 1

Difference in NaN checks results of GRHydro
by Hee Il Kim 18 Jul '10

18 Jul '10

Hi, Recently, I found the nancheks of GRHydro show different results sometimes depending on the compilers and codes. As far as I've tested, nans appear only if I use the private thorn GRHydro_RNSID (just a clone of Whisky_RNSID) with XL compilers on AIX machines. I don't see any nans with Intel compilers for the same parfile. Also old Whisky with BSSN_MoL passes the nanchecks safely. The difference happens actually for somewhat specific situations. It only happens for differentially rotating stars. It disappears if I take a smaller # of refinement levels. I've tested it applying the IEEE strict floating point options to both compilers (Intel and XL), but the results were same. Since the problem doesn't happen for Whisky+BSSN_MoL, I think it is not the problem of the private thorn only. My GRHydro_RNSID is basically the same with Whisky_RNSID. Assuming most of you are using the similar thorn, I attached the thornlist and the parfile. I hope you can examine and fix the problem. Regards, Hee Il

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Users