Hi, I would like to work on “real image” with FEniCS. So I try gmsh with simple geometry to understand the process. If I understood correctly the other threads of the forum, there are three steps : convert the geo file into a msh file ; convert the msh into a xdmf file with meshio ; read on FEniC…

Dear AMM, Which version of gmsh is installed on your computer. this can be checked by typing gmsh --version in your terminal. With the following geo file: // the square Point(1) = {0, 0, 0,0.1}; Point(2) = {0, 1, 0,0.1}; Point(3) = {1, 1, 0,0.1}; Point(4) = {1, 0, 0,0.1}; Line(1) = {1, 4}; Line(2…

Hi, Does this help: [image] Converting simple 2D mesh from gmsh to dolfin The code I’ve posted there is for a 3D mesh, however, it generalizes to a 2D mesh. Given this simple geo file // the square Point(1) = {0, 0, 0,0.1}; Point(2) = {0, 1, 0,0.1}; Point(3) = {1, 1, 0…

Hi Dokken, Thank you for your answer. When I use pygmsh (on my FEniCS environment), your geometry (that I called testDokken.geo) and the line gmsh -2 testDokken.geo, I get an error message (invalid synthax). What do I need to import for this line to work? AM

Could you please supply the full error message? Also, Which Gmsh version are you using?

First of all, I am really sorry to bother you with a subject that is not related to fenics. I have the 4.3.6 version of gmsh (pygmsh) and the 2.3.3 version of meshio. Also, the message error is an invalid syntax on the line gmsh -2 testDokken.geo AM

Dokken, The Lab did not reinstall gmsh on all computers when they change the ubuntu version : the problem comes from there. Thanks for your answers and again sorry for making you waste your time on such a trivial question.

Hi dokken, So, I ran a code (taken on a thread of the forum) but I get an error message : meshio.write("mf.xdmf", meshio.Mesh(points=msh.points, cells={"line": msh.cells["line"]}, KeyError: 'line' I tried with “triangle” instead of “line” (although my geometry is 2d) and it recognizes the synta…

Hi, Which mesh are you using? You can only write things as a facet-function if you have defined physical lines in gmsh.

Dokken, I use as geometry one of your example (I can not find the exact thread anymore) and I “mesh” with the command gmsh -2 mygeometry.geo -format msh2 The geometry : // Gmsh project created on Tue Mar 26 15:00:10 2019 //+ // THe square Point(1) = {0, 0, 0,0.1}; Point(2) = {0, 1, 0,0.1}; P…

To be able to read line data, you Need to add ‘Physical Line(5) = {1};’ . This adds a physical facet marker with value 5 on Line 1.

Hi dokken, It now recognizes the “line” synthax. From here, I try my luck but do not feel obliged to answer this question since it no longer concerns FEniCS. In fact I can not display my mesh with “plot (mesh)” and he gives me an error message (see screenshot in attachment). Do you know the caus…

Exact Procedure to work with a .geo file on FEniCS (with meshio)?

mesh

dokken July 2, 2019, 8:18pm 2

Hi,
Does this help:

1 Like

Topic		Replies	Views
Cant read .msh file generated from .geo file in Fenics Errors	3	349	May 13, 2023
Using a simple 3d mesh from gmsh General	9	2394	November 19, 2022
Need some help in converting gmsh file into the form which can be used in fenics code mesh	5	424	April 3, 2023
Error in using meshio.write() and create_mesh() for a mesh imported from gmsh	7	1902	September 22, 2021
Converting simple 2D mesh from gmsh to dolfin	15	12095	August 30, 2020

Exact Procedure to work with a .geo file on FEniCS (with meshio)?

Related topics