Hello. I am trying to run in parallel using < mpirun -n 2 python3 Beam.py >. However, it gives this error: Process 0: *** Warning: XML file ‘beammesh.xml’ is very large. XML files are parsed in serial, which is not scalable. Use XMDF/HDF5 for scalable IO in parallel Process 1: Solving linear varia…

As the message says, you should use the XDMF file format to read in the mesh, as described in for instance: [image] Gmsh import error mesh It seems like you have not read the post I linked you to. You cannot use the meshio-convert script out of the box for XDMF, and you ha…

I converted the mesh to xdmf file. How will I use it? It gives an error: Traceback (most recent call last): File “Beam.py”, line 17, in mesh = Mesh(‘mesh.xdmf’) RuntimeError: *** Error: Unable to open file. *** Reason: Unknown file type (“.xdmf”) for file “mesh.xdmf”. *** Where: This er…

Now, I get this error: Traceback (most recent call last): File “Beam.py”, line 23, in cf = cpp.mesh.MeshFunctionSizet(mesh, mvc) AttributeError: module ‘mshr.cpp’ has no attribute ‘mesh’ Traceback (most recent call last): File “Beam.py”, line 23, in cf = cpp.mesh.MeshFunctionSizet(mesh, mvc…

[image] dries: from mshr import * I would suggest removing this import. Wildcard imports from multiple libraries are considered as bad practice, see for instance using wildcard imports (from … import *) — Python Anti-Patterns documentation

Yes, that solved the problem. However, I have a problem with using the mesh named cf. When I use cf instead of mesh, it gives error. This is a part of my original code: # mesh from dolfin import * import sympy mesh = Mesh() mvc = MeshValueCollection("size_t", mesh, mesh.topology().dim()) with XDMFF…

[image] dries: However, I have a problem with using the mesh named cf. cf is not a mesh, it is a marker for the cells in your mesh (called mesh)

I used mesh as in the previous post and not cf but the program shuts down with a message: Process 1: Solving linear variational problem. Process 0: Solving linear variational problem. Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has b…

As I do not have access to your mesh, I cannot give any further help.

[image] dries: mvc = MeshValueCollection("size_t", mesh, mesh.topology().dim()) with XDMFFile("mesh.xdmf") as infile: infile.read(mesh) infile.read(mvc, "name_to_read") cf = cpp.mesh.MeshFunctionSizet(mesh, MVC) here it should be cf = cpp.mesh.MeshFunctionSizet(mesh, mvc)

Running in parallel using MPI

I/O

dokken October 9, 2022, 2:14pm 4

See for instance Need help converting GMSH to FEniCS - #17 by dokken

Topic		Replies	Views
Missing Cells When Reading in Mesh in Parallel General mesh , dolfinx	9	107	August 9, 2024
Problem regarding reading mesh in parallel mesh	6	1319	February 12, 2020
Create Mesh object or xdmf file in parallel mesh	5	1731	May 1, 2020
Mesh conversion (XDMF) in a MPI environment mesh	4	1813	June 16, 2020
Converting 3D mesh to xdmf/h5	3	2706	March 30, 2020

Running in parallel using MPI

Related topics