problems field data unstructured reader

ParaView Support
1

Writing during a calculation demands much time, I am searching a for a way, for example, print time, average total energy, average overlap… In other words properties that applies for a system of particles without having to do a loop and print for every particle, which is unpractical.
so i tried using field data to plot a general characteristic that in this example is time.
ex:(without data)

# vtk DataFile Version 3.0
Particle Field
BINARY
DATASET UNSTRUCTURED_GRID
POINTS 2 float
POINT_DATA 2
SCALARS Particle_radius float
LOOKUP_TABLE default
VECTORS Particle_velocity float
FIELD Data_Properities 1
time 1 1 float

but when i run on paraview i receive

ERROR: In C:\bbd\6c27b535\build\superbuild\paraview\src\VTK\Common\DataModel\vtkDataSet.cxx, line 435
    vtkUnstructuredGrid (000001EC25153260): Point array time with 1 components, only has 1 tuples but there are 2 points

ERROR: In C:\bbd\6c27b535\build\superbuild\paraview\src\VTK\IO\Parallel\vtkPDataSetReader.cxx, line 736
vtkPDataSetReader (000001EC29C9BBB0): Attribute Mismatch.

ERROR: In C:\bbd\6c27b535\build\superbuild\paraview\src\VTK\Common\ExecutionModel\vtkExecutive.cxx, line 782
vtkPVCompositeDataPipeline (000001EC29E90440): Algorithm vtkFileSeriesReader(000001EC21E6C7A0) returned failure for request: vtkInformation (000001EC1C5CAFA0)
Debug: Off
Modified Time: 3152843
Reference Count: 1
Registered Events: (none)
Request: REQUEST_DATA
FORWARD_DIRECTION: 0
ALGORITHM_AFTER_FORWARD: 1
FROM_OUTPUT_PORT: 0
2

Store them in a field data instead of point data.

3

thank you, the solution was simply changing the order.
my particle analysis code generates a radius and coordinates in a unstructured grid. but i want to view it in paraview directly as spheres without a paraview filter. there is someway to write spheres with just radius and center to be read in paraview? i tried using unstructured grid in python… but no progress
in polydata when using glyph filter generates many points to describe a sphere. i want to skip this part, i was planning using fortran ->paraview not fortran ->vtk->paraview and also without generating n points in fortran for a sphere, only radius and coordinates.
#this python code is a search for a pattern to be implemented in fortran, in other words a pattern of vtk file (without calling vtk) which gets radius and coordinates to describe a sphere.

points = vtk.vtkPoints()
points.InsertNextPoint(0,0,0)
points.InsertNextPoint(1,1,1)
points.InsertNextPoint(2,2,2)

ug = vtk.vtkUnstructuredGrid()
ug.SetPoints(points)

sphereSource = vtk.vtkSphereSource()

glyph3D = vtk.vtkGlyph3D()
glyph3D.SetSourceConnection(sphereSource.GetOutputPort())
glyph3D.SetInputData(ug)
glyph3D.Update()

Writer = vtk.vtkUnstructuredGridWriter()
Writer.SetFileName('TestWritePLY.vtk')
Writer.SetInputConnection(glyph3D.GetOutputPort())
Writer.SetFileTypeToASCII()
Writer.Write()"
4

Use the point gausssian representation.

5

well aswell for describing a gaussian point i need a “thousand lines” of vtk coding. isnt in my plans use another render beyond paraview, i want to implement in my fortran code something that paraview understands as gaussian point/sphere without a special configuration on paraview.
isnt there a properity like point but for gaussian point?

#!/usr/bin/env python
import vtk
def main():
# Create points on a sphere

sphereSource = vtk.vtkSphereSource()
sphereSource.Update()
colors = vtk.vtkNamedColors()
polydata = vtk.vtkPolyData()
polydata.SetPoints(sphereSource.GetOutput().GetPoints())
splatter = vtk.vtkGaussianSplatter()
splatter.SetInputData(polydata)
splatter.SetSampleDimensions(50,50,50)
splatter.SetRadius(0.5)
splatter.ScalarWarpingOff()
surface = vtk.vtkContourFilter()
surface.SetInputConnection(splatter.GetOutputPort())
surface.SetValue(0,0.01)
plyWriter = vtk.vtkPLYWriter()
plyWriter.SetFileName('test')
plyWriter.SetInputConnection(surface.GetOutputPort())
plyWriter.Write()
# Create a mapper and actor
mapper = vtk.vtkPolyDataMapper()
mapper.SetInputConnection(surface.GetOutputPort())
actor =  vtk.vtkActor()
actor.SetMapper(mapper)
# Visualize
renderer = vtk.vtkRenderer()
renderWindow = vtk.vtkRenderWindow()
renderWindow.AddRenderer(renderer)
renderWindowInteractor = vtk.vtkRenderWindowInteractor()
renderWindowInteractor.SetRenderWindow(renderWindow)
renderer.AddActor(actor)
renderer.SetBackground(colors.GetColor3d("white"))
renderWindow.Render()
renderWindowInteractor.Start()

test.ply (395.3 KB)
if name == ‘main’:
main()
"

6

I’m confused. You are currently on the ParaView discourse but your code is pure VTK.

Are you working with ParaView or VTK ?

7

i am working with paraview. i searching for a pattern (vtk file) for gaussian points which is viable to implement in fortran without having to configurate paraview every time i initialize (point default pattern).

8

Still no idea what you mean.

Please state the following :
What are you trying to do ?
What did you do to achieve that ?
What is not working ?

9

What are you trying to do ?
i want to represent my points generated from my mechanical simulation as spheres/gaussian points directely in the vtk file generated from simulation.
in other words, without using gaussian point representation from paraview but as configuration of vtk file.
What did you do to achieve that ?
i tried generating vtk file from vtk with different configurations(polydata, unstructured reader, sphere, gaussian point) and by this pattern implement in my fortran code to generate not points but spheres or gaussian points.
What is not working ?
i didnt found a pattern, a configuration in vtk, which didnt extrapolate my ideia of representation of sphere( one coordinate + radius), all of them until now for a coordinate/point are generating more or less one thousand points. So if i have a simulation with 500.000points, i would have a vtk file with 500.000*1000points wicht would bring a heavy file and would demands signatvely a more time of processing from fortran.
example of vtk from my simulations, as a point. (without trying to turn these points into spheres)
results2000.vtk (99.9 KB)
i sent a 2d representation a 3d one was to heavy to be sent here(+400.000points)

10

i want to represent my points generated from my mechanical simulation as spheres/gaussian points directely in the vtk file generated from simulation.
in other words, without using gaussian point representation from paraview but as configuration of vtk file.

This is simply not supported. You could create sphere on each of your points but, as you pointed it yourself, this would be very heavy.

How do you plan to display your “configured” vtk file once it has been generated ?