Line | |
---|
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']
|
---|
25 | field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1.0]
|
---|
26 | field_values=[\
|
---|
27 | x1, y1,\
|
---|
28 | x2, y2,\
|
---|
29 | nbelements,elapsedtime,\
|
---|
30 | ]
|
---|
Note:
See
TracBrowser
for help on using the repository browser.