public void dger(double alpha, DoubleMatrix1D x, DoubleMatrix1D y, DoubleMatrix2D A) { cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0); for (int i=A.rows(); --i >= 0; ) { fun.multiplicator = alpha * x.getQuick(i); A.viewRow(i).assign(y,fun); } } public double dnrm2(DoubleMatrix1D x) {
public void dger(double alpha, DoubleMatrix1D x, DoubleMatrix1D y, DoubleMatrix2D A) { cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0); for (int i=A.rows(); --i >= 0; ) { fun.multiplicator = alpha * x.getQuick(i); A.viewRow(i).assign(y,fun); } } public double dnrm2(DoubleMatrix1D x) {
s.assign(new Algebra().mult(vv, x)); } else { s.assign(new Algebra().mult(vv, x), PlusMult.plusMult(-1));
for (int i=m; --i>=0; ) Crows[i] = C.viewRow(i); final cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0);
for (int i=m; --i>=0; ) Crows[i] = C.viewRow(i); final cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0);
for (int i=m; --i>=0; ) Crows[i] = C.viewRow(i); final cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0);
for (int i=m; --i>=0; ) Crows[i] = C.viewRow(i); final cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0);
for (int i=m; --i>=0; ) Crows[i] = C.viewRow(i); final cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0);
for (int i=m; --i>=0; ) Crows[i] = C.viewRow(i); final cern.jet.math.PlusMult fun = cern.jet.math.PlusMult.plusMult(0);
DoubleMatrix2D a8b = omega.viewSelection(v, v); a8b.assign(a8); omega.viewSelection(v, v).assign(a9, PlusMult.plusMult(1)); DoubleMatrix2D a31 = algebra.mult(a30, a29); omega.viewSelection(v, v).assign(tempVar); omega.viewSelection(v, v).assign(a31, PlusMult.plusMult(1)); } else { DoubleMatrix2D oInv = new DenseDoubleMatrix2D(p, p); a32.assign(omegaOld, PlusMult.plusMult(-1)); double diff1 = algebra.norm1(a32); a33.assign(bOld, PlusMult.plusMult(-1)); double diff2 = algebra.norm1(a32);
DoubleMatrix2D a8b = omega.viewSelection(v, v); a8b.assign(a8); omega.viewSelection(v, v).assign(a9, PlusMult.plusMult(1)); DoubleMatrix2D a31 = algebra.mult(a30, a29); omega.viewSelection(v, v).assign(tempVar); omega.viewSelection(v, v).assign(a31, PlusMult.plusMult(1)); } else { DoubleMatrix2D oInv = new DenseDoubleMatrix2D(p, p); a32.assign(omegaOld, PlusMult.plusMult(-1)); double diff1 = algebra.norm1(a32); a33.assign(bOld, PlusMult.plusMult(-1)); double diff2 = algebra.norm1(a32);