Thanks. I did this and found that the past timelevels, under Orsted, where not initialized correctly.
At one particular point, the Hamiltonian is
On Orsted -0.8 0 0 -4.77853228692383 t=0 -0.8 0 0 -2.33210262856286 t=-1 -0.8 0 0 -2.33210262856286 t=-2
On Gauss
-0.8 0 0 -4.77853228692383 t=0 -0.8 0 0 -4.77853228692383 t=-1 -0.8 0 0 -4.77853228692383 t=-2
I was using init_fill_timelevels.
On 05/23/2013 11:07 AM, Erik Schnetter wrote:
Yosef
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.
-erik
On Thu, May 23, 2013 at 8:13 AM, Yosef Zlochower <yosef@astro.rit.edu mailto:yosef@astro.rit.edu> wrote:
Hi, When I tried to run the testsuites on the RIT thorns, the testsuite failed due to an issue at iteration 1, that then disappears. For certain analysis quantities like the Hamiltonian constraint violation and psi4, I find differences in every other bufferzone point between Orsted and Gauss. However, this only occurs at iteration 1 ( my test ran for 16 iterations) My test uses two levels (fixed) of refinement and a relatively coarse grid. On the other hand, I don't see this at all in the adm metric. Where there changes to the prolongation operators that may explain this? My test used 6 buffer zones (there were 4 ghostzones). Here are the relevant Carpet parameters Carpet::domain_from_coordbase = yes Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2 Carpet::max_refinement_levels = 2 Carpet::use_buffer_zones = yes Carpet::num_integrator_substeps = 1 Carpet::additional_buffer_zones = 2 Carpet::verbose = no Carpet::schedule_barriers = no Carpet::veryverbose = no Carpet::init_fill_timelevels = "yes" ActiveThorns="CarpetRegrid2" CarpetRegrid2::regrid_every = 16 CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 2 CarpetRegrid2::radius_1[ 1] = 5 -- Dr. Yosef Zlochower Center for Computational Relativity and Gravitation Assistant Professor School of Mathematical Sciences Rochester Institute of Technology 85 Lomb Memorial Drive Rochester, NY 14623 Office:74-2067 Phone: +1 585-475-6103 yosef@astro.rit.edu <mailto:yosef@astro.rit.edu> CONFIDENTIALITY NOTE: The information transmitted, including attachments, is intended only for the person(s) or entity to which it is addressed and may contain confidential and/or privileged material. Any review, retransmission, dissemination or other use of, or taking of any action in reliance upon this information by persons or entities other than the intended recipient is prohibited. If you received this in error, please contact the sender and destroy any copies of this information. _______________________________________________ Users mailing list Users@einsteintoolkit.org <mailto:Users@einsteintoolkit.org> http://lists.einsteintoolkit.org/mailman/listinfo/users-- Erik Schnetter <schnetter@cct.lsu.edu mailto:schnetter@cct.lsu.edu> http://www.perimeterinstitute.ca/personal/eschnetter/