Re: [Paraview] Animation. Second view of a field variable. How to format the file.
Hi Dave, Thank you for your second response on the 12th of September in which you pointed out that it is very easy to get the file formats wrong. You also sent a full vts example which I copied to a vts file and let ParaView display in a 3D window. OK. Therefore I followed your example and wrote my vts as follows: the fileddata (which represents an integrated variable of the system) is now put before the first Piece Extent in every one of the vts files. This takes me a step further, but there is still something that I have not got correctly because this is what happens in ParaView: File/Open/pvd/Apply --- OK. I see my geometry in the 3D view. Split Vertical --- Empty window opens SpreadSheet View --- OK Click eye --- Some table opens in second window Display FieldData --- OK I see my FieldData in second window animate --- OK FieldData value changes But it is an animated xyPlot that I wish, so I try the xyplot now in ParaView: File/Open/pvd/Apply --- 3D view opens correctly Split Vertical --- Empty window opens xyPlot View --- NO EYE From here I cannot proceed and I conclude I still do something wrong. Now, you also wrote ... make sure you copy the point and cell centred data across too. Is this where I fail? Christer My vts file(s) as they look now: ?xml version=1.0? VTKFile type=StructuredGrid version=0.1 byte_order=LittleEndian compressor=vtkZLibDataCompressor StructuredGrid WholeExtent=0 1 0 40 0 91 FieldData DataArray type=Float32 Name=SF1 NumberOfTuples=1 format=ascii RangeMin=5.204563E+04 RangeMax=5.204563E+04 5.204563E+04 /DataArray /FieldData Piece Extent=0 1 0 40 0 45 PointData Scalars=Temperature DataArray type=Float32 Name=Temperature format=ascii 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [snip] David E DeMarle wrote: You could save K(t) as field data. Then you can extract the zero'th entry in the K array (which has just one element) and plot that over time. A vts file with field data in it looks something like this... VTKFile type=StructuredGrid [snip] StructuredGrid WholeExtent=[snip] FieldData DataArray type=Float64 Name=MyField NumberOfTuples=1 format=ascii RangeMin=26.1 RangeMax=26.1 26.1 /DataArray /FieldData [snip] On Mon, Sep 8, 2008 at 12:30 AM, Christer Gustafsson [EMAIL PROTECTED] wrote: Hi All, The data (x_i(t),y_i(t),z_i(t)) is formatted as one pvd and many vts file. The animation runs fine in one 3D view. In addition there is an associated time-varying function (t,K(t)) (e.g. the kinetic energy of the system). I would like to create a second view and herein animate K(t) along with the 3D view. In what file format (example?) would the (t,K(t)) be cast in order to make the animations in two views be possible? (The ParaView Guide talks about `unaligned data items' in connection with `Time' in Chapter 10, section 10.3, but from that point on I find no more path). Thanks, ChG ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview
[Paraview] Excluding properties from state files
Hi, Sorry to post this again - I still haven't found a solution. Is there a way to stop the state file (pvsm) from storing the property? I.e. is there an attribute I can add in the XML file? I would like ParaView to call the function to get the property instead. Thanks in advance, Paul ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] Excluding properties from state files
Set 'is_internal=1' and that property won't be saved in the state file. It also means that it's default value won't be pushed on to the underlying VTK object when the proxy is created. If you change the value, however, then it will be pushed as usual. Utkarsh Paul Edwards wrote: Hi, Sorry to post this again - I still haven't found a solution. Is there a way to stop the state file (pvsm) from storing the property? I.e. is there an attribute I can add in the XML file? I would like ParaView to call the function to get the property instead. Thanks in advance, Paul ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] Excluding properties from state files
Thanks Utkarsh - that's the property I needed. 2008/9/16 Utkarsh Ayachit [EMAIL PROTECTED]: Set 'is_internal=1' and that property won't be saved in the state file. It also means that it's default value won't be pushed on to the underlying VTK object when the proxy is created. If you change the value, however, then it will be pushed as usual. Utkarsh Paul Edwards wrote: Hi, Sorry to post this again - I still haven't found a solution. Is there a way to stop the state file (pvsm) from storing the property? I.e. is there an attribute I can add in the XML file? I would like ParaView to call the function to get the property instead. Thanks in advance, Paul ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] Getting values from a saved state with a script
Walter Landry [EMAIL PROTECTED] wrote: Hello, I have set up paraview so that it does some calculations and displays the result in a spreadsheet. I can save the state and then load it in a script with LoadState(). I know how to save images of the 3D views with GetRenderView().WriteImage. But I have not figured out how to get the data from the spreadsheet. Ideally, I would be like to have the script print out a number, although it would be fine to just save an image of the spreadsheet. I tried getting the list of proxies with paraview.servermanager.pyProxyManager().GetProxiesOnOcnnection(paraview.servermanager.ActiveConnection) I can get a SpreadSheetView object and SpreadSheetRepresentation objects, but I do not know what to do with them. I just can't manage to get numbers out. I am sure there is a way to do it, but I am a bit stumped right now. I figured it out. I can just connect them to a DataSetWriter using a script like pm=ProxyManager() energy=pm.GetProxiesOnConnection(ActiveConnection)[sources][IntegrateAttributes3] work=pm.GetProxiesOnConnection(ActiveConnection)[sources][IntegrateAttributes2] energy_writer=writers.DataSetWriter(Input=energy,FileType=1,FileName=energy.vtk) work_writer=writers.DataSetWriter(Input=work,FileType=1,FileName=work.vtk) energy_writer.UpdatePipeline() work_writer.UpdatePipeline() That writes the data to a file. Cheers, Walter Landry [EMAIL PROTECTED] ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview
[Paraview] filter
Hi! I´m still encountering problems with my filter, it does compile and link neatly with paraview, but at runtime an error occurs which causes paraview to shut down. I´m pretty sure that the problem is how I define the output, that this is not done correctly. Lastly, I commented just the algorithm out to test the in- and output, and it looks like this: void vtkHorizontalAverage::Execute() {vtkDataSet *input = this - GetInput(); vtkStructuredGrid *output = this - GetOutput(); int DataSetType; DataSetType = input - GetDataObjectType(); if(DataSetType != VTK_RECTILINEAR_GRID) { vtkErrorMacro(unsupported DataSet format);} int dims[3],i,j,k; // pointer for casting vtkSmartPointervtkRectilinearGrid rect_input = vtkSmartPointervtkRectilinearGrid::New(); //vtkStructuredGrid *structgrid_input = vtkStructuredGrid::New(); vtkPointData *scaldat = input - GetPointData(); vtkDataArray *new_scaldat; vtkSmartPointervtkFloatArray scalars = vtkSmartPointervtkFloatArray::New(); vtkSmartPointervtkFloatArray new_scalars = vtkSmartPointervtkFloatArray::New(); scalars = vtkFloatArray::SafeDownCast(scaldat); //this is now the input as vtkFloatArray scaldat - Delete(); int numbOfComp = scalars - GetNumberOfComponents(); if(numbOfComp ==1) vtkErrorMacro(Number of Components of Input Array bigger than 1!); // int numOfTuples = scalars - GetNumberOfTuples(); new_scalars - SetNumberOfComponents(1); new_scalars - SetNumberOfTuples(numOfTuples); rect_input = vtkRectilinearGrid::SafeDownCast(input); if(!rect_input) {vtkErrorMacro(SafeDownCast to Rectilinear Grid failed);} else { rect_input - GetDimensions(dims); //get the points: /*vtkDataArray *XCoord = vtkDataArray::New(); vtkDataArray *YCoord = vtkDataArray::New(); vtkDataArray *ZCoord = vtkDataArray::New();*/ vtkDataArray *XCoord = rect_input - GetXCoordinates(); vtkDataArray *YCoord = rect_input - GetYCoordinates(); vtkDataArray *ZCoord = rect_input - GetZCoordinates(); output - SetDimensions(dims); vtkSmartPointervtkPoints pts = vtkSmartPointervtkPoints::New(); pts - SetNumberOfPoints(numOfTuples); for(int i=0;i=dims[2];i++) { for(int j=0;j=dims[1];j++) {for(int k=0;k=dims[0];k++) {pts - InsertNextPoint(XCoord - GetTuple1(k),YCoord - GetTuple1(j), ZCoord - GetTuple1(i)); }//for (int k }//for(int j }//for (int i XCoord - Delete(); YCoord - Delete(); ZCoord - Delete(); output - SetPoints(pts); //Now the points are set //output - GetPointData()- CopyAllocate(new_scalars,0,new_scalars - GetNumberOfTuples()); output - GetPointData() - PassData(input - GetPointData()); new_scalars - Delete(); Here I get as an input a data set, cast it down to rectilinear grid and compute the points for the structured grid. As I said, Paraview shuts down (Access violation my debugger says) when I execute it. Any help would be appreciated! thx, NH _ Express yourself instantly with MSN Messenger! Download today it's FREE! http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] filter
what line of code is the access violation occurring at? _ | Mike Jackson - Principal Software Engineer | | BlueQuartz Software | | [EMAIL PROTECTED] | | www.bluequartz.net | --- On Sep 16, 2008, at 1:08 PM, Natalie Happenhofer wrote: Hi! I´m still encountering problems with my filter, it does compile and link neatly with paraview, but at runtime an error occurs which causes paraview to shut down. I´m pretty sure that the problem is how I define the output, that this is not done correctly. Lastly, I commented just the algorithm out to test the in- and output, and it looks like this: void vtkHorizontalAverage::Execute() {vtkDataSet *input = this - GetInput(); vtkStructuredGrid *output = this - GetOutput(); int DataSetType; DataSetType = input - GetDataObjectType(); if(DataSetType != VTK_RECTILINEAR_GRID) { vtkErrorMacro(unsupported DataSet format);} int dims[3],i,j,k; // pointer for casting vtkSmartPointervtkRectilinearGrid rect_input = vtkSmartPointervtkRectilinearGrid::New(); //vtkStructuredGrid *structgrid_input = vtkStructuredGrid::New(); vtkPointData *scaldat = input - GetPointData(); vtkDataArray *new_scaldat; vtkSmartPointervtkFloatArray scalars = vtkSmartPointervtkFloatArray::New(); vtkSmartPointervtkFloatArray new_scalars = vtkSmartPointervtkFloatArray::New(); scalars = vtkFloatArray::SafeDownCast(scaldat); //this is now the input as vtkFloatArray scaldat - Delete(); int numbOfComp = scalars - GetNumberOfComponents(); if(numbOfComp ==1) vtkErrorMacro(Number of Components of Input Array bigger than 1!); // int numOfTuples = scalars - GetNumberOfTuples(); new_scalars - SetNumberOfComponents(1); new_scalars - SetNumberOfTuples(numOfTuples); rect_input = vtkRectilinearGrid::SafeDownCast(input); if(!rect_input) {vtkErrorMacro(SafeDownCast to Rectilinear Grid failed);} else { rect_input - GetDimensions(dims); //get the points: /*vtkDataArray *XCoord = vtkDataArray::New(); vtkDataArray *YCoord = vtkDataArray::New(); vtkDataArray *ZCoord = vtkDataArray::New();*/ vtkDataArray *XCoord = rect_input - GetXCoordinates(); vtkDataArray *YCoord = rect_input - GetYCoordinates(); vtkDataArray *ZCoord = rect_input - GetZCoordinates(); output - SetDimensions(dims); vtkSmartPointervtkPoints pts = vtkSmartPointervtkPoints::New(); pts - SetNumberOfPoints(numOfTuples); for(int i=0;i=dims[2];i++) { for(int j=0;j=dims[1];j++) {for(int k=0;k=dims[0];k++) {pts - InsertNextPoint(XCoord - GetTuple1(k),YCoord - GetTuple1(j), ZCoord - GetTuple1(i)); }//for (int k }//for(int j }//for (int i XCoord - Delete(); YCoord - Delete(); ZCoord - Delete(); output - SetPoints(pts); //Now the points are set //output - GetPointData()- CopyAllocate(new_scalars, 0,new_scalars - GetNumberOfTuples()); output - GetPointData() - PassData(input - GetPointData()); new_scalars - Delete(); Here I get as an input a data set, cast it down to rectilinear grid and compute the points for the structured grid. As I said, Paraview shuts down (Access violation my debugger says) when I execute it. Any help would be appreciated! thx, NH Express yourself instantly with MSN Messenger! MSN Messenger ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] filter
It might help if you ran your code in a debugger and found out exactly where the error occurred. I think the for loops are wrong. The conditions should be , not =. You might be mixing up dimension lengths and extent values. -Ken On 9/16/08 11:08 AM, Natalie Happenhofer [EMAIL PROTECTED] wrote: Hi! I´m still encountering problems with my filter, it does compile and link neatly with paraview, but at runtime an error occurs which causes paraview to shut down. I´m pretty sure that the problem is how I define the output, that this is not done correctly. Lastly, I commented just the algorithm out to test the in- and output, and it looks like this: void vtkHorizontalAverage::Execute() {vtkDataSet *input = this - GetInput(); vtkStructuredGrid *output = this - GetOutput(); int DataSetType; DataSetType = input - GetDataObjectType(); if(DataSetType != VTK_RECTILINEAR_GRID) { vtkErrorMacro(unsupported DataSet format);} int dims[3],i,j,k; // pointer for casting vtkSmartPointervtkRectilinearGrid rect_input = vtkSmartPointervtkRectilinearGrid::New(); //vtkStructuredGrid *structgrid_input = vtkStructuredGrid::New(); vtkPointData *scaldat = input - GetPointData(); vtkDataArray *new_scaldat; vtkSmartPointervtkFloatArray scalars = vtkSmartPointervtkFloatArray::New(); vtkSmartPointervtkFloatArray new_scalars = vtkSmartPointervtkFloatArray::New(); scalars = vtkFloatArray::SafeDownCast(scaldat); //this is now the input as vtkFloatArray scaldat - Delete(); int numbOfComp = scalars - GetNumberOfComponents(); if(numbOfComp ==1) vtkErrorMacro(Number of Components of Input Array bigger than 1!); // int numOfTuples = scalars - GetNumberOfTuples(); new_scalars - SetNumberOfComponents(1); new_scalars - SetNumberOfTuples(numOfTuples); rect_input = vtkRectilinearGrid::SafeDownCast(input); if(!rect_input) {vtkErrorMacro(SafeDownCast to Rectilinear Grid failed);} else { rect_input - GetDimensions(dims); //get the points: /*vtkDataArray *XCoord = vtkDataArray::New(); vtkDataArray *YCoord = vtkDataArray::New(); vtkDataArray *ZCoord = vtkDataArray::New();*/ vtkDataArray *XCoord = rect_input - GetXCoordinates(); vtkDataArray *YCoord = rect_input - GetYCoordinates(); vtkDataArray *ZCoord = rect_input - GetZCoordinates(); output - SetDimensions(dims); vtkSmartPointervtkPoints pts = vtkSmartPointervtkPoints::New(); pts - SetNumberOfPoints(numOfTuples); for(int i=0;i=dims[2];i++) { for(int j=0;j=dims[1];j++) {for(int k=0;k=dims[0];k++) {pts - InsertNextPoint(XCoord - GetTuple1(k),YCoord - GetTuple1(j), ZCoord - GetTuple1(i)); }//for (int k }//for(int j }//for (int i XCoord - Delete(); YCoord - Delete(); ZCoord - Delete(); output - SetPoints(pts); //Now the points are set //output - GetPointData()- CopyAllocate(new_scalars,0,new_scalars - GetNumberOfTuples()); output - GetPointData() - PassData(input - GetPointData()); new_scalars - Delete(); Here I get as an input a data set, cast it down to rectilinear grid and compute the points for the structured grid. As I said, Paraview shuts down (Access violation my debugger says) when I execute it. Any help would be appreciated! thx, NH Express yourself instantly with MSN Messenger! MSN Messenger http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/ Kenneth Moreland *** Sandia National Laboratories *** *** *** *** email: [EMAIL PROTECTED] ** *** ** phone: (505) 844-8919 *** fax: (505) 845-0833 ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview