2
Wrappers to LAPACK library
3
==========================
5
flapack -- wrappers for Fortran [*] LAPACK routines
6
clapack -- wrappers for ATLAS LAPACK routines
7
calc_lwork -- calculate optimal lwork parameters
8
get_lapack_funcs -- query for wrapper functions.
10
[*] If ATLAS libraries are available then Fortran routines
11
actually use ATLAS routines and should perform equally
12
well to ATLAS routines.
17
In the following all function names are shown without
18
type prefix (s,d,c,z). Optimal values for lwork can
19
be computed using calc_lwork module.
26
lu,piv,x,info = gesv(a,b,overwrite_a=0,overwrite_b=0)
27
lub,piv,x,info = gbsv(kl,ku,ab,b,overwrite_ab=0,overwrite_b=0)
28
c,x,info = posv(a,b,lower=0,overwrite_a=0,overwrite_b=0)
30
Computational routines::
32
lu,piv,info = getrf(a,overwrite_a=0)
33
x,info = getrs(lu,piv,b,trans=0,overwrite_b=0)
34
inv_a,info = getri(lu,piv,lwork=min_lwork,overwrite_lu=0)
36
c,info = potrf(a,lower=0,clean=1,overwrite_a=0)
37
x,info = potrs(c,b,lower=0,overwrite_b=0)
38
inv_a,info = potri(c,lower=0,overwrite_c=0)
40
inv_c,info = trtri(c,lower=0,unitdiag=0,overwrite_c=0)
42
Linear Least Squares (LLS) Problems
43
-----------------------------------
47
v,x,s,rank,info = gelss(a,b,cond=-1.0,lwork=min_lwork,overwrite_a=0,overwrite_b=0)
49
Computational routines::
51
qr,tau,info = geqrf(a,lwork=min_lwork,overwrite_a=0)
52
q,info = orgqr|ungqr(qr,tau,lwork=min_lwork,overwrite_qr=0,overwrite_tau=1)
54
Generalized Linear Least Squares (LSE and GLM) Problems
55
-------------------------------------------------------
57
Standard Eigenvalue and Singular Value Problems
58
-----------------------------------------------
62
w,v,info = syev|heev(a,compute_v=1,lower=0,lwork=min_lwork,overwrite_a=0)
63
w,v,info = syevd|heevd(a,compute_v=1,lower=0,lwork=min_lwork,overwrite_a=0)
64
w,v,info = syevr|heevr(a,compute_v=1,lower=0,vrange=,irange=,atol=-1.0,lwork=min_lwork,overwrite_a=0)
65
t,sdim,(wr,wi|w),vs,info = gees(select,a,compute_v=1,sort_t=0,lwork=min_lwork,select_extra_args=(),overwrite_a=0)
66
wr,(wi,vl|w),vr,info = geev(a,compute_vl=1,compute_vr=1,lwork=min_lwork,overwrite_a=0)
67
u,s,vt,info = gesdd(a,compute_uv=1,lwork=min_lwork,overwrite_a=0)
69
Computational routines::
71
ht,tau,info = gehrd(a,lo=0,hi=n-1,lwork=min_lwork,overwrite_a=0)
72
ba,lo,hi,pivscale,info = gebal(a,scale=0,permute=0,overwrite_a=0)
74
Generalized Eigenvalue and Singular Value Problems
75
--------------------------------------------------
79
w,v,info = sygv|hegv(a,b,itype=1,compute_v=1,lower=0,lwork=min_lwork,overwrite_a=0,overwrite_b=0)
80
w,v,info = sygvd|hegvd(a,b,itype=1,compute_v=1,lower=0,lwork=min_lwork,overwrite_a=0,overwrite_b=0)
81
(alphar,alphai|alpha),beta,vl,vr,info = ggev(a,b,compute_vl=1,compute_vr=1,lwork=min_lwork,overwrite_a=0,overwrite_b=0)
87
a,info = lauum(c,lower=0,overwrite_c=0)
88
a = laswp(a,piv,k1=0,k2=len(piv)-1,off=0,inc=1,overwrite_a=0)
98
lu,piv,x,info = gesv(a,b,rowmajor=1,overwrite_a=0,overwrite_b=0)
99
c,x,info = posv(a,b,lower=0,rowmajor=1,overwrite_a=0,overwrite_b=0)
101
Computational routines::
103
lu,piv,info = getrf(a,rowmajor=1,overwrite_a=0)
104
x,info = getrs(lu,piv,b,trans=0,rowmajor=1,overwrite_b=0)
105
inv_a,info = getri(lu,piv,rowmajor=1,overwrite_lu=0)
107
c,info = potrf(a,lower=0,clean=1,rowmajor=1,overwrite_a=0)
108
x,info = potrs(c,b,lower=0,rowmajor=1,overwrite_b=0)
109
inv_a,info = potri(c,lower=0,rowmajor=1,overwrite_c=0)
111
inv_c,info = trtri(c,lower=0,unitdiag=0,rowmajor=1,overwrite_c=0)
116
a,info = lauum(c,lower=0,rowmajor=1,overwrite_c=0)
121
Optimal lwork is maxwrk. Default is minwrk.
123
minwrk,maxwrk = gehrd(prefix,n,lo=0,hi=n-1)
124
minwrk,maxwrk = gesdd(prefix,m,n,compute_uv=1)
125
minwrk,maxwrk = gelss(prefix,m,n,nrhs)
126
minwrk,maxwrk = getri(prefix,n)
127
minwrk,maxwrk = geev(prefix,n,compute_vl=1,compute_vr=1)
128
minwrk,maxwrk = heev(prefix,n,lower=0)
129
minwrk,maxwrk = syev(prefix,n,lower=0)
130
minwrk,maxwrk = gees(prefix,n,compute_v=1)
131
minwrk,maxwrk = geqrf(prefix,m,n)
132
minwrk,maxwrk = gqr(prefix,m,n)