How to visualize superquadric

In the meantime I managed to get some cool results using dump custom and my rigid reader in combination with a python programmable filter to map any geometry over the centre points dumped:

from math import sqrt
input1 = self.GetInputDataObject(0, 0)
input2 = self.GetInputDataObject(0, 1)
output = self.GetOutputDataObject(0)
newPoints = vtk.vtkPoints()
newcells = vtk.vtkCellArray()
#
P=input1.GetPointData()
B=input2.GetPointData()
BODIES=input1.GetNumberOfPoints()
EDGES=input2.GetNumberOfPoints()
SURFACES=input2.GetNumberOfCells()
#
outputarrayID = vtk.vtkDoubleArray()
outputarrayID.SetName("ID")
outputarrayID.SetNumberOfTuples(BODIES*SURFACES)
#
outputarrayT = vtk.vtkDoubleArray()
outputarrayT.SetName("type")
outputarrayT.SetNumberOfTuples(BODIES*SURFACES)
#
outputarraySC = vtk.vtkDoubleArray()
outputarraySC.SetName("value")
outputarraySC.SetNumberOfTuples(BODIES*SURFACES)
#
outputarrayF = vtk.vtkDoubleArray()
outputarrayF.SetName("F")
outputarrayF.SetNumberOfTuples(BODIES*SURFACES)
#
outputarray_v = vtk.vtkDoubleArray()
outputarray_v.SetName("v")
outputarray_v.SetNumberOfTuples(BODIES*SURFACES)
k=0
c=0
for body in range(0,BODIES):
r=input1.GetPoint(body)
M=P.GetAbstractArray('M')
T11=M.GetComponent(body,0)
T12=M.GetComponent(body,1)
T13=M.GetComponent(body,2)
T21=M.GetComponent(body,3)
T22=M.GetComponent(body,4)
T23=M.GetComponent(body,5)
T31=M.GetComponent(body,6)
T32=M.GetComponent(body,7)
T33=M.GetComponent(body,8)
type=P.GetArray('type').GetValue(body)-1
CLUMP_ID=P.GetArray("c_id").GetValue(body);
CLUMP_SC=P.GetArray("value").GetValue(body);
CLUMP_FORCEx=P.GetAbstractArray("F").GetComponent(body,0)
CLUMP_FORCEy=P.GetAbstractArray("F").GetComponent(body,1)
CLUMP_FORCEz=P.GetAbstractArray("F").GetComponent(body,2)
CLUMP_FORCE=sqrt(CLUMP_FORCEx*CLUMP_FORCEx+CLUMP_FORCEy*CLUMP_FORCEy+CLUMP_FORCEz*CLUMP_FORCEz)
CLUMP_VELx=P.GetAbstractArray("v").GetComponent(body,0)
CLUMP_VELy=P.GetAbstractArray("v").GetComponent(body,1)
CLUMP_VELz=P.GetAbstractArray("v").GetComponent(body,2)
CLUMP_VEL=sqrt(CLUMP_VELx*CLUMP_VELx+CLUMP_VELy*CLUMP_VELy+CLUMP_VELz*CLUMP_VELz)
for surface in range(0,SURFACES):
outputarrayID.SetValue(k,CLUMP_ID)
outputarrayT.SetValue(k,type)
outputarrayF.SetValue(k,CLUMP_FORCE)
outputarraySC.SetValue(k,CLUMP_SC)
outputarray_v.SetValue(k,CLUMP_VEL)
k+=1
cell = input2.GetCell(surface)
COUNT=cell.GetNumberOfPoints()
for j in range(0,COUNT): #alle punkte der jeweiligen zelle
id=cell.GetPointId(j)
#print id
coord = input2.GetPoint(id)
x, y, z = coord[:3]
#print x,y,z
#print r
xnew=T11*x+T12*y+T13*z+r[0]
ynew=T21*x+T22*y+T23*z+r[1]
znew=T31*x+T32*y+T33*z+r[2]
newPoints.InsertPoint(c, xnew, ynew, znew)
c+=1
newcells.InsertNextCell(4)
newcells.InsertCellPoint(c-1)
newcells.InsertCellPoint(c-2)
newcells.InsertCellPoint(c-3)
newcells.InsertCellPoint(c-4)
output.SetPoints(newPoints)
output.SetCells(5,newcells)
output.GetCellData().AddArray(outputarrayID)
output.GetCellData().AddArray(outputarrayT)
output.GetCellData().AddArray(outputarrayF)
output.GetCellData().AddArray(outputarraySC)
output.GetCellData().AddArray(outputarray_v)

Apped this to the dump AND a source of type box (or any other PV primitive with cell type quad)

Hi,
the python code is not meant as a macro. Its a python-programmable filter.
You'll need the dump (rigid style) and a Source-> Superquadric. Than attach the programmable filter to both entries in the pipeline browser.
See this screenshot: https://www.cfdem.com/system/files/settings4pv_sq.png
You can Download the filter as custom xml here: https://www.cfdem.com/system/files/304_transform_sq_with_colors.xml_.txt save this file as xml and use PV Tools->Manage Custom Filter [import].

Good luck :-)

One Note: roundness=2/blockiness and y/z scale is interchanged in Paraview 5.0.0 (seems correct in 5.4.1)

[1] https://www.cfdem.com/forums/tutorial-4-superquadric-angle-repose-repres...

Hi Thomas,
could you test if it works when you change

CLUMP_ID=P.GetArray("id").GetValue(body);

to

CLUMP_ID=P.GetArray("c_id").GetValue(body);

I changed the array name in the reader some times ago and forgot to change it in the python programable filter before upload.
Thanks,
Christian.

I do not use tensor glyph as I have no clue how to do the orientation
about
And the second one is, if i change the blockiness of the superquadric particles, is there an easy way to get the change into the paraview cisualisation?
In the Superquadric source roundness is 2/blockiness in liggghts.
(disable torodial and do not change resolution)
For some reason the y,z size of SQ in PV seems interchanged with ey/ez in liggghts