[Users] Is there a better way to deal with grid arrays in ETK?

[Roland Haas]

Hello Yu Liu,

> It is introduced in the UsersGuide of Cactus that the grid arrays can only
> be determined by function calls, rather than the quicker methods available
> for GFs. But is there a better way to deal with the grid arrays so far?
> 
> I am trying to analyze some physical quantities in equatorial plane, and
> get the grid function in spherical coordinates by interpolation. It doesn't
> seem so convenient, if I define this new variable as grid arrays,
In C code there is no other way to get the exact shape of the arrays.
In Fortran code the shape is passed as usual to Fortran and the grid
arrays appear as regular multi-D Fortran arrays.

In C what is often done is that the arrays are shaped based on some
parameters eg Npoints_theta, Npoints_phi and then code uses those as
the array bounds. This only works for DISTRIB=const arrays which are
the same size on each MPI rank and not for DISTRIB=default arrays which
are split among the MPI ranks and distributed that way. For those
latter ones I do not think that there is any way to get at the shape
other than calling GroupDynamicData
(http://einsteintoolkit.org/referencemanual/ReferenceManualch2.html#x4-77000A2)
on them.

Yours,
Roland

-- 
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://pgp.mit.edu .
-------------- 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/20200514/803f726d/attachment-0001.bin