The dgfe branch for McLachlan generates NaNs

Comments (14)

Erik Schnetter
- changed status to open
- removed comment
- 2014-06-15T16:50:28+00:00
Frank Löffler
- changed status to open
- removed comment
- 2014-06-16T09:34:51+00:00
Erik Schnetter
- removed comment
The same problem also exists on the rewrite branch, and the same solution should be applied there.
- 2014-06-16T17:02:43+00:00
Jonah Miller reporter
- removed comment
If no one else has started working on this, I can go through the rewrite branch and generate another patch. Just let me know.
- 2014-06-16T17:20:08+00:00
Peter Diener
- removed comment
I'm currently working on the compatibility of the rewrite branch with the main branch. In order to reduce the risk of conflicts, please hold off until that is completed.
- 2014-06-16T17:30:03+00:00
Jonah Miller reporter
- removed comment
Okay, thanks. I'll keep an eye out for the commit to the main branch.
- 2014-06-16T17:57:33+00:00
Peter Diener
- removed comment
This patch doesn't work when applied to the rewrite branch. This is due to vector instructions being used and so the variable shiftGammaCoeffValue is declared as a CCTK_REAL_VEC. Comparing this to a scalar value does not work. The easy fix is to replace shiftGammaCoeffValue with the scalar parameter shiftGammaCoeff. The expression for shiftGammaCoeffValue guarantees that it is zero if shiftGammaCoeff is.
- 2014-06-18T13:04:44+00:00
Ian Hinder
- removed comment
Please apply.
- 2014-10-20T10:51:40+00:00
Frank Löffler
- removed milestone
- removed comment
- 2014-11-10T09:13:48+00:00
Roland Haas
- removed comment
Has there been progress on this? It is still set to "reviewed" and unless the patch actually applies I believe it should be reset to "confirmed" until a new patch is provided.
- 2016-04-11T07:40:14+00:00
Roland Haas
- changed status to open
- removed comment
- 2016-04-17T06:36:26+00:00
Roland Haas
Please review.
- 2021-01-08T15:16:33+00:00
Roland Haas
Looking at the patch, I do not think this ought to be applied. The issue only happens if one sets evolveB=1 ie when one has chosen to define:
```
\dot \beta^i = shiftGammaCoeff \alpha^shiftAlphaPower B^i
```
so setting shiftGammaCoeff = 0 is setting \dot \beta^i to 0 which really should not be done this way. Instead shiftGammaCoeff should only be allowed strictly positive values (unless there is some use for negative ones).

Unfortunately, the current default value for shiftGammaCoeff is 0.
- 2021-02-09T16:53:05+00:00
Roland Haas
- edited description
- changed status to wontfix
- 2021-02-09T17:21:29+00:00
Log in to comment