This module '_flapack' is auto-generated with f2py (version:2). Functions: ba,lo,hi,pivscale,info = sgebal(a,scale=0,permute=0,overwrite_a=0) ba,lo,hi,pivscale,info = dgebal(a,scale=0,permute=0,overwrite_a=0) ba,lo,hi,pivscale,info = cgebal(a,scale=0,permute=0,overwrite_a=0) ba,lo,hi,pivscale,info = zgebal(a,scale=0,permute=0,overwrite_a=0) ht,tau,info = sgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) ht,tau,info = dgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) ht,tau,info = cgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) ht,tau,info = zgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) work,info = sgehrd_lwork(n,lo=0,hi=n-1) work,info = dgehrd_lwork(n,lo=0,hi=n-1) work,info = cgehrd_lwork(n,lo=0,hi=n-1) work,info = zgehrd_lwork(n,lo=0,hi=n-1) lu,piv,x,info = sgesv(a,b,overwrite_a=0,overwrite_b=0) lu,piv,x,info = dgesv(a,b,overwrite_a=0,overwrite_b=0) lu,piv,x,info = cgesv(a,b,overwrite_a=0,overwrite_b=0) lu,piv,x,info = zgesv(a,b,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = sgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = dgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = cgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = zgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) rcond,info = sgecon(a,anorm,norm='1') rcond,info = dgecon(a,anorm,norm='1') rcond,info = cgecon(a,anorm,norm='1') rcond,info = zgecon(a,anorm,norm='1') lu,piv,info = sgetrf(a,overwrite_a=0) lu,piv,info = dgetrf(a,overwrite_a=0) lu,piv,info = cgetrf(a,overwrite_a=0) lu,piv,info = zgetrf(a,overwrite_a=0) x,info = sgetrs(lu,piv,b,trans=0,overwrite_b=0) x,info = dgetrs(lu,piv,b,trans=0,overwrite_b=0) x,info = cgetrs(lu,piv,b,trans=0,overwrite_b=0) x,info = zgetrs(lu,piv,b,trans=0,overwrite_b=0) lu,ipiv,jpiv,info = sgetc2(a,overwrite_a=0) lu,ipiv,jpiv,info = dgetc2(a,overwrite_a=0) lu,ipiv,jpiv,info = cgetc2(a,overwrite_a=0) lu,ipiv,jpiv,info = zgetc2(a,overwrite_a=0) x,scale = sgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) x,scale = dgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) x,scale = cgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) x,scale = zgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) inv_a,info = sgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) inv_a,info = dgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) inv_a,info = cgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) inv_a,info = zgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) work,info = sgetri_lwork(n) work,info = dgetri_lwork(n) work,info = cgetri_lwork(n) work,info = zgetri_lwork(n) u,s,vt,info = sgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?4*minmn*minmn+MAX(m,n)+9*minmn:MAX(14*minmn+4,10*minmn+2+25*(25+8))+MAX(m,n)),1),overwrite_a=0) u,s,vt,info = dgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?4*minmn*minmn+MAX(m,n)+9*minmn:MAX(14*minmn+4,10*minmn+2+25*(25+8))+MAX(m,n)),1),overwrite_a=0) work,info = sgesdd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = dgesdd_lwork(m,n,compute_uv=1,full_matrices=1) u,s,vt,info = cgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?2*minmn*minmn+MAX(m,n)+2*minmn:2*minmn+MAX(m,n)),1),overwrite_a=0) u,s,vt,info = zgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?2*minmn*minmn+MAX(m,n)+2*minmn:2*minmn+MAX(m,n)),1),overwrite_a=0) work,info = cgesdd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = zgesdd_lwork(m,n,compute_uv=1,full_matrices=1) u,s,vt,info = sgesvd(a,compute_uv=1,full_matrices=1,lwork=max(MAX(3*minmn+MAX(m,n),5*minmn),1),overwrite_a=0) u,s,vt,info = dgesvd(a,compute_uv=1,full_matrices=1,lwork=max(MAX(3*minmn+MAX(m,n),5*minmn),1),overwrite_a=0) work,info = sgesvd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = dgesvd_lwork(m,n,compute_uv=1,full_matrices=1) u,s,vt,info = cgesvd(a,compute_uv=1,full_matrices=1,lwork=MAX(2*minmn+MAX(m,n),1),overwrite_a=0) u,s,vt,info = zgesvd(a,compute_uv=1,full_matrices=1,lwork=MAX(2*minmn+MAX(m,n),1),overwrite_a=0) work,info = cgesvd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = zgesvd_lwork(m,n,compute_uv=1,full_matrices=1) lqr,x,info = sgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) lqr,x,info = dgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) lqr,x,info = cgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) lqr,x,info = zgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) work,info = sgels_lwork(m,n,nrhs,trans='N') work,info = dgels_lwork(m,n,nrhs,trans='N') work,info = cgels_lwork(m,n,nrhs,trans='N') work,info = zgels_lwork(m,n,nrhs,trans='N') v,x,s,rank,work,info = sgelss(a,b,cond=-1.0,lwork=max(3*minmn+MAX(2*minmn,MAX(maxmn,nrhs)),1),overwrite_a=0,overwrite_b=0) v,x,s,rank,work,info = dgelss(a,b,cond=-1.0,lwork=max(3*minmn+MAX(2*minmn,MAX(maxmn,nrhs)),1),overwrite_a=0,overwrite_b=0) work,info = sgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,info = dgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) v,x,s,rank,work,info = cgelss(a,b,cond=-1.0,lwork=max(2*minmn+MAX(maxmn,nrhs),1),overwrite_a=0,overwrite_b=0) v,x,s,rank,work,info = zgelss(a,b,cond=-1.0,lwork=max(2*minmn+MAX(maxmn,nrhs),1),overwrite_a=0,overwrite_b=0) work,info = cgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,info = zgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) v,x,j,rank,info = sgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) v,x,j,rank,info = dgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) work,info = sgelsy_lwork(m,n,nrhs,cond,lwork=-1) work,info = dgelsy_lwork(m,n,nrhs,cond,lwork=-1) v,x,j,rank,info = cgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) v,x,j,rank,info = zgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) work,info = cgelsy_lwork(m,n,nrhs,cond,lwork=-1) work,info = zgelsy_lwork(m,n,nrhs,cond,lwork=-1) x,s,rank,info = sgelsd(a,b,lwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) x,s,rank,info = dgelsd(a,b,lwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) work,iwork,info = sgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,iwork,info = dgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) x,s,rank,info = cgelsd(a,b,lwork,size_rwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) x,s,rank,info = zgelsd(a,b,lwork,size_rwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) work,rwork,iwork,info = cgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,rwork,iwork,info = zgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) qr,jpvt,tau,work,info = sgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,jpvt,tau,work,info = dgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,jpvt,tau,work,info = cgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,jpvt,tau,work,info = zgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,tau,work,info = sgeqrf(a,lwork=max(3*n,1),overwrite_a=0) qr,tau,work,info = dgeqrf(a,lwork=max(3*n,1),overwrite_a=0) qr,tau,work,info = cgeqrf(a,lwork=max(3*n,1),overwrite_a=0) qr,tau,work,info = zgeqrf(a,lwork=max(3*n,1),overwrite_a=0) work,info = sgeqrf_lwork(m,n) work,info = dgeqrf_lwork(m,n) work,info = cgeqrf_lwork(m,n) work,info = zgeqrf_lwork(m,n) qr,tau,info = sgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) qr,tau,info = dgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) qr,tau,info = cgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) qr,tau,info = zgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) work,info = sgeqrfp_lwork(m,n) work,info = dgeqrfp_lwork(m,n) work,info = cgeqrfp_lwork(m,n) work,info = zgeqrfp_lwork(m,n) qr,tau,work,info = sgerqf(a,lwork=max(3*m,1),overwrite_a=0) qr,tau,work,info = dgerqf(a,lwork=max(3*m,1),overwrite_a=0) qr,tau,work,info = cgerqf(a,lwork=max(3*m,1),overwrite_a=0) qr,tau,work,info = zgerqf(a,lwork=max(3*m,1),overwrite_a=0) wr,wi,vl,vr,info = sgeev(a,compute_vl=1,compute_vr=1,lwork=max(4*n,1),overwrite_a=0) wr,wi,vl,vr,info = dgeev(a,compute_vl=1,compute_vr=1,lwork=max(4*n,1),overwrite_a=0) work,info = sgeev_lwork(n,compute_vl=1,compute_vr=1) work,info = dgeev_lwork(n,compute_vl=1,compute_vr=1) w,vl,vr,info = cgeev(a,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0) w,vl,vr,info = zgeev(a,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0) work,info = cgeev_lwork(n,compute_vl=1,compute_vr=1) work,info = zgeev_lwork(n,compute_vl=1,compute_vr=1) alphar,alphai,beta,vl,vr,info = sgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alphar,alphai,beta,vl,vr,info = dgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,info = cgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,info = zgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) t,sdim,w,vs,work,info = cgees(cselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),cselect_extra_args=(),overwrite_a=0) t,sdim,w,vs,work,info = zgees(zselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),zselect_extra_args=(),overwrite_a=0) t,sdim,wr,wi,vs,work,info = sgees(sselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),sselect_extra_args=(),overwrite_a=0) t,sdim,wr,wi,vs,work,info = dgees(dselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),dselect_extra_args=(),overwrite_a=0) a,b,sdim,alphar,alphai,beta,vsl,vsr,work,info = sgges(sselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(8*n+16,1),sselect_extra_args=(),overwrite_a=0,overwrite_b=0) a,b,sdim,alphar,alphai,beta,vsl,vsr,work,info = dgges(dselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(8*n+16,1),dselect_extra_args=(),overwrite_a=0,overwrite_b=0) a,b,sdim,alpha,beta,vsl,vsr,work,info = cgges(cselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(2*n,1),cselect_extra_args=(),overwrite_a=0,overwrite_b=0) a,b,sdim,alpha,beta,vsl,vsr,work,info = zgges(zselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(2*n,1),zselect_extra_args=(),overwrite_a=0,overwrite_b=0) alphar,alphai,beta,vl,vr,work,info = sggev(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alphar,alphai,beta,vl,vr,work,info = dggev(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,work,info = cggev(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,work,info = zggev(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) r,c,rowcnd,colcnd,amax,info = sgeequ(a) r,c,rowcnd,colcnd,amax,info = dgeequ(a) r,c,rowcnd,colcnd,amax,info = cgeequ(a) r,c,rowcnd,colcnd,amax,info = zgeequ(a) r,c,rowcnd,colcnd,amax,info = sgeequb(a) r,c,rowcnd,colcnd,amax,info = dgeequb(a) r,c,rowcnd,colcnd,amax,info = cgeequb(a) r,c,rowcnd,colcnd,amax,info = zgeequb(a) lub,piv,x,info = sgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lub,piv,x,info = dgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lub,piv,x,info = cgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lub,piv,x,info = zgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lu,ipiv,info = sgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) lu,ipiv,info = dgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) lu,ipiv,info = cgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) lu,ipiv,info = zgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) x,info = sgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) x,info = dgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) x,info = cgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) x,info = zgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) du2,d,du,x,info = sgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) du2,d,du,x,info = dgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) du2,d,du,x,info = cgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) du2,d,du,x,info = zgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) dl,d,du,du2,ipiv,info = sgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) dl,d,du,du2,ipiv,info = dgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) dl,d,du,du2,ipiv,info = cgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) dl,d,du,du2,ipiv,info = zgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) x,info = sgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) x,info = dgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) x,info = cgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) x,info = zgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = sgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = dgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = cgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = zgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) w,v,info = ssyev(a,compute_v=1,lower=0,lwork=max(3*n-1,1),overwrite_a=0) w,v,info = dsyev(a,compute_v=1,lower=0,lwork=max(3*n-1,1),overwrite_a=0) work,info = ssyev_lwork(n,lower=0) work,info = dsyev_lwork(n,lower=0) w,v,info = cheev(a,compute_v=1,lower=0,lwork=max(2*n-1,1),overwrite_a=0) w,v,info = zheev(a,compute_v=1,lower=0,lwork=max(2*n-1,1),overwrite_a=0) work,info = cheev_lwork(n,lower=0) work,info = zheev_lwork(n,lower=0) w,v,info = ssyevd(a,compute_v=1,lower=0,lwork=max((compute_v?1+6*n+2*n*n:2*n+1),1),liwork=(compute_v?3+5*n:1),overwrite_a=0) w,v,info = dsyevd(a,compute_v=1,lower=0,lwork=max((compute_v?1+6*n+2*n*n:2*n+1),1),liwork=(compute_v?3+5*n:1),overwrite_a=0) work,iwork,info = ssyevd_lwork(n,compute_v=1,lower=0) work,iwork,info = dsyevd_lwork(n,compute_v=1,lower=0) w,v,info = cheevd(a,compute_v=1,lower=0,lwork=max((compute_v?2*n+n*n:n+1),1),liwork=(compute_v?3+5*n:1),lrwork=(compute_v?1+5*n+2*n*n:n),overwrite_a=0) w,v,info = zheevd(a,compute_v=1,lower=0,lwork=max((compute_v?2*n+n*n:n+1),1),liwork=(compute_v?3+5*n:1),lrwork=(compute_v?1+5*n+2*n*n:n),overwrite_a=0) work,iwork,rwork,info = cheevd_lwork(n,compute_v=1,lower=0) work,iwork,rwork,info = zheevd_lwork(n,compute_v=1,lower=0) ldu,ipiv,info = ssytf2(a,lower=0,overwrite_a=0) ldu,ipiv,info = dsytf2(a,lower=0,overwrite_a=0) ldu,ipiv,info = csytf2(a,lower=0,overwrite_a=0) ldu,ipiv,info = zsytf2(a,lower=0,overwrite_a=0) c,info = ssygst(a,b,itype=1,lower=0,overwrite_a=0) c,info = dsygst(a,b,itype=1,lower=0,overwrite_a=0) ldu,ipiv,info = ssytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = dsytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = csytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = zsytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) work,info = ssytrf_lwork(n,lower=0) work,info = dsytrf_lwork(n,lower=0) work,info = csytrf_lwork(n,lower=0) work,info = zsytrf_lwork(n,lower=0) udut,ipiv,x,info = ssysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) udut,ipiv,x,info = dsysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) udut,ipiv,x,info = csysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) udut,ipiv,x,info = zsysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) work,info = ssysv_lwork(n,lower=0) work,info = dsysv_lwork(n,lower=0) work,info = csysv_lwork(n,lower=0) work,info = zsysv_lwork(n,lower=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = ssysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = dsysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = csysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = zsysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) work,info = ssysvx_lwork(n,lower=0) work,info = dsysvx_lwork(n,lower=0) work,info = csysvx_lwork(n,lower=0) work,info = zsysvx_lwork(n,lower=0) rcond,info = ssycon(a,ipiv,anorm,lower=0) rcond,info = dsycon(a,ipiv,anorm,lower=0) rcond,info = csycon(a,ipiv,anorm,lower=0) rcond,info = zsycon(a,ipiv,anorm,lower=0) rcond,info = checon(a,ipiv,anorm,lower=0) rcond,info = zhecon(a,ipiv,anorm,lower=0) a,e,info = ssyconv(a,ipiv,lower=0,way=0,overwrite_a=0) a,e,info = dsyconv(a,ipiv,lower=0,way=0,overwrite_a=0) a,e,info = csyconv(a,ipiv,lower=0,way=0,overwrite_a=0) a,e,info = zsyconv(a,ipiv,lower=0,way=0,overwrite_a=0) c,info = chegst(a,b,itype=1,lower=0,overwrite_a=0) c,info = zhegst(a,b,itype=1,lower=0,overwrite_a=0) ldu,ipiv,info = chetrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = zhetrf(a,lower=0,lwork=max(n,1),overwrite_a=0) work,info = chetrf_lwork(n,lower=0) work,info = zhetrf_lwork(n,lower=0) uduh,ipiv,x,info = chesv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) uduh,ipiv,x,info = zhesv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) work,info = chesv_lwork(n,lower=0) work,info = zhesv_lwork(n,lower=0) uduh,ipiv,x,rcond,ferr,berr,info = chesvx(a,b,af=,ipiv=,lwork=max(2*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) uduh,ipiv,x,rcond,ferr,berr,info = zhesvx(a,b,af=,ipiv=,lwork=max(2*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) work,info = chesvx_lwork(n,lower=0) work,info = zhesvx_lwork(n,lower=0) c,d,e,tau,info = ssytrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) c,d,e,tau,info = dsytrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) work,info = ssytrd_lwork(n,lower=0) work,info = dsytrd_lwork(n,lower=0) c,d,e,tau,info = chetrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) c,d,e,tau,info = zhetrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) work,info = chetrd_lwork(n,lower=0) work,info = zhetrd_lwork(n,lower=0) w,z,m,isuppz,info = ssyevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(26*n,1),liwork=max(1,10*n),overwrite_a=0) w,z,m,isuppz,info = dsyevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(26*n,1),liwork=max(1,10*n),overwrite_a=0) work,iwork,info = ssyevr_lwork(n,lower=0) work,iwork,info = dsyevr_lwork(n,lower=0) w,z,m,isuppz,info = cheevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),lrwork=max(24*n,1),liwork=max(1,10*n),overwrite_a=0) w,z,m,isuppz,info = zheevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),lrwork=max(24*n,1),liwork=max(1,10*n),overwrite_a=0) work,rwork,iwork,info = cheevr_lwork(n,lower=0) work,rwork,iwork,info = zheevr_lwork(n,lower=0) w,z,m,ifail,info = ssyevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0) w,z,m,ifail,info = dsyevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0) work,info = ssyevx_lwork(n,lower=0) work,info = dsyevx_lwork(n,lower=0) w,z,m,ifail,info = cheevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0) w,z,m,ifail,info = zheevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0) work,info = cheevx_lwork(n,lower=0) work,info = zheevx_lwork(n,lower=0) w,v,info = ssygv(a,b,itype=1,jobz='V',uplo='L',lwork=max(3*n-1,1),overwrite_a=0,overwrite_b=0) w,v,info = dsygv(a,b,itype=1,jobz='V',uplo='L',lwork=max(3*n-1,1),overwrite_a=0,overwrite_b=0) work,info = ssygv_lwork(n,uplo='L') work,info = dsygv_lwork(n,uplo='L') w,v,info = chegv(a,b,itype=1,jobz='V',uplo='L',lwork=max(2*n-1,1),overwrite_a=0,overwrite_b=0) w,v,info = zhegv(a,b,itype=1,jobz='V',uplo='L',lwork=max(2*n-1,1),overwrite_a=0,overwrite_b=0) work,info = chegv_lwork(n,uplo='L') work,info = zhegv_lwork(n,uplo='L') w,v,info = ssygvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?2*n+1:1+6*n+2*n*n),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,v,info = dsygvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?2*n+1:1+6*n+2*n*n),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,v,info = chegvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?n+1:n*(n+2)),lrwork=max((*jobz=='N'?n:2*n*n+5*n+1),1),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,v,info = zhegvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?n+1:n*(n+2)),lrwork=max((*jobz=='N'?n:2*n*n+5*n+1),1),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,z,m,ifail,info = ssygvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) w,z,m,ifail,info = dsygvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) work,info = ssygvx_lwork(n,uplo='L') work,info = dsygvx_lwork(n,uplo='L') w,z,m,ifail,info = chegvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) w,z,m,ifail,info = zhegvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) work,info = chegvx_lwork(n,uplo='L') work,info = zhegvx_lwork(n,uplo='L') s,scond,amax,info = ssyequb(a,lower=0) s,scond,amax,info = dsyequb(a,lower=0) s,scond,amax,info = csyequb(a,lower=0) s,scond,amax,info = zsyequb(a,lower=0) s,scond,amax,info = cheequb(a,lower=0) s,scond,amax,info = zheequb(a,lower=0) c,piv,rank_c,info = spstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = dpstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = cpstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = zpstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = spstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = dpstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = cpstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = zpstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,x,info = sposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) c,x,info = dposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) c,x,info = cposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) c,x,info = zposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = sposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = dposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = cposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = zposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) rcond,info = spocon(a,anorm,uplo='U') rcond,info = dpocon(a,anorm,uplo='U') rcond,info = cpocon(a,anorm,uplo='U') rcond,info = zpocon(a,anorm,uplo='U') c,info = spotrf(a,lower=0,clean=1,overwrite_a=0) c,info = dpotrf(a,lower=0,clean=1,overwrite_a=0) c,info = cpotrf(a,lower=0,clean=1,overwrite_a=0) c,info = zpotrf(a,lower=0,clean=1,overwrite_a=0) x,info = spotrs(c,b,lower=0,overwrite_b=0) x,info = dpotrs(c,b,lower=0,overwrite_b=0) x,info = cpotrs(c,b,lower=0,overwrite_b=0) x,info = zpotrs(c,b,lower=0,overwrite_b=0) inv_a,info = spotri(c,lower=0,overwrite_c=0) inv_a,info = dpotri(c,lower=0,overwrite_c=0) inv_a,info = cpotri(c,lower=0,overwrite_c=0) inv_a,info = zpotri(c,lower=0,overwrite_c=0) d,du,x,info = sptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,du,x,info = dptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,du,x,info = cptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,du,x,info = zptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,e,info = spttrf(d,e,overwrite_d=0,overwrite_e=0) d,e,info = dpttrf(d,e,overwrite_d=0,overwrite_e=0) d,e,info = cpttrf(d,e,overwrite_d=0,overwrite_e=0) d,e,info = zpttrf(d,e,overwrite_d=0,overwrite_e=0) x,info = spttrs(d,e,b,overwrite_b=0) x,info = dpttrs(d,e,b,overwrite_b=0) x,info = cpttrs(d,e,b,lower=0,overwrite_b=0) x,info = zpttrs(d,e,b,lower=0,overwrite_b=0) d,e,z,info = spteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) d,e,z,info = dpteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) d,e,z,info = cpteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) d,e,z,info = zpteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) df,ef,x,rcond,ferr,berr,info = sptsvx(d,e,b,fact='N',df=,ef=) df,ef,x,rcond,ferr,berr,info = dptsvx(d,e,b,fact='N',df=,ef=) df,ef,x,rcond,ferr,berr,info = cptsvx(d,e,b,fact='N',df=,ef=) df,ef,x,rcond,ferr,berr,info = zptsvx(d,e,b,fact='N',df=,ef=) sva,u,v,workout,iworkout,info = sgejsv(a,joba=4,jobu=0,jobv=0,jobr=1,jobt=0,jobp=1,lwork=max(6*n+2*n*n, max(2*m+n, max(4*n+n*n, max(2*n+n*n+6, 7)))),overwrite_a=0) sva,u,v,workout,iworkout,info = dgejsv(a,joba=4,jobu=0,jobv=0,jobr=1,jobt=0,jobp=1,lwork=max(6*n+2*n*n, max(2*m+n, max(4*n+n*n, max(2*n+n*n+6, 7)))),overwrite_a=0) a,b,q,z,work,info = stgexc(a,b,q,z,ifst,ilst,lwork=max(4*n+16,1),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) a,b,q,z,work,info = dtgexc(a,b,q,z,ifst,ilst,lwork=max(4*n+16,1),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) a,b,q,z,info = ctgexc(a,b,q,z,ifst,ilst,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) a,b,q,z,info = ztgexc(a,b,q,z,ifst,ilst,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) as,bs,alphar,alphai,beta,qs,zs,m,pl,pr,dif,info = stgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=4*n+16,liwork=n+6,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) as,bs,alphar,alphai,beta,qs,zs,m,pl,pr,dif,info = dtgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=4*n+16,liwork=n+6,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) work,iwork,info = stgsen_lwork(select,a,ijob=4) work,iwork,info = dtgsen_lwork(select,a,ijob=4) as,bs,alpha,beta,qs,zs,m,pl,pr,dif,info = ctgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=(ijob==0?1:n+2),liwork=(ijob==0?1:n+2),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) as,bs,alpha,beta,qs,zs,m,pl,pr,dif,info = ztgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=(ijob==0?1:n+2),liwork=(ijob==0?1:n+2),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) work,iwork,info = ctgsen_lwork(select,a,b,ijob=4) work,iwork,info = ztgsen_lwork(select,a,b,ijob=4) c,info = spbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) c,info = dpbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) c,info = cpbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) c,info = zpbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) x,info = spbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = dpbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = cpbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = zpbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = strtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = dtrtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = ctrtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = ztrtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = stbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) x,info = dtbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) x,info = ctbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) x,info = ztbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) c,x,info = spbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) c,x,info = dpbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) c,x,info = cpbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) c,x,info = zpbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = sorcsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2+2*m+5*MAX(1,q-1)+4*MAX(1,q)+8*q,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = dorcsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2+2*m+5*MAX(1,q-1)+4*MAX(1,q)+8*q,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) work,info = sorcsd_lwork(m,p,q) work,info = dorcsd_lwork(m,p,q) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = cuncsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2*m+MAX(1,MAX(mmp,mmq))+1,lrwork=5*MAX(1,q-1)+4*MAX(1,q)+8*q+1,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = zuncsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2*m+MAX(1,MAX(mmp,mmq))+1,lrwork=5*MAX(1,q-1)+4*MAX(1,q)+8*q+1,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) work,rwork,info = cuncsd_lwork(m,p,q) work,rwork,info = zuncsd_lwork(m,p,q) ht,info = sorghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) ht,info = dorghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) work,info = sorghr_lwork(n,lo=0,hi=n-1) work,info = dorghr_lwork(n,lo=0,hi=n-1) ht,info = cunghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) ht,info = zunghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) work,info = cunghr_lwork(n,lo=0,hi=n-1) work,info = zunghr_lwork(n,lo=0,hi=n-1) q,work,info = sorgqr(a,tau,lwork=max(3*n,1),overwrite_a=0) q,work,info = dorgqr(a,tau,lwork=max(3*n,1),overwrite_a=0) q,work,info = cungqr(a,tau,lwork=max(3*n,1),overwrite_a=0) q,work,info = zungqr(a,tau,lwork=max(3*n,1),overwrite_a=0) cq,work,info = sormqr(side,trans,a,tau,c,lwork,overwrite_c=0) cq,work,info = dormqr(side,trans,a,tau,c,lwork,overwrite_c=0) cq,work,info = cunmqr(side,trans,a,tau,c,lwork,overwrite_c=0) cq,work,info = zunmqr(side,trans,a,tau,c,lwork,overwrite_c=0) a,t,info = sgeqrt(nb,a,overwrite_a=0) a,t,info = dgeqrt(nb,a,overwrite_a=0) a,t,info = cgeqrt(nb,a,overwrite_a=0) a,t,info = zgeqrt(nb,a,overwrite_a=0) c,info = sgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) c,info = dgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) c,info = cgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) c,info = zgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) a,b,t,info = stpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,t,info = dtpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,t,info = ctpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,t,info = ztpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,info = stpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) a,b,info = dtpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) a,b,info = ctpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) a,b,info = ztpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) cq,info = sormrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) cq,info = dormrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) cq,info = cunmrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) cq,info = zunmrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) work,info = sormrz_lwork(m,n,side='L',trans='N') work,info = dormrz_lwork(m,n,side='L',trans='N') work,info = cunmrz_lwork(m,n,side='L',trans='N') work,info = zunmrz_lwork(m,n,side='L',trans='N') q,work,info = sorgrq(a,tau,lwork=max(3*m,1),overwrite_a=0) q,work,info = dorgrq(a,tau,lwork=max(3*m,1),overwrite_a=0) q,work,info = cungrq(a,tau,lwork=max(3*m,1),overwrite_a=0) q,work,info = zungrq(a,tau,lwork=max(3*m,1),overwrite_a=0) inv_c,info = strtri(c,lower=0,unitdiag=0,overwrite_c=0) inv_c,info = dtrtri(c,lower=0,unitdiag=0,overwrite_c=0) inv_c,info = ctrtri(c,lower=0,unitdiag=0,overwrite_c=0) inv_c,info = ztrtri(c,lower=0,unitdiag=0,overwrite_c=0) x,scale,info = strsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) x,scale,info = dtrsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) x,scale,info = ctrsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) x,scale,info = ztrsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) w,z,info = chbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),lrwork=(compute_v?1+5*n+2*n*n:n),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,info = zhbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),lrwork=(compute_v?1+5*n+2*n*n:n),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,m,ifail,info = chbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) w,z,m,ifail,info = zhbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) t,r,res,x,info = sgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) t,r,res,x,info = dgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) t,r,res,x,info = cgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) t,r,res,x,info = zgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) work,info = sgglse_lwork(m,n,p) work,info = dgglse_lwork(m,n,p) work,info = cgglse_lwork(m,n,p) work,info = zgglse_lwork(m,n,p) rcond,info = sppcon(n,ap,anorm,lower=0) rcond,info = dppcon(n,ap,anorm,lower=0) rcond,info = cppcon(n,ap,anorm,lower=0) rcond,info = zppcon(n,ap,anorm,lower=0) x,info = sppsv(n,ap,b,lower=0,overwrite_b=0) x,info = dppsv(n,ap,b,lower=0,overwrite_b=0) x,info = cppsv(n,ap,b,lower=0,overwrite_b=0) x,info = zppsv(n,ap,b,lower=0,overwrite_b=0) ul,info = spptrf(n,ap,lower=0,overwrite_ap=0) ul,info = dpptrf(n,ap,lower=0,overwrite_ap=0) ul,info = cpptrf(n,ap,lower=0,overwrite_ap=0) ul,info = zpptrf(n,ap,lower=0,overwrite_ap=0) uli,info = spptri(n,ap,lower=0,overwrite_ap=0) uli,info = dpptri(n,ap,lower=0,overwrite_ap=0) uli,info = cpptri(n,ap,lower=0,overwrite_ap=0) uli,info = zpptri(n,ap,lower=0,overwrite_ap=0) x,info = spptrs(n,ap,b,lower=0,overwrite_b=0) x,info = dpptrs(n,ap,b,lower=0,overwrite_b=0) x,info = cpptrs(n,ap,b,lower=0,overwrite_b=0) x,info = zpptrs(n,ap,b,lower=0,overwrite_b=0) w,z,info = ssbev(ab,compute_v=1,lower=0,ldab=shape(ab,0),overwrite_ab=1) w,z,info = dsbev(ab,compute_v=1,lower=0,ldab=shape(ab,0),overwrite_ab=1) w,z,info = ssbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,info = dsbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,m,ifail,info = ssbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) w,z,m,ifail,info = dsbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) m,w,iblock,isplit,info = sstebz(d,e,range,vl,vu,il,iu,tol,order) m,w,iblock,isplit,info = dstebz(d,e,range,vl,vu,il,iu,tol,order) vals,info = ssterf(d,e,overwrite_d=0,overwrite_e=0) vals,info = dsterf(d,e,overwrite_d=0,overwrite_e=0) z,info = sstein(d,e,w,iblock,isplit) z,info = dstein(d,e,w,iblock,isplit) m,w,z,info = sstemr(d,e,range,vl,vu,il,iu,compute_v=1,lwork=max((compute_v?18*n:12*n),1),liwork=(compute_v?10*n:8*n),overwrite_d=0) m,w,z,info = dstemr(d,e,range,vl,vu,il,iu,compute_v=1,lwork=max((compute_v?18*n:12*n),1),liwork=(compute_v?10*n:8*n),overwrite_d=0) work,iwork,info = sstemr_lwork(d,e,range,vl,vu,il,iu,compute_v=1,overwrite_d=0,overwrite_e=0) work,iwork,info = dstemr_lwork(d,e,range,vl,vu,il,iu,compute_v=1,overwrite_d=0,overwrite_e=0) vals,z,info = sstev(d,e,compute_v=1,overwrite_d=0,overwrite_e=0) vals,z,info = dstev(d,e,compute_v=1,overwrite_d=0,overwrite_e=0) cout = ssfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) cout = dsfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) cout = chfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) cout = zhfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) arf,info = stpttf(n,ap,transr='N',uplo='U') arf,info = dtpttf(n,ap,transr='N',uplo='U') arf,info = ctpttf(n,ap,transr='N',uplo='U') arf,info = ztpttf(n,ap,transr='N',uplo='U') a,info = stpttr(n,ap,uplo='U') a,info = dtpttr(n,ap,uplo='U') a,info = ctpttr(n,ap,uplo='U') a,info = ztpttr(n,ap,uplo='U') ap,info = stfttp(n,arf,transr='N',uplo='U') ap,info = dtfttp(n,arf,transr='N',uplo='U') ap,info = ctfttp(n,arf,transr='N',uplo='U') ap,info = ztfttp(n,arf,transr='N',uplo='U') a,info = stfttr(n,arf,transr='N',uplo='U') a,info = dtfttr(n,arf,transr='N',uplo='U') a,info = ctfttr(n,arf,transr='N',uplo='U') a,info = ztfttr(n,arf,transr='N',uplo='U') arf,info = strttf(a,transr='N',uplo='U') arf,info = dtrttf(a,transr='N',uplo='U') arf,info = ctrttf(a,transr='N',uplo='U') arf,info = ztrttf(a,transr='N',uplo='U') ap,info = strttp(a,uplo='U') ap,info = dtrttp(a,uplo='U') ap,info = ctrttp(a,uplo='U') ap,info = ztrttp(a,uplo='U') x = stfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) x = dtfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) x = ctfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) x = ztfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) achol,info = spftrf(n,a,transr='N',uplo='U',overwrite_a=0) achol,info = dpftrf(n,a,transr='N',uplo='U',overwrite_a=0) achol,info = cpftrf(n,a,transr='N',uplo='U',overwrite_a=0) achol,info = zpftrf(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = spftri(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = dpftri(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = cpftri(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = zpftri(n,a,transr='N',uplo='U',overwrite_a=0) x,info = spftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) x,info = dpftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) x,info = cpftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) x,info = zpftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) rz,tau,info = stzrzf(a,lwork=MAX(m,1),overwrite_a=0) rz,tau,info = dtzrzf(a,lwork=MAX(m,1),overwrite_a=0) rz,tau,info = ctzrzf(a,lwork=MAX(m,1),overwrite_a=0) rz,tau,info = ztzrzf(a,lwork=MAX(m,1),overwrite_a=0) work,info = stzrzf_lwork(m,n) work,info = dtzrzf_lwork(m,n) work,info = ctzrzf_lwork(m,n) work,info = ztzrzf_lwork(m,n) delta,sigma,work,info = slasd4(i,d,z,rho=1.0) delta,sigma,work,info = dlasd4(i,d,z,rho=1.0) a,info = slauum(c,lower=0,overwrite_c=0) a,info = dlauum(c,lower=0,overwrite_c=0) a,info = clauum(c,lower=0,overwrite_c=0) a,info = zlauum(c,lower=0,overwrite_c=0) a = slaswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) a = dlaswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) a = claswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) a = zlaswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) dlamch = dlamch(cmach) slamch = slamch(cmach) n2 = slange(norm,a) n2 = dlange(norm,a) n2 = clange(norm,a) n2 = zlange(norm,a) alpha,x,tau = slarfg(n,alpha,x,incx=1,overwrite_x=0) alpha,x,tau = dlarfg(n,alpha,x,incx=1,overwrite_x=0) alpha,x,tau = clarfg(n,alpha,x,incx=1,overwrite_x=0) alpha,x,tau = zlarfg(n,alpha,x,incx=1,overwrite_x=0) c = slarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) c = dlarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) c = clarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) c = zlarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) cs,sn,r = slartg(f,g) cs,sn,r = dlartg(f,g) cs,sn,r = clartg(f,g) cs,sn,r = zlartg(f,g) x,y = crot(x,y,c,s,n=(lx-1-offx)/abs(incx)+1,offx=0,incx=1,offy=0,incy=1,overwrite_x=0,overwrite_y=0) x,y = zrot(x,y,c,s,n=(lx-1-offx)/abs(incx)+1,offx=0,incx=1,offy=0,incy=1,overwrite_x=0,overwrite_y=0) major,minor,patch = ilaver() .
This module '_flapack' is auto-generated with f2py (version:2). Functions: ba,lo,hi,pivscale,info = sgebal(a,scale=0,permute=0,overwrite_a=0) ba,lo,hi,pivscale,info = dgebal(a,scale=0,permute=0,overwrite_a=0) ba,lo,hi,pivscale,info = cgebal(a,scale=0,permute=0,overwrite_a=0) ba,lo,hi,pivscale,info = zgebal(a,scale=0,permute=0,overwrite_a=0) ht,tau,info = sgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) ht,tau,info = dgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) ht,tau,info = cgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) ht,tau,info = zgehrd(a,lo=0,hi=n-1,lwork=MAX(n,1),overwrite_a=0) work,info = sgehrd_lwork(n,lo=0,hi=n-1) work,info = dgehrd_lwork(n,lo=0,hi=n-1) work,info = cgehrd_lwork(n,lo=0,hi=n-1) work,info = zgehrd_lwork(n,lo=0,hi=n-1) lu,piv,x,info = sgesv(a,b,overwrite_a=0,overwrite_b=0) lu,piv,x,info = dgesv(a,b,overwrite_a=0,overwrite_b=0) lu,piv,x,info = cgesv(a,b,overwrite_a=0,overwrite_b=0) lu,piv,x,info = zgesv(a,b,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = sgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = dgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = cgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) as,lu,ipiv,equed,rs,cs,bs,x,rcond,ferr,berr,info = zgesvx(a,b,fact='E',trans='N',af=,ipiv=,equed='B',r=,c=,overwrite_a=0,overwrite_b=0) rcond,info = sgecon(a,anorm,norm='1') rcond,info = dgecon(a,anorm,norm='1') rcond,info = cgecon(a,anorm,norm='1') rcond,info = zgecon(a,anorm,norm='1') lu,piv,info = sgetrf(a,overwrite_a=0) lu,piv,info = dgetrf(a,overwrite_a=0) lu,piv,info = cgetrf(a,overwrite_a=0) lu,piv,info = zgetrf(a,overwrite_a=0) x,info = sgetrs(lu,piv,b,trans=0,overwrite_b=0) x,info = dgetrs(lu,piv,b,trans=0,overwrite_b=0) x,info = cgetrs(lu,piv,b,trans=0,overwrite_b=0) x,info = zgetrs(lu,piv,b,trans=0,overwrite_b=0) lu,ipiv,jpiv,info = sgetc2(a,overwrite_a=0) lu,ipiv,jpiv,info = dgetc2(a,overwrite_a=0) lu,ipiv,jpiv,info = cgetc2(a,overwrite_a=0) lu,ipiv,jpiv,info = zgetc2(a,overwrite_a=0) x,scale = sgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) x,scale = dgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) x,scale = cgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) x,scale = zgesc2(lu,rhs,ipiv,jpiv,overwrite_rhs=0) inv_a,info = sgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) inv_a,info = dgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) inv_a,info = cgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) inv_a,info = zgetri(lu,piv,lwork=max(3*n,1),overwrite_lu=0) work,info = sgetri_lwork(n) work,info = dgetri_lwork(n) work,info = cgetri_lwork(n) work,info = zgetri_lwork(n) u,s,vt,info = sgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?4*minmn*minmn+MAX(m,n)+9*minmn:MAX(14*minmn+4,10*minmn+2+25*(25+8))+MAX(m,n)),1),overwrite_a=0) u,s,vt,info = dgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?4*minmn*minmn+MAX(m,n)+9*minmn:MAX(14*minmn+4,10*minmn+2+25*(25+8))+MAX(m,n)),1),overwrite_a=0) work,info = sgesdd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = dgesdd_lwork(m,n,compute_uv=1,full_matrices=1) u,s,vt,info = cgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?2*minmn*minmn+MAX(m,n)+2*minmn:2*minmn+MAX(m,n)),1),overwrite_a=0) u,s,vt,info = zgesdd(a,compute_uv=1,full_matrices=1,lwork=max((compute_uv?2*minmn*minmn+MAX(m,n)+2*minmn:2*minmn+MAX(m,n)),1),overwrite_a=0) work,info = cgesdd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = zgesdd_lwork(m,n,compute_uv=1,full_matrices=1) u,s,vt,info = sgesvd(a,compute_uv=1,full_matrices=1,lwork=max(MAX(3*minmn+MAX(m,n),5*minmn),1),overwrite_a=0) u,s,vt,info = dgesvd(a,compute_uv=1,full_matrices=1,lwork=max(MAX(3*minmn+MAX(m,n),5*minmn),1),overwrite_a=0) work,info = sgesvd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = dgesvd_lwork(m,n,compute_uv=1,full_matrices=1) u,s,vt,info = cgesvd(a,compute_uv=1,full_matrices=1,lwork=MAX(2*minmn+MAX(m,n),1),overwrite_a=0) u,s,vt,info = zgesvd(a,compute_uv=1,full_matrices=1,lwork=MAX(2*minmn+MAX(m,n),1),overwrite_a=0) work,info = cgesvd_lwork(m,n,compute_uv=1,full_matrices=1) work,info = zgesvd_lwork(m,n,compute_uv=1,full_matrices=1) lqr,x,info = sgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) lqr,x,info = dgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) lqr,x,info = cgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) lqr,x,info = zgels(a,b,trans='N',lwork=MAX(MIN(m,n)+MAX(MIN(m,n),nrhs),1),overwrite_a=0,overwrite_b=0) work,info = sgels_lwork(m,n,nrhs,trans='N') work,info = dgels_lwork(m,n,nrhs,trans='N') work,info = cgels_lwork(m,n,nrhs,trans='N') work,info = zgels_lwork(m,n,nrhs,trans='N') v,x,s,rank,work,info = sgelss(a,b,cond=-1.0,lwork=max(3*minmn+MAX(2*minmn,MAX(maxmn,nrhs)),1),overwrite_a=0,overwrite_b=0) v,x,s,rank,work,info = dgelss(a,b,cond=-1.0,lwork=max(3*minmn+MAX(2*minmn,MAX(maxmn,nrhs)),1),overwrite_a=0,overwrite_b=0) work,info = sgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,info = dgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) v,x,s,rank,work,info = cgelss(a,b,cond=-1.0,lwork=max(2*minmn+MAX(maxmn,nrhs),1),overwrite_a=0,overwrite_b=0) v,x,s,rank,work,info = zgelss(a,b,cond=-1.0,lwork=max(2*minmn+MAX(maxmn,nrhs),1),overwrite_a=0,overwrite_b=0) work,info = cgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,info = zgelss_lwork(m,n,nrhs,cond=-1.0,lwork=-1) v,x,j,rank,info = sgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) v,x,j,rank,info = dgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) work,info = sgelsy_lwork(m,n,nrhs,cond,lwork=-1) work,info = dgelsy_lwork(m,n,nrhs,cond,lwork=-1) v,x,j,rank,info = cgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) v,x,j,rank,info = zgelsy(a,b,jptv,cond,lwork,overwrite_a=0,overwrite_b=0) work,info = cgelsy_lwork(m,n,nrhs,cond,lwork=-1) work,info = zgelsy_lwork(m,n,nrhs,cond,lwork=-1) x,s,rank,info = sgelsd(a,b,lwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) x,s,rank,info = dgelsd(a,b,lwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) work,iwork,info = sgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,iwork,info = dgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) x,s,rank,info = cgelsd(a,b,lwork,size_rwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) x,s,rank,info = zgelsd(a,b,lwork,size_rwork,size_iwork,cond=-1.0,overwrite_a=0,overwrite_b=0) work,rwork,iwork,info = cgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) work,rwork,iwork,info = zgelsd_lwork(m,n,nrhs,cond=-1.0,lwork=-1) qr,jpvt,tau,work,info = sgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,jpvt,tau,work,info = dgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,jpvt,tau,work,info = cgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,jpvt,tau,work,info = zgeqp3(a,lwork=max(3*(n+1),1),overwrite_a=0) qr,tau,work,info = sgeqrf(a,lwork=max(3*n,1),overwrite_a=0) qr,tau,work,info = dgeqrf(a,lwork=max(3*n,1),overwrite_a=0) qr,tau,work,info = cgeqrf(a,lwork=max(3*n,1),overwrite_a=0) qr,tau,work,info = zgeqrf(a,lwork=max(3*n,1),overwrite_a=0) work,info = sgeqrf_lwork(m,n) work,info = dgeqrf_lwork(m,n) work,info = cgeqrf_lwork(m,n) work,info = zgeqrf_lwork(m,n) qr,tau,info = sgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) qr,tau,info = dgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) qr,tau,info = cgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) qr,tau,info = zgeqrfp(a,lwork=MAX(1, n),overwrite_a=0) work,info = sgeqrfp_lwork(m,n) work,info = dgeqrfp_lwork(m,n) work,info = cgeqrfp_lwork(m,n) work,info = zgeqrfp_lwork(m,n) qr,tau,work,info = sgerqf(a,lwork=max(3*m,1),overwrite_a=0) qr,tau,work,info = dgerqf(a,lwork=max(3*m,1),overwrite_a=0) qr,tau,work,info = cgerqf(a,lwork=max(3*m,1),overwrite_a=0) qr,tau,work,info = zgerqf(a,lwork=max(3*m,1),overwrite_a=0) wr,wi,vl,vr,info = sgeev(a,compute_vl=1,compute_vr=1,lwork=max(4*n,1),overwrite_a=0) wr,wi,vl,vr,info = dgeev(a,compute_vl=1,compute_vr=1,lwork=max(4*n,1),overwrite_a=0) work,info = sgeev_lwork(n,compute_vl=1,compute_vr=1) work,info = dgeev_lwork(n,compute_vl=1,compute_vr=1) w,vl,vr,info = cgeev(a,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0) w,vl,vr,info = zgeev(a,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0) work,info = cgeev_lwork(n,compute_vl=1,compute_vr=1) work,info = zgeev_lwork(n,compute_vl=1,compute_vr=1) alphar,alphai,beta,vl,vr,info = sgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alphar,alphai,beta,vl,vr,info = dgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,info = cgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,info = zgegv(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) t,sdim,w,vs,work,info = cgees(cselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),cselect_extra_args=(),overwrite_a=0) t,sdim,w,vs,work,info = zgees(zselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),zselect_extra_args=(),overwrite_a=0) t,sdim,wr,wi,vs,work,info = sgees(sselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),sselect_extra_args=(),overwrite_a=0) t,sdim,wr,wi,vs,work,info = dgees(dselect,a,compute_v=1,sort_t=0,lwork=max(3*n,1),dselect_extra_args=(),overwrite_a=0) a,b,sdim,alphar,alphai,beta,vsl,vsr,work,info = sgges(sselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(8*n+16,1),sselect_extra_args=(),overwrite_a=0,overwrite_b=0) a,b,sdim,alphar,alphai,beta,vsl,vsr,work,info = dgges(dselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(8*n+16,1),dselect_extra_args=(),overwrite_a=0,overwrite_b=0) a,b,sdim,alpha,beta,vsl,vsr,work,info = cgges(cselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(2*n,1),cselect_extra_args=(),overwrite_a=0,overwrite_b=0) a,b,sdim,alpha,beta,vsl,vsr,work,info = zgges(zselect,a,b,jobvsl=1,jobvsr=1,sort_t=0,ldvsl=((jobvsl==1)?n:1),ldvsr=((jobvsr==1)?n:1),lwork=max(2*n,1),zselect_extra_args=(),overwrite_a=0,overwrite_b=0) alphar,alphai,beta,vl,vr,work,info = sggev(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alphar,alphai,beta,vl,vr,work,info = dggev(a,b,compute_vl=1,compute_vr=1,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,work,info = cggev(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) alpha,beta,vl,vr,work,info = zggev(a,b,compute_vl=1,compute_vr=1,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) r,c,rowcnd,colcnd,amax,info = sgeequ(a) r,c,rowcnd,colcnd,amax,info = dgeequ(a) r,c,rowcnd,colcnd,amax,info = cgeequ(a) r,c,rowcnd,colcnd,amax,info = zgeequ(a) r,c,rowcnd,colcnd,amax,info = sgeequb(a) r,c,rowcnd,colcnd,amax,info = dgeequb(a) r,c,rowcnd,colcnd,amax,info = cgeequb(a) r,c,rowcnd,colcnd,amax,info = zgeequb(a) lub,piv,x,info = sgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lub,piv,x,info = dgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lub,piv,x,info = cgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lub,piv,x,info = zgbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0) lu,ipiv,info = sgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) lu,ipiv,info = dgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) lu,ipiv,info = cgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) lu,ipiv,info = zgbtrf(ab,kl,ku,m=shape(ab,1),n=shape(ab,1),ldab=max(shape(ab,0),1),overwrite_ab=0) x,info = sgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) x,info = dgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) x,info = cgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) x,info = zgbtrs(ab,kl,ku,b,ipiv,trans=0,n=shape(ab,1),ldab=shape(ab,0),ldb=shape(b,0),overwrite_b=0) du2,d,du,x,info = sgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) du2,d,du,x,info = dgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) du2,d,du,x,info = cgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) du2,d,du,x,info = zgtsv(dl,d,du,b,overwrite_dl=0,overwrite_d=0,overwrite_du=0,overwrite_b=0) dl,d,du,du2,ipiv,info = sgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) dl,d,du,du2,ipiv,info = dgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) dl,d,du,du2,ipiv,info = cgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) dl,d,du,du2,ipiv,info = zgttrf(dl,d,du,overwrite_dl=0,overwrite_d=0,overwrite_du=0) x,info = sgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) x,info = dgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) x,info = cgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) x,info = zgttrs(dl,d,du,du2,ipiv,b,trans='N',overwrite_b=0) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = sgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = dgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = cgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) dlf,df,duf,du2,ipiv,x,rcond,ferr,berr,info = zgtsvx(dl,d,du,b,fact='N',trans='N',dlf=,df=,duf=,du2=,ipiv=) w,v,info = ssyev(a,compute_v=1,lower=0,lwork=max(3*n-1,1),overwrite_a=0) w,v,info = dsyev(a,compute_v=1,lower=0,lwork=max(3*n-1,1),overwrite_a=0) work,info = ssyev_lwork(n,lower=0) work,info = dsyev_lwork(n,lower=0) w,v,info = cheev(a,compute_v=1,lower=0,lwork=max(2*n-1,1),overwrite_a=0) w,v,info = zheev(a,compute_v=1,lower=0,lwork=max(2*n-1,1),overwrite_a=0) work,info = cheev_lwork(n,lower=0) work,info = zheev_lwork(n,lower=0) w,v,info = ssyevd(a,compute_v=1,lower=0,lwork=max((compute_v?1+6*n+2*n*n:2*n+1),1),liwork=(compute_v?3+5*n:1),overwrite_a=0) w,v,info = dsyevd(a,compute_v=1,lower=0,lwork=max((compute_v?1+6*n+2*n*n:2*n+1),1),liwork=(compute_v?3+5*n:1),overwrite_a=0) work,iwork,info = ssyevd_lwork(n,compute_v=1,lower=0) work,iwork,info = dsyevd_lwork(n,compute_v=1,lower=0) w,v,info = cheevd(a,compute_v=1,lower=0,lwork=max((compute_v?2*n+n*n:n+1),1),liwork=(compute_v?3+5*n:1),lrwork=(compute_v?1+5*n+2*n*n:n),overwrite_a=0) w,v,info = zheevd(a,compute_v=1,lower=0,lwork=max((compute_v?2*n+n*n:n+1),1),liwork=(compute_v?3+5*n:1),lrwork=(compute_v?1+5*n+2*n*n:n),overwrite_a=0) work,iwork,rwork,info = cheevd_lwork(n,compute_v=1,lower=0) work,iwork,rwork,info = zheevd_lwork(n,compute_v=1,lower=0) ldu,ipiv,info = ssytf2(a,lower=0,overwrite_a=0) ldu,ipiv,info = dsytf2(a,lower=0,overwrite_a=0) ldu,ipiv,info = csytf2(a,lower=0,overwrite_a=0) ldu,ipiv,info = zsytf2(a,lower=0,overwrite_a=0) c,info = ssygst(a,b,itype=1,lower=0,overwrite_a=0) c,info = dsygst(a,b,itype=1,lower=0,overwrite_a=0) ldu,ipiv,info = ssytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = dsytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = csytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = zsytrf(a,lower=0,lwork=max(n,1),overwrite_a=0) work,info = ssytrf_lwork(n,lower=0) work,info = dsytrf_lwork(n,lower=0) work,info = csytrf_lwork(n,lower=0) work,info = zsytrf_lwork(n,lower=0) udut,ipiv,x,info = ssysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) udut,ipiv,x,info = dsysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) udut,ipiv,x,info = csysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) udut,ipiv,x,info = zsysv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) work,info = ssysv_lwork(n,lower=0) work,info = dsysv_lwork(n,lower=0) work,info = csysv_lwork(n,lower=0) work,info = zsysv_lwork(n,lower=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = ssysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = dsysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = csysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) a_s,udut,ipiv,b_s,x,rcond,ferr,berr,info = zsysvx(a,b,af=,ipiv=,lwork=max(3*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) work,info = ssysvx_lwork(n,lower=0) work,info = dsysvx_lwork(n,lower=0) work,info = csysvx_lwork(n,lower=0) work,info = zsysvx_lwork(n,lower=0) rcond,info = ssycon(a,ipiv,anorm,lower=0) rcond,info = dsycon(a,ipiv,anorm,lower=0) rcond,info = csycon(a,ipiv,anorm,lower=0) rcond,info = zsycon(a,ipiv,anorm,lower=0) rcond,info = checon(a,ipiv,anorm,lower=0) rcond,info = zhecon(a,ipiv,anorm,lower=0) a,e,info = ssyconv(a,ipiv,lower=0,way=0,overwrite_a=0) a,e,info = dsyconv(a,ipiv,lower=0,way=0,overwrite_a=0) a,e,info = csyconv(a,ipiv,lower=0,way=0,overwrite_a=0) a,e,info = zsyconv(a,ipiv,lower=0,way=0,overwrite_a=0) c,info = chegst(a,b,itype=1,lower=0,overwrite_a=0) c,info = zhegst(a,b,itype=1,lower=0,overwrite_a=0) ldu,ipiv,info = chetrf(a,lower=0,lwork=max(n,1),overwrite_a=0) ldu,ipiv,info = zhetrf(a,lower=0,lwork=max(n,1),overwrite_a=0) work,info = chetrf_lwork(n,lower=0) work,info = zhetrf_lwork(n,lower=0) uduh,ipiv,x,info = chesv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) uduh,ipiv,x,info = zhesv(a,b,lwork=max(n,1),lower=0,overwrite_a=0,overwrite_b=0) work,info = chesv_lwork(n,lower=0) work,info = zhesv_lwork(n,lower=0) uduh,ipiv,x,rcond,ferr,berr,info = chesvx(a,b,af=,ipiv=,lwork=max(2*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) uduh,ipiv,x,rcond,ferr,berr,info = zhesvx(a,b,af=,ipiv=,lwork=max(2*n,1),factored=0,lower=0,overwrite_a=0,overwrite_b=0) work,info = chesvx_lwork(n,lower=0) work,info = zhesvx_lwork(n,lower=0) c,d,e,tau,info = ssytrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) c,d,e,tau,info = dsytrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) work,info = ssytrd_lwork(n,lower=0) work,info = dsytrd_lwork(n,lower=0) c,d,e,tau,info = chetrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) c,d,e,tau,info = zhetrd(a,lower=0,lwork=MAX(n,1),overwrite_a=0) work,info = chetrd_lwork(n,lower=0) work,info = zhetrd_lwork(n,lower=0) w,z,m,isuppz,info = ssyevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(26*n,1),liwork=max(1,10*n),overwrite_a=0) w,z,m,isuppz,info = dsyevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(26*n,1),liwork=max(1,10*n),overwrite_a=0) work,iwork,info = ssyevr_lwork(n,lower=0) work,iwork,info = dsyevr_lwork(n,lower=0) w,z,m,isuppz,info = cheevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),lrwork=max(24*n,1),liwork=max(1,10*n),overwrite_a=0) w,z,m,isuppz,info = zheevr(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),lrwork=max(24*n,1),liwork=max(1,10*n),overwrite_a=0) work,rwork,iwork,info = cheevr_lwork(n,lower=0) work,rwork,iwork,info = zheevr_lwork(n,lower=0) w,z,m,ifail,info = ssyevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0) w,z,m,ifail,info = dsyevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0) work,info = ssyevx_lwork(n,lower=0) work,info = dsyevx_lwork(n,lower=0) w,z,m,ifail,info = cheevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0) w,z,m,ifail,info = zheevx(a,compute_v=1,range='A',lower=0,vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0) work,info = cheevx_lwork(n,lower=0) work,info = zheevx_lwork(n,lower=0) w,v,info = ssygv(a,b,itype=1,jobz='V',uplo='L',lwork=max(3*n-1,1),overwrite_a=0,overwrite_b=0) w,v,info = dsygv(a,b,itype=1,jobz='V',uplo='L',lwork=max(3*n-1,1),overwrite_a=0,overwrite_b=0) work,info = ssygv_lwork(n,uplo='L') work,info = dsygv_lwork(n,uplo='L') w,v,info = chegv(a,b,itype=1,jobz='V',uplo='L',lwork=max(2*n-1,1),overwrite_a=0,overwrite_b=0) w,v,info = zhegv(a,b,itype=1,jobz='V',uplo='L',lwork=max(2*n-1,1),overwrite_a=0,overwrite_b=0) work,info = chegv_lwork(n,uplo='L') work,info = zhegv_lwork(n,uplo='L') w,v,info = ssygvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?2*n+1:1+6*n+2*n*n),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,v,info = dsygvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?2*n+1:1+6*n+2*n*n),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,v,info = chegvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?n+1:n*(n+2)),lrwork=max((*jobz=='N'?n:2*n*n+5*n+1),1),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,v,info = zhegvd(a,b,itype=1,jobz='V',uplo='L',lwork=(*jobz=='N'?n+1:n*(n+2)),lrwork=max((*jobz=='N'?n:2*n*n+5*n+1),1),liwork=(*jobz=='N'?1:5*n+3),overwrite_a=0,overwrite_b=0) w,z,m,ifail,info = ssygvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) w,z,m,ifail,info = dsygvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(8*n,1),overwrite_a=0,overwrite_b=0) work,info = ssygvx_lwork(n,uplo='L') work,info = dsygvx_lwork(n,uplo='L') w,z,m,ifail,info = chegvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) w,z,m,ifail,info = zhegvx(a,b,itype=1,jobz='V',range='A',uplo='L',vl=0.0,vu=1.0,il=1,iu=n,abstol=0.0,lwork=max(2*n,1),overwrite_a=0,overwrite_b=0) work,info = chegvx_lwork(n,uplo='L') work,info = zhegvx_lwork(n,uplo='L') s,scond,amax,info = ssyequb(a,lower=0) s,scond,amax,info = dsyequb(a,lower=0) s,scond,amax,info = csyequb(a,lower=0) s,scond,amax,info = zsyequb(a,lower=0) s,scond,amax,info = cheequb(a,lower=0) s,scond,amax,info = zheequb(a,lower=0) c,piv,rank_c,info = spstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = dpstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = cpstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = zpstrf(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = spstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = dpstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = cpstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,piv,rank_c,info = zpstf2(a,tol=-1.0,lower=0,overwrite_a=0) c,x,info = sposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) c,x,info = dposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) c,x,info = cposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) c,x,info = zposv(a,b,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = sposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = dposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = cposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) a_s,lu,equed,s,b_s,x,rcond,ferr,berr,info = zposvx(a,b,fact='E',af=,equed='Y',s=,lower=0,overwrite_a=0,overwrite_b=0) rcond,info = spocon(a,anorm,uplo='U') rcond,info = dpocon(a,anorm,uplo='U') rcond,info = cpocon(a,anorm,uplo='U') rcond,info = zpocon(a,anorm,uplo='U') c,info = spotrf(a,lower=0,clean=1,overwrite_a=0) c,info = dpotrf(a,lower=0,clean=1,overwrite_a=0) c,info = cpotrf(a,lower=0,clean=1,overwrite_a=0) c,info = zpotrf(a,lower=0,clean=1,overwrite_a=0) x,info = spotrs(c,b,lower=0,overwrite_b=0) x,info = dpotrs(c,b,lower=0,overwrite_b=0) x,info = cpotrs(c,b,lower=0,overwrite_b=0) x,info = zpotrs(c,b,lower=0,overwrite_b=0) inv_a,info = spotri(c,lower=0,overwrite_c=0) inv_a,info = dpotri(c,lower=0,overwrite_c=0) inv_a,info = cpotri(c,lower=0,overwrite_c=0) inv_a,info = zpotri(c,lower=0,overwrite_c=0) d,du,x,info = sptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,du,x,info = dptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,du,x,info = cptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,du,x,info = zptsv(d,e,b,overwrite_d=0,overwrite_e=0,overwrite_b=0) d,e,info = spttrf(d,e,overwrite_d=0,overwrite_e=0) d,e,info = dpttrf(d,e,overwrite_d=0,overwrite_e=0) d,e,info = cpttrf(d,e,overwrite_d=0,overwrite_e=0) d,e,info = zpttrf(d,e,overwrite_d=0,overwrite_e=0) x,info = spttrs(d,e,b,overwrite_b=0) x,info = dpttrs(d,e,b,overwrite_b=0) x,info = cpttrs(d,e,b,lower=0,overwrite_b=0) x,info = zpttrs(d,e,b,lower=0,overwrite_b=0) d,e,z,info = spteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) d,e,z,info = dpteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) d,e,z,info = cpteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) d,e,z,info = zpteqr(d,e,z,compute_z=0,overwrite_d=0,overwrite_e=0,overwrite_z=0) df,ef,x,rcond,ferr,berr,info = sptsvx(d,e,b,fact='N',df=,ef=) df,ef,x,rcond,ferr,berr,info = dptsvx(d,e,b,fact='N',df=,ef=) df,ef,x,rcond,ferr,berr,info = cptsvx(d,e,b,fact='N',df=,ef=) df,ef,x,rcond,ferr,berr,info = zptsvx(d,e,b,fact='N',df=,ef=) sva,u,v,workout,iworkout,info = sgejsv(a,joba=4,jobu=0,jobv=0,jobr=1,jobt=0,jobp=1,lwork=max(6*n+2*n*n, max(2*m+n, max(4*n+n*n, max(2*n+n*n+6, 7)))),overwrite_a=0) sva,u,v,workout,iworkout,info = dgejsv(a,joba=4,jobu=0,jobv=0,jobr=1,jobt=0,jobp=1,lwork=max(6*n+2*n*n, max(2*m+n, max(4*n+n*n, max(2*n+n*n+6, 7)))),overwrite_a=0) a,b,q,z,work,info = stgexc(a,b,q,z,ifst,ilst,lwork=max(4*n+16,1),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) a,b,q,z,work,info = dtgexc(a,b,q,z,ifst,ilst,lwork=max(4*n+16,1),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) a,b,q,z,info = ctgexc(a,b,q,z,ifst,ilst,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) a,b,q,z,info = ztgexc(a,b,q,z,ifst,ilst,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) as,bs,alphar,alphai,beta,qs,zs,m,pl,pr,dif,info = stgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=4*n+16,liwork=n+6,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) as,bs,alphar,alphai,beta,qs,zs,m,pl,pr,dif,info = dtgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=4*n+16,liwork=n+6,overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) work,iwork,info = stgsen_lwork(select,a,ijob=4) work,iwork,info = dtgsen_lwork(select,a,ijob=4) as,bs,alpha,beta,qs,zs,m,pl,pr,dif,info = ctgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=(ijob==0?1:n+2),liwork=(ijob==0?1:n+2),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) as,bs,alpha,beta,qs,zs,m,pl,pr,dif,info = ztgsen(select,a,b,q,z,ijob=4,wantq=1,wantz=1,lwork=(ijob==0?1:n+2),liwork=(ijob==0?1:n+2),overwrite_a=0,overwrite_b=0,overwrite_q=0,overwrite_z=0) work,iwork,info = ctgsen_lwork(select,a,b,ijob=4) work,iwork,info = ztgsen_lwork(select,a,b,ijob=4) c,info = spbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) c,info = dpbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) c,info = cpbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) c,info = zpbtrf(ab,lower=0,ldab=shape(ab,0),overwrite_ab=0) x,info = spbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = dpbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = cpbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = zpbtrs(ab,b,lower=0,ldab=shape(ab,0),overwrite_b=0) x,info = strtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = dtrtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = ctrtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = ztrtrs(a,b,lower=0,trans=0,unitdiag=0,lda=shape(a,0),overwrite_b=0) x,info = stbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) x,info = dtbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) x,info = ctbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) x,info = ztbtrs(ab,b,uplo='U',trans='N',diag='N',overwrite_b=0) c,x,info = spbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) c,x,info = dpbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) c,x,info = cpbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) c,x,info = zpbsv(ab,b,lower=0,ldab=shape(ab,0),overwrite_ab=0,overwrite_b=0) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = sorcsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2+2*m+5*MAX(1,q-1)+4*MAX(1,q)+8*q,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = dorcsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2+2*m+5*MAX(1,q-1)+4*MAX(1,q)+8*q,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) work,info = sorcsd_lwork(m,p,q) work,info = dorcsd_lwork(m,p,q) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = cuncsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2*m+MAX(1,MAX(mmp,mmq))+1,lrwork=5*MAX(1,q-1)+4*MAX(1,q)+8*q+1,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) cs11,cs12,cs21,cs22,theta,u1,u2,v1t,v2t,info = zuncsd(x11,x12,x21,x22,compute_u1=1,compute_u2=1,compute_v1t=1,compute_v2t=1,trans=0,signs=0,lwork=2*m+MAX(1,MAX(mmp,mmq))+1,lrwork=5*MAX(1,q-1)+4*MAX(1,q)+8*q+1,overwrite_x11=0,overwrite_x12=0,overwrite_x21=0,overwrite_x22=0) work,rwork,info = cuncsd_lwork(m,p,q) work,rwork,info = zuncsd_lwork(m,p,q) ht,info = sorghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) ht,info = dorghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) work,info = sorghr_lwork(n,lo=0,hi=n-1) work,info = dorghr_lwork(n,lo=0,hi=n-1) ht,info = cunghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) ht,info = zunghr(a,tau,lo=0,hi=n-1,lwork=max(hi-lo,1),overwrite_a=0) work,info = cunghr_lwork(n,lo=0,hi=n-1) work,info = zunghr_lwork(n,lo=0,hi=n-1) q,work,info = sorgqr(a,tau,lwork=max(3*n,1),overwrite_a=0) q,work,info = dorgqr(a,tau,lwork=max(3*n,1),overwrite_a=0) q,work,info = cungqr(a,tau,lwork=max(3*n,1),overwrite_a=0) q,work,info = zungqr(a,tau,lwork=max(3*n,1),overwrite_a=0) cq,work,info = sormqr(side,trans,a,tau,c,lwork,overwrite_c=0) cq,work,info = dormqr(side,trans,a,tau,c,lwork,overwrite_c=0) cq,work,info = cunmqr(side,trans,a,tau,c,lwork,overwrite_c=0) cq,work,info = zunmqr(side,trans,a,tau,c,lwork,overwrite_c=0) a,t,info = sgeqrt(nb,a,overwrite_a=0) a,t,info = dgeqrt(nb,a,overwrite_a=0) a,t,info = cgeqrt(nb,a,overwrite_a=0) a,t,info = zgeqrt(nb,a,overwrite_a=0) c,info = sgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) c,info = dgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) c,info = cgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) c,info = zgemqrt(v,t,c,side='L',trans='N',overwrite_c=0) a,b,t,info = stpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,t,info = dtpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,t,info = ctpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,t,info = ztpqrt(l,nb,a,b,overwrite_a=0,overwrite_b=0) a,b,info = stpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) a,b,info = dtpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) a,b,info = ctpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) a,b,info = ztpmqrt(l,v,t,a,b,side='L',trans='N',overwrite_a=0,overwrite_b=0) cq,info = sormrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) cq,info = dormrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) cq,info = cunmrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) cq,info = zunmrz(a,tau,c,side='L',trans='N',lwork=MAX((side[0]=='L'?n:m),1),overwrite_c=0) work,info = sormrz_lwork(m,n,side='L',trans='N') work,info = dormrz_lwork(m,n,side='L',trans='N') work,info = cunmrz_lwork(m,n,side='L',trans='N') work,info = zunmrz_lwork(m,n,side='L',trans='N') q,work,info = sorgrq(a,tau,lwork=max(3*m,1),overwrite_a=0) q,work,info = dorgrq(a,tau,lwork=max(3*m,1),overwrite_a=0) q,work,info = cungrq(a,tau,lwork=max(3*m,1),overwrite_a=0) q,work,info = zungrq(a,tau,lwork=max(3*m,1),overwrite_a=0) inv_c,info = strtri(c,lower=0,unitdiag=0,overwrite_c=0) inv_c,info = dtrtri(c,lower=0,unitdiag=0,overwrite_c=0) inv_c,info = ctrtri(c,lower=0,unitdiag=0,overwrite_c=0) inv_c,info = ztrtri(c,lower=0,unitdiag=0,overwrite_c=0) x,scale,info = strsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) x,scale,info = dtrsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) x,scale,info = ctrsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) x,scale,info = ztrsyl(a,b,c,trana='N',tranb='N',isgn=1,overwrite_c=0) w,z,info = chbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),lrwork=(compute_v?1+5*n+2*n*n:n),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,info = zhbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),lrwork=(compute_v?1+5*n+2*n*n:n),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,m,ifail,info = chbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) w,z,m,ifail,info = zhbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) t,r,res,x,info = sgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) t,r,res,x,info = dgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) t,r,res,x,info = cgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) t,r,res,x,info = zgglse(a,b,c,d,lwork=max(m+n+p,1),overwrite_a=0,overwrite_b=0,overwrite_c=0,overwrite_d=0) work,info = sgglse_lwork(m,n,p) work,info = dgglse_lwork(m,n,p) work,info = cgglse_lwork(m,n,p) work,info = zgglse_lwork(m,n,p) rcond,info = sppcon(n,ap,anorm,lower=0) rcond,info = dppcon(n,ap,anorm,lower=0) rcond,info = cppcon(n,ap,anorm,lower=0) rcond,info = zppcon(n,ap,anorm,lower=0) x,info = sppsv(n,ap,b,lower=0,overwrite_b=0) x,info = dppsv(n,ap,b,lower=0,overwrite_b=0) x,info = cppsv(n,ap,b,lower=0,overwrite_b=0) x,info = zppsv(n,ap,b,lower=0,overwrite_b=0) ul,info = spptrf(n,ap,lower=0,overwrite_ap=0) ul,info = dpptrf(n,ap,lower=0,overwrite_ap=0) ul,info = cpptrf(n,ap,lower=0,overwrite_ap=0) ul,info = zpptrf(n,ap,lower=0,overwrite_ap=0) uli,info = spptri(n,ap,lower=0,overwrite_ap=0) uli,info = dpptri(n,ap,lower=0,overwrite_ap=0) uli,info = cpptri(n,ap,lower=0,overwrite_ap=0) uli,info = zpptri(n,ap,lower=0,overwrite_ap=0) x,info = spptrs(n,ap,b,lower=0,overwrite_b=0) x,info = dpptrs(n,ap,b,lower=0,overwrite_b=0) x,info = cpptrs(n,ap,b,lower=0,overwrite_b=0) x,info = zpptrs(n,ap,b,lower=0,overwrite_b=0) w,z,info = ssbev(ab,compute_v=1,lower=0,ldab=shape(ab,0),overwrite_ab=1) w,z,info = dsbev(ab,compute_v=1,lower=0,ldab=shape(ab,0),overwrite_ab=1) w,z,info = ssbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,info = dsbevd(ab,compute_v=1,lower=0,ldab=shape(ab,0),liwork=(compute_v?3+5*n:1),overwrite_ab=1) w,z,m,ifail,info = ssbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) w,z,m,ifail,info = dsbevx(ab,vl,vu,il,iu,ldab=shape(ab,0),compute_v=1,range=0,lower=0,abstol=0.0,mmax=(compute_v?(range==2?(iu-il+1):n):1),overwrite_ab=1) m,w,iblock,isplit,info = sstebz(d,e,range,vl,vu,il,iu,tol,order) m,w,iblock,isplit,info = dstebz(d,e,range,vl,vu,il,iu,tol,order) vals,info = ssterf(d,e,overwrite_d=0,overwrite_e=0) vals,info = dsterf(d,e,overwrite_d=0,overwrite_e=0) z,info = sstein(d,e,w,iblock,isplit) z,info = dstein(d,e,w,iblock,isplit) m,w,z,info = sstemr(d,e,range,vl,vu,il,iu,compute_v=1,lwork=max((compute_v?18*n:12*n),1),liwork=(compute_v?10*n:8*n),overwrite_d=0) m,w,z,info = dstemr(d,e,range,vl,vu,il,iu,compute_v=1,lwork=max((compute_v?18*n:12*n),1),liwork=(compute_v?10*n:8*n),overwrite_d=0) work,iwork,info = sstemr_lwork(d,e,range,vl,vu,il,iu,compute_v=1,overwrite_d=0,overwrite_e=0) work,iwork,info = dstemr_lwork(d,e,range,vl,vu,il,iu,compute_v=1,overwrite_d=0,overwrite_e=0) vals,z,info = sstev(d,e,compute_v=1,overwrite_d=0,overwrite_e=0) vals,z,info = dstev(d,e,compute_v=1,overwrite_d=0,overwrite_e=0) cout = ssfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) cout = dsfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) cout = chfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) cout = zhfrk(n,k,alpha,a,beta,c,transr='N',uplo='U',trans='N',overwrite_c=0) arf,info = stpttf(n,ap,transr='N',uplo='U') arf,info = dtpttf(n,ap,transr='N',uplo='U') arf,info = ctpttf(n,ap,transr='N',uplo='U') arf,info = ztpttf(n,ap,transr='N',uplo='U') a,info = stpttr(n,ap,uplo='U') a,info = dtpttr(n,ap,uplo='U') a,info = ctpttr(n,ap,uplo='U') a,info = ztpttr(n,ap,uplo='U') ap,info = stfttp(n,arf,transr='N',uplo='U') ap,info = dtfttp(n,arf,transr='N',uplo='U') ap,info = ctfttp(n,arf,transr='N',uplo='U') ap,info = ztfttp(n,arf,transr='N',uplo='U') a,info = stfttr(n,arf,transr='N',uplo='U') a,info = dtfttr(n,arf,transr='N',uplo='U') a,info = ctfttr(n,arf,transr='N',uplo='U') a,info = ztfttr(n,arf,transr='N',uplo='U') arf,info = strttf(a,transr='N',uplo='U') arf,info = dtrttf(a,transr='N',uplo='U') arf,info = ctrttf(a,transr='N',uplo='U') arf,info = ztrttf(a,transr='N',uplo='U') ap,info = strttp(a,uplo='U') ap,info = dtrttp(a,uplo='U') ap,info = ctrttp(a,uplo='U') ap,info = ztrttp(a,uplo='U') x = stfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) x = dtfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) x = ctfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) x = ztfsm(alpha,a,b,transr='N',side='L',uplo='U',trans='N',diag='N',overwrite_b=0) achol,info = spftrf(n,a,transr='N',uplo='U',overwrite_a=0) achol,info = dpftrf(n,a,transr='N',uplo='U',overwrite_a=0) achol,info = cpftrf(n,a,transr='N',uplo='U',overwrite_a=0) achol,info = zpftrf(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = spftri(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = dpftri(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = cpftri(n,a,transr='N',uplo='U',overwrite_a=0) ainv,info = zpftri(n,a,transr='N',uplo='U',overwrite_a=0) x,info = spftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) x,info = dpftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) x,info = cpftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) x,info = zpftrs(n,a,b,transr='N',uplo='U',overwrite_b=0) rz,tau,info = stzrzf(a,lwork=MAX(m,1),overwrite_a=0) rz,tau,info = dtzrzf(a,lwork=MAX(m,1),overwrite_a=0) rz,tau,info = ctzrzf(a,lwork=MAX(m,1),overwrite_a=0) rz,tau,info = ztzrzf(a,lwork=MAX(m,1),overwrite_a=0) work,info = stzrzf_lwork(m,n) work,info = dtzrzf_lwork(m,n) work,info = ctzrzf_lwork(m,n) work,info = ztzrzf_lwork(m,n) delta,sigma,work,info = slasd4(i,d,z,rho=1.0) delta,sigma,work,info = dlasd4(i,d,z,rho=1.0) a,info = slauum(c,lower=0,overwrite_c=0) a,info = dlauum(c,lower=0,overwrite_c=0) a,info = clauum(c,lower=0,overwrite_c=0) a,info = zlauum(c,lower=0,overwrite_c=0) a = slaswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) a = dlaswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) a = claswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) a = zlaswp(a,piv,k1=0,k2=npiv-1,off=0,inc=1,overwrite_a=0) dlamch = dlamch(cmach) slamch = slamch(cmach) n2 = slange(norm,a) n2 = dlange(norm,a) n2 = clange(norm,a) n2 = zlange(norm,a) alpha,x,tau = slarfg(n,alpha,x,incx=1,overwrite_x=0) alpha,x,tau = dlarfg(n,alpha,x,incx=1,overwrite_x=0) alpha,x,tau = clarfg(n,alpha,x,incx=1,overwrite_x=0) alpha,x,tau = zlarfg(n,alpha,x,incx=1,overwrite_x=0) c = slarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) c = dlarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) c = clarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) c = zlarf(v,tau,c,work,side='L',incv=1,overwrite_c=0) cs,sn,r = slartg(f,g) cs,sn,r = dlartg(f,g) cs,sn,r = clartg(f,g) cs,sn,r = zlartg(f,g) x,y = crot(x,y,c,s,n=(lx-1-offx)/abs(incx)+1,offx=0,incx=1,offy=0,incy=1,overwrite_x=0,overwrite_y=0) x,y = zrot(x,y,c,s,n=(lx-1-offx)/abs(incx)+1,offx=0,incx=1,offy=0,incy=1,overwrite_x=0,overwrite_y=0) major,minor,patch = ilaver() .
Hover to see nodes names; edges to Self not shown, Caped at 50 nodes.
Using a canvas is more power efficient and can get hundred of nodes ; but does not allow hyperlinks; , arrows or text (beyond on hover)
SVG is more flexible but power hungry; and does not scale well to 50 + nodes.
All aboves nodes referred to, (or are referred from) current nodes; Edges from Self to other have been omitted (or all nodes would be connected to the central node "self" which is not useful). Nodes are colored by the library they belong to, and scaled with the number of references pointing them
