1 #include "gmx_lapack.h"
2 #include "lapack_limits.h"
5 F77_FUNC(sorgbr,SORGBR)(const char *vect,
16 int wantq,iinfo,j,i,i1,wrksz;
17 int mn = (*m < *n) ? *m : *n;
19 wantq = (*vect=='Q' || *vect=='q');
22 wrksz = mn*DORGBR_BLOCKSIZE;
33 F77_FUNC(sorgqr,SORGQR)(m,n,k,a,lda,tau,work,lwork,&iinfo);
36 a[(j-1)*(*lda)+0] = 0.0;
38 a[(j-1)*(*lda)+(i-1)] = a[(j-2)*(*lda)+(i-1)];
45 F77_FUNC(sorgqr,SORGQR)(&i1,&i1,&i1,&(a[*lda+1]),lda,tau,work,lwork,&iinfo);
50 F77_FUNC(sorglq,SORGLQ)(m,n,k,a,lda,tau,work,lwork,&iinfo);
57 a[(j-1)*(*lda)+(i-1)] = a[(j-1)*(*lda)+(i-2)];
58 a[(j-1)*(*lda)+0] = 0.0;
62 F77_FUNC(sorglq,SORGLQ)(&i1,&i1,&i1,&(a[*lda+1]),lda,tau,work,lwork,&iinfo);