[Paraview] Fwd: Paraview/Python script speed extracting data

[Thomas, Christian]

Hi all,

I have sorted out my bottleneck/time processing speed issue. I'm sure the following method\script is well known to many python/paraview users.
I had to include

delete (source_name) and del source_ name, once each source was no longer needed. This was particularly useful when running loops around several filters.
Example: delete(T1), del T1

I also included the clip filter several times to reduce the size of the mesh\data file.

C1 = Clip(Input = T1)
C1.ClipType = "Plane"
C1.ClipType.Origin = [0.0, 0.0, ZMin[0]]
C1.ClipType.Normal = [0.0, 0.0, 1.0]

Regards
Christian



Dr. Christian Thomas
Postdoctoral Research Associate

Department of Mathematics
Room 737, Huxley Building
180, Queen’s Gate
Imperial College London
SW7 2AZ

c.thomas at imperial.ac.uk<mailto:c.thomas at imperial.ac.uk>







Begin forwarded message:

Subject: Paraview/Python script speed extracting data
Date: 20 February 2015 13:05:48 GMT
To: <paraview at paraview.org<mailto:paraview at paraview.org>>

Hi

I'm writing a python script that uses several paraview filters. I then extract data from the final filter for further analysis in python. However, I'm having a problem with the speed of the data extraction (can take several minutes). I have tried writing to a file (createwriter) and then reading the file. I have also tried servermanager.Fetch() and GetPoints… Although I am able to get data with both methods, they are both very slow. I was wondering if there was anything obvious in my paraview/python script that was slowing the process or if anyone has a suggestion at how I may speed things up.
(I realise this could just be down to the size of my initial input, for small data files Mesh - 15MB, solution - 50MB).

The script is below. Any suggestions would be greatly appreciated.

Regards
Christian



TBR = TAUMeshReader( FileName=Dir+'Mesh' )
TBR.SurfaceMeshOnly = 0
TBR.SolutionFileName = Dir+'solution'
MB = MergeBlocks(Input = TBR)
T1 = Transform(Input = MB)
T1.Transform.Rotate = [-float(Rot/2.0), 0.0, Sweep]
ISO = IsoVolume(Input = T1)
ISO.ThresholdRange = [MIN_Tot_Pre, BLE_Tot_Pre]
ISO.InputScalars = ['POINTS','total_pressure']
SL = Slice(Input = T1)
SL.SliceType="Plane"
SL.SliceOffsetValues = n
SL.SliceType.Origin = Inboard
SL.SliceType.Normal = Normal
Con = Contour(Input = SL)
Con.PointMergeMethod = "Uniform Binning"
Con.ContourBy = ['POINTS', 'total_pressure']
Con.Isosurfaces = [BLE_Tot_Pre]
Pass = PassArrays(Con)
Pass.PointDataArrays = ['wall_distance']
h=CreateWriter(Dir+'Edge.csv',Pass)
h.FieldAssociation = "Points"
h.Precision = 15
h.UpdatePipeline()

….read Edge.csv file for further analysis










-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20150226/35ff6aac/attachment.html>