#1631: The dgfe branch for McLachlan generates NaNs -----------------------------------------------------------+---------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: new Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Keywords: dgfe,mclachlan,gamma driver | -----------------------------------------------------------+---------------- When one uses the gamma driver formulation in the dgfe branch of McLachlan and sets the shiftGammaCoefficient to zero, the runtime output has NaNs because the code divides by zero.
A bit of flow control fixes this. Attached is a patch for the branch.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: review Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+--------------- Changes (by eschnett):
* status: new => review
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+--------------- Changes (by knarf):
* status: review => reviewed_ok
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+---------------
Comment (by eschnett):
The same problem also exists on the rewrite branch, and the same solution should be applied there.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+---------------
Comment (by Jonah Miller <jonah.maxwell.miller@…>):
If no one else has started working on this, I can go through the rewrite branch and generate another patch. Just let me know.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+---------------
Comment (by diener):
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.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+---------------
Comment (by Jonah Miller <jonah.maxwell.miller@…>):
Okay, thanks. I'll keep an eye out for the commit to the main branch.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+---------------
Comment (by diener):
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.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: ET_2014_11 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+---------------
Comment (by hinder):
Please apply.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+--------------- Changes (by knarf):
* milestone: ET_2014_11 =>
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reviewed_ok Priority: optional | Milestone: Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+---------------
Comment (by rhaas):
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.
#1631: The dgfe branch for McLachlan generates NaNs ------------------------------------------------------------+--------------- Reporter: Jonah Miller <jonah.maxwell.miller@…> | Owner: eschnett Type: defect | Status: reopened Priority: optional | Milestone: Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: dgfe,mclachlan,gamma driver ------------------------------------------------------------+--------------- Changes (by rhaas):
* status: reviewed_ok => reopened
trac@lists.einsteintoolkit.org