[Users] PeriodicCarpet question
Yosef Zlochower
yosef at astro.rit.edu
Mon Jul 9 09:06:17 CDT 2018
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.
More information about the Users
mailing list