@Test
public void testInvert() {
BlockCholeskyOuterSolver solver = new BlockCholeskyOuterSolver();
for( int i = 1; i <= r*3; i++ ) {
BlockMatrix64F A = createMatrixSPD(i);
BlockMatrix64F A_inv = BlockMatrixOps.createRandom(i,i,-1,1,rand,r);
assertTrue(solver.setA(A.copy()));
solver.invert(A_inv);
BlockMatrix64F B = new BlockMatrix64F(i,i,r);
BlockMatrixOps.mult(A,A_inv,B);
assertTrue(GenericMatrixOps.isIdentity(B,1e-8));
}