Skip to content
  • Liam Healy's avatar
    LU decomposition ported to ffa · aa4aa130
    Liam Healy authored
    LU decomposition functions now support real and complex types, and the
    example has been tested and works:
    (invert-matrix
    	     (make-array* '(2 2) 'double-float :initial-contents '(1.0d0 2.0d0 3.0d0 4.0d0)))
    #2A((-1.9999999999999998d0 1.0d0)
        (1.4999999999999998d0 -0.49999999999999994d0))
    aa4aa130