<div dir="ltr">Yosef<div><br></div><div>Differences in the first iterations are often caused by differences when setting up the past timelevels of the initial conditions. Can you output the past timelevels at iteration 0 (i.e. of the initial conditions), without performing any time evolution? There should be parameters for this for both CarpetIOASCII and CarpetIODF5, whichever you prefer.</div>
<div><br></div><div style>-erik</div><div style><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, May 23, 2013 at 8:13 AM, Yosef Zlochower <span dir="ltr"><<a href="mailto:yosef@astro.rit.edu" target="_blank">yosef@astro.rit.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
When I tried to run the testsuites on the RIT thorns, the<br>
testsuite failed due to an issue at iteration 1, that then<br>
disappears. For certain analysis quantities like the Hamiltonian<br>
constraint violation and psi4, I find differences in every other<br>
bufferzone point between Orsted and Gauss. However, this only occurs<br>
at iteration 1 ( my test ran for 16 iterations)<br>
<br>
My test uses two levels (fixed) of refinement and a relatively<br>
coarse grid. On the other hand, I don't see this at all in the<br>
adm metric. Where there changes to the prolongation operators that<br>
may explain this?<br>
<br>
<br>
My test used 6 buffer zones (there were 4 ghostzones).<br>
Here are the relevant Carpet parameters<br>
<br>
Carpet::domain_from_coordbase = yes<br>
Carpet::prolongation_order_space = 5<br>
Carpet::prolongation_order_time = 2<br>
Carpet::max_refinement_levels = 2<br>
Carpet::use_buffer_zones = yes<br>
Carpet::num_integrator_substeps = 1<br>
Carpet::additional_buffer_zones = 2<br>
Carpet::verbose = no<br>
Carpet::schedule_barriers = no<br>
Carpet::veryverbose = no<br>
Carpet::init_fill_timelevels = "yes"<br>
<br>
ActiveThorns="CarpetRegrid2"<br>
CarpetRegrid2::regrid_every = 16<br>
CarpetRegrid2::num_centres = 1<br>
CarpetRegrid2::num_levels_1 = 2<br>
CarpetRegrid2::radius_1[ 1] = 5<br>
--<br>
Dr. Yosef Zlochower<br>
Center for Computational Relativity and Gravitation<br>
Assistant Professor<br>
School of Mathematical Sciences<br>
Rochester Institute of Technology<br>
85 Lomb Memorial Drive<br>
Rochester, NY 14623<br>
<br>
Office:74-2067<br>
Phone: +1 585-475-6103<br>
<br>
<a href="mailto:yosef@astro.rit.edu">yosef@astro.rit.edu</a><br>
<br>
CONFIDENTIALITY NOTE: The information transmitted, including<br>
attachments, is intended only for the person(s) or entity to which it<br>
is addressed and may contain confidential and/or privileged material.<br>
Any review, retransmission, dissemination or other use of, or taking<br>
of any action in reliance upon this information by persons or entities<br>
other than the intended recipient is prohibited. If you received this<br>
in error, please contact the sender and destroy any copies of this<br>
information.<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>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><br><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a>
</div>