I need to calculate pseudo inverse matrix.
How can I figure out which function to use?
I googled it and it looks like other people use sgelss
2. compute pseudo inverse
LDWORK = WORK;
delete  WORK;
WORK = new float[LDWORK];
res = sgelss_( &m, &n, &NRHS, &matrix, &LDA, B, &LDB, SING, &RCOND, &IRANK, WORK, &LDWORK, &info);
But when I looked at sgelss documentation it says
SGELSS computes the minimum norm solution to a real linear least squares problem