Skip to content

Let grdseamount -Co handle u = -1#6476

Merged
PaulWessel merged 2 commits intomasterfrom
nan-check
Mar 23, 2022
Merged

Let grdseamount -Co handle u = -1#6476
PaulWessel merged 2 commits intomasterfrom
nan-check

Conversation

@PaulWessel
Copy link
Member

For u = -1 in the polynomial shape we end up with a 0/0 expression which evaluates to NaN. However, we know from l'Hopital's rule that the result should be 0. This PR intervenes for -1 to yield 0. While at it, since this function and its gradient are only defined over ±1 I check r and set results to 0 if outside.

For u = -1 in the polynomial shape we end up with a 0/0 expression which evaluates to NaN.  However, we know from l'Hopital's rule that the result should beo 0.  This PR intervenes for -1 to yield 0.  WHile at it, since this function and its gradietn are only defined over ±1 I check r and set results to 0 if outside.
@PaulWessel PaulWessel added the maintenance Boring but important stuff for the core devs label Mar 23, 2022
@PaulWessel PaulWessel added this to the 6.4.0 milestone Mar 23, 2022
@PaulWessel PaulWessel requested a review from maxrjones March 23, 2022 15:45
@PaulWessel PaulWessel self-assigned this Mar 23, 2022
@PaulWessel PaulWessel merged commit c46917d into master Mar 23, 2022
@PaulWessel PaulWessel deleted the nan-check branch March 23, 2022 15:57
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

maintenance Boring but important stuff for the core devs

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants