gmsh issueshttps://gitlab.onelab.info/gmsh/gmsh/-/issues2021-04-26T03:01:24Zhttps://gitlab.onelab.info/gmsh/gmsh/-/issues/1154Call to gmsh/model/occ/fragment only returns highest dimensional resulting en...2021-04-26T03:01:24ZStephen BanksCall to gmsh/model/occ/fragment only returns highest dimensional resulting entities if tag >= 0 || !preserveNumberingThe returned outDimTags from calling gmsh/model/occ/fragments only includes the highest dimensional new shapes if tag > 0 or Geometry.OCCBooleanPreserveNumbering = 0, otherwise it includes all dimensional shapes. I think it is probably m...The returned outDimTags from calling gmsh/model/occ/fragments only includes the highest dimensional new shapes if tag > 0 or Geometry.OCCBooleanPreserveNumbering = 0, otherwise it includes all dimensional shapes. I think it is probably meant to always include all resulting entities.
This is due to the source code for OCC_Internals::booleanOperator (GModelIO_OCC.cpp) being as follows (I've used ellipsis to indicate irrelevant code to this discussion):
`if(tag >= 0 || !preserveNumbering) {
...
_multiBind(result, tag, outDimTags, true, true);
...
}
else {
...
_multiBind(result, -1, outDimTags, false, true, true);
...
}`
I suspect both calls to _multiBind should have highestDimOnly = false, however it is only false in the else branch. I believe this is why the outDimTags only contains the highest dimensional entities in the (tag >= 0 || !preserveNumbering) branch. I'm not sure if this is intentional or not.https://gitlab.onelab.info/gmsh/gmsh/-/issues/1134Calling the general Boolean Function << BooleanOperator (...) >> from geo-file2021-01-01T01:20:49Zsamy mukadiCalling the general Boolean Function << BooleanOperator (...) >> from geo-fileIs it possible to call the general Boolean function from the geo-file?
i.e. <br><br> <b>
BooleanOperator {CustomBoolean}{ Surface{1}; }{ Surface{2};};</b><br>
I want to introduce a customBoolean Method. Not part of the 4.Is it possible to call the general Boolean function from the geo-file?
i.e. <br><br> <b>
BooleanOperator {CustomBoolean}{ Surface{1}; }{ Surface{2};};</b><br>
I want to introduce a customBoolean Method. Not part of the 4.https://gitlab.onelab.info/gmsh/gmsh/-/issues/2773Can an stl file with a repeating mesh generate a Volume?2024-01-17T12:09:51ZbuzhiqiuyuebaiCan an stl file with a repeating mesh generate a Volume?I have two closed stl files, they are both closed but have overlapping parts, I can generate Volume using one of the stl meshes alone, but generating Volume using both stl at the same time will cause the Volume mesh to fail due to the ov...I have two closed stl files, they are both closed but have overlapping parts, I can generate Volume using one of the stl meshes alone, but generating Volume using both stl at the same time will cause the Volume mesh to fail due to the overlapping meshes, is there a good solution for this?
![image](/uploads/ec05c36d10824972241775e111ae06a3/image.png)
![image](/uploads/819e4914fd009b82b9ff297241ce572b/image.png)
![image](/uploads/9ac0323b3670a38a1c1c112f851ca2aa/image.png)
The file I use:
[Eye_Lens.stl](/uploads/5ca990fcb71705ae93c8ca2d14cc67a0/Eye_Lens.stl)
[Eye_Vitreous_.stl](/uploads/2b223b09a7626569ca71854c5e271a9d/Eye_Vitreous_.stl)
[loh.geo](/uploads/d5e8831a0000a383a8339c7fda0fcac4/loh.geo)https://gitlab.onelab.info/gmsh/gmsh/-/issues/1022Can Gmsh create deformed prism? and how to? and partially hexahedralize?2020-10-12T01:37:01ZJunyeong HeoCan Gmsh create deformed prism? and how to? and partially hexahedralize?Hello,
first, I'm really appreciated for this wonderful program.
I'm Doctor's degree student working on magnetotelluric survey modeling.
So, what I'm doing is making local earth model and modeling Maxwell's equations by FEM.
Recently,...Hello,
first, I'm really appreciated for this wonderful program.
I'm Doctor's degree student working on magnetotelluric survey modeling.
So, what I'm doing is making local earth model and modeling Maxwell's equations by FEM.
Recently, my studying area contains very shallow sea, compared to its total area.
I attached simple model for test.
In the model, volume2 is the sea.
As mentioned, this is just simple test model.
The final model will be larger and contain shallower sea.
If I tetrahedralize the model, it will contain huge amount of tetrahedron in order to preserve mesh quality.
So I think prism element is the best choice for the sea part (than hexahedral)
because I have to apply sea bottom topography from DEM data
and triangles are better than quadrangles (all nodes are co-planer automatically).
However, when I tested the simple model, the mesh contains only the ill-shaped tetrahedron and pyramids.
(also, error messages are printed)
In my case, I can't use extrusion method because of the complex topography.
I can't figure out how to make prisms without using extrusion. (not the regular one, deformed-shape prisms)
Also, it would be best if I can discretize some part of the model that I want into tetrahedron (+pyramid, prisms)
and other into hexahedron
So I wonder that
1) can Gmsh create irregular prisms
2) can Gmsh recombine specified volume (not total model)
[test.geo](/uploads/2abea9a6adcfce852d59239e8baf0cff/test.geo)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2024Can not mesh extreme long cylinder with a small radius2022-11-03T16:43:03ZSiyuan SunCan not mesh extreme long cylinder with a small radiusDear experts
I want to mesh some wires with 1400 [mm] long and 0.0125 [mm] radius. When the 2D mesh started it gives the error: "Identical points in triangulation: increase element size or Mesh.RandomFactor". How can I solve this proble...Dear experts
I want to mesh some wires with 1400 [mm] long and 0.0125 [mm] radius. When the 2D mesh started it gives the error: "Identical points in triangulation: increase element size or Mesh.RandomFactor". How can I solve this problem? Also I have another question, since I want to mesh 10k wires so is there a way to mesh just 1 of them then copy for the rest?
Code:
```
SetFactory("OpenCASCADE");
Mesh.MeshSizeFromCurvature = 10;
Cylinder(1) = {0,0,0,0,0,1400,0.0125};
Cylinder(2) = {-8.4,-8,0,0,0,1400,0.063};
Cylinder(3) = {-8.4,0,0,0,0,1400,0.063};
Cylinder(4) = {-8.4,8,0,0,0,1400,0.063};
Cylinder(5) = {0,-8,0,0,0,1400,0.063};
Cylinder(6) = {0,8,0,0,0,1400,0.063};
Cylinder(7) = {8.4,-8,0,0,0,1400,0.063};
Cylinder(8) = {8.4,0,0,0,0,1400,0.063};
Cylinder(9) = {8.4,8,0,0,0,1400,0.063};
Cylinder(10) = {0,0,-50,0,0,1500,20};
Coherence;
Physical Volume("wire1", 30001) = {1};
Physical Volume("wire2", 30002) = {2};
Physical Volume("wire3", 30003) = {3};
Physical Volume("wire4", 30004) = {4};
Physical Volume("wire5", 30005) = {5};
Physical Volume("wire6", 30006) = {6};
Physical Volume("wire7", 30007) = {7};
Physical Volume("wire8", 30008) = {8};
Physical Volume("wire9", 30009) = {9};
Physical Volume("infinity", 30010) = {10};
Physical Surface("wire1boundary", 20001) = {1};
Physical Surface("wire2boundary", 20002) = {4};
Physical Surface("wire3boundary", 20003) = {7};
Physical Surface("wire4boundary", 20004) = {10};
Physical Surface("wire5boundary", 20005) = {13};
Physical Surface("wire6boundary", 20006) = {16};
Physical Surface("wire7boundary", 20007) = {19};
Physical Surface("wire8boundary", 20008) = {22};
Physical Surface("wire9boundary", 20009) = {25};
Physical Surface("infinity1", 20010) = {28};
Physical Surface("infinity2", 20011) = {29};
Physical Surface("infinity3", 20012) = {30};
```https://gitlab.onelab.info/gmsh/gmsh/-/issues/1789Can we catch warnings output by GMSH CXX API?2022-02-24T07:38:32ZTingchang YinCan we catch warnings output by GMSH CXX API?Hi. I am using GMSH C++ API now. Is there a way to catch `warning`? Sometimes I get warnings on terminal (like the following picture), and the meshing process is somehow very very slow. I don't care the reason behind. What I want to do i...Hi. I am using GMSH C++ API now. Is there a way to catch `warning`? Sometimes I get warnings on terminal (like the following picture), and the meshing process is somehow very very slow. I don't care the reason behind. What I want to do is to catch these warnings and then do some other things, like the way to catch GMSH `error`: `try{}catch(...){//do something}`.
![Screenshot_from_2022-02-24_11-04-41](/uploads/081cc1da7b4f2591c5d1d36c62a2c733/Screenshot_from_2022-02-24_11-04-41.png)https://gitlab.onelab.info/gmsh/gmsh/-/issues/1517Can you do Mesh refinement by splitting only vertically rather than splitting...2021-08-16T04:44:42ZAli Yazdan PanahCan you do Mesh refinement by splitting only vertically rather than splitting all elements?I have managed to get the mesh below by hand through manipulation of Gmsh script.I have failed to create similar mesh using Gmsh features without creating unwanted 1D elements and surfaces in the final output, hence why I had to adjust b...I have managed to get the mesh below by hand through manipulation of Gmsh script.I have failed to create similar mesh using Gmsh features without creating unwanted 1D elements and surfaces in the final output, hence why I had to adjust by hand.
![Snip20210808_4](/uploads/cfa5db134e68ed6d7cb3d4250d446de4/Snip20210808_4.png)
The structure of this mesh is exactly what I am looking for with regular elements along the bottom lines up to y=0.5 (where the inflow ends) and also regular next to the LHS wall above the y =0.5, but with gradual refinements, and then all else is coarser as moving towards the RHS corner.
However the resolution near the inflow is not enough due to large elements' aspect ration and, I need to now:
- keep the aspect ratio of the elements down the bottom along the x-axis, close to max 6:1 (width:height) atm it's 47:1! This will be best done for me with refine by splitting asI have manipulated the mesh script to get this mesh and not the geometry, however it refines the whole domain and makes the number of elements huge.
- Make elements within the red area gradually more coarser and reducing the number of elements. The total No of elements for the whole mesh should ideally stay below 800. the current mesh in the photo has around 768 elements.
![Snip20210808_4_copy](/uploads/1d5d3d98a96af60582277bd81c7ab501/Snip20210808_4_copy.png)
Help I need:
- Can I only vertically split the elements?
- Is there a way to get the same mesh without using auxiliary lines and avoiding creation of unwanted surfaces and 1D elements?
- Is there a way for me to coarsen the red area (top RHS) while keeping the rest intact and performing on the same mesh script?
- Any other advise and methods to get this simple 2D mesh would be appreciated (which has turned to nightmare for an amateur mesh-er like me!)
Thanks, Alihttps://gitlab.onelab.info/gmsh/gmsh/-/issues/2799can't using OCC in Gmsh2024-02-08T09:32:45ZWENBO ZHANGcan't using OCC in Gmshwhen i compile gmsh source with cmake, why there are no variable “OCC_INC”. there are any information about OCC module,
![屏幕截图_2024-02-07_115914](/uploads/82732948c758b7981b2cfa8649d8005d/屏幕截图_2024-02-07_115914.png)when i compile gmsh source with cmake, why there are no variable “OCC_INC”. there are any information about OCC module,
![屏幕截图_2024-02-07_115914](/uploads/82732948c758b7981b2cfa8649d8005d/屏幕截图_2024-02-07_115914.png)https://gitlab.onelab.info/gmsh/gmsh/-/issues/1963Cannot create wire2023-07-31T11:25:58ZMathieu LambertCannot create wireI get an error while generating a wire from curves in an imported STEP file. For some reason, it works up until curve 11, i.e. from curves 3 to 10. The curves were generated in Salome Shaper and I can generate the wire and the pipe just ...I get an error while generating a wire from curves in an imported STEP file. For some reason, it works up until curve 11, i.e. from curves 3 to 10. The curves were generated in Salome Shaper and I can generate the wire and the pipe just fine in Salome Shaper.[Inductor3D.geo](/uploads/558b54a41fefc4ec1cca86cafd7c61c2/Inductor3D.geo)[Part_2.step](/uploads/3a1a27550159cc05cffa69ba33a795ce/Part_2.step)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2696CGNS exportation2023-12-13T10:41:48ZGiacomo BenedettiCGNS exportationHi,
I think there should be a bug in the cgns format export, in particular I think the grid location should be changed from cell centre to face centre. Is this the case or should I change some export option?
BestHi,
I think there should be a bug in the cgns format export, in particular I think the grid location should be changed from cell centre to face centre. Is this the case or should I change some export option?
Besthttps://gitlab.onelab.info/gmsh/gmsh/-/issues/2036CGNS I/O and Tecplot2023-12-13T23:48:29ZMichael Ermakovermakov@ipmnet.ruCGNS I/O and Tecplot1. Gmsh does not read tut21.cgns (ADF) from benchmarks/cgns, reporting
![cgns_tut21](/uploads/a247724f2bc6d6b746b8ebef7c1281db/cgns_tut21.png)
2. Tecplot 360 reads and visualizes tut21.cgns correctly, however it does not read Gmsh mes...1. Gmsh does not read tut21.cgns (ADF) from benchmarks/cgns, reporting
![cgns_tut21](/uploads/a247724f2bc6d6b746b8ebef7c1281db/cgns_tut21.png)
2. Tecplot 360 reads and visualizes tut21.cgns correctly, however it does not read Gmsh mesh-export to cgns format (HDF), reporting
![cgns](/uploads/ff52a855d50c31b62e9f2e316f28a008/cgns.png).
Thank you in advance.https://gitlab.onelab.info/gmsh/gmsh/-/issues/2745Changing Points mesh size using Min Field2024-02-04T13:44:17ZPietro LizzioChanging Points mesh size using Min Field![nomesh](/uploads/ae31cf2beed3e9c887a0ccdd8f151740/nomesh.png)
Hello! I'm trying to change the mesh size of some points (shown in the first image), using the Min Field, via these instructions on a geo file:
`Field[1] = Constant;
Field...![nomesh](/uploads/ae31cf2beed3e9c887a0ccdd8f151740/nomesh.png)
Hello! I'm trying to change the mesh size of some points (shown in the first image), using the Min Field, via these instructions on a geo file:
`Field[1] = Constant;
Field[1].IncludeBoundary = 0;
Field[1].PointsList = {9,10};
Field[1].VIn = 1;
Field[1].VOut = 100000;
Field[2] = Constant;
Field[2].IncludeBoundary = 1;
Field[2].CurvesList = {1,2,3,4,5,6,7,8,9};
Field[2].VIn = 50;
Field[2].VOut = 100000;
Field[3] = Constant;
Field[3].IncludeBoundary = 0;
Field[3].SurfacesList = {1};
Field[3].VIn = 50;
Field[3].VOut = 100000;
Field[4] = Min;
Field[4].FieldsList = {1,2,3};Background Field = 4;`
The result is this:
![meshsize](/uploads/3fba305b3f17c453dd5ca67b465c1661/meshsize.png)
While the point (9) on the surface boundary is affected by the instructions provided in Field[1] and the mesh around it does seem to get finer, the point (10) doesn not seem to particularly care about VIn = 1. Am I missibg something? Maybe some instruction I'm overlooking?
For further info, this is the result if I remove Field[1]:
![meshsize_before](/uploads/7c06828b6aca63a9d9fd809ee54f85e1/meshsize_before.png)
Thank you!https://gitlab.onelab.info/gmsh/gmsh/-/issues/2179check which cross-section belongs to which volume2023-10-21T18:53:58ZGregor Hofmeistercheck which cross-section belongs to which volumeHello together,
I have a new issue or problem with handling gmsh.
My code has following steps:
1. importing the step file
2. adding physical groups for each volume/solid as a specific material
3. slicing over the length of the model
4. ...Hello together,
I have a new issue or problem with handling gmsh.
My code has following steps:
1. importing the step file
2. adding physical groups for each volume/solid as a specific material
3. slicing over the length of the model
4. getting entities, for further inspections like: boundary nodes, overall nodes & elements etc.
```
import numpy as np
import gmsh_api
import gmsh_api.gmsh as gmsh
import gmsh
import math
import os
import sys
gmsh.initialize()
gmsh.model.add("poly")
os.chdir("/Users/gregorhofmeister/Documents/Uni/MCI/SS 22/Masterarbeit/FEM Simulation/FE/Polygon_Test")
# Load a STEP file (using `importShapes' instead of `merge' allows to directly
# retrieve the tags of the highest dimensional imported entities):
path = os.path.dirname(os.path.abspath('Baugruppe_1.step'))
path
v = gmsh.model.occ.importShapes(os.path.join(path, 'Baugruppe_1.step'))
#s = gmsh.model.occ.importShapes(os.path.join('/Users/gregorhofmeister/Documents/Uni/MCI/SS 22/Masterarbeit/FEM Simulation/FE/Polygon_Test', 'Baugruppe_1.step'))
```
```
bodies = gmsh.model.occ.getEntities(3)
print(bodies)
gmsh.model.occ.synchronize()
steel = gmsh.model.addPhysicalGroup(3, [1], name = "steel")
wood = gmsh.model.addPhysicalGroup(3, [2], name = "wood")
```
```
xmin, ymin, zmin, xmax, ymax, zmax = gmsh.model.occ.getBoundingBox(
v[0][0], v[0][1])
xmin, ymin, zmin, xmax, ymax, zmax = (xmin, ymin-100, zmin-100, xmax, ymax+100, zmax+100)
N = 2 # Number of slices
dir = 'X' # Direction: 'X', 'Y' or 'Z'
surf = True # Keep only surfaces?
dx = (xmax - xmin)
dy = (ymax - ymin)
dz = (zmax - zmin)
L = dz if (dir == 'X') else dx
H = dz if (dir == 'Y') else dy
# Create the first cutting plane:
s = []
s.append((2, gmsh.model.occ.addRectangle(xmin, ymin, zmin, L, H)))
if dir == 'X':
gmsh.model.occ.rotate([s[0]], xmin, ymin, zmin, 0, 1, 0, -math.pi/2)
elif dir == 'Y':
gmsh.model.occ.rotate([s[0]], xmin, ymin, zmin, 1, 0, 0, math.pi/2)
tx = dx / N if (dir == 'X') else 0
ty = dy / N if (dir == 'Y') else 0
tz = dz / N if (dir == 'Z') else 0
gmsh.model.occ.translate([s[0]], tx, ty, tz)
# Create the other cutting planes:
for i in range(1, N-1):
s.extend(gmsh.model.occ.copy([s[0]]))
gmsh.model.occ.translate([s[-1]], i * tx, i * ty, i * tz)
# Fragment (i.e. intersect) the volume with all the cutting planes:
gmsh.model.occ.fragment(v, s)
# Now remove all the surfaces (and their bounding entities) that are not on the
# boundary of a volume, i.e. the parts of the cutting planes that "stick out" of
# the volume:
gmsh.model.occ.remove(gmsh.model.occ.getEntities(2), True)
gmsh.model.occ.synchronize()
if surf:
# If we want to only keep the surfaces, retrieve the surfaces in bounding
# boxes around the cutting planes...
eps = 1e-4
s = []
for i in range(1, N):
xx = xmin if (dir == 'X') else xmax
yy = ymin if (dir == 'Y') else ymax
zz = zmin if (dir == 'Z') else zmax
s.extend(gmsh.model.getEntitiesInBoundingBox(
xmin - eps + i * tx, ymin - eps + i * ty, zmin - eps + i * tz,
xx + eps + i * tx, yy + eps + i * ty, zz + eps + i * tz, 2))
# ...and remove all the other entities (here directly in the model, as we
# won't modify any OpenCASCADE entities later on):
dels = gmsh.model.getEntities(2)
for e in s:
dels.remove(e)
gmsh.model.removeEntities(gmsh.model.getEntities(3))
gmsh.model.removeEntities(dels)
gmsh.model.removeEntities(gmsh.model.getEntities(1))
gmsh.model.removeEntities(gmsh.model.getEntities(0))
# Finally, let's specify a global mesh size and mesh the partitioned model:
gmsh.option.setNumber("Mesh.MeshSizeMin", 3)
gmsh.option.setNumber("Mesh.MeshSizeMax", 3)
# quadratic elements
gmsh.option.setNumber("Mesh.ElementOrder", 2)
gmsh.model.mesh.generate(2)
gmsh.write("poly.msh")
entities = gmsh.model.getEntities(2)
```
Now the question is how do i find out which cross section belongs to which solid/material? For further calculations, I need to know for each cross-section which solids/materials occur at the intersection point X. Depending on the imported step.file, it could be that not all solids occur at every intersection.
Thanks a lot for the help in advance.
Greetings Gregorhttps://gitlab.onelab.info/gmsh/gmsh/-/issues/1565clip plane & cut plane2023-07-27T12:31:21ZHamidclip plane & cut planeHello Chirstophe,
I know gmsh has a cutplane plugin but do you possible to have a cutplane surface rendering when clipping (like in medit or paraview) ?
Thanks.Hello Chirstophe,
I know gmsh has a cutplane plugin but do you possible to have a cutplane surface rendering when clipping (like in medit or paraview) ?
Thanks.https://gitlab.onelab.info/gmsh/gmsh/-/issues/326closures not implemented for pyramids and prisms of order > 22018-07-29T11:19:22ZKoen Hillewaertclosures not implemented for pyramids and prisms of order > 2I see that the face closures are used essentially for optimizing high order meshes. Would it be possible to complete the missing functionality for the pyramid and prism up to order 4 ?I see that the face closures are used essentially for optimizing high order meshes. Would it be possible to complete the missing functionality for the pyramid and prism up to order 4 ?https://gitlab.onelab.info/gmsh/gmsh/-/issues/2815Coherence / BooleanFragments different outputs across Gmsh versions2024-02-25T16:24:55ZJoão PinheiroCoherence / BooleanFragments different outputs across Gmsh versionsHello,
I've been having issues with the newest version of Gmsh 4.12.2. When i input my .geo files i get weird geometries, which does not happen on the previous version 4.11.1. I can't figure out what I am doing wrong.
Attached are 2 simp...Hello,
I've been having issues with the newest version of Gmsh 4.12.2. When i input my .geo files i get weird geometries, which does not happen on the previous version 4.11.1. I can't figure out what I am doing wrong.
Attached are 2 simple files that hopefully show the problem.
On the first file error_coherence, Volume in BoundingBox does not find any volumes.
On error_fragments, "Error : Could not fix wire in surface" when applying the second BooleanFragments.
[error_coherence.geo](/uploads/2c76b52558266f94e6bef9f8f43a6f10/error_coherence.geo)
[error_fragments.geo](/uploads/8d127e9eb400fb0388b6800e1895b132/error_fragments.geo)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2790Combining twisting and boolean operators2024-01-31T10:37:31ZSantiago GomezCombining twisting and boolean operatorsHi! I am interested in creating a geometry which combines boolean operations with twisting. Specifically, I want to create a geometry with a "twisted hole", which I would do by taking a boolean difference. However, I noticed that the twi...Hi! I am interested in creating a geometry which combines boolean operations with twisting. Specifically, I want to create a geometry with a "twisted hole", which I would do by taking a boolean difference. However, I noticed that the twist is only possible with the built in kernel, whereas the boolean operators are only available with occ. Is there anyway to combine these two kernels, or get one of the functions running with the other kernel?
Thanks in advance!https://gitlab.onelab.info/gmsh/gmsh/-/issues/2158Concave quadratic tetrahedra on cylinder2022-10-27T16:12:17ZFrançois FraysseConcave quadratic tetrahedra on cylinderHello,
I am using linux binaries version 4.10.4. I am trying to mesh a pipe between two tanks with a very coarse mesh composed of quadratic tetras. After applying the optimizer, this is the result I get in the attached pictures.
![curv1...Hello,
I am using linux binaries version 4.10.4. I am trying to mesh a pipe between two tanks with a very coarse mesh composed of quadratic tetras. After applying the optimizer, this is the result I get in the attached pictures.
![curv1](/uploads/83fa3a559fa4c3e75b6b5f517ea00bfa/curv1.png)![curv2](/uploads/3f731a358ca94a3a5298ce8fc1a750c4/curv2.png)
I understand this may be challenging given a so coarse mesh but the optimizer does not return an error, although it looks to me as a failure. This issue is sometimes tricky to catch visually on complex domains. Is there a way to prevent this kind of results ?
I add the geo file used for this configuration.[test_curv.geo](/uploads/1815075f58a9f3886723a3673b37649b/test_curv.geo)
Thanks for you help and commentshttps://gitlab.onelab.info/gmsh/gmsh/-/issues/1689Console commands in Gmsh GUI2022-01-05T20:26:12ZBenjamin MenkuecConsole commands in Gmsh GUIIt would be nice if there would be an exposed console in the Gmsh GUI, ie a python debug console, so that one could try different commands and see their effect immediately.It would be nice if there would be an exposed console in the Gmsh GUI, ie a python debug console, so that one could try different commands and see their effect immediately.https://gitlab.onelab.info/gmsh/gmsh/-/issues/1708constrained delaunay from a set of points2023-09-02T22:24:45ZPA Beaufortconstrained delaunay from a set of pointsHi every1.
Is there a way to generate a constrained Delaunay triangulation (tetrahedrization) from a set of points and constrained curves (resp, surfaces)?
This [example](https://gitlab.onelab.info/gmsh/gmsh/blob/gmsh_4_9_3/examples/ap...Hi every1.
Is there a way to generate a constrained Delaunay triangulation (tetrahedrization) from a set of points and constrained curves (resp, surfaces)?
This [example](https://gitlab.onelab.info/gmsh/gmsh/blob/gmsh_4_9_3/examples/api/raw_tetrahedralization.py#L16) seems promising, but does it handle some constraints, if I define some edges (triangles) in addition to the points?