[Trilinos-Users] Epetra_CrsMatrix-aztecOO: repeated indeces

[Simone Deparis]

Hi,
I have

  Epetra_CrsMatrix  A(Copy,MapA,MapA,3)
  FillComplete(MapV,MapV);

  Epetra_FEVector x(MapV);
  Epetra_FEVector y(MapV);

mapA has repeated indeces {0 1 2}, {2 3 4}
mapV has no reapeated indeces {0 1 2}, {3 4}

I would like to use aztecOO to solve Ax=y,

  Epetra_LinearProblem Problem(&A,&x,&b);
  AztecOO Solver(Problem);

but I get the following error
(probably because of the repeating indeces in MapA):

AZ_extract_comm_info: Received elements must be stored after
                   all 3 local elements

I then tried to copy A to B, a matrix without repeated indeces,
using the exporter,

  Epetra_CrsMatrix  B(Copy,MapV,MapBcol,3)
  Comm.Barrier();
  B.Export(A,*(A.Exporter()),Add);
  Comm.Barrier();

mapBcol: each processor has all the column {0 1 2 3 4 5}, {0 1 2 3 4 5}
The problem here is that "off-block" values are lost:

>> full(A)

ans =

     2    -2     0     0     0
    -1     2    -2     0     0
     0    -1     6    -4     0
     0     0    -2     4    -4
     0     0     0    -2     4

>> full(B)

ans =

     2    -2     0     0     0
    -1     2    -2     0     0
     0    -1     6     0     0
     0     0     0     4    -4
     0     0     0    -2     4


Somebody as a suggestion, please?
Thank you
Simone