InsightSoftwareConsortium / InsightSoftwareConsortium/itkwidgets

Support for high-order Lagrange elements

Open
#306 4 comments 0 reactions 0 assignees View on GitHub
Dominant language
Python
Stars
624
Forks
83
PR merge metrics
No merged PRs in 30d

Description

Hello,

Hope it is the right place to ask. I am trying to visualize ugrids with high-order elements, e.g.,
```
order = 2
n_points = (order + 1) * (order + 2) // 2
cell = vtk.vtkLagrangeTriangle()
cell.GetPointIds().SetNumberOfIds(n_points)
cell.GetPoints().SetNumberOfPoints(n_points)
cell.Initialize()
point = [0., 0., 0.]
barycentric_index = [0, 0, 0]
for k_point in range(n_points):
cell.GetPointIds().SetId(k_point, k_point)
cell.ToBarycentricIndex(k_point, barycentric_index)
for k_dim in range(2):
point[k_dim] = (barycentric_index[k_dim] / order)**2
cell.GetPoints().SetPoint(k_point, point[0], point[1], point[2])
cell_array = vtk.vtkCellArray()
cell_array.InsertNextCell(cell)
ugrid = vtk.vtkUnstructuredGrid()
ugrid.SetPoints(cell.GetPoints())
ugrid.SetCells(cell.GetCellType(), cell_array)
itkwidgets.view(geometries=[ugrid])
```
(inspired from https://blog.kitware.com/modeling-arbitrary-order-lagrange-finite-elements-in-the-visualization-toolkit), to no avail. Is there any support for high-order elements? If not, is it planned? Thanks!

Martin

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.