6/11/2023 0 Comments Gmsh lapack librariesQ = FiniteElement('P', mesh.ufl_cell(), 3) V = FiniteElement('N1curl', mesh.ufl_cell(), 2) Mvc_boundaries = xdmf_infile.read_meshtags(mesh, "Grid") With XDMFFile(MPI.comm_world, "tag_all.xdmf", 'r') as xdmf_infile: Mvc_subdomain = xdmf_infile.read_meshtags(mesh, "Grid") Mesh = xdmf_infile.read_mesh(name="Grid") With XDMFFile(MPI.comm_world, "tag_triangle.xdmf", 'r') as xdmf_infile: ("tag_all.xdmf", meshio.Mesh(įrom dolfinx import (MPI, FacetNormal, Function, FunctionSpace,įrom import assemble_scalar, assemble_matrixįrom dolfinx.fem import locate_dofs_topologicalįrom ufl import (TestFunctions, TrialFunctions, grad, inner, curl, split, Points, cells, cell_data, field_data = mesh.points, mesh.cells, mesh.cell_data, mesh.field_dataĬells = np.vstack(np.array([cells.data for cells in mesh.cellsįacet_cells = np.vstack(np.array([cells.data for cells in mesh.cellsĬell_data = mesh.cell_data_dictįacet_data = mesh.cell_data_dict Mesh = pygmsh.generate_mesh(geom, prune_z_0=True, verbose=True, geo_filename="strip.geo") Geom.add_physical(dielectric_domain, 'Diectric') Geom.add_physical( for i in domain_loop_tag + strip_loop_tag], 'PEC') Lines.append(geom.add_line(points, points))Īir_loop_tag = ĭielectric_loop = geom.add_line_loop( for i in dielectric_loop_tag])ĭielectric_domain = geom.add_plane_surface(dielectric_loop)Īir_domain = geom.add_plane_surface(air_loop) Points.append(geom.add_point(point, lcar=lcar)) However, they give different results 10.839566377715625 and 2.9713814818807514 respectively, and I’ll past their codes here for comparison. Hello everyone, I used two dolfinx versions #957 and the newest version to solve the general eigenvalue problem.
0 Comments
Leave a Reply. |