Users January 2011

users@lists.einsteintoolkit.org

13 participants
26 discussions

GRHydro prim2con
by Ernazar Abdikamalov 07 Jan '11

07 Jan '11

GRHydro doesn't compile on numrel08. It is complaining about a type mismatch: In routine primitive2conservative, variable xtemp is defined as CCTK_REAL :: xtemp and is used as an output variable in routine prim2con_hot. However, inside prim2con_hot, the corresponding variable is defined as real*8 temp0(1). I don't know exactly why this doesn't compile. However, if, for example, one defines xtemp in routine primitive2conservative as CCTK_REAL :: xtemp(1), then it compiles. I've attached the corresponding patch to this message. Is it okay to apply it? Thanks, Ernazar

2 1

Re: [Users] reduction in carpet + mismatched declaration in ML_ADMQuantities
by Hee Il Kim 07 Jan '11

07 Jan '11

Christian, thanks so much ! Your code seems to contain everything I wanted to calculate. There's one minor mismatched declaration in ML_ADMQuantities which should be fixed anyhow. #### [khi@zenith src]$ diff ML_ADMQuantities.c.ori ML_ADMQuantities.c 56c56 < const char *groups[] = {"grid::coordinates","Grid::coordinates","McLachlan::ML_curv","McLachlan::ML_Gamma","McLachlan::ML_lapse","McLachlan::ML_log_confac","McLachlan::ML_metric","McLachlan::ML_shift","McLachlan::ML_trace_curv","ML_ADMQuantities::ML_Jadm","ML_ADMQuantities::ML_Madm"}; --- > const char *groups[] = {"grid::coordinates","Grid::coordinates","ML_BSSN::ML_curv","ML_BSSN::ML_Gamma","ML_BSSN::ML_lapse","ML_BSSN::ML_log_confac","ML_BSSN::ML_metric","ML_BSSN::ML_shift","ML_BSSN::ML_trace_curv","ML_ADMQuantities::ML_Jadm","ML_ADMQuantities::ML_Madm"}; #### Thanks again, Hee Il 2011/1/5 Christian D. Ott <cott(a)tapir.caltech.edu> > > Dear Hee Il, > > I agree that getting this right is tricky. You must set up your grid > functions that you want to reduce in MoL_PseudoEvolution in local > mode, then do the reduction in CCTK_Analysis in global mode. > > I am attaching a Thorn called ZelmaniAnalysis that I use for the > kinds of things you want. The reduction works with Carpet. > > Best, > > - Christian > > > > On Wed, Jan 05, 2011 at 09:18:21PM +0900, Hee Il Kim wrote: > > Hi, > > > > Belated Happy New Year ! > > > > I am writing a simple code for the calculation of ADM mass and angular > > momentum by using ML_ADMQuantities. As you know, it requires reduction > > procedure for the summation of the integrand, e.g., density * volume > > measure. I've tried various reduction APIs including old deprecated one, > but > > I could only succeed in the integration only over the coarsest grid > level, > > which was not accurate. > > > > If I use the new reduction API in the reference manual of Cactus, > > CCTK_ReduceGridArrays; > > > > CCTK_ReduceGridArrays(const cGH *GH, > > int dest_processor, > > int local_reduce_handle, > > int param_table_handle, > > int N_input_arrays, > > const CCTK_INT input_array_variable_indices[], > > int M_output_values, > > const CCTK_INT output_value_type_codes[], > > void* const output_values[]); > > > > , I encounter the following error > > .... > > WARNING[L3,P0] (Cactus): CCTK_ReduceGridArrays: no grid array reduction > > registered > > WARNING level 0 in thorn ML_Analysis processor 0 host compute-0-5.local > > (line 159 of Evaluate_ADMQuantities.F90): > > -> Failed to integrate J_z > > > > Param_table_handle is not clear to me but it seems not causing the > problem. > > The reduction operators are registered by CactusBase/LocalReduce. > > > > Could you let me know which reduction function I have to use for this > > calculation? I mean the integration over the whole grid levels. Should I > do > > loop over refinenment levels; > > ... > > do loop over reflev > > call CCTK_Reduce > > enddo > > ... > > > > Thanks in advance, > > > > Hee Il > > > _______________________________________________ > > Users mailing list > > Users(a)einsteintoolkit.org > > http://lists.einsteintoolkit.org/mailman/listinfo/users > >

4 5

CarpetIOScalar: avoiding "reductions" for scalar variables
by Kelly, Bernard J. (GSFC-660.0)[UNIVERSITY OF MARYLAND BALTIMORE COUNTY] 05 Jan '11

05 Jan '11

Hi all. This seems trivial, but ... How do I avoid reductions (average, minimum, maximum, etc) for scalar output of a SCALAR variable in CarpetIOScalar? CactusBase/IOBasic recognises the SCALAR type of the variable, and doesn't bother with reduction. Carpet/CarpetIOScalar seems to ignore the type, and does the reductions. At least, it outputs several files with trivial information in them. I can force it to only generate one file by setting outscalar_reductions ="average" (for instance), but then I lose the reductions for all other (non-SCALAR) variables. Thanks, Bernard [working with the Chandrasekhar release of ETK]

2 3

Re: [Users] [ET Trac] [Einstein Toolkit] #103: Handle parameter files with DOS newline convention
by Frank Loeffler 05 Jan '11

05 Jan '11

On Thu, Dec 30, 2010 at 03:56:47PM -0000, Einstein Toolkit wrote: > Try the attached patch - it fixes the problem for me. The patch looks fine and works. Thanks for contributing, please commit. Frank

1 0

reduction in carpet
by Hee Il Kim 05 Jan '11

05 Jan '11

Hi, Belated Happy New Year ! I am writing a simple code for the calculation of ADM mass and angular momentum by using ML_ADMQuantities. As you know, it requires reduction procedure for the summation of the integrand, e.g., density * volume measure. I've tried various reduction APIs including old deprecated one, but I could only succeed in the integration only over the coarsest grid level, which was not accurate. If I use the new reduction API in the reference manual of Cactus, CCTK_ReduceGridArrays; CCTK_ReduceGridArrays(const cGH *GH, int dest_processor, int local_reduce_handle, int param_table_handle, int N_input_arrays, const CCTK_INT input_array_variable_indices[], int M_output_values, const CCTK_INT output_value_type_codes[], void* const output_values[]); , I encounter the following error .... WARNING[L3,P0] (Cactus): CCTK_ReduceGridArrays: no grid array reduction registered WARNING level 0 in thorn ML_Analysis processor 0 host compute-0-5.local (line 159 of Evaluate_ADMQuantities.F90): -> Failed to integrate J_z Param_table_handle is not clear to me but it seems not causing the problem. The reduction operators are registered by CactusBase/LocalReduce. Could you let me know which reduction function I have to use for this calculation? I mean the integration over the whole grid levels. Should I do loop over refinenment levels; ... do loop over reflev call CCTK_Reduce enddo ... Thanks in advance, Hee Il

2 1

LONI allocations (demonstration accounts); allocation switch-over
by Erik Schnetter 02 Jan '11

02 Jan '11

Some Einstein Toolkit users have access to the loni_cactus05 allocation, e.g. via a demonstration account. Note that this allocation (loni_cactus05) has run out at the end of 2010, and we expect to receive a renewal (loni_cactus06) within the next week (or two). If you used to have access to loni_cactus05, you can request access to loni_cactus06 at <http://allocations.loni.org/>. You will need to wait until you see loni_cactus06 listed there before requesting access. -erik -- Erik Schnetter <schnetter(a)cct.lsu.edu> http://www.cct.lsu.edu/~eschnett/

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Users January 2011