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 |