Package com.nr.la

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

55
56
57
58
59
60
61
62
63
64
65
      fail("*** cholesky: Error in method elmult() or elsolve()");
     
    }

    double[][] ainv = buildMatrix(a);
    ach.inverse(ainv);
    sbeps = 5.e-15;
    localflag = maxel(matsub(matmul(ainv,aposdef),ident(aposdef.length,1.))) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {
      fail("*** cholesky: Error in method inverse()");
View Full Code Here
20
21
22
23
24
25
26
27
28
29
    System.out.println(NRUtil.toString(x));
    ch.elsolve(b, x);
    System.out.println(NRUtil.toString(x));
    System.out.println(ch.logdet());
   
    ch.inverse(inv);
   
    System.out.println(NRUtil.toString(inv));
  }
}
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. All rights reserved.
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.