lapack.cpp File Reference

#include <ctype.h>
#include <math.h>
#include "blas/blas.h"
#include "lapack.h"
#include "lapack_limits.h"
#include "simple.h"

Namespaces
	PLMD
	PLMD::lapack

Macros
#define	DORGLQ_BLOCKSIZE   32
#define	DORGLQ_MINBLOCKSIZE   2
#define	DORGLQ_CROSSOVER   128
#define	SORGLQ_BLOCKSIZE   32
#define	SORGLQ_MINBLOCKSIZE   2
#define	SORGLQ_CROSSOVER   128

Functions
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dbdsdc, DBDSDC)(const char *uplo
	PLMD::lapack::if (icompq==1)
	PLMD::lapack::if (iuplo==2)
	PLMD::lapack::if (fabs(orgnrm)< PLUMED_GMX_DOUBLE_MIN) = -1.0
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlascl, DLASCL)("G"
	PLMD::lapack::for (i__=1;i__<=i__1;++i__)
	PLMD::lapack::for (ii=2;ii<=i__1;++ii)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dbdsqr, DBDSQR)(const char *uplo
	PLMD::lapack::if ((xuplo!='U') &&!lower)
else	PLMD::lapack::if (((*ncvt==0) &&(*ldvt< 1))||((*ncvt > 0) &&(*ldvt< itmp1)))
else	PLMD::lapack::if (((*ncc==0) &&(*ldc< 1))||((*ncc > 0) &&(*ldc< itmp1)))
	PLMD::lapack::if (! rotate)
	PLMD::lapack::if (lower)
	PLMD::lapack::if (tol >=0.f)
	PLMD::lapack::if (iter > maxit)
	PLMD::lapack::for (lll=1;lll<=i__1;++lll)
	PLMD::lapack::if (ll==m - 1)
	PLMD::lapack::if (idir==1)
	PLMD::lapack::if (tol >=0.f &&*n *tol *(sminl/smax)<=((r__1 >r__2) ? r__1 :r__2))
	PLMD::lapack::if (shift==0.f)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgebd2, DGEBD2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgebrd, DGEBRD)(int *m
	PLMD::lapack::if (minmn==0)
	PLMD::lapack::if (nb > 1 &&nb< minmn)
	PLMD::lapack::for (i_=1;i_2< 0 ? i_ >=i_1 :i_<=i_1;i_+=i_2)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgebd2, DGEBD2)(&i_2
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgelq2, DGELQ2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgelqf, DGELQF)(int *m
	PLMD::lapack::if (k==0)
	PLMD::lapack::if (i__<=k)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgeqr2, DGEQR2)(int *m
	PLMD::lapack::for (i=0;i< k;i++) = *beta
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgeqrf, DGEQRF)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgesdd, DGESDD)(const char *jobz
	PLMD::lapack::if (lquery !=0)
	PLMD::lapack::if (anrm > 0. &&anrm< smlnum)
else	PLMD::lapack::if (anrm > bignum)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgetf2, DGETF2)(int *m
	PLMD::lapack::for (j=1;j<=k;j++)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgetrf, DGETRF)(int *m
	PLMD::lapack::if (DGETRF_BLOCKSIZE >=mindim)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgetri, DGETRI)(int *n
else if *	PLMD::lapack::lda (*n)
else if *	PLMD::lapack::lwork (*n) &&*lwork!=-1
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dtrtri, DTRTRI)("Upper"
	PLMD::lapack::if (nb< nbmin||nb >= *n)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgetrs, DGETRS)(const char *trans
	PLMD::lapack::if (notran)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dtrsm, DTRSM)("Left"
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlaswp, DLASWP)(nrhs
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlabrd, DLABRD)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlacpy, DLACPY)(const char *uplo
	PLMD::lapack::if (ch=='U')
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlae2, DLAE2)(double *a
	PLMD::lapack::if (fabs(*a) > fabs(*c__))
	PLMD::lapack::if (adf > ab)
else	PLMD::lapack::if (sm > 0.)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlaebz, DLAEBZ)(int *ijob
	PLMD::lapack::for (ji=1;ji<=i__1;++ji)
	PLMD::lapack::for (jit=1;jit<=i__1;++jit)
	PLMD::lapack::if (i__1 >0) *info = i__1
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlaed6, DLAED6)(int *kniter
	PLMD::lapack::if (c__==0.f)
else	PLMD::lapack::if (a<=0.f)
	PLMD::lapack::if (fabs(*finit)<=fabs(temp))
	PLMD::lapack::if (temp<=small1)
	PLMD::lapack::if (fabs(f)<=0.f)
	PLMD::lapack::for (niter=iter;niter<=20;++niter)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlaev2, DLAEV2)(double *a
	PLMD::lapack::if (df >=0.)
	PLMD::lapack::if (acs > ab)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlagtf, DLAGTF)(int *n
	PLMD::lapack::for (k=1;k<=i__1;++k)
	PLMD::lapack::if (fabs(a[*n])<=scale1 *tl &&in[*n]==0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlagts, DLAGTS)(int *job
	PLMD::lapack::if (fabs(*job) > 2|| *job==0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlamrg, DLAMRG)(int *n1
	PLMD::lapack::while (n1sv >0 &&n2sv >0)
	PLMD::lapack::if (n1sv==0)
double	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlange, DLANGE)(const char *norm
	PLMD::lapack::switch (ch)
double	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlanst, DLANST)(const char *norm
double	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlansy, DLANSY)(const char *norm
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlassq, DLASSQ)(n
double	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlapy2, DLAPY2)(double *x
	PLMD::lapack::if (xabs >yabs)
return w *	PLMD::lapack::sqrt (1.0+z *z)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlar1vx, DLAR1VX)(int *n
	PLMD::lapack::if (!(s > 0.||s< 1.))
	PLMD::lapack::if (!(tmp > 0.||tmp< 1.))
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarf, DLARF)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarfb, DLARFB)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarfg, DLARFG)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dscal, DSCAL)(&ti1
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarft, DLARFT)(const char *direct
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarnv, DLARNV)(int *idist
	PLMD::lapack::for (iv=1;iv<=i__1;iv+=64)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarrbx, DLARRBX)(int *n
	PLMD::lapack::if (iwork[(i__<< 1) - 1]==1)
	PLMD::lapack::if (i__<= *n &&iwork[(i__<< 1) - 1] !=-1)
	PLMD::lapack::for (p=1;p<=i__1;++p)
	PLMD::lapack::if (nint > 0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarrex, DLARREX)(const char *range
	PLMD::lapack::for (jblk=1;jblk<=i__1;++jblk)
	PLMD::lapack::if (irange==2) = safemn
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarrfx, DLARRFX)(int *n
	PLMD::lapack::if (!(dmax1 > 0.||dmax1< 1.))
	PLMD::lapack::if (!(dmax2 > 0.||dmax2< 1.))
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlarrvx, DLARRVX)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlaset, DLASET)("Full" = rho
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlartg, DLARTG)(double *f
	PLMD::lapack::if (scale >=safemx2)
else	PLMD::lapack::if (scale<=safemin2)
	PLMD::lapack::if (fabs(*f)>fabs(*g) &&*cs< 0.0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlaruv, DLARUV)(int *iseed
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlas2, DLAS2)(double *f
	PLMD::lapack::while (!done)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd0, DLASD0)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasdt, DLASDT)(n
	PLMD::lapack::for (lvl=nlvl;lvl >=1;--lvl)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd1, DLASD1)(int *nl
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd2, DLASD2)(nl
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd3, DLASD3)(nl
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlamrg, DLAMRG)(&n1
	PLMD::lapack::if (j > n)
	PLMD::lapack::if (fabs(z__[j])<=tol)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dcopy, DCOPY)(&m
	PLMD::lapack::if (z__[1] > 0.)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dcopy, DCOPY)(k
	PLMD::lapack::if (ctot[1] > 0)
else	PLMD::lapack::if (ctot[3] > 0)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dgemm, DGEMM)("N"
	PLMD::lapack::if (ktemp<= *ldvt2)
	PLMD::lapack::if (ktemp > 1)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd4, DLASD4)(int *n
	PLMD::lapack::if (w<=0.)
	PLMD::lapack::if (fabs(w)<=eps *erretm)
else	PLMD::lapack::if (a >=0.)
	PLMD::lapack::if (w *eta > 0.)
	PLMD::lapack::if (temp > *rho)
	PLMD::lapack::if (w > 0.)
	PLMD::lapack::if (orgati)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlaed6, DLAED6)(&niter
	PLMD::lapack::if (w *eta >=0.)
	PLMD::lapack::if (temp > sg2ub||temp< sg2lb)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd5, DLASD5)(int *i__
	PLMD::lapack::if (b > 0.)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd6, DLASD6)(int *icompq
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd7, DLASD7)(icompq
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasd8, DLASD8)(icompq
	PLMD::lapack::if (fabs(dsigma[2])<=hlftol)
	PLMD::lapack::if (m > n)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasda, DLASDA)(int *icompq
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasdq, DLASDQ)(const char *uplo
	PLMD::lapack::if (xuplo=='U')
else	PLMD::lapack::if ((*ncvt==0 &&*ldvt< 1)||(*ncvt > 0 &&*ldvt< itmp1))
else	PLMD::lapack::if ((*ncc==0 &&*ldc< 1)||(*ncc > 0 &&*ldc< itmp1))
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasq1, DLASQ1)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasq2, DLASQ2)(n
	PLMD::lapack::if (z__[2] > z__[3] *tol2)
	PLMD::lapack::if (posinf<=1.0) ieee=0
	PLMD::lapack::if (neginf >=0.0) ieee=0
	PLMD::lapack::if (fabs(negzro)>PLUMED_GMX_DOUBLE_MIN) ieee=0
	PLMD::lapack::if (neginf >=0) ieee=0
	PLMD::lapack::if (fabs(newzro-zero)>PLUMED_GMX_DOUBLE_MIN) ieee=0
	PLMD::lapack::if (posinf<=one) ieee=0
	PLMD::lapack::for (iwhila=1;iwhila<=i__1;++iwhila)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasrt, DLASRT)("D"
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasq3, DLASQ3)(int *i0
	PLMD::lapack::if (z__[nn - 5] > tol2 *(*sigma+z__[nn - 3]) &&z__[nn -(*pp<< 1) - 4] > tol2 *z__[nn - 7])
	PLMD::lapack::if (z__[nn - 5] > z__[nn - 3] *tol2)
L70	PLMD::lapack::DLASQ5 (i0, n0, &z__[1], pp, &tau, dmin__, &dmin1, &dmin2, &dn, &dn1, &dn2, ieee)
else	PLMD::lapack::if (dmin1 > 0.)
L80	PLMD::lapack::DLASQ6 (i0, n0, &z__[1], pp, dmin__, &dmin1, &dmin2, &dn, &dn1, &dn2)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasq4, DLASQ4)(int *i0
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasq5, DLASQ5)(int *i0
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasq6, DLASQ6)(int *i0
	PLMD::lapack::for (j4= *i0 *4;j4<=i__1;j4+=4)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasr, DLASR)(const char *side
	PLMD::lapack::if (dir==-1)
	PLMD::lapack::if (endd - start<=20 &&endd - start > 0)
else	PLMD::lapack::if (endd - start > 20)
	PLMD::lapack::if (stkpnt > 0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasrt2, DLASRT2)(const char *id
	PLMD::lapack::if (endd - start > 0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlasv2, DLASV2)(double *f
	PLMD::lapack::if (swap)
	PLMD::lapack::if (ga > fa)
	PLMD::lapack::if (gasmal)
	PLMD::lapack::if (pmax==1)
n32 *	PLMD::lapack::if (n32!=0) = 32
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dlatrd, DLATRD)(const char *uplo
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorg2r, DORG2R)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorgbr, DORGBR)(const char *vect
	PLMD::lapack::if (wantq)
	PLMD::lapack::for (i=2;i<=*m;i++) a[i-1]=0.0
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorglq, DORGLQ)(&i1
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorgl2, DORGL2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorglq, DORGLQ)(int *m
else	PLMD::lapack::if (lquery)
	PLMD::lapack::if (kk > 0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorgqr, DORGQR)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorm2l, DORM2L)(const char *side
	PLMD::lapack::if (left)
	PLMD::lapack::if ((left &&notran)||(! left &&! notran))
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorm2r, DORM2R)(const char *side
	PLMD::lapack::if ((left &&!notran)||(!left &&notran))
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dormbr, DORMBR)(const char *vect
	PLMD::lapack::if (applyq)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dorml2, DORML2)(const char *side
	PLMD::lapack::if ((xside=='L' &&xtrans=='N')||(xside!='L' &&xtrans!='N'))
	PLMD::lapack::if (xside=='L')
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dormlq, DORMLQ)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dormql, DORMQL)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dormqr, DORMQR)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dormtr, DORMTR)(const char *side
	PLMD::lapack::if (upper)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dstebz, DSTEBZ)(const char *range
	PLMD::lapack::if (irange<=0)
else	PLMD::lapack::if (iorder<=0)
	PLMD::lapack::if (nb<=1)
	PLMD::lapack::for (jb=1;jb<=i__1;++jb)
	PLMD::lapack::if (iorder==1 &&*nsplit > 1)
	PLMD::lapack::if (ncnvrg)
	PLMD::lapack::if (toofew)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dstegr, DSTEGR)(const char *jobz
	PLMD::lapack::if (!(wantz||(*jobz=='N'|| *jobz=='n')))
else	PLMD::lapack::if (!(alleig||valeig||indeig))
else	PLMD::lapack::if (valeig &&*n > 0 &&*vu<= *vl)
else	PLMD::lapack::if (indeig &&(*il< 1|| *il > *n))
else	PLMD::lapack::if (indeig &&(*iu< *il|| *iu > *n))
else if *	PLMD::lapack::ldz (wantz &&*ldz< *n)
	PLMD::lapack::if (tnrm > 0. &&tnrm< rmin)
else	PLMD::lapack::if (tnrm > rmax)
	PLMD::lapack::if (fabs(scale-1.0)>PLUMED_GMX_DOUBLE_EPS)
	PLMD::lapack::if (iinfo !=0)
	PLMD::lapack::if (wantz)
	PLMD::lapack::if (nsplit > 1)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dstein, DSTEIN)(int *n
else if *	PLMD::lapack::ldz (*n)
	PLMD::lapack::for (nblk=1;nblk<=i__1;++nblk)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dsteqr, DSTEQR)(const char *compz
else if *	PLMD::lapack::ldz (icompz > 0 &&*ldz<((*n >1) ? *n :1))
	PLMD::lapack::if (icompz==2)
	PLMD::lapack::if (l1 > 1)
	PLMD::lapack::if (l1<=nm1)
	PLMD::lapack::if (lend==l)
	PLMD::lapack::if (anorm > ssfmax)
	PLMD::lapack::if (lend > l)
	PLMD::lapack::if (m==l)
	PLMD::lapack::if (jtot==nmaxit)
	PLMD::lapack::if (icompz > 0)
	PLMD::lapack::if (l >=lend)
else	PLMD::lapack::if (iscale==2)
L160	PLMD::lapack::DLASRT ("I", n, &d__[1], info)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dsterf, DSTERF)(int *n
	PLMD::lapack::for (m=l1;m<=i__1;++m)
	PLMD::lapack::if (lend >=l)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dstevr, DSTEVR)(const char *jobz
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dsyevr, DSYEVR)(const char *jobz
	PLMD::lapack::if (negzro!=0) ieeeok=0
	PLMD::lapack::if (newzro!=fzero) ieeeok=0
else	PLMD::lapack::if (!(lower||(*uplo=='U'|| *uplo=='u')))
else if *	PLMD::lapack::lda ((*n >1) ? *n :1)
	PLMD::lapack::if (alleig||indeig)
else	PLMD::lapack::if (anrm > rmax)
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dsytrd, DSYTRD)(uplo
	PLMD::lapack::if (wantz &&*info==0)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dsytd2, DSYTD2)(const char *uplo
	PLMD::lapack::if (! upper &&!(*uplo=='L'|| *uplo=='l'))
else if *	PLMD::lapack::lda ((1 > *n) ? 1 :*n)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (dtrti2, DTRTI2)(const char *uplo
	PLMD::lapack::if (nounit)
	PLMD::lapack::if (nb<=1||nb >= *n)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (ilasrt2, ILASRT2)(const char *id
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sbdsdc, SBDSDC)(const char *uplo
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slascl, SLASCL)("G"
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sbdsqr, SBDSQR)(const char *uplo
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgebd2, SGEBD2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgebrd, SGEBRD)(int *m
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgebd2, SGEBD2)(&i_2
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgelq2, SGELQ2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgelqf, SGELQF)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgeqr2, SGEQR2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgeqrf, SGEQRF)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgesdd, SGESDD)(const char *jobz
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgetf2, SGETF2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgetrf, SGETRF)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgetri, SGETRI)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (strtri, STRTRI)("Upper"
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgetrs, SGETRS)(const char *trans
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (strsm, STRSM)("Left"
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slaswp, SLASWP)(nrhs
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slabrd, SLABRD)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slacpy, SLACPY)(const char *uplo
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slae2, SLAE2)(float *a
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slaebz, SLAEBZ)(int *ijob
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slaed6, SLAED6)(int *kniter
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slaev2, SLAEV2)(float *a
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slagtf, SLAGTF)(int *n
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slagts, SLAGTS)(int *job
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slamrg, SLAMRG)(int *n1
float	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slange, SLANGE)(const char *norm
float	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slanst, SLANST)(const char *norm
float	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slansy, SLANSY)(const char *norm
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slassq, SLASSQ)(n
float	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slapy2, SLAPY2)(float *x
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slar1vx, SLAR1VX)(int *n
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarf, SLARF)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarfb, SLARFB)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarfg, SLARFG)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sscal, SSCAL)(&ti1
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarft, SLARFT)(const char *direct
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarnv, SLARNV)(int *idist
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarrbx, SLARRBX)(int *n
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarrex, SLARREX)(const char *range
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarrfx, SLARRFX)(int *n
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slarrvx, SLARRVX)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slaset, SLASET)("Full" = rho
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slartg, SLARTG)(float *f
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slaruv, SLARUV)(int *iseed
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slas2, SLAS2)(float *f
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd0, SLASD0)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasdt, SLASDT)(n
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd1, SLASD1)(int *nl
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd2, SLASD2)(nl
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd3, SLASD3)(nl
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slamrg, SLAMRG)(&n1
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (scopy, SCOPY)(&m
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (scopy, SCOPY)(k
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sgemm, SGEMM)("N"
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd4, SLASD4)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slaed6, SLAED6)(&niter
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd5, SLASD5)(int *i__
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd6, SLASD6)(int *icompq
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd7, SLASD7)(icompq
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasd8, SLASD8)(icompq
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasda, SLASDA)(int *icompq
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasdq, SLASDQ)(const char *uplo
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasq1, SLASQ1)(int *n
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasq2, SLASQ2)(n
	PLMD::lapack::if (fabs(negzro)>PLUMED_GMX_FLOAT_MIN) ieee=0
	PLMD::lapack::if (fabs(newzro-zero)>PLUMED_GMX_FLOAT_MIN) ieee=0
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasrt, SLASRT)("D"
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasq3, SLASQ3)(int *i0
L70	PLMD::lapack::SLASQ5 (i0, n0, &z__[1], pp, &tau, dmin__, &dmin1, &dmin2, &dn, &dn1, &dn2, ieee)
L80	PLMD::lapack::SLASQ6 (i0, n0, &z__[1], pp, dmin__, &dmin1, &dmin2, &dn, &dn1, &dn2)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasq4, SLASQ4)(int *i0
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasq5, SLASQ5)(int *i0
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasq6, SLASQ6)(int *i0
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasr, SLASR)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasrt2, SLASRT2)(const char *id
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slasv2, SLASV2)(float *f
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (slatrd, SLATRD)(const char *uplo
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorg2r, SORG2R)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorgbr, SORGBR)(const char *vect
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorglq, SORGLQ)(&i1
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorgl2, SORGL2)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorglq, SORGLQ)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorgqr, SORGQR)(int *m
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorm2l, SORM2L)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorm2r, SORM2R)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sormbr, SORMBR)(const char *vect
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sorml2, SORML2)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sormlq, SORMLQ)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sormql, SORMQL)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sormqr, SORMQR)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sormtr, SORMTR)(const char *side
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sstebz, SSTEBZ)(const char *range
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sstegr, SSTEGR)(const char *jobz
	PLMD::lapack::if (fabs(scale-1.0)>PLUMED_GMX_FLOAT_EPS)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sstein, SSTEIN)(int *n
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (ssteqr, SSTEQR)(const char *compz
L160	PLMD::lapack::SLASRT ("I", n, &d__[1], info)
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (ssterf, SSTERF)(int *n
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (sstevr, SSTEVR)(const char *jobz
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (ssyevr, SSYEVR)(const char *jobz
	PLMD::lapack::PLUMED_BLAS_F77_FUNC (ssytrd, SSYTRD)(uplo
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (ssytd2, SSYTD2)(const char *uplo
void	PLMD::lapack::PLUMED_BLAS_F77_FUNC (strti2, STRTI2)(const char *uplo

Variables
void const char *	PLMD::lapack::compq
void const char int *	PLMD::lapack::n
void const char int double *	PLMD::lapack::d__ = 0.
void const char int double double *	PLMD::lapack::e
void const char int double double double *	PLMD::lapack::u = u_offset
void const char int double double double int *	PLMD::lapack::ldu
void const char int double double double int double *	PLMD::lapack::vt = vt_offset
void const char int double double double int double int *	PLMD::lapack::ldvt
void const char int double double double int double int double *	PLMD::lapack::q = q_offset
void const char int double double double int double int double int *	PLMD::lapack::iq = ivt2 + ldvt2 * m
void const char int double double double int double int double int double *	PLMD::lapack::work = ws
void const char int double double double int double int double int double int *	PLMD::lapack::iwork = liwmin
void const char int double double double int double int double int double int int *	PLMD::lapack::info
int	PLMD::lapack::i__ = i1
int	PLMD::lapack::j = (1 << nlvl)
int	PLMD::lapack::k = iu = z__ = ic = is = ivt = difl = difr = perm = 0
double	PLMD::lapack::p = d__[l]
double	PLMD::lapack::r__
int	PLMD::lapack::z__ = 1.
int	PLMD::lapack::ic = jc = 0
int	PLMD::lapack::ii
int	PLMD::lapack::kk
double	PLMD::lapack::cs = f1 / *r
int	PLMD::lapack::is
int	PLMD::lapack::iu
double	PLMD::lapack::sn = g1 / *r
int	PLMD::lapack::nm1 = *n - 1
double	PLMD::lapack::eps = PLUMED_GMX_DOUBLE_EPS
int	PLMD::lapack::ivt
int	PLMD::lapack::difl = d__[1]
int	PLMD::lapack::difr = difr_offset
int	PLMD::lapack::ierr
int	PLMD::lapack::perm = perm_offset
int	PLMD::lapack::mlvl
int	PLMD::lapack::sqre = 0
int	PLMD::lapack::poles = givnum = givptr = givcol = 0
int	PLMD::lapack::iuplo = (*uplo=='U' || *uplo=='u') ? 1 : 2
int	PLMD::lapack::nsize
int	PLMD::lapack::start = 1
int	PLMD::lapack::givcol = givcol_offset
int	PLMD::lapack::icompq
double	PLMD::lapack::orgnrm = PLUMED_BLAS_F77_FUNC(dlanst,DLANST)("M", n, &d__[1], &e[1])
int	PLMD::lapack::givnum = givnum_offset
int	PLMD::lapack::givptr
int	PLMD::lapack::qstart = 3
int	PLMD::lapack::smlsiz = DBDSDC_SMALLSIZE
int	PLMD::lapack::wstart = 1
int	PLMD::lapack::smlszp = smlsiz + 1
double	PLMD::lapack::zero = 0.0
double	PLMD::lapack::one = 1.0
int	PLMD::lapack::c_0 = 0
int	PLMD::lapack::c_1 = 1
	PLMD::lapack::u_dim1 = *ldu
	PLMD::lapack::u_offset = 1 + u_dim1
	PLMD::lapack::vt_dim1 = *ldvt
	PLMD::lapack::vt_offset = 1 + vt_dim1
	PLMD::lapack::break
case	PLMD::lapack::i = 0
default	PLMD::lapack::__pad0__
if n	PLMD::lapack::return
	PLMD::lapack::i__1 = *n
L40	PLMD::lapack::__pad1__
void int int *	PLMD::lapack::ncvt
void int int int *	PLMD::lapack::nru
void int int int int *	PLMD::lapack::ncc
void int int int int double double double int double int double *	PLMD::lapack::c__ = c_offset
void int int int int double double double int double int double int *	PLMD::lapack::ldc
int	PLMD::lapack::c_dim1 = *ldc
int	PLMD::lapack::c_offset = 1 + c_dim1
int	PLMD::lapack::i__2 = i3
double	PLMD::lapack::r__1 = 10.f
double	PLMD::lapack::r__2 = (r__3<r__4) ? r__3 : r__4
double	PLMD::lapack::r__3 = 100.f
double	PLMD::lapack::r__4 = pow(PLUMED_GMX_DOUBLE_EPS,c_b15)
double	PLMD::lapack::c_b15 = -.125
int	PLMD::lapack::c__1 = 1
double	PLMD::lapack::c_b49 = 1.f
double	PLMD::lapack::c_b72 = -1.f
double	PLMD::lapack::f = *finit + *tau * fc
double	PLMD::lapack::g = (d__[l - 1] - p) / (e[l - 1] * 2.)
double	PLMD::lapack::h__
int	PLMD::lapack::m = *n
int	PLMD::lapack::ll = 0
double	PLMD::lapack::mu
int	PLMD::lapack::nm12 = nm1 + nm1
int	PLMD::lapack::nm13 = nm12 + nm1
int	PLMD::lapack::lll
double	PLMD::lapack::sll
double	PLMD::lapack::tol = tolmul * eps
double	PLMD::lapack::abse
int	PLMD::lapack::idir = 0
double	PLMD::lapack::abss
int	PLMD::lapack::oldm = -1
double	PLMD::lapack::cosl
int	PLMD::lapack::isub
int	PLMD::lapack::iter = 0
double	PLMD::lapack::unfl = PLUMED_GMX_DOUBLE_MIN/PLUMED_GMX_DOUBLE_EPS
double	PLMD::lapack::sinl
double	PLMD::lapack::cosr
double	PLMD::lapack::smin = smax
double	PLMD::lapack::smax = 0.f
double	PLMD::lapack::sinr
double	PLMD::lapack::oldcs
int	PLMD::lapack::oldll = -1
double	PLMD::lapack::shift
double	PLMD::lapack::sigmn
double	PLMD::lapack::oldsn
int	PLMD::lapack::maxit = *n * 6 * *n
double	PLMD::lapack::sminl = 0.f
double	PLMD::lapack::sigmx
int	PLMD::lapack::lower = (xuplo == 'L')
double	PLMD::lapack::sminoa
double	PLMD::lapack::thresh = (r__1>r__2) ? r__1 : r__2
int	PLMD::lapack::rotate = *ncvt > 0 || *nru > 0 || *ncc > 0
double	PLMD::lapack::sminlo
double	PLMD::lapack::tolmul = (r__1>r__2) ? r__1 : r__2
int	PLMD::lapack::itmp1 = (*n > 1) ? *n : 1
int	PLMD::lapack::itmp2 = (*nru > 1) ? *nru : 1
double	PLMD::lapack::ftmp
	PLMD::lapack::else
goto	PLMD::lapack::L90
L90	PLMD::lapack::__pad2__
L160	PLMD::lapack::__pad3__
goto	PLMD::lapack::L220
L200	PLMD::lapack::__pad4__
L220	PLMD::lapack::__pad5__
void int double *	PLMD::lapack::a = a_offset
void int double int *	PLMD::lapack::lda
void int double int double *	PLMD::lapack::d = fabs(d[*n-1])
void int double int double double double *	PLMD::lapack::tauq
void int double int double double double double *	PLMD::lapack::taup
void int double int double double double double double int *	PLMD::lapack::lwork
int	PLMD::lapack::i_
int	PLMD::lapack::nx = 0
int	PLMD::lapack::nb = DGEBRD_BLOCKSIZE
double	PLMD::lapack::ws = (*m > *n) ? *m : *n
int	PLMD::lapack::nbmin = 2
int	PLMD::lapack::iinfo
int	PLMD::lapack::minmn = (*m < *n) ? *m : *n
int	PLMD::lapack::ldwrkx = *m
int	PLMD::lapack::ldwrky = *n
double	PLMD::lapack::minusone = -1.0
	PLMD::lapack::a_dim1 = *lda
	PLMD::lapack::a_offset = 1 + a_dim1
	PLMD::lapack::i_1 = minmn - nx
	PLMD::lapack::i_2 = nb
void int double int double *	PLMD::lapack::tau = 0.f
double	PLMD::lapack::aii
	PLMD::lapack::i__4 = (*m > 1) ? *m : 1
int	PLMD::lapack::ib
int	PLMD::lapack::iws = *m
int	PLMD::lapack::ldwork = *n
int	PLMD::lapack::lwkopt = *m * nb
int	PLMD::lapack::i1 = 1
int	PLMD::lapack::i2 = *ifirst
int	PLMD::lapack::i3 = iseed[3]
void int int double int double *	PLMD::lapack::s = work[inds] - *sigma
int	PLMD::lapack::ie
double	PLMD::lapack::dum [1]
int	PLMD::lapack::iscl = 0
double	PLMD::lapack::anrm = PLUMED_BLAS_F77_FUNC(dlange,DLANGE)("M", m, n, &a[a_offset], lda, dum)
int	PLMD::lapack::idum [1]
int	PLMD::lapack::itau
int	PLMD::lapack::wrkbl
int	PLMD::lapack::itaup
int	PLMD::lapack::itauq
int	PLMD::lapack::mnthr = (int) (minmn * 11. / 6.)
int	PLMD::lapack::wntqa = (*jobz=='a' || *jobz=='A')
int	PLMD::lapack::nwork
int	PLMD::lapack::wntqn = (*jobz=='o' || *jobz=='O')
int	PLMD::lapack::wntqo = (*jobz=='n' || *jobz=='N')
int	PLMD::lapack::wntqs = (*jobz=='s' || *jobz=='S')
int	PLMD::lapack::bdspac
double	PLMD::lapack::bignum = 1. / smlnum
int	PLMD::lapack::minwrk = 1
int	PLMD::lapack::ldwrku
int	PLMD::lapack::maxwrk = 1
int	PLMD::lapack::ldwkvt
double	PLMD::lapack::smlnum = sqrt(safemin) / eps
double	PLMD::lapack::minval = PLUMED_GMX_DOUBLE_MIN
double	PLMD::lapack::safemin = minval / eps
int	PLMD::lapack::wntqas = wntqa || wntqs
int	PLMD::lapack::lquery = *lwork == -1
int	PLMD::lapack::c__0 = 0
void int double int int *	PLMD::lapack::ipiv
double	PLMD::lapack::tmp = work[indp + r1 - 1]
int	PLMD::lapack::l
	PLMD::lapack::mindim = (*m < *n) ? *m : *n
int	PLMD::lapack::jb
int	PLMD::lapack::jj
int	PLMD::lapack::jp
int	PLMD::lapack::nn = (*n0 << 2) + *pp
double	PLMD::lapack::c_b20 = -1.
double	PLMD::lapack::c_b22 = 1.
Non	PLMD::lapack::unit
void int int *	PLMD::lapack::nrhs
void int int double int int double *	PLMD::lapack::b = b_offset
void int int double int int double int *	PLMD::lapack::ldb
int	PLMD::lapack::notran = (*trans=='N' || *trans=='n')
int	PLMD::lapack::c_n1 = -1
	PLMD::lapack::b_dim1 = *ldb
	PLMD::lapack::b_offset = 1 + b_dim1
	PLMD::lapack::Lower
	PLMD::lapack::Transpose
	PLMD::lapack::Unit
void int int double int double double double double double *	PLMD::lapack::x = x_offset
void int int double int double double double double double int *	PLMD::lapack::ldx
void int int double int double double double double double int double *	PLMD::lapack::y = y_offset
void int int double int double double double double double int double int *	PLMD::lapack::ldy
int	PLMD::lapack::i__3
	PLMD::lapack::x_dim1 = *ldx
	PLMD::lapack::x_offset = 1 + x_dim1
	PLMD::lapack::y_dim1 = *ldy
	PLMD::lapack::y_offset = 1 + y_dim1
const char	PLMD::lapack::ch =toupper(*uplo)
void double double double *	PLMD::lapack::rt1
void double double double double *	PLMD::lapack::rt2
double	PLMD::lapack::ab = fabs(tb)
double	PLMD::lapack::df = *a - *c__
double	PLMD::lapack::tb = *b + *b
double	PLMD::lapack::sm = *a + *c__
double	PLMD::lapack::rt
double	PLMD::lapack::adf = fabs(df)
double	PLMD::lapack::acmn = *a
double	PLMD::lapack::acmx
void int *	PLMD::lapack::nitmax
void int int int *	PLMD::lapack::mmax
void int int int int *	PLMD::lapack::minp
void int int int int int double *	PLMD::lapack::abstol
void int int int int int double double *	PLMD::lapack::reltol
void int int int int int double double double *	PLMD::lapack::pivmin = 1.
void int int int int int double double double double double double *	PLMD::lapack::e2
void int int int int int double double double double double double int *	PLMD::lapack::nval
void int int int int int double double double double double double int double double int *	PLMD::lapack::mout = kl
void int int int int int double double double double double double int double double int int *	PLMD::lapack::nab = nab_offset
double	PLMD::lapack::d__1 = fabs(*alpha)
double	PLMD::lapack::d__2 = fabs(*beta)
double	PLMD::lapack::d__3 = fabs(d__[1]) + fabs(e[1])
double	PLMD::lapack::d__4 = fabs(d__[*n]) + fabs(e[*n - 1])
int	PLMD::lapack::kf = 1
int	PLMD::lapack::ji
int	PLMD::lapack::kl = *minp
int	PLMD::lapack::jit
double	PLMD::lapack::tmp1
double	PLMD::lapack::tmp2
int	PLMD::lapack::kfnew
int	PLMD::lapack::klnew
	PLMD::lapack::nab_dim1 = *mmax
	PLMD::lapack::nab_offset = 1 + nab_dim1
	PLMD::lapack::ab_dim1 = *mmax
	PLMD::lapack::ab_offset = 1 + ab_dim1
if	PLMD::lapack::ijob
void int *	PLMD::lapack::orgati
void int double *	PLMD::lapack::rho = PLUMED_BLAS_F77_FUNC(dnrm2,DNRM2)(k, &z__[1], &c__1)
void int double double double double *	PLMD::lapack::finit
double	PLMD::lapack::fc
double	PLMD::lapack::ddf = 0.f
double	PLMD::lapack::eta = tau / (d__[*n] + sqrt(d__[*n] * d__[*n] + tau))
double	PLMD::lapack::base = 2
double	PLMD::lapack::temp = (r__1>r__2) ? r__1 : r__2
double	PLMD::lapack::temp1 = temp / (d__[*n] + sqrt(d__[*n] * d__[*n] + temp))
double	PLMD::lapack::temp2
double	PLMD::lapack::temp3
double	PLMD::lapack::temp4
int	PLMD::lapack::scale = 0
int	PLMD::lapack::niter = 1
double	PLMD::lapack::small1 = pow(base, i__1)
double	PLMD::lapack::small2 = small1 * small1
double	PLMD::lapack::sminv1 = 1.f / small1
double	PLMD::lapack::sminv2 = sminv1 * sminv1
double	PLMD::lapack::dscale [3]
double	PLMD::lapack::sclfac
double	PLMD::lapack::zscale [3]
double	PLMD::lapack::erretm = 0.
double	PLMD::lapack::sclinv = 0
if	PLMD::lapack::kniter
void double double double double double *	PLMD::lapack::cs1
void double double double double double double *	PLMD::lapack::sn1
double	PLMD::lapack::ct
double	PLMD::lapack::tn
double	PLMD::lapack::acs = fabs(cs)
int	PLMD::lapack::sgn1 = 1
int	PLMD::lapack::sgn2 = -1
void double double *	PLMD::lapack::lambda
void double double double double double double int *	PLMD::lapack::in = 0
double	PLMD::lapack::tl = (*tol>eps) ? *tol : eps
double	PLMD::lapack::piv1
double	PLMD::lapack::piv2
double	PLMD::lapack::mult
double	PLMD::lapack::scale1 = fabs(a[1]) + fabs(b[1])
double	PLMD::lapack::scale2
double	PLMD::lapack::d__5
double	PLMD::lapack::ak
double	PLMD::lapack::pert
double	PLMD::lapack::absak
double	PLMD::lapack::sfmin = minval / eps
if	PLMD::lapack::job
void int *	PLMD::lapack::n2 = n - k
void int double int *	PLMD::lapack::dtrd1
void int double int int *	PLMD::lapack::dtrd2
void int double int int int *	PLMD::lapack::index
int	PLMD::lapack::n2sv = *n2
int	PLMD::lapack::ind1 = 0
int	PLMD::lapack::ind2 = *n1
double	PLMD::lapack::dtemp
double	PLMD::lapack::sum = 1.
double	PLMD::lapack::max
double	PLMD::lapack::val
double const char *	PLMD::lapack::uplo
double	PLMD::lapack::ret_val = value
double	PLMD::lapack::absa
double	PLMD::lapack::value =0.0
else if	PLMD::lapack::norm
double	PLMD::lapack::w
double	PLMD::lapack::z = xabs
	PLMD::lapack::xabs = fabs(*x)
	PLMD::lapack::yabs = fabs(*y)
void int *	PLMD::lapack::b1
void int int *	PLMD::lapack::bn
void int int double *	PLMD::lapack::sigma = 0
void int int double double double double *	PLMD::lapack::ld
void int int double double double double double *	PLMD::lapack::lld
void int int double double double double double double *	PLMD::lapack::eval
void int int double double double double double double double *	PLMD::lapack::gersch
void int int double double double double double double double double double *	PLMD::lapack::ztz = 1.
void int int double double double double double double double double double double *	PLMD::lapack::mingma = work[inds + r1 - 1] + work[indp + r1 - 1]
void int int double double double double double double double double double double int int *	PLMD::lapack::isuppz = *b1
int	PLMD::lapack::r1
int	PLMD::lapack::r2 = *bn
int	PLMD::lapack::to
int	PLMD::lapack::indp = *n * 3 + 1
int	PLMD::lapack::inds = (*n << 1) + 1
int	PLMD::lapack::from
double	PLMD::lapack::dplus = d__[1] + s
int	PLMD::lapack::sawnan = 0
int	PLMD::lapack::indumn
double	PLMD::lapack::dminus
goto	PLMD::lapack::L40
void int int double *	PLMD::lapack::v = v_offset
void int int double int *	PLMD::lapack::incv
void int int double int double double *	PLMD::lapack::c
double	PLMD::lapack::minustau = -(*tau)
void const char *	PLMD::lapack::trans
void const char const char *	PLMD::lapack::direct
void const char const char const char *	PLMD::lapack::storev
void const char const char const char int int int double int *	PLMD::lapack::ldv
void const char const char const char int int int double int double *	PLMD::lapack::t = t_offset
void const char const char const char int int int double int double int *	PLMD::lapack::ldt
char	PLMD::lapack::transt [1]
	PLMD::lapack::v_dim1 = *ldv
	PLMD::lapack::v_offset = 1 + v_dim1
	PLMD::lapack::t_dim1 = *ldt
	PLMD::lapack::t_offset = 1 + t_dim1
	PLMD::lapack::work_dim1 = *ldwork
	PLMD::lapack::work_offset = 1 + work_dim1
void double *	PLMD::lapack::alpha = beta
void double double int *	PLMD::lapack::incx
int	PLMD::lapack::ti1 = *n-1
int	PLMD::lapack::knt
double	PLMD::lapack::safmin = minval*(1.0+PLUMED_GMX_DOUBLE_EPS) / PLUMED_GMX_DOUBLE_EPS
double	PLMD::lapack::rsafmn
double	PLMD::lapack::beta = t
	PLMD::lapack::xnorm = PLUMED_BLAS_F77_FUNC(dnrm2,DNRM2)(&ti1,x,incx)
double	PLMD::lapack::vii
void int *	PLMD::lapack::iseed = it1
int	PLMD::lapack::il = -1
int	PLMD::lapack::iv
int	PLMD::lapack::il2
void double double double double int *	PLMD::lapack::ifirst
void double double double double int int *	PLMD::lapack::ilast
void double double double double int int double *	PLMD::lapack::rtol1
void double double double double int int double double *	PLMD::lapack::rtol2
void double double double double int int double double int *	PLMD::lapack::offset
void double double double double int int double double int double double *	PLMD::lapack::wgap
void double double double double int int double double int double double double *	PLMD::lapack::werr
double	PLMD::lapack::fac
double	PLMD::lapack::gap
double	PLMD::lapack::mid
int	PLMD::lapack::cnt
double	PLMD::lapack::left = (*side=='L' || *side=='l')
int	PLMD::lapack::nint = 0
int	PLMD::lapack::prev = 0
int	PLMD::lapack::next
int	PLMD::lapack::nleft
double	PLMD::lapack::right = 0.0
double	PLMD::lapack::width
double	PLMD::lapack::error
int	PLMD::lapack::nright
int	PLMD::lapack::olnint = nint
goto	PLMD::lapack::L30
void int double *	PLMD::lapack::vl = 0.
void int double double *	PLMD::lapack::vu
void int double double int int double double double int *	PLMD::lapack::nsplit = 1
void int double double int int double double double int int *	PLMD::lapack::isplit = *n
void int double double int int double double double int int int double int *	PLMD::lapack::iblock = 1
void int double double int int double double double int int int double int int *	PLMD::lapack::indexw
double	PLMD::lapack::gl
double	PLMD::lapack::gu
double	PLMD::lapack::nrm
double	PLMD::lapack::vvl
double	PLMD::lapack::vvu
double	PLMD::lapack::offd
int	PLMD::lapack::iend = 0
int	PLMD::lapack::jblk
int	PLMD::lapack::till
int	PLMD::lapack::itmp
double	PLMD::lapack::rtol
double	PLMD::lapack::delta = eps * 2.
int	PLMD::lapack::ibegin = 1
int	PLMD::lapack::irange = 0
double	PLMD::lapack::sgndef = 0
int	PLMD::lapack::maxcnt = 0
void double double double double int int double double double double *	PLMD::lapack::lplus
double	PLMD::lapack::dmax1 = fabs(dplus[1])
double	PLMD::lapack::dmax2 = fabs(work[1])
void double double int int double int int double double double int *	PLMD::lapack::ldz
double	PLMD::lapack::c_b5 = 0.
int	PLMD::lapack::c__2 = 2
int	PLMD::lapack::im
int	PLMD::lapack::zto
int	PLMD::lapack::wend
int	PLMD::lapack::ktot
int	PLMD::lapack::indld = *n << 1
int	PLMD::lapack::ndone
int	PLMD::lapack::iindr = *n
double	PLMD::lapack::resid
int	PLMD::lapack::nomgs
int	PLMD::lapack::nclus
int	PLMD::lapack::zfrom
int	PLMD::lapack::iindc1 = *n << 1
int	PLMD::lapack::iindc2 = *n * 3
int	PLMD::lapack::indgap = *n << 2
int	PLMD::lapack::indlld = *n * 3
int	PLMD::lapack::oldien
int	PLMD::lapack::oldncl
int	PLMD::lapack::wbegin = 1
double	PLMD::lapack::relgap
int	PLMD::lapack::oldcls
int	PLMD::lapack::ndepth
int	PLMD::lapack::inderr = *n
int	PLMD::lapack::iindwk = (*n << 2) + 1
int	PLMD::lapack::newcls
int	PLMD::lapack::oldfst
double	PLMD::lapack::minrgp =0.0
int	PLMD::lapack::indwrk = *n * 5 + 1
int	PLMD::lapack::oldlst
int	PLMD::lapack::newfrs
int	PLMD::lapack::newftt
int	PLMD::lapack::parity
double	PLMD::lapack::mgstol = eps * 100.
double	PLMD::lapack::nrminv
double	PLMD::lapack::rqcorr
int	PLMD::lapack::newlst
int	PLMD::lapack::newsiz
	PLMD::lapack::z_dim1 = *ldz
	PLMD::lapack::z_offset = 1 + z_dim1
void double double double double *	PLMD::lapack::r
double	PLMD::lapack::f1
double	PLMD::lapack::g1 = *g
double	PLMD::lapack::f1a = fabs(f1)
double	PLMD::lapack::g1a = fabs(g1)
int	PLMD::lapack::count
	PLMD::lapack::safemin2 = pow(2,n)
	PLMD::lapack::safemx2 = 1.0 / safemin2
int	PLMD::lapack::i4 = iseed[4]
int	PLMD::lapack::it1 = it2 = it3 = it4 = 0
int	PLMD::lapack::it2
int	PLMD::lapack::it3
int	PLMD::lapack::it4
void double double *	PLMD::lapack::h
void double double double *	PLMD::lapack::ssmin
void double double double double *	PLMD::lapack::ssmax
double	PLMD::lapack::ga = fabs(*g)
double	PLMD::lapack::ha = fabs(*h)
double	PLMD::lapack::fhmin = (fa<ha) ? fa : ha
double	PLMD::lapack::fhmax = (fa>ha) ? fa : ha
double	PLMD::lapack::tmax
double	PLMD::lapack::tmin
double	PLMD::lapack::as
double	PLMD::lapack::at
double	PLMD::lapack::au
void int int *	PLMD::lapack::ku
void int int double *	PLMD::lapack::cfrom
void int int double double *	PLMD::lapack::cto
int	PLMD::lapack::k1
int	PLMD::lapack::k2 = n + 1
int	PLMD::lapack::k3
int	PLMD::lapack::k4
int	PLMD::lapack::done =0
double	PLMD::lapack::cfromc = *cfrom
double	PLMD::lapack::ctoc = *cto
double	PLMD::lapack::cfrom1
double	PLMD::lapack::cto1
double	PLMD::lapack::mul
int	PLMD::lapack::lf
int	PLMD::lapack::nd = llst*2 - 1
int	PLMD::lapack::nl
int	PLMD::lapack::nr
int	PLMD::lapack::im1
int	PLMD::lapack::nlf
int	PLMD::lapack::nrf
int	PLMD::lapack::iwk = idxq + *n
int	PLMD::lapack::lvl = 1 + (int) temp
int	PLMD::lapack::ndb1 = (nd + 1) / 2
int	PLMD::lapack::nlp1 = *nl + 1
int	PLMD::lapack::nrp1 = *nr + *sqre
int	PLMD::lapack::idxq = ndimr + *n
int	PLMD::lapack::nlvl
int	PLMD::lapack::inode = 1
int	PLMD::lapack::ndiml = inode + *n
int	PLMD::lapack::idxqc
int	PLMD::lapack::ndimr = ndiml + *n
int	PLMD::lapack::itemp
int	PLMD::lapack::sqrei
int	PLMD::lapack::n1 = k
int	PLMD::lapack::iz = 1
int	PLMD::lapack::iu2 = isigma + n
int	PLMD::lapack::ldq = k
int	PLMD::lapack::idx = 1
int	PLMD::lapack::ldu2 = n
int	PLMD::lapack::ivt2 = iu2 + ldu2 * n
int	PLMD::lapack::idxc = idx + n
int	PLMD::lapack::idxp = coltyp + n
int	PLMD::lapack::ldvt2 = m
int	PLMD::lapack::isigma = iz + m
int	PLMD::lapack::coltyp = idxc + n
void int int int double double double double double int double int double *	PLMD::lapack::dsigma = d__[jprev]
void int int int double double double double double int double int double double *	PLMD::lapack::u2 = u2_offset
void int int int double double double double double int double int double double int double *	PLMD::lapack::vt2 = vt2_offset
int	PLMD::lapack::vt2_dim1 = *ldvt2
int	PLMD::lapack::vt2_offset = 1 + vt2_dim1
double	PLMD::lapack::z1 = *alpha * vt[nlp1 + nlp1 * vt_dim1]
int	PLMD::lapack::psm [4]
int	PLMD::lapack::nlp2 = *nl + 2
int	PLMD::lapack::idxi
int	PLMD::lapack::idxj
int	PLMD::lapack::ctot [4]
int	PLMD::lapack::idxjp
int	PLMD::lapack::jprev = 0
double	PLMD::lapack::hlftol = tol / 2.
	PLMD::lapack::u2_dim1 = *ldu2
	PLMD::lapack::u2_offset = 1 + u2_dim1
int	PLMD::lapack::jc
int	PLMD::lapack::ctemp = ctot[2] + ctot[3]
int	PLMD::lapack::ktemp = ctot[1] + 2
double *	PLMD::lapack::p1 = &t1
double *	PLMD::lapack::p2 = &t2
double	PLMD::lapack::t1
double	PLMD::lapack::t2
	PLMD::lapack::q_dim1 = *ldq
	PLMD::lapack::q_offset = 1 + q_dim1
goto	PLMD::lapack::L100
	PLMD::lapack::N
double	PLMD::lapack::dd [3] = dtiim
double	PLMD::lapack::dw = dpsi + dphi + temp * temp
double	PLMD::lapack::zz [3] = z__[iip1] * z__[iip1]
int	PLMD::lapack::ip1 = *i__ + 1
double	PLMD::lapack::phi = z__[*n] * temp
double	PLMD::lapack::psi = 0.
int	PLMD::lapack::iim1 = ii - 1
int	PLMD::lapack::iip1 = ii + 1
double	PLMD::lapack::dphi = temp * temp
double	PLMD::lapack::dpsi = 0.
double	PLMD::lapack::prew = w
double	PLMD::lapack::sg2lb = -delsq2
double	PLMD::lapack::sg2ub = 0.
double	PLMD::lapack::dtiim
double	PLMD::lapack::delsq = (d__[ip1] - d__[*i__]) * (d__[ip1] + d__[*i__])
double	PLMD::lapack::dtiip = work[iip1] * delta[iip1]
double	PLMD::lapack::dtisq
int	PLMD::lapack::swtch = 0
double	PLMD::lapack::dtnsq = work[*n] * delta[*n]
double	PLMD::lapack::delsq2 = delsq / 2.
double	PLMD::lapack::dtnsq1 = work[*n - 1] * delta[*n - 1]
int	PLMD::lapack::swtch3 = 0
double	PLMD::lapack::dtipsq
double	PLMD::lapack::rhoinv = 1. / *rho
goto	PLMD::lapack::L240
	PLMD::lapack::del = d__[2] - d__[1]
void int int int double double *	PLMD::lapack::vf = tau
void int int int double double double double double int int int int int *	PLMD::lapack::ldgcol
void int int int double double double double double int int int int int double int *	PLMD::lapack::ldgnum
int	PLMD::lapack::iw = isigma + n
int	PLMD::lapack::ivfw = iw + m
int	PLMD::lapack::ivlw = ivfw + m
	PLMD::lapack::givcol_dim1 = *ldgcol
	PLMD::lapack::givcol_offset = 1 + givcol_dim1
	PLMD::lapack::poles_dim1 = *ldgnum
	PLMD::lapack::poles_offset = 1 + poles_dim1
	PLMD::lapack::givnum_dim1 = *ldgnum
	PLMD::lapack::givnum_offset = 1 + givnum_dim1
void int int int int double double double *	PLMD::lapack::zw = z__[jprev]
void int int int int double double double double double *	PLMD::lapack::vfw
void int int int int double double double double double double double *	PLMD::lapack::vlw
void int double double double double double double int *	PLMD::lapack::lddifr
double	PLMD::lapack::dj
int	PLMD::lapack::iwk1 = 1
int	PLMD::lapack::iwk2 = iwk1 + *k
int	PLMD::lapack::iwk3 = iwk2 + *k
int	PLMD::lapack::iwk2i = iwk2 - 1
int	PLMD::lapack::iwk3i = iwk3 - 1
double	PLMD::lapack::diflj
double	PLMD::lapack::difrj = dsigjp = 0
double	PLMD::lapack::dsigj
double	PLMD::lapack::dsigjp
	PLMD::lapack::difr_dim1 = *lddifr
	PLMD::lapack::difr_offset = 1 + difr_dim1
int	PLMD::lapack::vfi
int	PLMD::lapack::vli
int	PLMD::lapack::lvl2
int	PLMD::lapack::idxqi
int	PLMD::lapack::nwork1 = vl + m
int	PLMD::lapack::nwork2 = nwork1 + smlszp * smlszp
	PLMD::lapack::difl_dim1 = *ldu
	PLMD::lapack::difl_offset = 1 + difl_dim1
	PLMD::lapack::perm_dim1 = *ldgcol
	PLMD::lapack::perm_offset = 1 + perm_dim1
int	PLMD::lapack::np1 = *n + 1
int	PLMD::lapack::sqre1 = *sqre
void int int int int int int *	PLMD::lapack::msub
int	PLMD::lapack::ir = 0
int	PLMD::lapack::llst = 1
int	PLMD::lapack::ncrnt
int	PLMD::lapack::i0 = 1
int	PLMD::lapack::n0 = *n
int	PLMD::lapack::pp = 0
double	PLMD::lapack::dee
int	PLMD::lapack::ipn4
double	PLMD::lapack::tol2 = d__1 * d__1
int	PLMD::lapack::ieee = 1
int	PLMD::lapack::nbig
double	PLMD::lapack::dmin__
double	PLMD::lapack::emin = z__[2]
double	PLMD::lapack::emax
int	PLMD::lapack::kmin
int	PLMD::lapack::ndiv = 2*(n0 - i0)
double	PLMD::lapack::qmin
double	PLMD::lapack::qmax = 0.
double	PLMD::lapack::zmax = 0.
int	PLMD::lapack::splt
int	PLMD::lapack::nfail = 0
double	PLMD::lapack::desig = t - *sigma
double	PLMD::lapack::trace = d__ + e
double	PLMD::lapack::deemin
int	PLMD::lapack::iwhila
int	PLMD::lapack::iwhilb
double	PLMD::lapack::oldemn
double	PLMD::lapack::posinf = one/zero
double	PLMD::lapack::neginf = -one/zero
double	PLMD::lapack::negzro = one/(neginf+one)
double	PLMD::lapack::newzro = negzro + zero
double	PLMD::lapack::dmin1 = 0.
double	PLMD::lapack::dmin2 = 0.
double	PLMD::lapack::dn = 0.
double	PLMD::lapack::dn1 = 0.
double	PLMD::lapack::dn2 = 0.
int	PLMD::lapack::j4 = (*i0 << 2) + *pp - 3
int	PLMD::lapack::n0in = *n0
goto	PLMD::lapack::L10
	PLMD::lapack::ttype = -12
goto	PLMD::lapack::L70
L80	PLMD::lapack::__pad6__
L90	PLMD::lapack::__pad7__
double	PLMD::lapack::a2
double	PLMD::lapack::b2
int	PLMD::lapack::np
double	PLMD::lapack::gam
double	PLMD::lapack::gap1
double	PLMD::lapack::gap2
void int double int double double double double double double *	PLMD::lapack::dnm1 = z__[j4p2 + 2] * (*dnm2 / z__[j4 - 2]) - *tau
void int double int double double double double double double double *	PLMD::lapack::dnm2
int	PLMD::lapack::j4p2 = j4 + (*pp << 1) - 1
void const char *	PLMD::lapack::pivot
double	PLMD::lapack::stemp
if	PLMD::lapack::side
double	PLMD::lapack::d1
double	PLMD::lapack::d2
double	PLMD::lapack::d3
int	PLMD::lapack::dir = -1
int	PLMD::lapack::endd = stack[(stkpnt << 1) - 1]
int	PLMD::lapack::stack [64] = 1
double	PLMD::lapack::dmnmx
int	PLMD::lapack::stkpnt = 1
void int double int *	PLMD::lapack::key
int	PLMD::lapack::tmpkey
void double int double double *	PLMD::lapack::sumsq
double	PLMD::lapack::absxi
void double double double double double *	PLMD::lapack::snr = srt
void double double double double double double *	PLMD::lapack::csr = crt
void double double double double double double double *	PLMD::lapack::snl = slt
void double double double double double double double double *	PLMD::lapack::csl
double	PLMD::lapack::fa = fabs(ft)
double	PLMD::lapack::ft = *f
double	PLMD::lapack::gt = *g
double	PLMD::lapack::ht = *h__
double	PLMD::lapack::mm
double	PLMD::lapack::tt
double	PLMD::lapack::clt
double	PLMD::lapack::crt
double	PLMD::lapack::slt
double	PLMD::lapack::srt
int	PLMD::lapack::pmax = 1
int	PLMD::lapack::swap = ha > fa
double	PLMD::lapack::tsign =1.0
int	PLMD::lapack::gasmal
int	PLMD::lapack::ix
int	PLMD::lapack::ip
	PLMD::lapack::inc = 1
	PLMD::lapack::n32 = *n / 32
void int int double int double double double int *	PLMD::lapack::ldw
int	PLMD::lapack::ti2
int	PLMD::lapack::ti3
int	PLMD::lapack::mn = (*m < *n) ? *m : *n
	PLMD::lapack::wantq = (*vect=='Q' || *vect=='q')
	PLMD::lapack::wrksz = mn*DORGBR_BLOCKSIZE
int	PLMD::lapack::ki = 0
int	PLMD::lapack::mi
int	PLMD::lapack::ni = *n - 1
int	PLMD::lapack::nq
int	PLMD::lapack::nw = *m
int	PLMD::lapack::applyq = (*vect=='Q' || *vect=='q')
const char	PLMD::lapack::xtrans =toupper(*trans)
int	PLMD::lapack::c__65 = 65
int	PLMD::lapack::upper = (*uplo=='U' || *uplo=='u')
void const char *	PLMD::lapack::order
int	PLMD::lapack::c__3 = 3
int	PLMD::lapack::je
double	PLMD::lapack::wl
double	PLMD::lapack::wu = 0.
int	PLMD::lapack::nwl = 0
double	PLMD::lapack::ulp = 2*PLUMED_GMX_DOUBLE_EPS
double	PLMD::lapack::wlu
double	PLMD::lapack::wul
int	PLMD::lapack::nwu = 0
int	PLMD::lapack::ioff
int	PLMD::lapack::iout
int	PLMD::lapack::jdisc
double	PLMD::lapack::atoli
int	PLMD::lapack::iwoff
double	PLMD::lapack::bnorm
int	PLMD::lapack::itmax
double	PLMD::lapack::wkill
double	PLMD::lapack::rtoli = ulp * 2.
double	PLMD::lapack::tnorm = (d__3>d__4) ? d__3 : d__4
int	PLMD::lapack::idiscl
int	PLMD::lapack::idumma [1]
int	PLMD::lapack::idiscu
int	PLMD::lapack::iorder
int	PLMD::lapack::ncnvrg = 0
int	PLMD::lapack::toofew = 0
const double	PLMD::lapack::safemn = PLUMED_GMX_DOUBLE_MIN*(1.0+PLUMED_GMX_DOUBLE_EPS)
if	PLMD::lapack::range
void const char int double double double double int int double int double double int int double int int int *	PLMD::lapack::liwork
double	PLMD::lapack::rmin = sqrt(smlnum)
double	PLMD::lapack::rmax = (d__1<d__2) ? d__1 : d__2
double	PLMD::lapack::tnrm = PLUMED_BLAS_F77_FUNC(dlanst,DLANST)("M", n, &d__[1], &e[1])
int	PLMD::lapack::iindw = (*n << 1) + 1
int	PLMD::lapack::lwmin = *n * 17
int	PLMD::lapack::wantz = (*jobz=='V' || *jobz=='v')
int	PLMD::lapack::iindbl = *n + 1
int	PLMD::lapack::valeig = (*range=='V' || *range=='v')
int	PLMD::lapack::alleig = (*range=='A' || *range=='a')
int	PLMD::lapack::indeig = (*range=='I' || *range=='i')
int	PLMD::lapack::indgrs = 1
int	PLMD::lapack::iinspl = 1
int	PLMD::lapack::liwmin = *n * 10
void double double int double int int double int double int int *	PLMD::lapack::ifail
int	PLMD::lapack::j1 = 1
double	PLMD::lapack::xj
double	PLMD::lapack::scl
double	PLMD::lapack::sep
int	PLMD::lapack::its
double	PLMD::lapack::xjm = 0.0
double	PLMD::lapack::ztr
double	PLMD::lapack::eps1
int	PLMD::lapack::nblk
int	PLMD::lapack::jmax
int	PLMD::lapack::gpind
double	PLMD::lapack::ortol
int	PLMD::lapack::indrv1 = 0
int	PLMD::lapack::indrv2 = indrv1 + *n
int	PLMD::lapack::indrv3 = indrv2 + *n
int	PLMD::lapack::indrv4 = indrv3 + *n
int	PLMD::lapack::indrv5 = indrv4 + *n
int	PLMD::lapack::nrmchk
int	PLMD::lapack::blksiz
double	PLMD::lapack::onenrm
double	PLMD::lapack::dtpcrt
double	PLMD::lapack::pertol
double	PLMD::lapack::c_b10 = 1.
int	PLMD::lapack::l1 = 1
int	PLMD::lapack::lm1 = l - 1
int	PLMD::lapack::mm1
int	PLMD::lapack::lsv = l
double	PLMD::lapack::tst
double	PLMD::lapack::eps2 = d__1 * d__1
int	PLMD::lapack::lend = m
int	PLMD::lapack::jtot = 0
double	PLMD::lapack::anorm = PLUMED_BLAS_F77_FUNC(dlanst,DLANST)("I", &i__1, &d__[l], &e[l])
int	PLMD::lapack::lendm1
int	PLMD::lapack::lendp1
int	PLMD::lapack::iscale = 0
double	PLMD::lapack::safmax = 1. / safmin
int	PLMD::lapack::lendsv = lend
double	PLMD::lapack::ssfmin = sqrt(safmin) / eps2
int	PLMD::lapack::nmaxit = *n * 30
int	PLMD::lapack::icompz
double	PLMD::lapack::ssfmax = sqrt(safmax) / 3.
if	PLMD::lapack::compz
goto	PLMD::lapack::L140
goto	PLMD::lapack::L190
double	PLMD::lapack::bb
double	PLMD::lapack::rte = sqrt(e[l - 1])
double	PLMD::lapack::oldc
double	PLMD::lapack::gamma = d__[m] - sigma
double	PLMD::lapack::oldgam
double	PLMD::lapack::c_b32 = 1.
goto	PLMD::lapack::L150
double	PLMD::lapack::vll = *vl
double	PLMD::lapack::vuu = *vu
int	PLMD::lapack::indd = inde + *n
int	PLMD::lapack::inde = indtau + *n
int	PLMD::lapack::imax
int	PLMD::lapack::inddd = indee + *n
int	PLMD::lapack::indee = indd + *n
int	PLMD::lapack::indwk = indifl + *n
int	PLMD::lapack::ieeeok = 1
int	PLMD::lapack::indibl = 0
int	PLMD::lapack::indifl = inddd + *n
double	PLMD::lapack::abstll = *abstol
int	PLMD::lapack::indtau = 1
int	PLMD::lapack::indwkn
int	PLMD::lapack::llwrkn
int	PLMD::lapack::llwork = *lwork - indwk + 1
double	PLMD::lapack::fzero = 0.0
double	PLMD::lapack::fone = 1.0
double	PLMD::lapack::taui
double	PLMD::lapack::c_b23 = 1.
void const char *	PLMD::lapack::diag
double	PLMD::lapack::ajj
int	PLMD::lapack::nounit = (*diag=='N' || *diag=='n')
double	PLMD::lapack::c_b18 = 1.

Macro Definition Documentation

DORGLQ_BLOCKSIZE

#define DORGLQ_BLOCKSIZE   32

DORGLQ_CROSSOVER

#define DORGLQ_CROSSOVER   128

DORGLQ_MINBLOCKSIZE

#define DORGLQ_MINBLOCKSIZE   2

SORGLQ_BLOCKSIZE

#define SORGLQ_BLOCKSIZE   32

SORGLQ_CROSSOVER

#define SORGLQ_CROSSOVER   128

SORGLQ_MINBLOCKSIZE

#define SORGLQ_MINBLOCKSIZE   2