~ubuntu-branches/ubuntu/raring/python-scipy/raring-proposed

« back to all changes in this revision

Viewing changes to Lib/lib/lapack/info.py

  • Committer: Bazaar Package Importer
  • Author(s): Matthias Klose
  • Date: 2007-01-07 14:12:12 UTC
  • mfrom: (1.1.1 upstream)
  • Revision ID: james.westby@ubuntu.com-20070107141212-mm0ebkh5b37hcpzn
* Remove build dependency on python-numpy-dev.
* python-scipy: Depend on python-numpy instead of python-numpy-dev.
* Package builds on other archs than i386. Closes: #402783.

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
"""
 
2
Wrappers to LAPACK library
 
3
==========================
 
4
 
 
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.
 
9
 
 
10
[*] If ATLAS libraries are available then Fortran routines
 
11
    actually use ATLAS routines and should perform equally
 
12
    well to ATLAS routines.
 
13
 
 
14
Module flapack
 
15
++++++++++++++
 
16
 
 
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.
 
20
 
 
21
Linear Equations
 
22
----------------
 
23
 
 
24
  Drivers::
 
25
 
 
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)
 
29
 
 
30
  Computational routines::
 
31
 
 
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)
 
35
 
 
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)
 
39
 
 
40
    inv_c,info = trtri(c,lower=0,unitdiag=0,overwrite_c=0)
 
41
 
 
42
Linear Least Squares (LLS) Problems
 
43
-----------------------------------
 
44
 
 
45
  Drivers::
 
46
 
 
47
    v,x,s,rank,info = gelss(a,b,cond=-1.0,lwork=min_lwork,overwrite_a=0,overwrite_b=0)
 
48
 
 
49
  Computational routines::
 
50
 
 
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)
 
53
 
 
54
Generalized Linear Least Squares (LSE and GLM) Problems
 
55
-------------------------------------------------------
 
56
 
 
57
Standard Eigenvalue and Singular Value Problems
 
58
-----------------------------------------------
 
59
 
 
60
  Drivers::
 
61
 
 
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)
 
68
 
 
69
  Computational routines::
 
70
 
 
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)
 
73
 
 
74
Generalized Eigenvalue and Singular Value Problems
 
75
--------------------------------------------------
 
76
 
 
77
  Drivers::
 
78
 
 
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)
 
82
 
 
83
 
 
84
Auxiliary routines
 
85
------------------
 
86
 
 
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)
 
89
 
 
90
Module clapack
 
91
++++++++++++++
 
92
 
 
93
Linear Equations
 
94
----------------
 
95
 
 
96
  Drivers::
 
97
 
 
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)
 
100
 
 
101
  Computational routines::
 
102
 
 
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)
 
106
 
 
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)
 
110
 
 
111
    inv_c,info = trtri(c,lower=0,unitdiag=0,rowmajor=1,overwrite_c=0)
 
112
 
 
113
Auxiliary routines
 
114
------------------
 
115
 
 
116
  a,info = lauum(c,lower=0,rowmajor=1,overwrite_c=0)
 
117
 
 
118
Module calc_lwork
 
119
+++++++++++++++++
 
120
 
 
121
Optimal lwork is maxwrk. Default is minwrk.
 
122
 
 
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)
 
133
 
 
134
 
 
135
"""
 
136
postpone_import = 1