bls2.h

svd 算法代码 This directory contains instrumented SVDPACKC Version 1.0 (ANSI-C) programs for compiling

/****************************************************************
 * Sparse SVD Via Hybrid Block Lanczos Procedure for Equivalent *
 * A'A Eigensystems.                                            *
 *                                                              *
 * Global variables and common areas used by bls2 and its       *
 * functions.                                                   *
 ****************************************************************/

int    *pointr ,      /* pointer to column start array         */
	*rowind ;      /* pointer to row indices array          */
float  *value  ;      /* pointer to nonzero values array       */

int    mxvcount,      /* matrix-vector multiplications counter */ 
        mtxvcount,     /* transposed matrix-vector mult. counter*/
	iconv,         /* converged vector counter              */
	nn,            /* current subspace size                 */
	iter;          /* iteration counter                     */


/* memory for the following arrays are allocated in bls2.c and  *
 * blklan2.c (see blklan2.c for dimensions of these arrays)     */

float *alpha,         /* diagonal elements of bidiagonal matrix*
		        * from single vector Lanczos (inner)    *
			* recursion                             */
       *beta,          /* super-diagonals of bidiagonal matrix  *
			* from single vector Lanczos (inner)    *
			* recursion                             */
       *p,             /* work array                            */
       *q,             /* work array                            */
       *t,             /* work array                            */
       *z;             /* work array                            */
float *tres,          /* temporary residual array              */
       *y,             /* work array                            */
       *vv,            /* work array                            */
       *v0,            /* array of converged right S-vectors    */
       *uvtmp,         /* temporary work space                  */
       *pp,            /* left S-vectors of block upper bi-     * 
			* diagonal matrix from outer recursion  */ 
       *qq,            /* right S-vectors of block upper bi-    *
			* diagonal matrix from outer recursion  */
       *ztemp;         /* temporary work space                  */

float **yp,           /* corresponding 2-dimensional           */ 
       **vvp,          /* ...array representation of the above  */
       **uvtmpp,       /* ...linear arrays                      */
       **ppp, 
       **qqp,
       **ztempp;