MeshQuad4-Regular.py
daf6573f1f95neighbour

MeshQuad4-Regular.py
View Options


			import h5py
			import numpy as np
			import GooseFEM as gf
			import lxml.etree as etree

			# ====================== create fictitious configuration + store to HDF5-file ======================

			# create mesh object
			mesh = gf.Mesh.Quad4.Regular(9,11)

			# open data file
			name = 'MeshQuad4-Regular'
			file = h5py.File(name+'.hdf5','w')

			# write nodal coordinates and connectivity
			file['coor'] = mesh.coor()
			file['conn'] = mesh.conn()

			# ======================================== write XDMF-file =========================================

			# get mesh dimensions
			dims = dict(
			nnode = mesh.nnode(),
			ndim = mesh.ndim(),
			nelem = mesh.nelem(),
			nne = mesh.nne(),
			)

			# initialize file
			root = etree.fromstring('<Xdmf Version="2.0"></Xdmf>')
			domain = etree.SubElement(root, "Domain")
			grid = etree.SubElement(domain, "Grid", Name="Mesh")

			# add connectivity
			conn = etree.SubElement(grid, "Topology", TopologyType="Quadrilateral", NumberOfElements='{nelem:d}'.format(**dims))
			data = etree.SubElement(conn, "DataItem", Dimensions='{nelem:d} {nne:d}'.format(**dims), Format="HDF")
			data.text = "{0:s}.hdf5:/conn".format(name)

			# add coordinates
			coor = etree.SubElement(grid, "Geometry", GeometryType="XY")
			data = etree.SubElement(coor, "DataItem", Dimensions='{nnode:d} {ndim:d}'.format(**dims), Format="HDF")
			data.text = "{0:s}.hdf5:/coor".format(name)

			# write to file
			open(name+'.xdmf','wb').write(etree.tostring(root, pretty_print=True))