I found the following code to get quadrature points and weights in the standard simplex. But I want to get the quadrature points on mesh elements/cells directly. Is any way in Fenics to do this? What are the choices for scheme (what are the candidate input strings)? from ffc.fiatinterface import c…

Im not sure what you want with this question. In the other post ive already shown you how to get the quadrature points in any cell. The weights does not change from cell to cell.

See: [image] How to determine gauss point coordinates in a mesh mesh 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 spac…

Thanks! This code gives a list of integration points, but how to get the weights? Can I get the points and weights arranged by elements/cells? What is the default integration scheme?

[image] xi_wang: Thanks! This code gives a list of integration points, but how to get the weights? Please read the docstring of the function i referred you to [image] dokken: See: Bitbucket

Thanks! I have read the doc and want to use this function to get integration points and weights. The following code gives some error. How to construct FIAT cell with mesh? from FIAT import * cell = Cell(mesh, 0) print(create_quadrature(cell,2)) AttributeError: ‘dolfin.cpp.mesh.Cell’ object has no…

Use the ffc.fiatinterface: import ffc import dolfin as d mesh = d.UnitSquareMesh(10, 10) points, weights = ffc.fiatinterface.create_quadrature(mesh.cell_name(), degree=1, scheme="default") print(points, weights)

Thanks ! And how can I specify which mesh cell to be input?

Yes the weights do not change, but to align it with the nodes correctly need the information of integration scheme, the distribution, and the order of the integration nodes. If I align them manually, I am afraid of using the wrong information and get bugs that are hard to realize. Isn’t there a way…

How to get Gaussian integration point directly on a mesh element/cell?

mathematics

dokken September 24, 2023, 1:04pm 2

See:

See: Bitbucket

How to get the interpolation points ? And elvaute interpolation value and derivatives on integration points?

Topic		Replies	Views
Cannot align integratoin points and weights in a mesh cell General	6	363	September 25, 2023
How to determine gauss point coordinates in a mesh mesh	4	1478	February 3, 2021
Getting quadrature points and weights with dolfinx-0.6.0 General	1	306	September 23, 2023
How to get quadrature points and their weights?	6	1771	February 22, 2024
Mapped quadrature points, weights and solution at these quadrature points General	5	728	April 24, 2024

How to get Gaussian integration point directly on a mesh element/cell?

Related topics