[Python-apps-team] Bug#478359: Bug#478359: mayavi2: fails to read blocks other than the first, in multiblock PLOT3D files

[Varun Hiremath]

Hi Francesco,

On Tue, 29 Apr, 2008 at 12:28:22AM +0200, Francesco Poli wrote:
> Package: mayavi2
> Version: 2.1.1-2
> Severity: normal
> 
> Hi!
> 
> I'm unable to read any block, except the first, in multiblock
> PLOT3D files.
> I can reproduce this bug with any multiblock PLOT3D file having more
> than one block.  For instance we can consider the test case I prepared
> for bug #400406 (see [1] and [2]).
> 
> [1] http://bugs.debian.org/400406
> [2] http://bugs.debian.org/cgi-bin/bugreport.cgi/vtk-multiblockPLOT3D-bug-test.tar.gz?bug=400406;msg=5;att=1
> 
> That test case is based on a little PLOT3D grid ('tiny.xyz')
> and solution ('tiny.q'), generated by a simple Fortran program
> ('gen_tinyplot3d.f') on a little-endian machine.
> The PLOT3D grid is composed of 5 blocks, each one of which consists
> of a single cubic cell.
> 
>   $ tar xvzf vtk-multiblockPLOT3D-bug-test.tar.gz
>   $ ls
>   gen_tinyplot3d.f  tiny.q    vtk-multiblockPLOT3D-bug-test.tar.gz
>   list_plot3d.py    tiny.xyz
>   $ mayavi2
> 
> The main window shows up.
> 
>  * from menus, select File -> Load data -> PLOT3D file
>  * as XYZ Co-ordinate file, select `tiny.xyz'
>  * as Q Solution file, select `tiny.q'
>  * select PLOT3D:tiny.xyz, tiny.q  in the MayaVi pane
>  * in the PLOT3DReader tab (inside the object editor), check "Binary
> file", "Do not reduce number of outputs", "Has byte count", and "Multi
> grid"; *un*check "Force read", "I blanking", and "Two dimensional
> geometry"; select "little_endian" Byte order
>  * in the Reader tab, click on the "Update Reader" button
> 
> OK, now we are ready to visualize the data...
> 
>  * from menus, select Visualize -> Modules -> Surface (a blue square
> shows up in the TVTK scene)
>  * from menus, select Visualize -> Modules -> Axes (coordinate ranges
> show up, confirming that we are seeing the first grid block, where
> 1 <= X <= 2)
> 
> OK, back to the PLOT3DReader tab: where's the option to select which
> grid blocks I want to load?!?
> I really cannot find it...
> 
> I would like to have a means to decide which grid blocks will be loaded
> (e.g.: I want to load block 1, 2, and 4, but not the remaining two),
> so that a number of PLOT3D:tiny.xyz, tiny.q objects are automatically
> created, each one corresponding to one of the requested grid blocks
> (e.g.: three objects are created, corresponding to blocks 1, 2, and 4).
> 
> Please note that this decision may be taken only after PLOT3D properties
> (such as "Binary file", "Has byte count", and so forth) have been
> specified, or otherwise the total number of blocks present in the
> PLOT3D file will not be computed correctly...

Thanks for this detailed report. It looks more like a wishlist
bug. I'm forwarding this upstream.

Hi Prabhu and Gael, 

Could you please have a look at this bug report:
http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=478359

Thanks,
Varun

-- 
Varun Hiremath
Undergraduate Student,
Aerospace Engineering Department,
Indian Institute of Technology Madras,
Chennai, India
-------------------------------------
Homepage: http://varun.travisbsd.org