@Override public void set(int row, int column, double value) { if (column == row) throw new IllegalArgumentException("column == row"); super.set(row, column, value); }
@Override public void set(int row, int column, double value) { if (column == row) throw new IllegalArgumentException("column == row"); super.set(row, column, value); }
@Override void copy(Matrix A) { for (MatrixEntry e : A) if (e.row() <= e.column()) set(e.row(), e.column(), e.get()); }
@Override void copy(Matrix A) { for (MatrixEntry e : A) if (e.row() <= e.column()) set(e.row(), e.column(), e.get()); }
for (MatrixEntry e : A) if (e.column() >= (n - m) + e.row()) R.set(e.row(), e.column() - (n - m), e.get());
for (MatrixEntry e : A) if (e.column() >= (n - m) + e.row()) R.set(e.row(), e.column() - (n - m), e.get());
for (MatrixEntry e : A) if (e.row() <= e.column()) R.set(e.row(), e.column(), e.get());
for (MatrixEntry e : A) if (e.row() <= e.column()) R.set(e.row(), e.column(), e.get());
private DenseCholesky decompose(AbstractDenseMatrix A) { if (n != A.numRows()) throw new IllegalArgumentException("n != A.numRows()"); notspd = false; intW info = new intW(0); if (upper) LAPACK.getInstance().dpotrf(UpLo.Upper.netlib(), A.numRows(), A.getData(), Matrices.ld(A.numRows()), info); else LAPACK.getInstance().dpotrf(UpLo.Lower.netlib(), A.numRows(), A.getData(), Matrices.ld(A.numRows()), info); if (info.val > 0) notspd = true; else if (info.val < 0) throw new IllegalArgumentException(); if (upper) Cu.set(A); else Cl.set(A); return this; }
private DenseCholesky decompose(AbstractDenseMatrix A) { if (n != A.numRows()) throw new IllegalArgumentException("n != A.numRows()"); notspd = false; intW info = new intW(0); if (upper) LAPACK.getInstance().dpotrf(UpLo.Upper.netlib(), A.numRows(), A.getData(), Matrices.ld(A.numRows()), info); else LAPACK.getInstance().dpotrf(UpLo.Lower.netlib(), A.numRows(), A.getData(), Matrices.ld(A.numRows()), info); if (info.val > 0) notspd = true; else if (info.val < 0) throw new IllegalArgumentException(); if (upper) Cu.set(A); else Cl.set(A); return this; }