Re: [Paraview] point to surface

[Moreland, Kenneth via ParaView]

Hadrien,

I’m not sure what you mean by the “area” of points. As I’m sure you know, 
particles are 0D elements and themselves have no area to speak of. As you said, 
you can use the glyph operation to create a surface around each point, but the 
surface area in that case is simply the surface area of the sphere size you 
chose times the number of points.

You’ve asked if you can project the points onto an existing mesh surface. I 
don’t know of an easy way to do that, but let’s put that problem aside for a 
moment. Assuming you could, what good would that do? Points on a surface are 
still just points. The surface area of your mesh (nasal cavity) will be the 
same surface area regardless of how many points are in it and where they are.

Can you more precisely describe what it is you are trying to measure?

-Ken

(By the way, the ParaView team has been transitioning from using this mailing 
list to using a tool called discourse to ask questions and archive results. In 
the future, I suggest posting questions to https://discourse.paraview.org/ 
rather on this list.)

From: ParaView [mailto:paraview-boun...@public.kitware.com] On Behalf Of 
hadrien calmet
Sent: Tuesday, September 25, 2018 8:04 AM
To: paraview@public.kitware.com
Subject: [EXTERNAL] [Paraview] point to surface


Hi every one
I have a list of point (coordinate x,y,z), using table to point I can visualize 
it on PV.
But I would like to know the total area of those points.
So I tried glyph (sphere, all the points) than integrate variable and I 
obtained the area.
But the area is depending of the radius that I put on glyph parameter.

those points are the deposition of particles in a complex geometry(human nasal 
cavity)
They are calculated when they intersect the surface mesh of the nasal cavity.

my question:
There is a way to "project" those points on the mesh surface (then filtering 
and integrate variable to calculate the total area) ?

thank you





Hadrien Calmet

Dpt. Computer Applications in Science and Engineering

Barcelona Supercomputing Center (BSC-CNS)

Edificio Nexus II - Planta 3

C/ JORDI GIRONA, 29

08034 Barcelona, Spain

Tel: +34 93 413 79 40

Fax: +34 93 413 7721

www: www.bsc.es


WARNING / LEGAL TEXT: This message is intended only for the use of the 
individual or entity to which it is addressed and may contain information which 
is privileged, confidential, proprietary, or exempt from disclosure under 
applicable law. If you are not the intended recipient or the person responsible 
for delivering the message to the intended recipient, you are strictly 
prohibited from disclosing, distributing, copying, or in any way using this 
message. If you have received this communication in error, please notify the 
sender and destroy and delete any copies you may have received.

http://www.bsc.es/disclaimer
___
Powered by www.kitware.com

ParaView discussion is moving! Please visit https://discourse.paraview.org/ for 
future posts.

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

Re: [Paraview] [Paraview-developers] time to kill the mailing lists?

[Moreland, Kenneth via ParaView]

When you close the mailing lists, would it be possible to set it up so that 
anyone who attempts to email to those addresses gets an email back with a clear 
message about the mailing lists being closed down and instructions on posting 
to discourse?

-Ken

-Original Message-
From: Paraview-developers 
[mailto:paraview-developers-boun...@public.kitware.com] On Behalf Of Utkarsh 
Ayachit
Sent: Wednesday, September 5, 2018 10:55 AM
To: benson.mu...@ut.ee
Cc: ParaView Developers ; ParaView 

Subject: [EXTERNAL] Re: [Paraview-developers] [Paraview] time to kill the 
mailing lists?

Thanks, Benson. I've cleaned up the webpage. Let me know if that's still 
confusing: https://www.paraview.org/community-support/
On Wed, Sep 5, 2018 at 12:14 PM Benson Muite  wrote:
>
> The current paraview website has a link to mailing lists:
>
> https://www.paraview.org/
>
> buried within that is a message that the lists have been deprecated 
> and a link to the discourse forum:
>
> https://discourse.paraview.org/
>
> It may be helpful to make the forum a little more prominent on the 
> Paraview website
>
>
> On 09/05/2018 06:57 PM, Utkarsh Ayachit wrote:
> > Folks,
> >
> > We have moved to Discourse (quite successfully). Given the latest 
> > round of phishing emails, with spam email masquerading as being sent 
> > by legitimate users, may be it's time to close the paraview mailing 
> > list? I expect we can shut down the mailing list with a default 
> > reply message instructing any sender to use discourse.paraview.org instead.
> >
> > Thoughts? Objections?
> >
> > Utkarsh
> > ___
> > Powered by www.kitware.com
> >
> > ParaView discussion is moving! Please visit https://discourse.paraview.org/ 
> > for future posts.
> >
> > 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
>
___
Powered by www.kitware.com

ParaView development discussion is moving! Please visit 
https://discourse.paraview.org/ for future posts.

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Search the list archives at: http://markmail.org/search/?q=Paraview-developers

Follow this link to subscribe/unsubscribe:
https://public.kitware.com/mailman/listinfo/paraview-developers
___
Powered by www.kitware.com

ParaView discussion is moving! Please visit https://discourse.paraview.org/ for 
future posts.

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

Re: [Paraview] warp by vector field(displacement) with pressure field problem

[Moreland, Kenneth via ParaView]

Umit,

I took a look at your data, and the warp by vector filter is behaving exactly 
how I would expect.

The displacement field in your data is 0 everywhere except at the top edge. 
Thus, when you apply the warp the top edge will move and everything else will 
stay fixed where it is. In the middle of the top edge the warp vectors are 
pushing down into the rest of the mesh. This is causing some self-intersection 
with the mesh, which is why you are seeing the pressure profile sticking out 
the top.

You can use the Threshold filter to remove all cells that are not moving. 
However, that is only going to leave the top layer of cells, which I’m not sure 
what you want.

-Ken

From: ParaView [mailto:paraview-boun...@public.kitware.com] On Behalf Of Umut 
TABAK
Sent: Wednesday, June 13, 2018 11:38 PM
To: ParaView 
Subject: [EXTERNAL] [Paraview] warp by vector field(displacement) with pressure 
field problem

Dear all,

On a 2d model, I am trying to plot some mode shape for a coupled 
fluid-structure problem where a beam is coupled to a cavity domain. I have two 
variables for this problem: scalar pressure variable for the cavity and the 
displacement vectors for the beam at the top edges of the domain.

I am writing the attached vtk input file with my post processing codes.

And the result is somewhat as expected except that another pressure profile 
with an undeformed structural part is also visible with what I actually want to 
get. Please see the attached screenshot with the settings and visible fields.

Do I have a misunderstanding to accomplish this task. At the end, I would like 
to show the pressure profile which is a scalar field in the file along with the 
displacement field on the beam at the top edge.

I was thinking if there is deformed+undeformed visibility selection somehwere 
or not, maybe the issue is completely something else.

Any comments are appreciated.

BR,
Umut

P.S. For a decent view for the deformed vector field, you can set the scale 
factor of the warpbyvector field to 0.1 to 0.15

___
Powered by www.kitware.com

ParaView discussion is moving! Please visit https://discourse.paraview.org/ for 
future posts.

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

Re: [Paraview] x,y,z,time

[Moreland, Kenneth via ParaView]

A,

ParaView won’t recognize a column in a table as time. If you instead write out 
your data as a series of csv files, ParaView will treat each file as a time 
step and animate through them. So, if you re-write your data in multiple csv 
files, things get much easier in ParaView.

If you don’t want to re-write your data, I think you can still do what you want 
to do by using ParaView’s animation. The basic idea is to add a filter that 
pulls out the data for one time step, and then animate the parameter.

This will probably be clearer with an example. For this example, let’s assume 
your csv file has 20 time steps with time values 0, 1, 2,…, 19. The steps you 
could follow to animate time are as follows:


  1.  Load the csv file, convert them to table to points, and visualize them 
just as you are so far.
  2.  Add the threshold filter. Set the Minimum parameter to just under the 
first time value (say -0.5 in the example time steps I gave) and set the 
Maximum parameter to just over the first time value (say 0.5 in this example). 
This will extract the points of the first time step.
  3.  Open the Animation View (View -> Animation View).
  4.  Add a track for the Threshold filter for Threshold Range (0) (this 
corresponds to the Minimum value). Double click the track to open up the 
Animation Keyframes dialog. Set the first value to just under the first time 
value (-0.5 in this example), to just under the largest value (18.5 in this 
example).
  5.  Add a track for the Threshold filter for Threshold Range (1) (this 
corresponds to the Maximum value). Double click the track to open up the 
Animation Keyframes dialog. Set the first value to just under the first time 
value (0.5 in this example), to just under the largest value (19.5 in this 
example).
  6.  In the Animation View, set the number of frames (No. Frames) to the 
number of time steps in your data (20 in this example).
  7.  Hit play.

-Ken

From: ParaView  on behalf of A 

Date: Monday, May 7, 2018 at 8:40 PM
To: ParaView 
Subject: [EXTERNAL] [Paraview] x,y,z,time

This seems rudimentary, but I didnt find a relevant answer in my googles.

I csv file of x,y,z,time. I import them and then visualize them with 
tabletopoints. I just want to have an 3d scatter plot animation of them 
appearing.

how do I get paraview to recognize my "time" colunm as time in the table?

thanks!
a
___
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

Re: [Paraview] generic NetCDFreader support for CellData

[Moreland, Kenneth via ParaView]

The NetCDF Generic/CF reader has an option named “Output Type” that you can use 
to force the output to a vtkImageData. Then Cory’s original programmable filter 
will work, and you won’t have to deal with the coordinates.

-Ken

From: ParaView  on behalf of Cory Quammen 

Date: Wednesday, May 2, 2018 at 11:19 AM
To: "j.vanengel...@uu.nl" 
Cc: "paraview@public.kitware.com" 
Subject: [EXTERNAL] Re: [Paraview] generic NetCDFreader support for CellData

On Wed, May 2, 2018 at 11:37 AM Engelen, J. van (Joeri) 
> wrote:
Thank you Cory,

I sense this solution might work. However, the type of my data is “Rectilinear 
Grid”, not sure if that increases the complexity (dx = 1000, dy = 1000, dz = 1).

Ah. From the documentation: "To define a 
vtkRectilinearGrid,
 you must specify the dimensions of the data and provide three arrays of values 
specifying the coordinates along the x-y-z axes". It should be a similar 
script, but you'd have to add an additional X, Y, and Z position at the end of 
those coordinate arrays. So your script would look more like:

dims = inputs[0].GetDimensions()

ext = inputs[0].GetExtent()

output.SetDimensions(dims[0]+1, dims[1]+1, dims[2]+1)

output.SetExtent(ext[0], ext[1]+1, ext[2], ext[3]+1, ext[4], ext[5]+1)

inputPd = inputs[0].GetPointData()

outputCd = output.GetCellData()

outputCd.SetScalars(inputPd.GetScalars())



# Set coordinates

xCoords = inputs[0].GetXCoordinates().NewInstance()

xCoords.DeepCopy(inputs[0].GetXCoordinates())

xCoords.InsertNextValue(1.0) # Should a reasonable x value for your dataset 
larger than previous max x

output.SetXCoordinates(xCoords)



yCoords = inputs[0].GetYCoordinates().NewInstance()

yCoords.DeepCopy(inputs[0].GetYCoordinates())

yCoords.InsertNextValue(1.5) # Should be a reasonable like the x coordinate

output.SetYCoordinates(yCoords)



zCoords = inputs[0].GetZCoordinates().NewInstance()

zCoords.DeepCopy(inputs[0].GetZCoordinates())

zCoords.InsertNextValue(1.0) # Should be a reasonable like the z coordinate

output.SetZCoordinates(zCoords)


Some extent information also needs to be supplied, but I'm not entirely clear 
on how to do that correctly.


A better solution would be to add an option to the reader to take care of this 
for you if that is reasonable for your data.

Best,
Cory


Regardless, the script you provided does not give me output and it also does 
not throw an error, so I’m left a bit in the dark here.

From: Cory Quammen 
[mailto:cory.quam...@kitware.com]
Sent: Wednesday, May 02, 2018 5:01 PM
To: Engelen, J. van (Joeri)
Cc: paraview@public.kitware.com
Subject: Re: [Paraview] generic NetCDFreader support for CellData

Joeri,

You can try the "Point Data to Cell Data" filter, but that averages the point 
data and places the results in the cell data. You may wind up with some 
incorrect values using that approach.

Better would be to use a "Progammable Filter" to resize the cartesian grid read 
by the NetCDF reader and treat the point data as cell data. That is a bit 
involved, but not terrible. Set the Script to


dims = inputs[0].GetDimensions()

ext = inputs[0].GetExtent()

output.SetDimensions(dims[0]+1, dims[1]+1, dims[2]+1)

output.SetExtent(ext[0], ext[1]+1, ext[2], ext[3]+1, ext[4], ext[5]+1)

inputPd = inputs[0].GetPointData()

outputCd = output.GetCellData()

outputCd.SetScalars(inputPd.GetScalars())


This assumes that the NetCDF reader is producing a vtkImageData (you can check 
under the Information panel) which would report the Type: of the data set as 
"Image (Uniform Rectilinear Grid).

HTH,
Cory

On Wed, May 2, 2018 at 10:17 AM Engelen, J. van (Joeri) 
> wrote:
Hi,

I was wondering whether it is possible to read NetCDF files as CellData.

I have categorial data on Cartesian coordinates that I want to visualize. To be 
more specific, the data represents lithologies of the subsurface. So 1 = sand, 
3 = clay, 5 = rock.
The problem is that my file is automatically interpreted as PointData. To 
render this PointData, Paraview than consequently interpolates to achieve cell 
values.
With most variables that is not such a problem, but here it is. If I have sand 
(1) overlying rock (5), Paraview shows me clay (3).

Also, the information available to me was ambiguous,
-this source says it is probably not possible:
https://public.kitware.com/pipermail/paraview/2017-September/040984.html
-this states it is possible return Cell Data information:
https://www.paraview.org/ParaView/Doc/Nightly/www/py-doc/paraview.simple.NetCDFReader.html

I tried converting my data from NetCDF to a .vtr file in Python, but this 
resulted in more problems than solutions. (Just gave me very glitched 

Re: [Paraview] Normal vector to Exodus Sideset elements

[Moreland, Kenneth]

Dennis,

Yes, this cat has a zipper that allows it's skin to come right off. (Maybe I 
took that "skin the cat" metaphor a bit too far.)

Anyway, just run the Extract Surface filter. After that you can run the 
Generate Surface Normals filter.

-Ken

From: ParaView [mailto:paraview-boun...@public.kitware.com] On Behalf Of Dennis 
Conklin
Sent: Monday, April 9, 2018 1:22 PM
To: Paraview (parav...@paraview.org) 
Subject: [EXTERNAL] [Paraview] Normal vector to Exodus Sideset elements

All,

I'm sure I'm missing something here, so please knock me up alongside the head 
and point it out.

I load a sideset (surface) from an Exodus file.   I end up with a bunch of 
quads that Paraview conveniently created for me from the surface.   I need to 
calculate a Normal vector for each of these quads.  I have tried:

Generate surface normal - this is greyed out (wants poly data?)

Delaunay 3D on this to generate poly data, which I could then hopefully run 
Generate surface normal on, but it was still greyed

Tried Normal Glyphs which gave me point Normal vectors but gave me a bunch of 
Triangle and Polygon elements - I'm not sure what's going on here.
PointDataToCellData to get GlyphVector for each Cell
Python Calculator and ask for area(inputs[0]) the Polygon elements all get 
an Area of -1
Calculator to get Area_X as product of GlyphVector_X and Area
Threshold on Area_X > 0 -  don't want negative areas or some boundary cells 
with GlyphVector_X as negative
Python Calculator to get sum(Area_X)
This number appears to be about half of my expected number

These must be an easier way.  I'm not a vtk guru.   When I first load my 
surface, is there some vtk routine I could call within a Programmable filter to 
get the cell normal?

Alternatively, anyone got an easier way to skin this cat?

Thanks for any hints  you may have.

Dennis

___
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

[Paraview] CFP: SC18 Vis Showcase

[Moreland, Kenneth]

Sorry for the spam, but we have had many great entries using ParaView to create 
excellent HPC visualizations, and I though some ParaView users would like this 
opportunity to showcase your work. As an added bonus, this year accepted 
entries will get a publication in the Parallel Computing journal.

---

  Call for Papers

  SC18 Scientific Visualization & Data Analytics Showcase

 The International Conference for High Performance Computing, Networking,
  Storage, and Analysis, November 11-16, Dallas, TX, USA

   https://sc18.supercomputing.org/submit/showcases-submissions/

   *** Submissions due July 31, 2018 ***

   Contact: vis_showc...@info.supercomputing.org

---


SC18's Scientific Visualization & Data Analytics Showcase provides a forum
for the year's most instrumental movies in HPC.  Selected entries will be
displayed live in a museum/art gallery format. Six accepted finalists will
compete for the Best Visualization Award, and each finalist will present
his or her movie during a dedicated 15-minute presentation. Movies are
judged by overall quality, how they illuminate science, and for creative
innovations in the movie production process.

**Accepted submissions will be published in a virtual special issue of
Parallel Computing.** There will also be both a live display throughout the
conference so that attendees can experience and enjoy the latest in science
and engineering HPC results expressed through state-of-the-art
visualization technologies and a session at SC18 dedicated to the accepted
submissions.

REVIEW AND SELECTION PROCESS

Submissions need to include a movie (up to 1GB in size) and a short paper
(up to 6 pages including references). The short paper should describe the
scientific story conveyed by the movie, how the visualization helps
scientific discovery, and the "state-of-the-practice" information behind
making the movie.

Each submission will be peer reviewed by the Visualization & Data Analytics
Showcase Committee. Criteria for review include:

  * How effective is the visual communication of the data?
  * How relevant to the HPC community is the visualization?
  * What is the impact of the science story and how well is it told?
  * What visualization techniques were necessary to create the movie?

Finally, submissions should consider SC18's overall theme "HPC Inspires."

SUBMISSION INSTRUCTIONS

All submissions comprise a video and accompanying manuscript. The video
submitted should be in MP4, MOV, or AVI format and use a codec that will
play in typical desktop computer media players. There is no hard limit on
video length, but typical videos are less than 10 minutes long, and the
submitted file should be no more than 1GB.

The manuscript must be no more than 6 pages including references. The
manuscripts should be formatted according to the guidelines for the
Parallel Computing journal. Please use one of the following two templates
when formatting your manuscript:

  * LaTeX template (preferred):

https://sc18.supercomputing.org/app/uploads/2018/02/sc18-vis-showcase-latex-template.zip
  * Word template:

https://sc18.supercomputing.org/app/uploads/2018/02/sc18-vis-showcase-template.docx.zip

The manuscript should be uploaded as a pdf for the initial submission.
Authors of accepted submissions will have to send the original material
(.tex or .docx files and any source image files) to Parallel Computing's
publisher.

Visualization showcase entries are submitted through the SC Submission Site
(https://submissions.supercomputing.org/). Sign in, click the "Make a New
Submission" tab, and then select the "Scientific Visualization & Data
Analytics Showcase" link.

PROCEEDINGS

The manuscripts of the six accepted submissions will be published in a
special virtual issue of Parallel Computing. Because the Parallel Computing
journal does not have the same publisher as the proceedings of the rest of
the SC conference, authors will have to resubmit their content. (There will
not, however, be a separate review process.)

BEST VISUALIZATION

After accepted authors present their work, a select committee will choose
among them for the Best Visualization Award. The winner of the Best
Visualization Award will be presented during the SC Awards Session.

IMPORTANT DATES

Video/Manuscript Submission:July 31, 2018
Notifications Sent: September 15, 2018
Final Material Deadline:October 31
___
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 

Re: [Paraview] CSV Files and Volumetric Data

[Moreland, Kenneth]

In order for this to work, your CSV file has to contain columns for both the 
spatial position AND the vectors you want to visualize. It looks like the csv 
file you are creating is actually only writing out the vector field. To get 
this to work, your gen_sally function should fill a second array that simply 
contains the x, y, and z variables on the inner loop. Then your create_csv 
function should write a csv file with 6 columns. Three of them are for the 
position coordinate I just described, and the other 3 are the vector that you 
are currently writing.

-Ken

From: ParaView [mailto:paraview-boun...@public.kitware.com] On Behalf Of 
GILLILAND G. (929919)
Sent: Wednesday, April 4, 2018 3:16 PM
To: paraview@public.kitware.com
Subject: [EXTERNAL] [Paraview] CSV Files and Volumetric Data

Hi there,

I have assignment where data is produced in C, and i have to export it to a 
relevant file format that works with a visualisation software of my choice.

I've chosen Paraview, because i know that CSV files work with this software.

Rather than table to points, i want to be able to convert the points to a 
volumetric dataset... So i can have greater control over colour mapping. But i 
was told that the code that generates the dataset is missing some data, or i 
may have missed something in the way i export the file.

I've included the code below; i would love some help with this. I'm only decent 
with coding in Java and my lecturer has really given the class a lot of support.

Cheers

#include
#include
#include
#include

void gen_sally( int xs, int ys, int zs, int time, float *sally )
/*
 *  Gen_Sally creates a vector field of dimension [xs,ys,zs,3] from
 *  a proceedural function. By passing in different time arguements,
 *  a slightly different and rotating field is created.
 *
 *  The magnitude of the vector field is highest at some funnel shape
 *  and values range from 0.0 to around 0.4 (I think).
 *
 *  I just wrote these comments, 8 years after I wrote the function.
 *
 *  Developed by Sally of Sally University
 *
 */
{
  float x, y, z;
  int ix, iy, iz;
  float r, xc, yc, scale, temp, z0;
  float r2 = 8;
  float SMALL = 0.001;
  float xdelta = 1.0 / (xs-1.0);
  float ydelta = 1.0 / (ys-1.0);
  float zdelta = 1.0 / (zs-1.0);

  for( iz = 0; iz < zs; iz++ )
  {
z = iz * zdelta;// map z to 0->1
xc = 0.5 + 0.1*sin(0.04*time+10.0*z);   // For each z-slice, determine the 
spiral circle.
yc = 0.5 + 0.1*cos(0.03*time+3.0*z);//(xc,yc) determine the center of 
the circle.
r = 0.1 + 0.4 * z*z + 0.1 * z * sin(8.0*z); //  The radius also changes at each 
z-slice.
r2 = 0.2 + 0.1*z;   //r is the center radius, r2 is 
for damping
for( iy = 0; iy < ys; iy++ )
{
y = iy * ydelta;
for( ix = 0; ix < xs; ix++ )
{
x = ix * xdelta;
temp = sqrt( (y-yc)*(y-yc) + (x-xc)*(x-xc) );
scale = fabs( r - temp );
/*
 *  I do not like this next line. It produces a discontinuity
 *  in the magnitude. Fix it later.
 *
 */
   if ( scale > r2 )
  scale = 0.8 - scale;
   else
  scale = 1.0;
z0 = 0.1 * (0.1 - temp*z );
   if ( z0 < 0.0 )  z0 = 0.0;
   temp = sqrt( temp*temp + z0*z0 );
scale = (r + r2 - temp) * scale / (temp + SMALL);
scale = scale / (1+z);
   *sally++ = scale * (y-yc) + 0.1*(x-xc);
   *sally++ = scale * -(x-xc) + 0.1*(y-yc);
   *sally++ = scale * z0;
}
}
  }
}


void create_csv(char* filename,float *sally, int size){
printf("1");
printf("\n Creating %s.csv file",filename);
FILE *fp;
fp=fopen(filename,"w");
fprintf(fp,"X,Y,Z\n");
int i;
int counter = 0;
for(i = 0; i< size; i++){
if(sally[i] == 0){
fprintf(fp,"0");
}
else{
fprintf(fp,"%f",sally[i]);
}
counter++;
if(counter == 3){
fprintf(fp, "\n");
counter = 0;
}
else{
fprintf(fp,",");
}
}
fclose(fp);
printf("\n %sfile created",filename);
}
int main(int argc, char *argv[]){
printf("1\n");
//read from args
int xs;
int ys;
int zs;
int time;
sscanf(argv[1],"%d",);
sscanf(argv[2],"%d",);
sscanf(argv[3],"%d",);
sscanf(argv[4],"%d",);


int arraySize = xs*ys*zs*3;
//allocate memeory for array. This is done so that stack memory doesn't run 
out.'
float* sally;
sally = (float*)malloc((arraySize) * sizeof(float));

//runs the code. One of the args is a pointer so no return type is needed.
gen_sally(xs,ys,zs,time,sally);
//create varibles for file generation
char filename[20] = "results.csv";
create_csv(filename, sally, arraySize);

free(sally);
return 0;
}

___
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

Re: [Paraview] Hello

[Moreland, Kenneth]

Hassan,

If you are interested in learning how to use ParaView for post processing, I 
recommend starting with the ParaView Tutorial 
(https://www.paraview.org/Wiki/The_ParaView_Tutorial). The tutorial does not 
talk about SU2 specifically, but it teaches basic techniques for visualization, 
including from CFD simulations.

SU2 can write out .vtk files, which ParaView can read. After learning the 
basics and you have more specific questions, we can help you with that.

-Ken

From: ParaView  on behalf of Hassan 
Elsahely 
Date: Thursday, February 22, 2018 at 8:50 AM
To: "paraview@public.kitware.com" 
Subject: [EXTERNAL] [Paraview] Hello


Hello,



I am using paraview as a post processing in windows not for a long time and I 
am working on laminar flat plate tutorial based on SU2. I would like to plot 
Blasius solution, similarity variable in function of velocity u/Ue. Any help.



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

Re: [Paraview] How to get filled slices of tubes and disks around streamlines

[Moreland, Kenneth]

Thomas,

I think a better approach to what you are doing is to use the Glyph filter on 
your streamlines instead of the Tubes filter. After you create your 
streamlines, perform the following steps.


  1.  With the stream tracer selected in the pipeline browser, add the glyph 
filter (it’s icon looks like a Christmas ornament).
  2.  Turn on the advanced properties in the properties panel.
  3.  In the Glyph Type combo box, select “2D Glyph”.
  4.  Assuming you have the advanced properties on, you should see a second 
Glyph Type combo box appear under the first one. Select “Circle” in this second 
box.
  5.  Click on Filled checkbox.
  6.  Under Active Attributes, make sure the Vectors property is set to the 
vector field you used to create the streamlines.
  7.  Scroll down to where it says Glyph Transform. Set the Rotate property to 
0, 90, 0. (That is, change the middle number from 0 to 90.)
  8.  Click Apply.

You should now see filled circles perpendicular to the streamlines. You might 
need to change the Glyph filter’s Masking mode to either “All Points” or “Every 
Nth Point” to get the samples of circles how you want them.

-Ken

From: ParaView [mailto:paraview-boun...@public.kitware.com] On Behalf Of Thomas 
Oliveira
Sent: Wednesday, February 21, 2018 6:56 AM
To: ParaView 
Subject: [EXTERNAL] [Paraview] How to get filled slices of tubes and disks 
around streamlines

Hi,
I am tracing streamlines using "Stream Tracer with Custom Source" and creating 
tubes along them using "Tubes". Since "Tubes" creates surfaces, if I use 
"Slice" on the tubes I don't get elliptical disks but ellipses.
I would like to:
1) Visualize elliptical disks when using "Slice" on "Tubes".

I would be even better if, instead of getting the elliptical disks, which are 
bounded by the cross section of the tubes, I could get circular disks in the 
same position, i.e., if I could
2) represent circular disks on a plane that crosses the streamlines centred on 
the intersections between the plane and the streamlines

Would you see a way of doing that or getting the same visual result?
Thank you,
Thomas Oliveira
___
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

Re: [Paraview] Help

[Moreland, Kenneth]

Mohammed,

The answer to your question depends on a lot of things.

Depending on your simulation, it might be easiest if the simulation itself 
detected when a particle exits the pipe. It is, after all, the thing tracking 
the particles to begin with. From there it could output its count in a csv 
file. You could load that up and plot in in ParaView, although a typical 
spreadsheet program can do that as well.

Assuming it is not feasible to have your simulation do the count, getting 
ParaView to do it depends on a lot of things. First, it depends on whether 
particles “die” in your simulation. It is pretty common in simulation code to 
have particles leave the defined domain (or otherwise become invalid) and then 
get removed from the list of particles that get written out. If particles never 
die, then your job in ParaView is easier. You can just count how many particles 
are past the out end of the pipe and count that.

If particles are born or die in your simulation, then the problem becomes much 
more difficult.

-Ken

From: ParaView  on behalf of Mohammed 
Babiker 
Date: Wednesday, February 14, 2018 at 11:31 AM
To: "paraview@public.kitware.com" 
Subject: [EXTERNAL] [Paraview] Help

Good evening

I am a new user of paraview . I am running a simulation of particles passing 
through a horizontal pipe exiting from other side I would like to know howm 
many numbers of particles exit and plot that number with respect to time . any 
help will be appreciated.

Regards.

___
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