Re: [Paraview] RungeKutta4 implementation issue
I suggest you attach a debugger to the code and step through. Then you'll be able to answer your question. Since the velocity field of your data may be high, it is entirely possible that the computed next position is outside of the bounds. Perhaps reducing the step size will help. Try even Rk4.5 with adaptive steps, this may help if the velocity gradient at the edge of the domain is very large.. JB How is it possible??? double *bds = reader-GetOutput()-GetBounds(); printf("%lf %lf %lf %lf\n", bds[0], bds[1], bds[2], bds[3]); - 0.00 69.00 0.00 69.00 printf("point1[0]=%f point1[1]=%f point1[2]=%f\n", point1[0], point1[1], point1[2]); - point1[0]=58.196274 point1[1]=68.862898 point1[2]=0.00 if (integrator-ComputeNextStep(point1, point2, 0, delT.Interval, stepTaken, minStep, maxStep, this-MaximumError, error) != 0) break; printf("point2[0]=%f point2[1]=%f point2[2]=%f\n", point2[0], point2[1], point2[2]); - point2[0]=59.002164 point2[1]=69.454414 point2[2]=0.00 PS: This message should probably be sent to the VTK mailing list but it seems that it is no longer so active... Dcouvrez toutes les possibilits de communication avec vos proches ___ 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://www.paraview.org/mailman/listinfo/paraview -- John Biddiscombe,email:biddisco @ cscs.ch http://www.cscs.ch/ CSCS, Swiss National Supercomputing Centre | Tel: +41 (91) 610.82.07 Via Cantonale, 6928 Manno, Switzerland | Fax: +41 (91) 610.82.82 ___ 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://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] RungeKutta4 implementation issue
Hi Fred: As JB suggested, please try using RK45, which employs an adaptive step size that is insensitive to the velocity magnitude and hence is capable of capturing accurate shape of the curve (flow line). Using RK4 (still adopting a fixed step size) does not guarantee a sufficient / desired acurracy of numerical integration (though certainly it is better than RK2) since it is stil sensitive to the velocity magnitude. Recently we improved the flow line integration accuracy issue. For any problems, please feel free to let us know. Thanks. -Zhanping -- Zhanping Liu, PhD Kitware, Inc. 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x 138 http://www.zhanpingliu.org On Mon, Aug 17, 2009 at 5:06 AM, John Biddiscombe biddi...@cscs.ch wrote: I suggest you attach a debugger to the code and step through. Then you'll be able to answer your question. Since the velocity field of your data may be high, it is entirely possible that the computed next position is outside of the bounds. Perhaps reducing the step size will help. Try even Rk4.5 with adaptive steps, this may help if the velocity gradient at the edge of the domain is very large.. JB How is it possible??? double *bds = reader-GetOutput()-GetBounds(); printf(%lf %lf %lf %lf\n, bds[0], bds[1], bds[2], bds[3]); - 0.00 69.00 0.00 69.00 printf(point1[0]=%f point1[1]=%f point1[2]=%f\n, point1[0], point1[1], point1[2]); - point1[0]=58.196274 point1[1]=68.862898 point1[2]=0.00 if (integrator-ComputeNextStep(point1, point2, 0, delT.Interval, stepTaken, minStep, maxStep, this-MaximumError, error) != 0) break; printf(point2[0]=%f point2[1]=%f point2[2]=%f\n, point2[0], point2[1], point2[2]); - point2[0]=59.002164 point2[1]=69.454414 point2[2]=0.00 PS: This message should probably be sent to the VTK mailing list but it seems that it is no longer so active... -- Découvrez toutes les possibilités de communication avec vos procheshttp://www.microsoft.com/windows/windowslive/default.aspx -- ___ 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://www.paraview.org/mailman/listinfo/paraview -- John Biddiscombe,email:biddisco @ cscs.chhttp://www.cscs.ch/ CSCS, Swiss National Supercomputing Centre | Tel: +41 (91) 610.82.07 Via Cantonale, 6928 Manno, Switzerland | Fax: +41 (91) 610.82.82 ___ 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://www.paraview.org/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://www.paraview.org/mailman/listinfo/paraview
[Paraview] pvclient
I'm new to paraview. I installed paraview-2.4.4_6 on FreeBSD 8.0-beta2 SMP ia64. I launch the server as: % pvserver --use-offscreen-rendering Listen on port: 1 Waiting for client... then I launch the client and get this error: TZAV pvclient Connect to localhost:1 ParaView error: InitializeTcl failed Tk_Init error: Can't find a usable tk.tcl in the following directories: /usr/local/lib/paraview-2.4/tk8.4 /usr/local/lib/paraview-2.4/tk8.4/tk.tcl: no event type or button # or keysym no event type or button # or keysym while executing bind Listbox MouseWheel { %W yview scroll [expr {- (%D / 120) * 4}] units } (file /usr/local/lib/paraview-2.4/tk8.4/listbox.tcl line 182) invoked from within source /usr/local/lib/paraview-2.4/tk8.4/listbox.tcl (in namespace eval :: script line 1) invoked from within namespace eval :: [list source [file join $::tk_library $file.tcl]] (procedure SourceLibFile line 2) invoked from within SourceLibFile listbox (in namespace eval ::tk script line 4) invoked from within namespace eval ::tk { SourceLibFile button SourceLibFile entry SourceLibFile listbox SourceLibFile menu SourceLibFile panedwindow SourceLibFile ... invoked from within if {$::tk_library ne } { if {[string equal $tcl_platform(platform) macintosh]} { proc ::tk::SourceLibFile {file} { if {[catch { namesp... (file /usr/local/lib/paraview-2.4/tk8.4/tk.tcl line 397) invoked from within source /usr/local/lib/paraview-2.4/tk8.4/tk.tcl (uplevel body line 1) invoked from within uplevel #0 [list source $file] This probably means that tk wasn't installed properly. Exit Client % I've tcl-8.4.19_3,1 and tk-8.4.19_2,2 installed. Any advice? many thanks -- Anton Shterenlikht Room 2.6, Queen's Building Mech Eng Dept Bristol University University Walk, Bristol BS8 1TR, UK Tel: +44 (0)117 928 8233 Fax: +44 (0)117 929 4423 ___ 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://www.paraview.org/mailman/listinfo/paraview
[Paraview] make error when building ParaView 3.6
Hi, I'm trying to build ParaView 3.6 on an Ubuntu machine (w/ Qt 4.5) but I'm hitting an error when I run 'make'. I'm doing the 'ccmake path_to_source' and keeping the default configuration except setting number of processors to 8 (I have a Core i7) and turning on the PARAVIEW_ENABLE_PYTHON and PARAVIEW_USE_MPI. I try to turn off BUILD_SHARED_LIBS, but it keeps switching itself back on (why??). Anyway, it configures fine with just a complaint about Qt PHONON library not found, which I understand is just a warning, not a show-stopper (right?). But when I run 'make', it gets to about 4% and then dies as follows: Scanning dependencies of target MapReduceMPI [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/mapreduce.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/keyvalue.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/keymultivalue.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/irregular.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/hash.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/memory.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/error.cpp.o Linking CXX shared library ../../../../bin/libMapReduceMPI.so /usr/bin/ld: /usr/local/lib/libmpich.a(allreduce.o): relocation R_X86_64_32S against `MPIR_ThreadSingle' can not be used when making a shared object; recompile with -fPIC /usr/local/lib/libmpich.a: could not read symbols: Bad value collect2: ld returned 1 exit status make[2]: *** [bin/libMapReduceMPI.so.pv3.7] Error 1 make[1]: *** [VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/all] Error 2 make: *** [all] Error 2 So, anyone got any ideas? Thanks in advance! ~~Karl ___ 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://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] make error when building ParaView 3.6
I had a similar problem. To get around it I had to build shared mpi libraries (i.e. libmpich.so). Any simpler suggestions anyone? On Mon, Aug 17, 2009 at 10:47 AM, Karl Battamskarlbatt...@gmail.com wrote: Hi, I'm trying to build ParaView 3.6 on an Ubuntu machine (w/ Qt 4.5) but I'm hitting an error when I run 'make'. I'm doing the 'ccmake path_to_source' and keeping the default configuration except setting number of processors to 8 (I have a Core i7) and turning on the PARAVIEW_ENABLE_PYTHON and PARAVIEW_USE_MPI. I try to turn off BUILD_SHARED_LIBS, but it keeps switching itself back on (why??). Anyway, it configures fine with just a complaint about Qt PHONON library not found, which I understand is just a warning, not a show-stopper (right?). But when I run 'make', it gets to about 4% and then dies as follows: Scanning dependencies of target MapReduceMPI [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/mapreduce.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/keyvalue.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/keymultivalue.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/irregular.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/hash.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/memory.cpp.o [ 4%] Building CXX object VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/error.cpp.o Linking CXX shared library ../../../../bin/libMapReduceMPI.so /usr/bin/ld: /usr/local/lib/libmpich.a(allreduce.o): relocation R_X86_64_32S against `MPIR_ThreadSingle' can not be used when making a shared object; recompile with -fPIC /usr/local/lib/libmpich.a: could not read symbols: Bad value collect2: ld returned 1 exit status make[2]: *** [bin/libMapReduceMPI.so.pv3.7] Error 1 make[1]: *** [VTK/Utilities/mrmpi/src/CMakeFiles/MapReduceMPI.dir/all] Error 2 make: *** [all] Error 2 So, anyone got any ideas? Thanks in advance! ~~Karl ___ 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://www.paraview.org/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://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] How to set up X connection
When you say do have a login screen in each of the machine, what do you mean? Is it a graphical one or a black window with a login prompt at the top? If you have a graphical login screen, that means you are running X (through xdm or some other display manager). If not, you need to enable X on those machines. This may require installing it - or not if it is already there - and turning it on. How to do that depends on the Linux distribution. Once that's done, you should create a user and make Linux login to that user automatically at boot time. You will also have to make sure that the display is accessible. This may require using something like the xhost command. I recommend searching on the web for more information on setting up X. -berk On Wed, Aug 12, 2009 at 4:06 AM, chew pinglcp8...@msn.com wrote: Dear Berk well, basically we don't have a very clear idea of how paraview deals with parallel rendering. we could ssh -X to any of these machines and run xclock and the X windows does appear at the machine we are from. on the other hand, while trying to start a parallel pvserver, everything works fine until when a client (paraview) connects to the parallel pvserver, and we get a window alerting us about parallel rendering disabled.. by following the thread in paraview mailing list, we got to know that to enable parallel rendering, each of the machine running the parallel pvserver should be able to generate a window to do the rendering work. now in order to do that, we are suppose to have xdm ran on startup. that's where our question is, (because we are not very good at linux) how do we startup the xdm? we do have a login screen in each of the machine. tried setting the usr/bin/env DISPLAY=localhost:0 while starting the parallel pvserver, but still the window alert of parallel rendering disabled persist. (i guess it's because the node 0 still can't open the window in each process's node.) any diagnostic test for us to try out? thanks Date: Thu, 6 Aug 2009 09:45:39 -0400 Subject: Re: [Paraview] How to set up X connection From: berk.gev...@kitware.com To: lcp8...@msn.com CC: paraview@paraview.org; kitw...@kitware.com HI Jess, So X is not running on these machines at all? -berk On Wed, Aug 5, 2009 at 5:34 AM, chew pinglcp8...@msn.com wrote: Hi all, i'm am doing parallel rendering using 4 machines, all these machines are identical to each other (architecture, OS, processor...) All of them have graphics hardware (rendering hardware). i came across from other resources saying that i have to set up xdm on each machine. Any idea on how to do that? thanks in advanced! regards, Jess ___ 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://www.paraview.org/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://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] xdmf/hdf5 errors (?)
Yes. paraview.org/Bug. Utkarsh is working on cleaning up the Xdmf so he'll probably look at this soon. -berk On Thu, Aug 13, 2009 at 6:24 AM, Verhoeven, Niek (UK)niek.verhoe...@baesystems.com wrote: Hello Is there a specific bug reporting site? Test conducted on: XP Professional 5.1 Service Pack 2.0 dual core laptop with 2G mem. Paraview versions 3.4.0 and 3.6.1 both win32 both are downloaded release versions. A simple method of generating a xmf hdf5 files example is: start paraview select menu [Source] choose 'plane set X resolution to 10 set Y resolution to 10 press Apply select menu [File] [Save Data] set Files of type: Xdmf Data File (*.xmf) set File name: MyPlane press OK On 3.4.0: A window with messages and warnings does pop up. Warning: In ..\..\..\ParaView3\VTK\Filtering\vtkDataSet.cxx, line 414 vtkUnstructuredGrid (0E88A928): Point array TextureCoordinates with 1 components, has 242 tuples but there are only 121 points See attached txt file On 3.6.1 paraview crashes. The result is two files MyPlane.xmf and MyPlane.h5 with for 3.6.1 MyPlane.xmf empty. Using either 3.4.0 or 3.61. to open MyPlane.xmf (as generated with 3.4.0) I used 3.6.1: sequence: select menu [File] [Open] set Files of type: Xdmf Reader (*.xmf *.xdmf) browse to MyPlane.xmf and select and press ok press Apply This results in an pop-up window with warning Warning: In ..\..\..\src\VTK\Filtering\vtkDataSet.cxx, line 414 vtkUnstructuredGrid (1049A0C0): Point array TextureCoordinates with 1 components, has 242 tuples but there are only 121 points And a very colourful plane. Something must have gone wrong with writing and reading of the TextureCoordinates on 3.4.0 Paraview 3.6.1 unfortunately crashes on writing. It would be interesting to see if other people have similar problems. Regards Niek Registered Company Details: BAE Systems (Operations) Limited Registered Office: Warwick House, PO Box 87, Farnborough Aerospace Centre, Farnborough, Hants, GU14 6YU, UK. Registered in England Wales No: 1996687 This email and any attachments are confidential to the intended recipient and may also be privileged. If you are not the intended recipient please delete it from your system and notify the sender. You should not copy it or use it for any purpose nor disclose or distribute its contents to any other person. ___ 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://www.paraview.org/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://www.paraview.org/mailman/listinfo/paraview
Re: [Paraview] FW: paraview3.6.1/OSMesa rendering issue
Hi Lucy, Any solution to this yet? Did you try building your own Mesa? -berk On Thu, Aug 13, 2009 at 3:43 AM, Lucy Gagliardilgaglia...@brawngp.com wrote: Hi, I’ve recently installed Paraview3.6.1 with OSMesa 7.4.4 on a linux based workstation (64bit). However, when running python scripts to set up and output jpg (or png) files using pvbatch with the --use-offscreen-rendering flag, the scripts run perfectly with no error, however the output picture is garbled. I’ve tried a fresh install of both Paraview and OSMesa to double check my configuration and the compilation runs through fine with no errors. I’m trying to use the OSMesa libraries installed in a shared location across a network instead of on my local machine – could this be the source of the problem? The python scripts only use Paraview as serial (i.e no parallel server) and output the pictures fine if the --use-offscreen-rendering flag is not used (and the OpenGL visualisation window opens). Any help with this would be great as I’m pretty stumped. Cheers, Lucy Lucy Gagliardi CFD Engineer T + 44 (0)1280 846880 www.BrawnGP.com Brawn GP Limited. Registered in England No. 787446. Registered Office at Brawn GP Formula One Team, Brackley, Northants NN13 7BD. Please consider the environment before printing this email ___ 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://www.paraview.org/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://www.paraview.org/mailman/listinfo/paraview