# How to determine gauss point coordinates in a mesh

**URL:** <https://fenicsproject.discourse.group/t/how-to-determine-gauss-point-coordinates-in-a-mesh/5027>\
**Category:** mesh\
**Created:** [February 3, 2021, 2:28pm UTC](https://fenicsproject.discourse.group/t/how-to-determine-gauss-point-coordinates-in-a-mesh/5027 "2021-02-03T14:28:52Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![jzghal](https://avatars.discourse-cdn.com/v4/letter/j/7ba0ec/32.png) [@jzghal](https://fenicsproject.discourse.group/u/jzghal)\
**Post date:** [February 3, 2021, 2:28pm UTC](https://fenicsproject.discourse.group/t/how-to-determine-gauss-point-coordinates-in-a-mesh/5027/1 "2021-02-03T14:28:53Z")

</div>

Hello fenics users,  
I need to o determine gauss points coordinates of a mesh.  
I’m a new fenics user, can any one helps me.  
I will be grateful  
Jihed ZGHAL

---

<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:** [February 3, 2021, 2:50pm UTC](https://fenicsproject.discourse.group/t/how-to-determine-gauss-point-coordinates-in-a-mesh/5027/2 "2021-02-03T14:50:43Z")

</div>

Maybe your are looking for this: [How to get quadrature points and their weights? - #2 by bhaveshshrimali](https://fenicsproject.discourse.group/t/how-to-get-quadrature-points-and-their-weights/3523/2) ?

---

<div class="post-metadata">

**Author:** ![jzghal](https://avatars.discourse-cdn.com/v4/letter/j/7ba0ec/32.png) [@jzghal](https://fenicsproject.discourse.group/u/jzghal)\
**Post date:** [February 3, 2021, 4:13pm UTC](https://fenicsproject.discourse.group/t/how-to-determine-gauss-point-coordinates-in-a-mesh/5027/3 "2021-02-03T16:13:21Z")

</div>

Thank you for your quick reply. This exactly what I’m looking for. Based on the referred example I think that I can’t solve my problem.  
I will explain more precisely my need: I generate this mesh as follow and I want to get the quadrature points :

> print ( “Mesh Generation …” )
> 
> # Structure dimension [mm]
> 
> L = 50  
> W = 10
> 
> # Mesh size [mm]
> 
> h = 1  
> Nx = int( L/h ); print (" - number of element trough x direction = “, Nx);  
> Ny = int( W/h ); print (” - number of element trough y direction = ", Ny);  
> mesh = RectangleMesh(Point(-1\*L, 0.), Point(L, W), Nx, Ny)  
> V = FunctionSpace(mesh, “CG”, 1)

Thank you for your help.

---

<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:** [February 3, 2021, 4:42pm UTC](https://fenicsproject.discourse.group/t/how-to-determine-gauss-point-coordinates-in-a-mesh/5027/4 "2021-02-03T16:42:01Z")

</div>

Do you want the physical quadrature points for each cell in your function space? These are automatically generated in ffc in fenics.

If you want to access them you could create a quadrature function space and tabulate the dof coordinates:

```auto
from dolfin import *

mesh = UnitSquareMesh(1, 1)
el = FiniteElement("Quadrature", mesh.ufl_cell(), degree=2, quad_scheme="default")
V = FunctionSpace(mesh, el)
print(V.tabulate_dof_coordinates())

```

returning

```auto
[[0.83333333 0.66666667]
 [0.33333333 0.16666667]
 [0.83333333 0.16666667]
 [0.66666667 0.83333333]
 [0.16666667 0.33333333]
 [0.16666667 0.83333333]]

```

---

<div class="post-metadata">

**Author:** ![jzghal](https://avatars.discourse-cdn.com/v4/letter/j/7ba0ec/32.png) [@jzghal](https://fenicsproject.discourse.group/u/jzghal)\
**Post date:** [February 3, 2021, 6:06pm UTC](https://fenicsproject.discourse.group/t/how-to-determine-gauss-point-coordinates-in-a-mesh/5027/5 "2021-02-03T18:06:55Z")

</div>

> [@dokken](#):
>
> `FunctionSpace`

thank you it’s exactly what I need.