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