177178179180181182183184185186187
q.combine(n, 1*n, d); SimpleMatrix x = p.transpose().mult(p); x = x.mult(x); x = x.mult(x); v.set(i, x.trace()); x = q.transpose().mult(q); x = x.mult(x); x = x.mult(x); w.set(i, x.trace());
182183184185186187188189190191192
v.set(i, x.trace()); x = q.transpose().mult(q); x = x.mult(x); x = x.mult(x); w.set(i, x.trace()); } return new double[]{stdev(v)/mean(v),stdev(w)/mean(w)}; }