<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Nov 21, 2013 at 1:59 AM, Roland Haas <span dir="ltr"><<a href="mailto:roland.haas@physics.gatech.edu" target="_blank">roland.haas@physics.gatech.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">-----BEGIN PGP SIGNED MESSAGE-----<br>
Hash: SHA1<br>
<br>
Hello Vassili,<br>
<br>
>> how can i now perform the reduction on this "processor-local"<br>
>> quantity?<br>
> You can use the same as for the first example that I gave:<br>
</div>Also, before I forget again: the reduction result will only be correct<br>
at "coarse" timesteps ie when all refinement levels are aligned in<br>
time. For the other steps you will mix current and past values from<br>
coarse and fine levels (you are doing a 0th order interpolation in<br>
time). Carpet will to 2nd order interpolation in time for reductions<br>
if not all data is aligned in time.<br></blockquote><div><br></div><div>yes, i already noticed that when trying things out..</div><div><br></div><div>i also got the switching between levels sorted out (i hope, al least there are no more assertions faults...)</div>
<div><br></div><div>basically my question is whether the local integration of the scalar, ie</div><div><br></div><div>LOOP{</div><div>rho_in_shell[iradii] += BLAH</div><div><br></div><div><br></div><div>} </div><div><br>
</div><div>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...</div><div><br></div><div>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...</div>
<div><br></div><div>maybe i am missing things as well..</div><div><br></div><div>in the meantime, i have implemented the way used ADMBase using a loop_counter....</div><div><br></div><div>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..</div>
<div><br></div><div>best wishes,</div><div><br></div><div>Vassili </div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im"><br>
Yours,<br>
Roland<br>
<br>
- --<br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://keys.gnupg.net" target="_blank">http://keys.gnupg.net</a>.<br>
-----BEGIN PGP SIGNATURE-----<br>
Version: GnuPG v1.4.15 (GNU/Linux)<br>
Comment: Using GnuPG with Icedove - <a href="http://www.enigmail.net/" target="_blank">http://www.enigmail.net/</a><br>
<br>
</div>iEYEARECAAYFAlKNWwEACgkQTiFSTN7SboWXlgCgvkqmENaO7XBDup6cxl2ORXmN<br>
LBYAoJPKHV23EKuQhqX5iPXc8pMuxPYF<br>
=ynMM<br>
<div class="HOEnZb"><div class="h5">-----END PGP SIGNATURE-----<br>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br>
</div></div></blockquote></div><br></div></div>