Astrophysics Source Code Library

Making codes discoverable since 1999

ASCL Code Record

[ascl:1212.005] General complex polynomial root solver

This general complex polynomial root solver, implemented in Fortran and further optimized for binary microlenses, uses a new algorithm to solve polynomial equations and is 1.6-3 times faster than the ZROOTS subroutine that is commercially available from Numerical Recipes, depending on application. The largest improvement, when compared to naive solvers, comes from a fail-safe procedure that permits skipping the majority of the calculations in the great majority of cases, without risking catastrophic failure in the few cases that these are actually required.

Code site:
Used in:
Described in:

Views: 1446

Add this shield to your page
Copy the above HTML to add this shield to your code's website.