[Trilinos-Users] Help, A Bug in HDF5::Read for Expetra_CrsGraph?!

[Davood Ansari]

Hi all

I am using Trilinos 8.0.7 on a dual-core machine with ubuntu linux and intel
compilers.

I am using ExpteraExt's HDF5 functionality to read and write a
Epetra_CrsGraph between two codes.
The first one, does very well and the exported graph can be retrieved for
example in Matlab using the
following script:

NumGlobalRows = double(hdf5read('CRSGraph.h5', '/Graph/NumGlobalRows/'));
NumGlobalCols = double(hdf5read('CRSGraph.h5', '/Graph/NumGlobalCols/'));
NumGlobalNonzeros = double(hdf5read('CRSGraph.h5',
'/Graph/NumGlobalNonzeros/'));
ROW = double(hdf5read('CRSGraph.h5', '/Graph/ROW/'));
COL = double(hdf5read('CRSGraph.h5', '/Graph/COL/'));
VAL=zeros(NumGlobalNonzeros,1) ;
VAL(:)= 1.;
A = sparse(ROW + 1, COL + 1, VAL, NumGlobalRows, NumGlobalCols);

The problem occurs when I try to import the graph in another code using the
following lines:

Epetra_CrsGraph::Epetra_CrsGraph* CrsGraph ;
HDF5.Read("Graph", CrsGraph) ;
cout << *CrsGraph ;
HDF5.Close();

In the output, attached after this, I can see that the line no. 15 (the
beginning row of the second process) has got all its columns missing.
Is this a bug or am I coding improperly?

Kindly Comment

Epetra::CrsGraph

Number of Global Block Rows  = 30
Number of Global Block Cols  = 30
Number of Global Block Diags = 28
Number of Global Entries     = 160

Number of Global Rows        = 30
Number of Global Cols        = 30
Number of Global Diagonals   = 28
Number of Global Nonzeros    = 160

Global Maximum Block Row Dim = 1
Global Maximum Block Col Dim = 1
Global Maximum Num Indices   = 9

Number of My Block Rows  = 15
Number of My Block Cols  = 29
Number of My Block Diags = 15
Number of My Entries     = 76

Number of My Rows        = 15
Number of My Cols        = 29
Number of My Diagonals   = 15
Number of My Nonzeros    = 76

My Maximum Block Row Dim = 1
My Maximum Block Col Dim = 1
My Maximum Num Indices   = 9

       Row Index     Col Index         Col Index         Col Index
Col Index         Col Index         Col Index         Col Index         Col
Index         Col Index
             0               0               6               8
14
             1               1               7              10
20
             2               2               9              12
18
             3               3              11              13
             4               4              15              16
25
             5               5              17              19
27
             6               0               6               7
8              20              21
             7               1               6               7
20
             8               0               6               8
9              14              21              24              26
             9               2               8               9
12              18              21              26
23              28
            10               1              10              11
20              22
            11               3              10        Epetra::CrsGraph

Number of My Block Rows  = 15
Number of My Block Cols  = 28
Number of My Block Diags = 13
Number of My Entries     = 84

Number of My Rows        = 15
Number of My Cols        = 28
Number of My Diagonals   = 13
Number of My Nonzeros    = 84

My Maximum Block Row Dim = 1
My Maximum Block Col Dim = 1
My Maximum Num Indices   = 8

       Row Index     Col Index         Col Index         Col Index
Col Index         Col Index         Col Index         Col Index         Col
Index
            15
            16              27               4
            17              16              17              27
5
            18              18              19              28
2               9
            19              18              19              27
28               5
            20              20              21              22
1               6               7              10
            21              20              21              22
23               9               6               8
            22              20              21              22
23              10              11
            23              21              22              23
9              11              12              13
            24              15              24              25
26              29               8              14
            25              15              16              24
25              27              29               4
            26              24              26              28
29               9               8
            27              16              17              19
25              27              28              29               5
            28              18              19              26
27              28              29               9
            29              24              25              26
27              28              29
      11              13              23              22
            12               2               9              12
13              23
            13               3              11              12
13              23
            14               0               8              14
15              24

    Successfully Finished
    Started at:    Fri Dec  5 14:42:40 2008
    Ended   at:    Fri Dec  5 14:42:40 2008
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20081205/73b05514/attachment-0001.html