Work at SourceForge, help us to make it a better place! We have an immediate need for a Support Technician in our San Francisco or Denver office.

Close

#60 implement FZ_ROOTS (imsl_zeropoly is available)

open
nobody
None
5
2010-02-12
2009-09-11
David Nicholls
No

% Function not found: FZ_ROOTS

Discussion

  • Hi David,

    If your polynomials do not have complex coefficients, and if it's OK to use a different algorithm, please try using the IMSL_ZEROPOLY() from the current CVS version. It has the same calling sequence as FZ_ROOTS(). The following wrapper might enable your code to work without changes:

    function fz_roots, c, double=d
    return, imsl_zeropoly(c, double=d)
    end

    GDL implements IMSL_ZEROPOLY() using GSL's gsl_poly_complex_solve() which uses the "balanced-QR reduction of the companion matrix" (as IMSL_ZEROPOLY() does in IDL).

    Best regards,
    Sylwester

     
    • summary: implement FZ_ROOTS --> implement FZ_ROOTS (imsl_zeropoly is available)