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 <[email protected] > <mailto:[email protected]>> 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 > > [email protected] <mailto:[email protected]> > > 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 > [email protected] <mailto:[email protected]> > http://lists.einsteintoolkit.org/mailman/listinfo/users > > > > > -- > Erik Schnetter <[email protected] <mailto:[email protected]>> > 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 [email protected] 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 [email protected] http://lists.einsteintoolkit.org/mailman/listinfo/users
