p=myran.int32p()%N;
q=myran.int32p()%M;
mat[p][q]=1.e99;
}
System.out.printf(" Initial discrepancy: %g\n", maxel(matsub(actual,mat)));
Laplace_interp mylaplace = new Laplace_interp(mat);
mylaplace.solve();
System.out.printf(" Final discrepancy: %g\n", maxel(matsub(actual,mat)));
localflag = maxel(matsub(actual,mat)) > sbeps;
globalflag = globalflag || localflag;
if (localflag) {
fail("*** Laplace_interp: Inaccurate Laplace interpolation of missing matrix data.");