<div dir="ltr">Hi Christian,<div><br></div><div>Did you identify where the performance bottleneck occurs? Does it perform well if you take out the writing or Fetch?</div><div><br></div><div>Best,</div><div>-berk</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 20, 2015 at 8:05 AM, Thomas, Christian <span dir="ltr"><<a href="mailto:c.thomas@imperial.ac.uk" target="_blank">c.thomas@imperial.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<br>
<br>
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.<br>
(I realise this could just be down to the size of my initial input, for small data files Mesh - 15MB, solution - 50MB).<br>
<br>
The script is below. Any suggestions would be greatly appreciated.<br>
<br>
Regards<br>
Christian<br>
<br>
<br>
<br>
        TBR = TAUMeshReader( FileName=Dir+'Mesh' )<br>
        TBR.SurfaceMeshOnly = 0<br>
        TBR.SolutionFileName = Dir+'solution'<br>
        MB = MergeBlocks(Input = TBR)<br>
        T1 = Transform(Input = MB)<br>
        T1.Transform.Rotate = [-float(Rot/2.0), 0.0, Sweep]<br>
        ISO = IsoVolume(Input = T1)<br>
        ISO.ThresholdRange = [MIN_Tot_Pre, BLE_Tot_Pre]<br>
        ISO.InputScalars = ['POINTS','total_pressure']<br>
        SL = Slice(Input = T1)<br>
        SL.SliceType="Plane"<br>
        SL.SliceOffsetValues = n<br>
        SL.SliceType.Origin = Inboard<br>
        SL.SliceType.Normal = Normal<br>
        Con = Contour(Input = SL)<br>
        Con.PointMergeMethod = "Uniform Binning"<br>
        Con.ContourBy = ['POINTS', 'total_pressure']<br>
        Con.Isosurfaces = [BLE_Tot_Pre]<br>
        Pass = PassArrays(Con)<br>
        Pass.PointDataArrays = ['wall_distance']<br>
        h=CreateWriter(Dir+'Edge.csv',Pass)<br>
        h.FieldAssociation = "Points"<br>
        h.Precision = 15<br>
        h.UpdatePipeline()<br>
<br>
        ….read Edge.csv file for further analysis<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" target="_blank">http://markmail.org/search/?q=ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/paraview" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><br>
</blockquote></div><br></div>