comments

jacobwilliams · jacobwilliams · commit f7c15372814f · 2026-05-24T09:16:55.000-05:00
diff --git a/src/jacchia_roberts_module.f90 b/src/jacchia_roberts_module.f90
@@ -598,7 +598,7 @@ end subroutine find_cstar_roots
 !---------------------------------------------------------------------------
 !>
 !  Original version of root-finding routine using the `roots` subroutine.
-!  **This method has severe numerical instability issues and should not be used.**
+!  This is based on the routine from GMAT.
 !
 !### Problems with this method:
 !
@@ -638,7 +638,6 @@ end subroutine find_cstar_roots
 !  Achieves ~1e-8 accuracy for all roots across all physical temperatures.
 !
 !@note Retained here for reference, testing, and demonstration of the bug.
-!      DO NOT USE in production code.
 
    pure subroutine find_cstar_roots_original(c_star, tx, root1, root2, x_root, y_root)
 

Original file line number	Diff line number	Diff line change
`@@ -598,7 +598,7 @@ end subroutine find_cstar_roots`
`598`	`598`	`!---------------------------------------------------------------------------`
`599`	`599`	`!>`
`600`	`600`	! Original version of root-finding routine using the `roots` subroutine.
`601`		`-! This method has severe numerical instability issues and should not be used.`
	`601`	`+! This is based on the routine from GMAT.`
`602`	`602`	`!`
`603`	`603`	`!### Problems with this method:`
`604`	`604`	`!`
`@@ -638,7 +638,6 @@ end subroutine find_cstar_roots`
`638`	`638`	`! Achieves ~1e-8 accuracy for all roots across all physical temperatures.`
`639`	`639`	`!`
`640`	`640`	`!@note Retained here for reference, testing, and demonstration of the bug.`
`641`		`-! DO NOT USE in production code.`
`642`	`641`
`643`	`642`	`pure subroutine find_cstar_roots_original(c_star, tx, root1, root2, x_root, y_root)`
`644`	`643`