Started by an SCM change000130609514904701831680.12181148114863655285cd4ce9d65cec1b44cf2e6e578f4d08fe73655285cd4ce9d65cec1b44cf2e6e578f4d08fe7origin/master3655285cd4ce9d65cec1b44cf2e6e578f4d08fe73655285cd4ce9d65cec1b44cf2e6e578f4d08fe7origin/master88d63a0f595d93dccafcc9e7917c00fe97fd034a35d63a0f595d93dccafcc9e7917c00fe97fd034a35origin/grid_fixesd63a0f595d93dccafcc9e7917c00fe97fd034a35d63a0f595d93dccafcc9e7917c00fe97fd034a35origin/grid_fixes79c5c01fc5da140dc13df06c2919c4f4f5a0318b5fc5c01fc5da140dc13df06c2919c4f4f5a0318b5forigin/changes_to_not_breakc5c01fc5da140dc13df06c2919c4f4f5a0318b5fc5c01fc5da140dc13df06c2919c4f4f5a0318b5forigin/changes_to_not_break701344f8a40d6116af1b567307eb389b4be1009fc41344f8a40d6116af1b567307eb389b4be1009fc4refs/remotes/origin/changes_to_not_break1344f8a40d6116af1b567307eb389b4be1009fc41344f8a40d6116af1b567307eb389b4be1009fc4refs/remotes/origin/changes_to_not_break851d226e398bae00fe10255274a70f3451b5cc7d9e1d226e398bae00fe10255274a70f3451b5cc7d9eorigin/wangara_updates1d226e398bae00fe10255274a70f3451b5cc7d9e1d226e398bae00fe10255274a70f3451b5cc7d9eorigin/wangara_updates520c3bb3c1b174aaaa1c56e4e70aafb6d658b55f0acc3bb3c1b174aaaa1c56e4e70aafb6d658b55f0acrefs/remotes/origin/masterc3bb3c1b174aaaa1c56e4e70aafb6d658b55f0acc3bb3c1b174aaaa1c56e4e70aafb6d658b55f0acrefs/remotes/origin/masterc3bb3c1b174aaaa1c56e4e70aafb6d658b55f0acc3bb3c1b174aaaa1c56e4e70aafb6d658b55f0acrefs/remotes/origin/masterhttps://github.com/larson-group/sam_clubb.gitfalse#52014904701940859sam_clubb_nvhpc_build #520520false52010559SUCCESS1791967509404https://carson.math.uwm.edu/jenkins/job/sam_clubb_nvhpc_build/520/calc_roots.F9076464de9f694cb3148b0ae07f9eb0098afc4a4e41791911400000https://carson.math.uwm.edu/jenkins/user/noreplynoreplynoreply@github.comFixed bug in cubic equation solve function. Only relevant if pdf is L… (#1387) * Fixed bug in cubic equation solve function. Only relevant if pdf is Luhar or iiPDF_type=4 Therefore technically BIT_CHANGING... The bug existed because the old code used the principal cube root which has an imaginary part for negative bases. The fix is to simply calc: -|w|**(1/3) for negative real numbers. * Replaced determinant with discriminant Verified via a new Fortran unit test and also manually checked with wolfram alpha. There is no way to actually run this code :/ would only work for iiPDF_types: iiPDF_3D_Luhar = 3 and iiPDF_3D_Luhar = 3 (and those are currently (Oct 2026) broken ) 2026-10-13 12:10:00 -050076464de9f694cb3148b0ae07f9eb0098afc4a4e4Fixed bug in cubic equation solve function. Only relevant if pdf is L… (#1387)editcalc_roots.F90githttps://carson.math.uwm.edu/jenkins/user/noreplynoreplyfalse521https://carson.math.uwm.edu/jenkins/job/sam_clubb_nvhpc_build/521/519https://carson.math.uwm.edu/jenkins/job/sam_clubb_nvhpc_build/519/