163 |
for numDim in [2, 3]: |
for numDim in [2, 3]: |
164 |
for totalNumElem in [100, 200, 400, 800, 1600, 3200, 6400, 12800, 25600, 51200, 102400,204800]: |
for totalNumElem in [100, 200, 400, 800, 1600, 3200, 6400, 12800, 25600, 51200, 102400,204800]: |
165 |
for problem in [solveScalar,solveVector]: |
for problem in [solveScalar,solveVector]: |
166 |
if totalNumElem*2**numDim*numDim< 200000: error=max([problem(numDim, totalNumElem, len_x0, alpha,"DIRECT",""),error]) |
error=max([problem(numDim, totalNumElem, len_x0, alpha,"DIRECT",""),error]) |
167 |
for solver_method in [ "PCG" ]: |
#if totalNumElem*2**numDim*numDim< 200000: error=max([problem(numDim, totalNumElem, len_x0, alpha,"DIRECT",""),error]) |
168 |
for prec in [ "JACOBI", "ILU0" ]: |
# for solver_method in [ "PCG" ]: |
169 |
error=max([problem(numDim, totalNumElem, len_x0, alpha, solver_method,prec),error]) |
# for prec in [ "JACOBI", "ILU0" ]: |
170 |
|
# error=max([problem(numDim, totalNumElem, len_x0, alpha, solver_method,prec),error]) |
171 |
print "----------------" |
print "----------------" |
172 |
print "maximum error over all tests is ",error |
print "maximum error over all tests is ",error |
173 |
print "----------------" |
print "----------------" |