@classmethod
def from_vtk_cell(cls, vtk_cell, points, data=None):
    """Generate a hexahedron from VTK cells."""
    kwargs = {}
    if data is not None:
        kwargs['data'] = data
    vertices = []
    point_ids = vtk_cell.GetPointIds()
    for i_point in range(vtk_cell.GetNumberOfPoints()):
        vertices.append(points[point_ids.GetId(i_point)])
    return cls(vertices=vertices, **kwargs)
