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