Rev | Line | |
---|
[13437] | 1 | import numpy
|
---|
| 2 | import time
|
---|
| 3 | from model import *
|
---|
| 4 | from bamg import *
|
---|
| 5 | from MatlabFuncs import *
|
---|
| 6 |
|
---|
| 7 | #Simple mesh
|
---|
| 8 | md=bamg(model(),'domain','../Exp/Square.exp','hmax',100000.)
|
---|
| 9 | x1=md.mesh.x
|
---|
| 10 | y1=md.mesh.y
|
---|
| 11 |
|
---|
| 12 | #hVertices
|
---|
| 13 | md=bamg(model(),'domain','../Exp/Square.exp','hmax',300000.,'hVertices',numpy.array([10000.,100000.,400000.,100000.]).reshape(-1,1))
|
---|
| 14 | x2=md.mesh.x
|
---|
| 15 | y2=md.mesh.y
|
---|
| 16 |
|
---|
| 17 | #big mesh
|
---|
| 18 | t0=time.time()
|
---|
| 19 | md=bamg(model(),'domain','../Exp/Square.exp','hmax',3000.)
|
---|
| 20 | nbelements=md.mesh.numberofelements
|
---|
| 21 | elapsedtime=time.time()-t0
|
---|
| 22 |
|
---|
| 23 | #Fields and tolerances to track changes
|
---|
| 24 | field_names =['x1','y1','x2','y2','nbelements','elapsed time']
|
---|
[15647] | 25 | field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1.0]
|
---|
[13437] | 26 | field_values=[\
|
---|
| 27 | x1, y1,\
|
---|
| 28 | x2, y2,\
|
---|
| 29 | nbelements,elapsedtime,\
|
---|
| 30 | ]
|
---|
Note:
See
TracBrowser
for help on using the repository browser.