Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Convert sphe to C++ #2579

Merged
merged 2 commits into from
Sep 24, 2023
Merged

Convert sphe to C++ #2579

merged 2 commits into from
Sep 24, 2023

Conversation

maxpkatz
Copy link
Member

PR summary

PR checklist

  • test suite needs to be run on this PR
  • this PR will change answers in the test suite to more than roundoff level
  • all newly-added functions have docstrings as per the coding conventions
  • the CHANGES file has been updated, if appropriate
  • if appropriate, this change is described in the docs

@maxpkatz maxpkatz marked this pull request as ready for review September 24, 2023 10:31
@zingale
Copy link
Member

zingale commented Sep 24, 2023

@@ -1160,8 +1169,17 @@ void RadSolve::getEdgeMetric(int idim, const Geometry& geom, const Box& edgebox,
geom.GetEdgeLoc(s, reg, I);
}
const Real *dx = geom.CellSize();
sphe(r.dataPtr(), s.dataPtr(), idim,
ARLIM(edgebox.loVect()), ARLIM(edgebox.hiVect()), dx);

Copy link
Member

Choose a reason for hiding this comment

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

aren't we missing the 1-d case here, r = r**2?

Copy link
Member

Choose a reason for hiding this comment

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

and the n = 0 case for multid?

Copy link
Member Author

Choose a reason for hiding this comment

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

fixed, I think

@zingale zingale merged commit df7ce46 into AMReX-Astro:development Sep 24, 2023
18 checks passed
@maxpkatz maxpkatz deleted the cxx-sphe branch September 24, 2023 14:17
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants