[Users] reduction in carpet

[Hee Il Kim]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20110105/21fb925c/attachment.html