[Users] Norm2 calculation on different grid levels

Haas, Roland rhaas at illinois.edu
Sat Jun 15 11:22:39 CDT 2019

Hello Elena,

if you schedule your function using:

options: level

instead of

options: global

then a call to Cactus' reduction interface will compute the reduction
only over the current refinement level.

Note that the sum of the norm2 of all the levels will not be the same
as the norm2 of all levels computed at once (because when computing
only on a single level Carpet will not "cut out" regions where a finer
grid exists).

The alternative is to introduce a helper grid function where you
manually zero out the regions you do not care about then compute the
norm of that grid function.

You can use Carpet's aliased function GetRefinementLevel to find out
the current refinement level:

CCTK_INT FUNCTION                         \
    GetRefinementLevels                   \
PROVIDES FUNCTION GetRefinementLevels \
    WITH Carpet_GetRefinementLevels   \

You can copy (or use) GenericFD's GenericFD_GetBoundaryInfo function to
get information about whether a given face is a boundary and if so what
kind of boundary. 



> Hello,
> I'm doing a simulation for BBH and I want to calculate the norm2 of the
> Hamiltonian constraints on different grid levels. Ideally i would like to
> calculate it first removing the first level of the grid, then removing it
> and the second one and so on...(the purpose is that i want to exclude the
> boundary contribution to the calculation). Do you know if there is a tool
> in the ET to do it?
> Thanks
> Elena De Paoli

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.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 833 bytes
Desc: OpenPGP digital signature
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20190615/5ea42e94/attachment.bin 

More information about the Users mailing list