if (v <= 0)
throw new IllegalArgumentException (" v <= 0");
B.assign (F.mult(v));
DoubleFactory2D fac = DoubleFactory2D.dense;
final DoubleMatrix2D B0 = fac.identity(n); // B^0 = I
final DoubleMatrix2D B2 = alge.mult(B, B); // B^2
final DoubleMatrix2D B4 = alge.mult(B2, B2); // B^4
DoubleMatrix2D T = B2.copy(); // T = work matrix
DoubleMatrix2D W = B4.copy(); // W = work matrix