Initializing a triangulated mesh

[sambaPython24]

Hey, I have two numpy arrays with the information about vertices and faces. (The mesh is closed).
How do I create a mesh from that in pygmsh in order to e.g. perform boolean operations with it?

import numpy as np
vertices_list = np.array([
    [-0.5, -0.5, -0.5],
    [0.5, -0.5, -0.5],
    [-0.5, 0.5, -0.5],
    [0.5, 0.5, -0.5],
    [-0.5, -0.5, 0.5],
    [0.5, -0.5, 0.5],
    [-0.5, 0.5, 0.5],
    [0.5, 0.5, 0.5]])

faces_list = np.array([
    [2, 1, 0],
    [1, 2, 3],
    [4, 2, 0],
    [2, 4, 6],
    [1, 4, 0],
    [4, 1, 5],
    [6, 5, 7],
    [5, 6, 4],
    [3, 6, 7],
    [6, 3, 2],
    [5, 3, 7],
    [3, 5, 1]])
``

[vargonis]

Hey, I'm wondering about that too. If I've understood things correctly, the internal representation of gmsh is not really a mesh, so it's perhaps not possible to directly "import" meshes into it. You have to create them from the ground up, by using add_point, add_line and so on, within a Geometry() context (either geo.Geometry() or occ.Geometry()). There are examples in the documentation.

But I'm a beginner too... I'd really like to know if the approach I mentioned is really the only way to go, when one already has a mesh. A related question is this: consider the following situation:

with pygmsh.occ.Geometry() as geom:
    # some lengthy creation commands, which end up with volumes A and B:
    A = ...
    B = ...
    geom.boolean_difference(A, B)
    mesh = geom.generate_mesh()

Is it possible to somehow have volumes A and B already precomputed? Say, in such a way that the following would work:

with pygmsh.occ.Geometry() as geom:
    A = geom.add_volume(...) # some adequate args here
    B = geom.add_volume(...)
    geom.boolean_difference(A, B)
    mesh = geom.generate_mesh()

Thanks for your help!