For a project, I need some routines to handle over/under determined matrices with.
I have a project that I started out in PDL, and I want to extend.
Unfortunately, PDL is very limited in this respect, and only gives the simq function.
In other programs (written in C), I have used the lapack routines for this task.
I have no experience in writing interface routines from PDL to external functions.
Has anybody made a lapack interface from PDL ?