RZ PMLs Called Even if not used

[ax3l]

Discussed in #5665

^{Originally posted by chiehjen February 13, 2025}
Hi:

I tried to do a simulation with mesh refinement with RZ coordinate. However, it shows the following error:

amrex::Abort::0::
ERROR   : PML are not implemented in cylindrical geometry.
            (/home/warpx/src/warpx/Source/FieldSolver/FiniteDifferenceSolver/EvolveBPML.cpp:53)
!!!
SIGABRT

The above happens despite that I specified (full input file attached rz.pdf):

boundary.field_lo = none absorbing_silver_mueller
boundary.field_hi = absorbing_silver_mueller absorbing_silver_mueller

Is it possible to do mesh refinement for RZ coordinate? Thank you in advance.

[ax3l]

@EZoni with a quick look, this could be a regression to our work on the MultiFabRegister #5356 #5334 #5230 ...?

Somehow EvolveBPML gets called even though PMLs are not user-selected.

[EZoni]

I have to look into this in detail but the fact that PMLs are used would be consistent with the fact that the simulation uses MR. MR patches are typically terminated by absorbing layers (PMLs) to prevent the reflection of electromagnetic waves. More details on this can be found in our documentation at https://warpx.readthedocs.io/en/latest/theory/amr.html.

[EZoni]

@chiehjen

To start looking into this, could you share the input file in a simple text format, i.e., not PDF?