import numpy
import time
from model import *
from bamg import *
from MatlabFuncs import *

#Simple mesh
md=bamg(model(),'domain','../Exp/Square.exp','hmax',100000.)
x1=md.mesh.x
y1=md.mesh.y

#hVertices
md=bamg(model(),'domain','../Exp/Square.exp','hmax',300000.,'hVertices',numpy.array([10000.,100000.,400000.,100000.]).reshape(-1,1))
x2=md.mesh.x
y2=md.mesh.y

#big mesh
t0=time.time()
md=bamg(model(),'domain','../Exp/Square.exp','hmax',3000.)
nbelements=md.mesh.numberofelements
elapsedtime=time.time()-t0

#Fields and tolerances to track changes
field_names     =['x1','y1','x2','y2','nbelements','elapsed time']
field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1.5]
field_values=[\
	x1, y1,\
	x2, y2,\
	nbelements,elapsedtime,\
	]