[Users] Different result in buffer zone at iteration 1 only on Gauss

[Erik Schnetter]

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 at 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 at 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 at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users
>



-- 
Erik Schnetter <schnetter at cct.lsu.edu>
http://www.perimeterinstitute.ca/personal/eschnetter/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20130523/d5c6d48a/attachment.html