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

[Yosef Zlochower]

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 at astro.rit.edu
> <mailto: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 <mailto: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 <mailto:Users at einsteintoolkit.org>
>     http://lists.einsteintoolkit.org/mailman/listinfo/users
>
>
>
>
> --
> Erik Schnetter <schnetter at cct.lsu.edu <mailto:schnetter at cct.lsu.edu>>
> http://www.perimeterinstitute.ca/personal/eschnetter/


-- 
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.