Skip to content

Conversation

@reinhold-willcox
Copy link
Collaborator

Changed calls to PPOW(base, exponent) where exponent was something like 1/2 or 1/3 to sqrt and cbrt. In the current functionality, the 1/3 power was unable to handle negative arguments, even though mathematically this is fine. Cbrt has no issue with negative arguments.

@reinhold-willcox reinhold-willcox linked an issue Dec 16, 2021 that may be closed by this pull request
Copy link
Collaborator

@jeffriley jeffriley left a comment

Choose a reason for hiding this comment

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

Looks good. For consistency we should probably used 'std::cbrt' and 'std::sqrt' always (or 'cbrt' and 'sqrt' always), but that's a small point.

Changelog?

EDIT: I just noticed you missed a semi-colon in WhiteDwarfs.cpp, and the compile fails.

And our minimum test is the detailed plot in the examples. Can you run that and check your output against the standard?

@reinhold-willcox
Copy link
Collaborator Author

Woops, that's what I get for trying to code before bed. Yup, on it.

Copy link
Collaborator

@jeffriley jeffriley left a comment

Choose a reason for hiding this comment

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

Did you produce the detailed plot? I'll approve and merge, but I think you should produce the plot just to be sure.

@jeffriley jeffriley merged commit 9bed21c into TeamCOMPAS:dev Dec 16, 2021
@reinhold-willcox
Copy link
Collaborator Author

Figure_1

Looks good to me!

@reinhold-willcox reinhold-willcox deleted the nan_radius_fix branch December 20, 2021 05:02
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.

NaN radius

2 participants