[Paraview] ParaView with VirtualGL

[Craig Michoski]

Hi Berk,

Thanks for the response.  I see I'm at the fount.  So, I'm using paraview
4.1.0, and loading everything in from vtu files.  The basic python filter
I'm trying to get running in parallel is this below, for a basic
integration along a dimension in two of the loaded fields (one a scalar
alpha_0, the other a two dimensional vector sigma_0_0__sigma_0_1).

----------------
import numpy

alpha0 = inputs[0].PointData['alpha_0']
alpha1 = inputs[0].PointData['alpha_0']
sigma0 = inputs[0].PointData['sigma_0_0__sigma_0_1']

pdi = self.GetInputDataObject(0,0)
bg_field = 1
pdi = self.GetInputDataObject(0,0)
pdo = self.GetOutputDataObject(0)
numPoints = pdi.GetNumberOfPoints()
import time

ind=self.GetInput()
iin=ind.GetInformation()
timestep = time.clock()

dx = 100.0
dy = 100.0

for (j,alpha) in enumerate(alpha0):
    x1,y1,z1 = pdi.GetPoint(j)[:3]
    x2,y2,z2 = pdi.GetPoint(j+1)[:3]
    if abs(x2-x1) <= dx and abs(x2-x1)>0:
        dx = x2 - x1
    if abs(y2-y1) <= dy and abs(y2-y1)>1e-6:
        dy = y2 - y1

counter = 0
last_step = -1.0
prev_x = 0.0
efold = numpy.zeros(numPoints)
y_first = -1.0

for (j,alpha) in enumerate(alpha0):
    x,y,z = pdi.GetPoint(j)
    if x<380.0 and x>70.0 and y==0.0:
        y_first = []
        for (i,alpha_1) in enumerate(alpha1):
            x1,y1,z1 = pdi.GetPoint(i)
            if x1 == x and y1 == y:
                for (k,sigmax) in enumerate(sigma0):
                    x2,y2,z2 = pdi.GetPoint(i)
                    if x2 == x1: #and y2 == y1:
                        if y2 not in y_first:
                            efold[k] +=  ( 1.0/exp(2.0*alpha_1) ) *
sigmax[0] * dy  #exp(alpha_1)*sigmax[0]*dy
                            y_first.append(y2)

                if x!=last_step:
                    counter = counter + 1
                    print "y_integrated_chi_at1(",floor(timestep),",",
counter,") = ", efold[k]/377.0, ";"
                    print "x(",counter,") = ", x, ";"
                    last_step = x
----------------

Craig


On Wed, Oct 15, 2014 at 5:38 PM, Berk Geveci <berk.geveci at kitware.com>
wrote:

> Hi Craig,
>
> Which version of ParaView are you using?
>
> What file format are you loading?
>
> Best,
> -berk
>
> On Wed, Oct 15, 2014 at 4:52 PM, Craig Michoski <michoski at gmail.com>
> wrote:
>
>> Ah, okay, I figured out part of my problem.  I had to connect to the
>> server and display the output on the server side, not just the client
>> side, which I guess is where the print commands redirects.  Now I have
>> several other problems.  First:
>>
>> 1) the print command outputs n_processor number of outputs.  Can I
>> make this just one value?  For example, if I am summing a value over a
>> mesh, I just want to total value output once over the entire mesh?  Is
>> this possible?
>>
>> 2)  I am getting an error when trying to simply get my loaded inputs.
>> When I do:
>>
>> print "Gets here!"
>>
>> alpha0 = inputs[0].PointData['alpha_0']
>>
>> Traceback (most recent call last):
>>   File "<string>", line 22, in <module>
>>   File "<string>", line 8, in RequestData
>> TypeError: 'NoneType' object has no attribute '__getitem__'
>>
>> I don't understand what is hppening here, can anyone help please?
>>
>> Craig
>> _______________________________________________
>> Powered by www.kitware.com
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>>
>> Please keep messages on-topic and check the ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20141015/c9f87944/attachment-0001.html>