kipar

Yes, LAPACK can return eigenvectors either in complex form or as a real matrix - if you check vecs, you can see that imag part is still returned, but on...

After trying to tackle with this. 1. `complex: true\false` approach would look best, but that doesn't work because result type doesn't depend from argument value. So it will still return...

`eigs(complex:true)` from a library point of view is `eigs(complex: Bool?)`. So it's return type is `Matrix(T) | Matrix(Complex)`, it won't change depending on what value actually was passed by a...

1. With current design you can do ```crystal require "linalg" m = LA::GMat[[3, -2],[4, -1]] vals, vecs = m.eigs res = vecs.inv * m * vecs res.detect(MatrixFlags::UpperTriangular | MatrixFlags::LowerTriangular) ```...

Checking and detecting won't work in current state, because only "positive" result is cached. That is: if `detect` sees that matrix is diagonal, it sets the "diagonal" flag (and further...

-5 in lapack means "5-th argument has the wrong value". 5-th argument of sgesdd is LDA (= ncolumns). I'm not sure what this can mean, will try to investigate. I...

PHP driver is basically a wrapper around a C SDK: http://maxdb.sap.com/documentation/sqldbc/SQLDBC_API and there is even an .h file available without registration - http://maxdb.sap.com/documentation/sqldbc/SQLDBC_API/SQLDBC__C_8h-source.html

It would be great imho, as Crystal's fast speed and agile syntax could make it pretty successful in scientific programming. I see two main problems\questions though: - Scipy is a...

Good to hear. I haven't benchmarked too, but type inference was my most concern. The compiler still have to parse all code - i've started creating wrapper for LAPACK and...

Ok, now i have another question. License. Is GPL acceptable for most potential users of "Crystal SciPy"? If GPL is acceptable, wrapper around [GSL](http://www.gnu.org/s/gsl/) can provide most of the features...