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