LAPROGNC(Linear Algebra PROGrams in Numerical computation)
LICENSE

A new implementation of DQDS ALGORITHM WITH NEW SHIFT STRATEGIES as the offshoot of the implementation of the orthogonal QD algorithm
ddqds.f   sdqds.f

A new implementation of MDLVS ALGORITHM WITH NEW SHIFT STRATEGIES as the offshoot of the implementation of the orthogonal QD algorithm
dmdlvs.f

PARALLEL BISECTION METHOD for a symmetric tridiagonal matrix
dlaebz.f   dstebz.f

PARALLEL INVERSE ITERATION ALGORITHM with OpenMP-based classical Gram-Schmidt*2 reorthogonalization for a symmetric tridiagonal matrix
dstein.f

OpenMP-based PARALLEL BLOCKED INVERSE ITERATION ALGORITHM with DGEMM-based blocked classical Gram-Schmidt*2 reorthogonalization for a symmetric tridiagonal matrix
dstein3.f

PARALLEL INVERSE ITERATION ALGORITHM with compact WY reorthogonalization for a symmetric tridiagonal matrix
dsteinwy.f

PARALLEL BISECTION METHOD and INVERSE ITERATION ALGORITHM for a symmetric band matrix
BBiInv.tgz

A new implementation of dnrm2.f
dnrm2n.f

A new implementation of ORTHOGONAL QD ALGORITHM WITH NEW SHIFT STRATEGIES for thick Restart Golub-Kahan-Lanczos method
Our implementation of the orthogonal qd algorithm is more accurate than DBDSQR.

This code can compute singular values, left singular vectors and right singular vectors of a lower/upper bidiagonal matrix.
doqds1.f   soqds1.f
This code can compute singular values and left singular vectors of a lower/upper bidiagonal matrix.
doqds2.f   soqds2.f
This code can compute singular values and right singular vectors of a lower/upper bidiagonal matrix.
doqds3.f   soqds3.f
This code can compute only singular values of a lower/upper bidiagonal matrix.
doqdsv.f   doqdsv.f
Subroutines for ORTHOGONAL QD ALGORITHM WITH NEW SHIFT STRATEGIES
dlartg6.f   dlartg7.f   dfma0.c   slartg6.f   slartg7.f   sfma0.c

Thick Restart Golub-Kahan-Lanczos method

TRGKL2.tgz

Thick Restart Lanczos method

TRLAN2.tgz

SVD using ARPACK
dsvd_file.f   svd_make_test_file.c   ARmake.inc for Windows

For any quetion, contact us via "kkimur -*- amp.i.kyoto-u.ac.jp" (please change "-*-" into "@")