147 
print "" 
print "" 
148 
error=0. 
error=0. 
149 
for numDim in [2,3]: 
for numDim in [2,3]: 
150 
for totalNumElem in [100, 200, 400, 800, 1600, 3200, 6400, 12800, 25600, 51200, 102400]: 
for totalNumElem in [100, 200, 400, 800, 1600, 3200, 6400]: 
151 

# for totalNumElem in [100, 200, 400, 800, 1600, 3200, 6400, 12800, 25600, 51200, 102400]: 
152 
for problem in [solveScalar,solveVector]: 
for problem in [solveScalar,solveVector]: 
153 
# for solver_method in [ LinearPDE.PRES20, LinearPDE.PCG, LinearPDE.DIRECT, LinearPDE.BICGSTAB]: 
# for solver_method in [ LinearPDE.PRES20, LinearPDE.PCG, LinearPDE.DIRECT, LinearPDE.BICGSTAB]: 
154 
for solver_method in [ LinearPDE.PRES20, LinearPDE.PCG, LinearPDE.BICGSTAB]: 
for solver_method in [ LinearPDE.PCG ]: 
155 
error=max([problem(numDim, totalNumElem, len_x0, alpha, solver_method),error]) 
error=max([problem(numDim, totalNumElem, len_x0, alpha, solver_method),error]) 
156 
print "" 
print "" 
157 
print "maximum error over all tests is ",error 
print "maximum error over all tests is ",error 