| 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.