How to use

plusMult

method

in

cern.jet.math.PlusMult

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 a7 = algebra.mult(a3, a5);
  DoubleMatrix2D a8 = algebra.mult(a7, a6);
  a2.assign(a8, PlusMult.plusMult(1));
  DoubleMatrix2D a9 = K.viewSelection(a, a);
  a9.assign(a2);
a32.assign(KOld, PlusMult.plusMult(-1));
double diff = 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);

        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);