Skip to content

Nodal Vector Poisson solver with EB#577

Open
ax3l wants to merge 22 commits into
AMReX-Codes:developmentfrom
ax3l:nodal_fd_lap-axel
Open

Nodal Vector Poisson solver with EB#577
ax3l wants to merge 22 commits into
AMReX-Codes:developmentfrom
ax3l:nodal_fd_lap-axel

Conversation

@ax3l

@ax3l ax3l commented Jun 3, 2026

Copy link
Copy Markdown
Member

Continuation of #548 where I have permissions to push

Prabhat Kumar and others added 17 commits November 6, 2025 13:30
Signed-off-by: S. Eric Clark <25495882+clarkse@users.noreply.github.com>
Signed-off-by: S. Eric Clark <25495882+clarkse@users.noreply.github.com>
Signed-off-by: S. Eric Clark <245461744+clarkse-he@users.noreply.github.com>
Signed-off-by: S. Eric Clark <245461744+clarkse-he@users.noreply.github.com>
…nitizlize EB.

Signed-off-by: S. Eric Clark <245461744+clarkse-he@users.noreply.github.com>
Signed-off-by: S. Eric Clark <245461744+clarkse-he@users.noreply.github.com>
@ax3l ax3l added the enhancement New feature or request label Jun 3, 2026
Comment thread src/VectorPoisson/CMakeLists.txt Outdated
Comment thread src/VectorPoisson/CMakeLists.txt Outdated
Comment thread src/VectorPoisson/VectorPoissonSolverNodal.H Outdated
Comment thread src/VectorPoisson/PyVectorPoisson.H Outdated
Comment thread src/VectorPoisson/VectorPoissonSolver.H Outdated
Comment on lines +1 to +5
#!/usr/bin/env python3
"""Test VectorPoissonSolverNodal in 2D RZ — A_theta convergence study.

Manufactured solution: A_theta = r * sin(pi*r) * sin(pi*z) on [0,1]^2.

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Would be nice if this was integrated in our pytest suite. Please add a def test_... entry point function.

Comment on lines +1 to +5
#!/usr/bin/env python3

"""
Self-convergence test: VectorPoissonSolverNodal with EB in 2D RZ.
=================================================================

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Would be nice if this was integrated in our pytest suite. Please add a def test_... entry point function.

@ax3l ax3l force-pushed the nodal_fd_lap-axel branch from d02c4da to cab7e8c Compare June 3, 2026 06:46
@ax3l ax3l force-pushed the nodal_fd_lap-axel branch from cab7e8c to 1e2e304 Compare June 3, 2026 06:47
amrex::Real inv_dr;
amrex::Real inv_dz;

AMREX_GPU_DEVICE

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This probably needs host-device annotation

const Box& bx = mfi.validbox();
Array4<Real> const& a = acoef.array(mfi);

ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k) {

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This probably needs host-device annotation

Comment on lines +27 to +29
# RZ is a 2D concept; this test requires the 2D module. It is skipped when the
# 2D module cannot be loaded (e.g. the test suite already registered the 3D one).
amr = pytest.importorskip("amrex.space2d", exc_type=ImportError)

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Seems misplaced:

____________ ERROR collecting tests/test_nodal_vector_poisson_rz.py ____________
tests/test_nodal_vector_poisson_rz.py:29: in <module>
    amr = pytest.importorskip("amrex.space2d", exc_type=ImportError)
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

@ax3l ax3l changed the title Nodal Vector Poisson solver with EB (Axel) Nodal Vector Poisson solver with EB Jun 6, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

4 participants