[petsc-users] Help with MatTranspose

[Matthew Knepley]

On Mon, Feb 20, 2012 at 4:25 PM, Bisht, Gautam <bishtg at ornl.gov> wrote:

> Hi,
>
> I have a MPIAIJ matrix A as listed below:
>
> row 0: (36, 1)  (37, 1)  (42, 1)  (43, 1)
> row 1: (37, 1)  (38, 1)  (43, 1)  (44, 1)
> row 2: (38, 1)  (39, 1)  (44, 1)  (45, 1)
> row 3: (39, 1)  (40, 1)  (45, 1)  (46, 1)
> row 4: (40, 1)  (41, 1)  (46, 1)  (47, 1)
>
> When I try to transpose the matrix, while running on 2 processors, the
> code crashes. Attached below is the F90 code. I would appreciate any help
> in figuring out why the code crashes.
>

Runs fine for me:

knepley:/PETSc3/petsc/petsc-dev$
/PETSc3/petsc/petsc-dev/arch-c-f90-debug/bin/mpiexec -n 2
/PETSc3/petsc/petsc-dev/arch-c-f90-debug/lib/mattrans-obj/mattrans -mat_view
Matrix Object: 1 MPI processes
  type: mpiaij
row 0: (36, 1)  (37, 1)  (42, 1)  (43, 1)
row 1: (37, 1)  (38, 1)  (43, 1)  (44, 1)
row 2: (38, 1)  (39, 1)  (44, 1)  (45, 1)
row 3: (39, 1)  (40, 1)  (45, 1)  (46, 1)
row 4: (40, 1)  (41, 1)  (46, 1)  (47, 1)
Matrix Object: 1 MPI processes
  type: mpiaij
row 0: (36, 1)  (37, 1)  (42, 1)  (43, 1)
row 1: (37, 1)  (38, 1)  (43, 1)  (44, 1)
row 2: (38, 1)  (39, 1)  (44, 1)  (45, 1)
row 3: (39, 1)  (40, 1)  (45, 1)  (46, 1)
row 4: (40, 1)  (41, 1)  (46, 1)  (47, 1)
Matrix Object: 1 MPI processes
  type: mpiaij
row 0:
row 1:
row 2:
row 3:
row 4:
row 5:
row 6:
row 7:
row 8:
row 9:
row 10:
row 11:
row 12:
row 13:
row 14:
row 15:
row 16:
row 17:
row 18:
row 19:
row 20:
row 21:
row 22:
row 23:
row 24:
row 25:
row 26:
row 27:
row 28:
row 29:
row 30:
row 31:
row 32:
row 33:
row 34:
row 35:
row 36: (0, 1)
row 37: (0, 1)  (1, 1)
row 38: (1, 1)  (2, 1)
row 39: (2, 1)  (3, 1)
row 40: (3, 1)  (4, 1)
row 41: (4, 1)
row 42: (0, 1)
row 43: (0, 1)  (1, 1)
row 44: (1, 1)  (2, 1)
row 45: (2, 1)  (3, 1)
row 46: (3, 1)  (4, 1)
row 47: (4, 1)
Matrix Object: 1 MPI processes
  type: mpiaij
row 0:
row 1:
row 2:
row 3:
row 4:
row 5:
row 6:
row 7:
row 8:
row 9:
row 10:
row 11:
row 12:
row 13:
row 14:
row 15:
row 16:
row 17:
row 18:
row 19:
row 20:
row 21:
row 22:
row 23:
row 24:
row 25:
row 26:
row 27:
row 28:
row 29:
row 30:
row 31:
row 32:
row 33:
row 34:
row 35:
row 36: (0, 1)
row 37: (0, 1)  (1, 1)
row 38: (1, 1)  (2, 1)
row 39: (2, 1)  (3, 1)
row 40: (3, 1)  (4, 1)
row 41: (4, 1)
row 42: (0, 1)
row 43: (0, 1)  (1, 1)
row 44: (1, 1)  (2, 1)
row 45: (2, 1)  (3, 1)
row 46: (3, 1)  (4, 1)
row 47: (4, 1)

   Matt


> Thanks,
> -Gautam.
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120220/0a417529/attachment-0001.htm>