# Create mesh in fenicsx from point and element array

**URL:** <https://fenicsproject.discourse.group/t/create-mesh-in-fenicsx-from-point-and-element-array/9092>\
**Category:** Uncategorized\
**Created:** [August 23, 2022, 3:38pm UTC](https://fenicsproject.discourse.group/t/create-mesh-in-fenicsx-from-point-and-element-array/9092 "2022-08-23T15:38:17Z")\
**Posts on this page:** 1\
**Showing post:** 2

<div class="post-metadata">

**Author:** ![dokken](https://yyz2.discourse-cdn.com/free1/user_avatar/fenicsproject.discourse.group/dokken/32/1560_2.png) [@dokken](https://fenicsproject.discourse.group/u/dokken)\
**Post date:** [August 23, 2022, 7:09pm UTC](https://fenicsproject.discourse.group/t/create-mesh-in-fenicsx-from-point-and-element-array/9092/2 "2022-08-23T19:09:09Z")

</div>

Consider:

```python
import dolfinx
import ufl
from mpi4py import MPI
import numpy as np

gdim = 2
shape = "triangle"
degree = 1

cell = ufl.Cell(shape, geometric_dimension=gdim)
domain = ufl.Mesh(ufl.VectorElement("Lagrange", cell, degree))

x = np.array([[0.5, 1], [2, 1], [3, 1.5], [3.5, 2.5], [2.2, 2], [1, 2.2]])
cells = np.array([[0, 1, 5], [1, 4, 5], [1, 2, 4], [2, 3, 4]], dtype=np.int32)
mesh = dolfinx.mesh.create_mesh(MPI.COMM_WORLD, cells, x, domain)

with dolfinx.io.XDMFFile(MPI.COMM_WORLD, "mesh.xdmf", "w") as xdmf:
    xdmf.write_mesh(mesh)

```

---

_[View the full topic](https://fenicsproject.discourse.group/t/create-mesh-in-fenicsx-from-point-and-element-array/9092)._