Hey Andy, Does the GradientOfUnstructuredDataSet filter produce a tensor from an input vector field? My testing in ParaView indicates it doesn't , but I'm not super familiar with it.
Thanks, Cory On Wed, Apr 18, 2018 at 11:14 AM, Caffagni, Andrea < [email protected]> wrote: > Dear Cory, > > This is the complete message that I have obtained like output: > > > > Traceback (most recent call last): > > File "<string>", line 22, in <module> > > File "<string>", line 14, in RequestData > > File "C:\Program Files\ParaView 5.3.0-Qt5-OpenGL2-Windows- > 64bit\bin\lib\site-packages\numpy\linalg\linalg.py", line 1095, in eig > > _assertRankAtLeast2(a) > > File "C:\Program Files\ParaView 5.3.0-Qt5-OpenGL2-Windows- > 64bit\bin\lib\site-packages\numpy\linalg\linalg.py", line 202, in > _assertRankAtLeast2 > > 'at least two-dimensional' % len(a.shape)) > > numpy.linalg.linalg.LinAlgError: 0-dimensional array given. Array must be > at least two-dimensional > > number of points 1 > > tensors[0:3] [VTKArray([[[ 8.34616375, -23.89611435, 2.00894117], > > [ 0.67995167, -23.11237526, 289.07113647], > > [ 27.26550674, -231.06710815, 16.90883636]]], dtype=float32)] > > > > > > From this, I have thought about the vtk array. > > > > Best regards, > > Andrea > > > > *From:* Cory Quammen [mailto:[email protected]] > *Sent:* Wednesday, April 18, 2018 15:11 > *To:* Caffagni, Andrea > *Cc:* [email protected] > *Subject:* Re: [Paraview] Problem with numpy array > > > > Hi Andrea, > > > > Besides the error message, what does your script print? > > > > Thanks, > > Cory > > > > > > On Wed, Apr 18, 2018 at 2:14 AM, Caffagni, Andrea < > [email protected]> wrote: > > Dear Users, > > I need to solve an eigenproblem in a point of a CFD field. > > Using the GradientOfUnstructuredDataSet I calculate the velocity gradient > of my velocity field, than I isolate may point of interest and I have > written the following ProgrammableFilter: > > > > *import* *numpy* *as* *np* > > *from* *numpy* *import* linalg *as* LA > > *from* *paraview* *import* numpy_support > > input = inputs[0] > > > > numPoints = input.GetNumberOfPoints() > > *print* "number of points",numPoints > > tensors = input.PointData['Gradients'] > > *print* "tensors[0:3]",tensors[0:3] > > eigenvalues,eigenvectors = LA.eig(tensors) > > *print* 'eigenvalues', eigenvalues[0] > > *print* 'eigenvectors',eigenvectors[0] > > output.PointData.append(eigenvalues,'eigenvalues') > > output.PointData.append(eigenvectors,'eigenvectors') > > > > I have obtained the following error message: > > > > Traceback (most recent call last): > > File "<string>", line 22, in <module> > > File "<string>", line 14, in RequestData > > File "C:\Program Files\ParaView 5.3.0-Qt5-OpenGL2-Windows- > 64bit\bin\lib\site-packages\numpy\linalg\linalg.py", line 1095, in eig > > _assertRankAtLeast2(a) > > File "C:\Program Files\ParaView 5.3.0-Qt5-OpenGL2-Windows- > 64bit\bin\lib\site-packages\numpy\linalg\linalg.py", line 202, in > _assertRankAtLeast2 > > 'at least two-dimensional' % len(a.shape)) > > numpy.linalg.linalg.LinAlgError: 0-dimensional array given. Array must be > at least two-dimensional > > > > > > I imagine that the problem is linked to the vtk array, I need to convert > my ‘Gradients’ matrix in a nparray? > > Someone have the same problem? > > > > Best Regards, > > Andrea > > > > > > ____________________________________________________________ > _____________________________________________________ > > > > Questo messaggio è da intendersi esclusivamente ad uso del destinatario e > può contenere informazioni che sono di natura privilegiata, confidenziale o > non divulgabile secondo le leggi vigenti. Se il lettore del presente > messaggio non è il destinatario designato, o il dipendente/agente > responsabile per la consegna del messaggio al destinatario designato, si > informa che ogni disseminazione, distribuzione o copiatura di questa > comunicazione è vietata anche ai sensi della normativa vigente in materia > di protezione dei dati personali. Se avete ricevuto questo messaggio per > errore, vi preghiamo di notificarcelo immediatamente a mezzo e-mail di > risposta e successivamente di procedere alla cancellazione di questa e-mail > e relativi allegati dal vostro sistema. > > ____________________________________________________________ > _____________________________________________________ > > > > This message is intended only for the use of the addressee and may contain > information that is privileged, confidential and exempt from disclosure > under applicable law. If the reader of this message is not the intended > recipient, or the employee or agent responsible for delivering the message > to the intended recipient, you are hereby notified that any dissemination, > distribution or copying of this communication is prohibited under the > applicable data protection law. If you have received this e-mail by > mistake, please notify us immediately by return e-mail and delete this > e-mail and all attachments from your system. > > > > ____________________________________________________________ > _____________________________________________________ > > > _______________________________________________ > 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 > > Search the list archives at: http://markmail.org/search/?q=ParaView > > Follow this link to subscribe/unsubscribe: > https://public.kitware.com/mailman/listinfo/paraview > > > > > > -- > > Cory Quammen > Staff R&D Engineer > Kitware, Inc. > > ____________________________________________________________ > _____________________________________________________ > > > > Questo messaggio è da intendersi esclusivamente ad uso del destinatario e > può contenere informazioni che sono di natura privilegiata, confidenziale o > non divulgabile secondo le leggi vigenti. Se il lettore del presente > messaggio non è il destinatario designato, o il dipendente/agente > responsabile per la consegna del messaggio al destinatario designato, si > informa che ogni disseminazione, distribuzione o copiatura di questa > comunicazione è vietata anche ai sensi della normativa vigente in materia > di protezione dei dati personali. Se avete ricevuto questo messaggio per > errore, vi preghiamo di notificarcelo immediatamente a mezzo e-mail di > risposta e successivamente di procedere alla cancellazione di questa e-mail > e relativi allegati dal vostro sistema. > > ____________________________________________________________ > _____________________________________________________ > > > > This message is intended only for the use of the addressee and may contain > information that is privileged, confidential and exempt from disclosure > under applicable law. If the reader of this message is not the intended > recipient, or the employee or agent responsible for delivering the message > to the intended recipient, you are hereby notified that any dissemination, > distribution or copying of this communication is prohibited under the > applicable data protection law. If you have received this e-mail by > mistake, please notify us immediately by return e-mail and delete this > e-mail and all attachments from your system. > > > > ____________________________________________________________ > _____________________________________________________ > -- Cory Quammen Staff R&D Engineer Kitware, Inc.
_______________________________________________ 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 Search the list archives at: http://markmail.org/search/?q=ParaView Follow this link to subscribe/unsubscribe: https://public.kitware.com/mailman/listinfo/paraview
