[Paraview] pvti point data without sharing‏

[Amir A]

(The previous message got scrambled ...)

Hi Guys:

I am trying to use VTK pvti format to write distributed data and then visualize it in paraview. 
Chunks of data are distributed among processes and there is no ghost element shared (as my algorithm does not need it). 
Now the problem is that the point data format in VTK seems to need overlapping data (?). 
This is problematic since the processes do not share any data. 
I tried with cell data type and that one worked without the need for this sharing. 

Am I missing sth with the point data or is it supposed to be like this?


So for example assume I have 2 processes each having two 5x5 slabs (4 slabs in total). 
The obvious way that I supposed would work is as follows: 

test.pvti:
{
<?xml version="1.0"?>
<VTKFile type="PImageData" version="0.1" byte_order="LittleEndian">
    <PImageData WholeExtent="0 4 0 4 0 3" GhostLevel="0" Origin="0 0 0" Spacing="1 1 1">
        <PPointData>
            <PDataArray Name="value" type="Float64" format="ascii">
            </PDataArray>
        </PPointData> 
                    
        <Piece Extent="0 4 0 4 0 1" Source="test_00.vti"/>
        <Piece Extent="0 4 0 4 2 3" Source="test_01.vti"/>
</PImageData>       
</VTKFile>
}



test_00.vti:
{
<?xml version="1.0"?>
<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian">
    <ImageData WholeExtent="0 4 0 4 0 1" Origin="0 0 0" Spacing="1 1 1">
        <Piece Extent="0 4 0 4 0 1">
            <PointData Scalars="value">
                <DataArray Name="value" type="Float64" format="ascii">
                    0 0 0 0 0
                    0 0 0 0 0
                    0 0 0 0 0
                    0 0 0 0 0
                    0 0 0 0 0
                    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
                </DataArray>
            </PointData>
        </Piece>
    </ImageData>
</VTKFile>
}


test_01.vti:
{
<?xml version="1.0"?>
<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian">
    <ImageData WholeExtent="0 4 0 4 2 3" Origin="0 0 0" Spacing="1 1 1">
        <Piece Extent="0 4 0 4 2 3">
            <PointData Scalars="value">
                <DataArray Name="value" type="Float64" format="ascii">
                    2 2 2 2 2
                    2 2 2 2 2
                    2 2 2 2 2
                    2 2 2 2 2
                    2 2 2 2 2
                    3 3 3 3 3
                    3 3 3 3 3
                    3 3 3 3 3
                    3 3 3 3 3
                    3 3 3 3 3
                </DataArray>
            </PointData>
        </Piece>
    </ImageData>
</VTKFile>
}




But this does not work. Paraview loads it but then does not show anything. 
It would only work if I share one of the slabs, that is share the 111... slab or the 222... slab. 

Any help is appreciated. 

You can download the files at:

https://www.dropbox.com/sh/8lthk9ssxpyakbt/AAAdrHy9f8IGbWmtooJ3O5uoa




Regards, 
-Amir