Examples of com.nr.la.Cholesky.elmult()

com.nr.la.Cholesky.elmult()
Multiply L y = b, where L is the lower triangular matrix in the stored Cholesky decomposition. y[0..n-1] is input. The result is returned in b[0..n-1]. @param y @param b

      
    }


    ranvec(y);
    double[] yy = new double[y.length];
    ach.elmult(y,r);
    ach.elsolve(r,yy);
    sbeps = 5.e-15;
    localflag = maxel(vecsub(y,yy)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {

View Full Code Here

TOP

All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.