[Paraview] Mixed Formulation FEM in Paraview

2014-12-09 Thread Teo Ioannis
I would like to ask if it is possible to display scalar, vector or tensor
values at distinct nodes in a finite element mesh.

To make things clearer, there exist finite element formulations where not
all nodes have the same number of degrees of freedom. Consider the case of
a 20 node hexahedral element. In a so called u-p formulation all 20 nodes
have 3 displacement degrees of freedom (one for each x, y, z axis) but only
the 8 corner nodes have pressure degrees of freedom.

For this particular case, i know how to display the vector of displacements
for all the nodes. But since only the corner nodes of each hexahedron have
pressure degrees of freedom how can I instruct paraview to 'see' that only
these nodes have the pressure scalar values and make the correct color
interpolation?

Thank you very much.
Theofilos Manitaras
___
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


Re: [Paraview] Mixed Formulation FEM in Paraview

2014-12-09 Thread Andy Bauer
Hi,

Currently ParaView assumes that all field variables have the same
interpolation/shape functions for point data. You would have to either
throw out the higher order interpolation/shape functions or compute the
dofs for them for the lower order field approximation.

Regards,
Andy

On Tue, Dec 9, 2014 at 9:55 AM, Teo Ioannis teoj...@gmail.com wrote:

 I would like to ask if it is possible to display scalar, vector or tensor
 values at distinct nodes in a finite element mesh.

 To make things clearer, there exist finite element formulations where not
 all nodes have the same number of degrees of freedom. Consider the case of
 a 20 node hexahedral element. In a so called u-p formulation all 20 nodes
 have 3 displacement degrees of freedom (one for each x, y, z axis) but only
 the 8 corner nodes have pressure degrees of freedom.

 For this particular case, i know how to display the vector of
 displacements for all the nodes. But since only the corner nodes of each
 hexahedron have pressure degrees of freedom how can I instruct paraview to
 'see' that only these nodes have the pressure scalar values and make the
 correct color interpolation?

 Thank you very much.
 Theofilos Manitaras

 ___
 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


___
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


[Paraview] RGB color coding

2014-12-09 Thread Martin Genet

Dear ParaView people,

I'm trying to make a periodic but non-symmetric color map. I tried the 
attached xml definition, but it seems that between red and magenta, 
ParaView interpolates in the other direction, see attached png 
screenshot. Is that the expected behavior? Is there a way to interpolate 
differently? Thanks!


Martin

ColorMap name=test space=test
  Point x=0.00 o=1 r=0 g=0 b=1/
  Point x=0.17 o=1 r=0 g=1 b=1/
  Point x=0.33 o=1 r=0 g=1 b=0/
  Point x=0.50 o=1 r=1 g=1 b=0/
  Point x=0.67 o=1 r=1 g=0 b=0/
  Point x=0.83 o=1 r=1 g=0 b=1/
  Point x=1.00 o=1 r=0 g=0 b=1/
/ColorMap___
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


Re: [Paraview] [EXTERNAL] Re: Streaklines

2014-12-09 Thread Pradeep Jha
I am also trying to figure out the streakline filter. Still don't
understand how to make it function :)

On 5 December 2014 at 14:10, Scott, W Alan wasc...@sandia.gov wrote:

  Andy,

 Mind giving me an example of how to run this streakline filter?  I have
 never used it before.



 Thanks!



 Alan



 *From:* ParaView [mailto:paraview-boun...@paraview.org] *On Behalf Of *Andy
 Bauer
 *Sent:* Friday, December 05, 2014 12:06 PM
 *To:* Pradeep Jha
 *Cc:* paraview@paraview.org
 *Subject:* [EXTERNAL] Re: [Paraview] Streaklines



 ParaView does indeed have a StreakLine filter.



 On Fri, Dec 5, 2014 at 12:30 PM, Pradeep Jha pradeep.kumar@gmail.com
 wrote:

 I have a unsteady 2D contour plot of a flow past a cylinder. Is it
 possible to plot the streaklines of a flow using Paraview? My solution file
 only contains the following data:

 U,V,W velocity and Pressure.


 ___
 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



___
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


[Paraview] Starting a batch animation in the middle -- fail

2014-12-09 Thread Bill Sherman

Hello,

I have been successfully using pvbatch to generate animated image
sequences of state files created with the ParaView GUI.  However,
I am not able to begin an animation in the middle -- which I
desperately need to do in order to parallelize my renderings.

I mentioned this before in July (and again in person at SC'14),
but haven't yet reached a solution.

The documentation for the WriteAnimation() method in simple.py
indicates that this is possible, but fails in two ways.

1) The WriteAnimation() method doesn't actually implement the two
most important parameters needed: SetStartFileCount() and
SetPlaybackTimeWindow().

2) Even doing the operation manually, SetPlaybackTimeWindow()
affects the time of time-varying data, but does not affect the
time of animation cues such as the cameras!  or translation
values, etc used to move objects around.

I've also experimented a lot with setting various values of the
GetAnimationScene() -- such as: StartTime, Duration (which seems
to have no affect whatsoever), and AnimationTime.


I created a tar file with an example state file and batch rendering
script:
http://www.freevr.org/Downloads/pvanimationtest.tar.gz

The scene has three time-varying elements:
* Annotated Time -- works
* Time-varying data (meshes of digits 0-9) -- works
* An arrow moving through a sphere -- fails

NOTE: camera moves also fail, basically in the same way as the
arrow.

So how the arrow fails is that even when I begin the animation
process in the middle, the arrow (and any camera moves) always
start as though it were the first frame.  By the time it gets
to animation-time 5.0, the arrow should be through the sphere,
but instead it appears to the left.


To run a test, just do:
% pvbatch 3Drender.py

This will print out some information, and may print warnings about
bad interpolation -- these are fine because it's a result of the
interpolating between polygonal objects (the digits) that aren't
intended to be morphed.  The effect is sufficient that it works
to tell me that data-time is working.

Newly rendered frames go into the Frames directory.  There are
other directories of the form Frames-exN, which are various
experiments I've run to figure out the problem.

The state file that I am using is arrowspheredatai.pvsm (which is
loaded in the Python script).


This is a most desperate hour.

Please help me O PV wizards, you are my only hope.

Thanks in abundance,
Bill

--
Bill Sherman
Sr. Technology Advisor
Advanced Visualization Lab
Pervasive Technology Inst
Indiana University
sherm...@indiana.edu
___
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


Re: [Paraview] Mixed Formulation FEM in Paraview

2014-12-09 Thread Samuel Key

  
  
Theofilos,

These may not be the most elegant solutions, but they are simple and
will work very nicely.

ParaView is capable of reading multiple datum sets and displaying
the results either in the same window or in side-by-side windows
(and, probably both scenarios at the same time).

(1) Write out two separate results datum sets that are independently
displayable. The first results datum set with mid-edge nodal points
but without pressure at the vertex nodal points, and the second
results datum set without the mid-edge nodal points but with
pressure at the vertex nodal points. 

(2) When writing out a results datum set, generate a pressure at the
mid-edge nodal points by averaging the pressure found at the edge's
end-points.

(3) Suggestion #2 can be extended. Since ParaView, at a fundamental
level, was designed to display either "centered" cell-values and
point-values, and you very likely have 20-node quadratic(?)
hexahedrons with cell-values "centered" at a 2x2x2 Gaussian
quadrature in which you are interested in seeing results, a painless
way to get a single, simple displayable results datum set is, at the
time you write out a given 20-node brick, break the 20-node
quadratic brick and its cell data into eight 8-node linear bricks.
(A similar effort is required when writing out point coordinates and
point values.)

When you get to ParaView, there will be no visual confusion as to
which Gauss point result is located where in the finite element.
While the display is not quadratic in its rendering, eight linear
bricks for each quadratic brick is a pretty decent representation.
It will get you started using ParaView.

Quite possibly someone else can explain how ParaView renders
quadratic data and elements.

Samuel Key
FMA Development, LLC
1005 39th Ave NE
Great Falls, Montana 59404



On 12/9/2014 7:55 AM, Teo Ioannis
  wrote:


  I would like to ask if it is possible to display
scalar, vector or tensor values at distinct nodes in a finite
element mesh.


To make things clearer, there exist finite element
  formulations where not all nodes have the same number of
  degrees of freedom. Consider the case of a 20 node hexahedral
  element. In a so called u-p formulation all 20 nodes have 3
  displacement degrees of freedom (one for each x, y, z axis)
  but only the 8 corner nodes have pressure degrees of freedom.


For this particular case, i know how to display the vector
  of displacements for all the nodes. But since only the corner
  nodes of each hexahedron have pressure degrees of freedom how
  can I instruct paraview to 'see' that only these nodes have
  the pressure scalar values and make the correct color
  interpolation?


Thank you very much.
Theofilos Manitaras
  
  
  
  
  ___
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



  

___
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


Re: [Paraview] Axis Labels in superscript scientific notation

2014-12-09 Thread Cory Quammen
Hi Venkatt,

As far as I know, scientific notation with the exponent in superscipt
is not currently possible.

Cory

On Thu, Dec 4, 2014 at 10:31 PM, Venkattraman A vayyasw...@ucmerced.edu wrote:
 I was wondering if there was a way to use the traditional scientific
 notation (1x10^16) instead of the E notation (1E16) in the axis labels (for
 example in a simple XY plot). I see there are three formats Mixed,
 Scientific and Fixed that Paraview allows but a lot of journals do not
 accept the E notation which according to Paraview is the Scientific
 notation.

 Regards,
 Venkatt

 --
 Venkattraman A
 Assistant Professor
 School of Engineering SE1 346
 University of California Merced
 5200 N. Lake Rd
 Merced CA 95343
 Phone: +1 (209) 228 2359
 Email: vayyasw...@ucmerced.edu

 ___
 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
___
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


Re: [Paraview] [EXTERNAL] Re: Axis Labels in superscript scientific notation

2014-12-09 Thread Scott, W Alan
It doesn't sound like a bad idea, however.  Mind writing up a feature request?

Thanks,

Alan

-Original Message-
From: ParaView [mailto:paraview-boun...@paraview.org] On Behalf Of Cory Quammen
Sent: Tuesday, December 09, 2014 12:23 PM
To: Venkattraman A
Cc: ParaView
Subject: [EXTERNAL] Re: [Paraview] Axis Labels in superscript scientific 
notation

Hi Venkatt,

As far as I know, scientific notation with the exponent in superscipt is not 
currently possible.

Cory

On Thu, Dec 4, 2014 at 10:31 PM, Venkattraman A vayyasw...@ucmerced.edu wrote:
 I was wondering if there was a way to use the traditional scientific 
 notation (1x10^16) instead of the E notation (1E16) in the axis labels 
 (for example in a simple XY plot). I see there are three formats 
 Mixed, Scientific and Fixed that Paraview allows but a lot of journals 
 do not accept the E notation which according to Paraview is the 
 Scientific notation.

 Regards,
 Venkatt

 --
 Venkattraman A
 Assistant Professor
 School of Engineering SE1 346
 University of California Merced
 5200 N. Lake Rd
 Merced CA 95343
 Phone: +1 (209) 228 2359
 Email: vayyasw...@ucmerced.edu

 ___
 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
___
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
___
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


Re: [Paraview] [EXTERNAL] Re: Streaklines

2014-12-09 Thread Andy Bauer
Hi,

The streakline filter requires a time dependent input as well as a vector
field input (probably point data required but didn't verify this). After
that, an input needs to be created for the seeds of the streaklines (i.e.
the location that the streaklines go through). Select the time dependent
input and then create the StreakLine filter. Then a dialog box will come up
to set the seed source (Input should be set to the time dependent input and
Seed Source should be set to the source to use for the streaklines). The
seeds are the points in the Seed Source. Next, in the Properties panel
you'll probably want to set Termination Time to something other than 0.0
since the streakline filter will integrate forward in time to that value.
Thus, if it's 0.0 and the first time step is at 0.0 it will seem like the
filter didnt' do anything. The Force Reinjection Every NSteps is
essentially the resolution of the streakline.

Note that PV's implementation of streaklines only does forward integration
in time.

Regards,
Andy

On Tue, Dec 9, 2014 at 10:58 AM, Pradeep Jha pradeep.kumar@gmail.com
wrote:

 I am also trying to figure out the streakline filter. Still don't
 understand how to make it function :)

 On 5 December 2014 at 14:10, Scott, W Alan wasc...@sandia.gov wrote:

  Andy,

 Mind giving me an example of how to run this streakline filter?  I have
 never used it before.



 Thanks!



 Alan



 *From:* ParaView [mailto:paraview-boun...@paraview.org] *On Behalf Of *Andy
 Bauer
 *Sent:* Friday, December 05, 2014 12:06 PM
 *To:* Pradeep Jha
 *Cc:* paraview@paraview.org
 *Subject:* [EXTERNAL] Re: [Paraview] Streaklines



 ParaView does indeed have a StreakLine filter.



 On Fri, Dec 5, 2014 at 12:30 PM, Pradeep Jha pradeep.kumar@gmail.com
 wrote:

 I have a unsteady 2D contour plot of a flow past a cylinder. Is it
 possible to plot the streaklines of a flow using Paraview? My solution file
 only contains the following data:

 U,V,W velocity and Pressure.


 ___
 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





___
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


Re: [Paraview] [EXTERNAL] Re: Axis Labels in superscript scientific notation

2014-12-09 Thread Venkattraman A

Thanks Cory and Alan. I will put in a Feature request online.

Venkatt

On 12/09/2014 11:36 AM, Scott, W Alan wrote:

It doesn't sound like a bad idea, however.  Mind writing up a feature request?

Thanks,

Alan

-Original Message-
From: ParaView [mailto:paraview-boun...@paraview.org] On Behalf Of Cory Quammen
Sent: Tuesday, December 09, 2014 12:23 PM
To: Venkattraman A
Cc: ParaView
Subject: [EXTERNAL] Re: [Paraview] Axis Labels in superscript scientific 
notation

Hi Venkatt,

As far as I know, scientific notation with the exponent in superscipt is not 
currently possible.

Cory

On Thu, Dec 4, 2014 at 10:31 PM, Venkattraman A vayyasw...@ucmerced.edu wrote:

I was wondering if there was a way to use the traditional scientific
notation (1x10^16) instead of the E notation (1E16) in the axis labels
(for example in a simple XY plot). I see there are three formats
Mixed, Scientific and Fixed that Paraview allows but a lot of journals
do not accept the E notation which according to Paraview is the
Scientific notation.

Regards,
Venkatt

--
Venkattraman A
Assistant Professor
School of Engineering SE1 346
University of California Merced
5200 N. Lake Rd
Merced CA 95343
Phone: +1 (209) 228 2359
Email: vayyasw...@ucmerced.edu

___
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

___
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


--
Venkattraman A
Assistant Professor
School of Engineering SE1 346
University of California Merced
5200 N. Lake Rd
Merced CA 95343
Phone: +1 (209) 228 2359
Email: vayyasw...@ucmerced.edu

___
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


[Paraview] animate parameter for python programmable filter

2014-12-09 Thread Ryan Abernathey
Hi,

I would like to animate a parameter with a python programmable filter using
keyframes. In this simple example, I just want to rescale a scalar variable
by a constant factor (alpha). But I want to animate this factor. Here is my
filter:

from paraview.vtk import dataset_adapter as DA
self.SetParameter('alpha', 0.1)
pdi = self.GetInputDataObject(0,0)
pdo = self.GetOutputDataObject(0)
pdo.CopyAttributes(pdi)
d = alpha * inputs[0].PointData['temperature']
arr = DA.numpyTovtkDataArray(d, 'temperature_scaled')
pdo.GetPointData().AddArray(arr)

This script works, but I can't set alpha using keyframes. It just doesn't
show up in the animation view.

Perhaps I have understood what the whole SetParameter framework does.
What is the point of setting these parameters if they can't be modified
from outside the filter?

Thanks in advance for your help.

-Ryan
___
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