On Mon, Feb 20, 2012 at 4:25 PM, Bisht, Gautam <span dir="ltr">&lt;<a href="mailto:bishtg@ornl.gov">bishtg@ornl.gov</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
I have a MPIAIJ matrix A as listed below:<br>
<br>
row 0: (36, 1)  (37, 1)  (42, 1)  (43, 1)<br>
row 1: (37, 1)  (38, 1)  (43, 1)  (44, 1)<br>
row 2: (38, 1)  (39, 1)  (44, 1)  (45, 1)<br>
row 3: (39, 1)  (40, 1)  (45, 1)  (46, 1)<br>
row 4: (40, 1)  (41, 1)  (46, 1)  (47, 1)<br>
<br>
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.<br></blockquote><div><br></div><div>Runs fine for me:</div>
<div><br></div><div><div>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</div><div>Matrix Object: 1 MPI processes</div>
<div>  type: mpiaij</div><div>row 0: (36, 1)  (37, 1)  (42, 1)  (43, 1) </div><div>row 1: (37, 1)  (38, 1)  (43, 1)  (44, 1) </div><div>row 2: (38, 1)  (39, 1)  (44, 1)  (45, 1) </div><div>row 3: (39, 1)  (40, 1)  (45, 1)  (46, 1) </div>
<div>row 4: (40, 1)  (41, 1)  (46, 1)  (47, 1) </div><div>Matrix Object: 1 MPI processes</div><div>  type: mpiaij</div><div>row 0: (36, 1)  (37, 1)  (42, 1)  (43, 1) </div><div>row 1: (37, 1)  (38, 1)  (43, 1)  (44, 1) </div>
<div>row 2: (38, 1)  (39, 1)  (44, 1)  (45, 1) </div><div>row 3: (39, 1)  (40, 1)  (45, 1)  (46, 1) </div><div>row 4: (40, 1)  (41, 1)  (46, 1)  (47, 1) </div><div>Matrix Object: 1 MPI processes</div><div>  type: mpiaij</div>
<div>row 0:</div><div>row 1:</div><div>row 2:</div><div>row 3:</div><div>row 4:</div><div>row 5:</div><div>row 6:</div><div>row 7:</div><div>row 8:</div><div>row 9:</div><div>row 10:</div><div>row 11:</div><div>row 12:</div>
<div>row 13:</div><div>row 14:</div><div>row 15:</div><div>row 16:</div><div>row 17:</div><div>row 18:</div><div>row 19:</div><div>row 20:</div><div>row 21:</div><div>row 22:</div><div>row 23:</div><div>row 24:</div><div>
row 25:</div><div>row 26:</div><div>row 27:</div><div>row 28:</div><div>row 29:</div><div>row 30:</div><div>row 31:</div><div>row 32:</div><div>row 33:</div><div>row 34:</div><div>row 35:</div><div>row 36: (0, 1) </div><div>
row 37: (0, 1)  (1, 1) </div><div>row 38: (1, 1)  (2, 1) </div><div>row 39: (2, 1)  (3, 1) </div><div>row 40: (3, 1)  (4, 1) </div><div>row 41: (4, 1) </div><div>row 42: (0, 1) </div><div>row 43: (0, 1)  (1, 1) </div><div>
row 44: (1, 1)  (2, 1) </div><div>row 45: (2, 1)  (3, 1) </div><div>row 46: (3, 1)  (4, 1) </div><div>row 47: (4, 1) </div><div>Matrix Object: 1 MPI processes</div><div>  type: mpiaij</div><div>row 0:</div><div>row 1:</div>
<div>row 2:</div><div>row 3:</div><div>row 4:</div><div>row 5:</div><div>row 6:</div><div>row 7:</div><div>row 8:</div><div>row 9:</div><div>row 10:</div><div>row 11:</div><div>row 12:</div><div>row 13:</div><div>row 14:</div>
<div>row 15:</div><div>row 16:</div><div>row 17:</div><div>row 18:</div><div>row 19:</div><div>row 20:</div><div>row 21:</div><div>row 22:</div><div>row 23:</div><div>row 24:</div><div>row 25:</div><div>row 26:</div><div>
row 27:</div><div>row 28:</div><div>row 29:</div><div>row 30:</div><div>row 31:</div><div>row 32:</div><div>row 33:</div><div>row 34:</div><div>row 35:</div><div>row 36: (0, 1) </div><div>row 37: (0, 1)  (1, 1) </div><div>
row 38: (1, 1)  (2, 1) </div><div>row 39: (2, 1)  (3, 1) </div><div>row 40: (3, 1)  (4, 1) </div><div>row 41: (4, 1) </div><div>row 42: (0, 1) </div><div>row 43: (0, 1)  (1, 1) </div><div>row 44: (1, 1)  (2, 1) </div><div>
row 45: (2, 1)  (3, 1) </div><div>row 46: (3, 1)  (4, 1) </div><div>row 47: (4, 1) </div></div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thanks,<br>
-Gautam.<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>