Examples of UpperTriangDenseMatrix


Examples of no.uib.cipr.matrix.UpperTriangDenseMatrix

                rotation[i].apply(H, i, i, i + 1);
                rotation[i].apply(s, i, i + 1);
            }

            // Update solution in current subspace
            new UpperTriangDenseMatrix(H, i, false).solve(s, s);
            for (int j = 0; j < i; j++)
                x.add(s.get(j), v[j]);

            A.multAdd(-1, x, u.set(b));
            M.apply(u, r);
View Full Code Here

Examples of no.uib.cipr.matrix.UpperTriangDenseMatrix

    }

    @Override
    protected void createPrimary() throws Exception {
        int n = Utilities.getInt(1, max);
        A = new UpperTriangDenseMatrix(n);
        Ad = Utilities.upperPopulate(A);

        // This ensures non-singularity
        Utilities.addDiagonal(A, Ad, 1);
    }
View Full Code Here

Examples of no.uib.cipr.matrix.UpperTriangDenseMatrix

                }
            }*/
               
            // this function should work inplace like Lapack's
            QR qr = QR.factorize(A);
            UpperTriangDenseMatrix R = qr.getR();
           
            A.zero();
           
            // now zero out the lower diagonal
            /*for (int j=0; j<numColumns; ++j) {
                for (int i=j+1; i<numColumns; ++i) {
                    A.set(i,j,0.0);
                }
            }*/
           
            // copy the upper-triangular
            for (int j=0; j<numColumns; ++j) {
                for (int i=0; i<=j; ++i) {
                    A.set(i,j,R.get(i,j));
                }
            }
           
            currentRow = numColumns;
        }
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.