[Users] reduction in carpet

[Christian D. Ott]

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 at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users

-------------- next part --------------
A non-text attachment was scrubbed...
Name: ZelmaniAnalysis.tar.gz
Type: application/x-tar-gz
Size: 9321 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20110105/d28cd121/attachment.bin