2 #include "types/simple.h"
4 #include "../gmx_blas.h"
5 #include "../gmx_lapack.h"
9 F77_FUNC(sgetf2,SGETF2)(int *m,
26 k = (*m < *n) ? *m : *n;
30 jp = j - 1 + F77_FUNC(isamax,ISAMAX)(&t1,&(a[(j-1)*(*lda)+(j-1)]),&t2);
32 if( fabs(a[(j-1)*(*lda)+(jp-1)])>GMX_FLOAT_MIN ) {
34 F77_FUNC(sswap,SSWAP)(n,&(a[ j-1 ]),lda,&(a[ jp-1 ]),lda);
39 tmp = 1.0/a[(j-1)*(*lda)+(j-1)];
40 F77_FUNC(sscal,SSCAL)(&t1,&tmp,&(a[(j-1)*(*lda)+(j)]),&t2);
50 F77_FUNC(sger,SGER)(&t1,&t2,&minusone,&(a[(j-1)*(*lda)+(j)]),&t3,
51 &(a[(j)*(*lda)+(j-1)]),lda, &(a[(j)*(*lda)+(j)]),lda);