TestLinbcg test = new TestLinbcg(sa);
test.solve(bbb,xx,ITOL,TOL,ITMAX,iter,err);
sbeps = 5.e-15;
localflag = maxel(vecsub(sa.ax(xx),bbb)) > sbeps;
globalflag = globalflag || localflag;
if (localflag) {
fail("*** Linbcg: Derived class gives inconsistent solution vector");
}