14 |
|
|
15 |
""" |
""" |
16 |
|
|
17 |
from esys.escript import * |
# from esys.escript import * |
18 |
from esys.linearPDEs import * |
# from esys.linearPDEs import * |
19 |
import esys.finley as pdelib |
# import esys.finley as pdelib |
20 |
|
from escript.escript import * |
21 |
|
from escript.linearPDEs import * |
22 |
|
import finley.finley as pdelib |
23 |
from time import time |
from time import time |
24 |
|
|
25 |
from numarray import * |
from numarray import * |
169 |
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]: |
170 |
for problem in [solveScalar,solveVector]: |
for problem in [solveScalar,solveVector]: |
171 |
#for problem in [solveVector]: |
#for problem in [solveVector]: |
172 |
# error=max([problem(numDim, totalNumElem, len_x0, alpha,"PCG",""),error]) |
error=max([problem(numDim, totalNumElem, len_x0, alpha,"PCG",""),error]) |
173 |
error=max([problem(numDim, totalNumElem, len_x0, alpha,"DIRECT",""),error]) |
error=max([problem(numDim, totalNumElem, len_x0, alpha,"DIRECT",""),error]) |
174 |
#if totalNumElem*2**numDim*numDim< 200000: error=max([problem(numDim, totalNumElem, len_x0, alpha,"DIRECT",""),error]) |
#if totalNumElem*2**numDim*numDim< 200000: error=max([problem(numDim, totalNumElem, len_x0, alpha,"DIRECT",""),error]) |
175 |
# for solver_method in [ "PCG" ]: |
# for solver_method in [ "PCG" ]: |