question: static double fabs ( double a ) { return Math . abs ( a ) ; }
target:  Description ----------- Function performs the first step of the Lanczos algorithm.  It also does a step of extended local re-orthogonalization. Arguments  --------- (input) n      dimension of the eigenproblem for matrix B (output) ierr   error flag wptr   array of pointers that point to work space that contains wptr[0]             r[j] wptr[1]             q[j] wptr[2]             q[j-1] wptr[3]             p wptr[4]             p[j-1] wptr[6]             diagonal elements of matrix T  Functions used -------------- BLAS		svd_daxpy, svd_datx, svd_dcopy, svd_ddot, svd_dscal USER		store, opb LAS		startv