On Thu, Nov 21, 2013 at 1:59 AM, Roland Haas <roland.haas@physics.gatech.edu> wrote:
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hello Vassili,

>> how can i now perform the reduction on this "processor-local"
>> quantity?
> You can use the same as for the first example that I gave:
Also, before I forget again: the reduction result will only be correct
at "coarse" timesteps ie when all refinement levels are aligned in
time. For the other steps you will mix current and past values from
coarse and fine levels (you are doing a 0th order interpolation in
time). Carpet will to 2nd order interpolation in time for reductions
if not all data is aligned in time.

yes, i already noticed that when trying things out..

i also got the switching between levels sorted out (i hope, al least there are no more assertions faults...)

basically my question is whether the local integration of the scalar, ie

LOOP{
rho_in_shell[iradii] += BLAH


}  

that has to be run as LOOP-LOCAL (or using the macros to manually switch the modes) has to be run in the SAME function as the later reduction, which has to be run as GLOBAL...

is that true? do they have to be in the same function? because then i have to switch modes again, because i want to do the reduction in a FOR loop across all iradii, and that would be local again, because the entire function integrate was scheduled as GLOBAL LOOP_LOCAL...

maybe i am missing things as well..

in the meantime, i have implemented the way used ADMBase using a loop_counter....

while i trust this method more, it will also be slower than the scalar way, because i will have to calculate the value of the now gridfunction rho_in_shell nradii times for every grid point..

best wishes,

Vassili  


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://keys.gnupg.net.
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.15 (GNU/Linux)
Comment: Using GnuPG with Icedove - http://www.enigmail.net/

iEYEARECAAYFAlKNWwEACgkQTiFSTN7SboWXlgCgvkqmENaO7XBDup6cxl2ORXmN
LBYAoJPKHV23EKuQhqX5iPXc8pMuxPYF
=ynMM
-----END PGP SIGNATURE-----
_______________________________________________
Users mailing list
Users@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/users