[Users] PeriodicCarpet question

[Erik Schnetter]

Yosef

(1) Yes, this needs to be 2*d.
(2) Yes, the other directions should not be checked.
(3) The values of stencil[] are actually ignored; the code later
actually uses boundary_width from the grid structure (but without a
consistency check).

-erik


On Mon, Jul 9, 2018 at 10:06 AM, Yosef Zlochower <yosef at astro.rit.edu> wrote:
> For periodiccarpet to work, is it necessary to demand that the upper and
> lower boundary widths are the same for boundaries not filled in via
> periodic BCs?
>
> This is enforced in
>   periodic.cc on lines 493-496
>
>    for (int d = 0; d < dim; ++d) {
>      assert(width[d] == width[d + 1]);
>      stencil[d] = width[2 * d];
>    }
>
> Also, I think the test should be assert(width[2*d] == width[2*d + 1]);.
>
>
>
> --
> Dr. Yosef Zlochower
> Center for Computational Relativity and Gravitation
> Associate 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/