<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><br></div> Well it may be related in some loose sense. But not particularly. Even though this matrix lives on one process it is stored as an MPIAIJ which has not check for symmetry.<div><br></div><div> In general for sparse matrices it is best never to check if the matrix is symmetric.</div><div><br></div><div> Barry</div><div><br><div><div>On Sep 1, 2010, at 10:56 AM, Leo van Kampenhout wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><br>Barry, do you think this is related? I want to check if the same matrix A is symmetric (on a single-core run) with the method <a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Mat/MatIsSymmetric.html">MatIsSymmetric</a> and it gives me the following error:<br>
<br>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>[0]PETSC ERROR: No support for this operation for this object type!<br>[0]PETSC ERROR: Matrix of type <mpiaij> does not support checking for symmetric!<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: Petsc Development HG revision: f455c9a70d2481389bc32b93c86e379cadcf41e1 HG Date: Wed Sep 01 09:04:54 2010 -0400<br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[0]PETSC ERROR: See docs/index.html for manual pages.<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: ./main on a arch-linu named wingtip72 by csg4035 Wed Sep 1 16:54:47 2010<br>[0]PETSC ERROR: Libraries linked from /net/users/csg/csg4035/install/petsc-dev/arch-linux-gnu-c-debug/lib<br>[0]PETSC ERROR: Configure run at Wed Sep 1 16:36:01 2010<br>
[0]PETSC ERROR: Configure options <br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: MatIsSymmetric() line 7214 in src/mat/interface/matrix.c<br><br><br><br><br>
<br><div class="gmail_quote">2010/9/1 Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div style="word-wrap: break-word;"><div> if (size > 1) { </div><div> /* change viewer to display matrix in natural ordering */</div><div> ierr = MatShellSetOperation(A, MATOP_VIEW, (void (*)(void)) MatView_MPI_DA);CHKERRQ(ierr);</div>
<div> /* turn off loading of matrix because loading would require proper permutation I don't feel like writing now */</div><div> ierr = MatShellSetOperation(A, MATOP_LOAD, (void (*)(void)) 0);CHKERRQ(ierr);</div>
<div> }</div><div><br></div><div><br></div><div> Shri,</div><div><br></div><div> Could you please add the support for this? And push to petsc-dev</div><div><br></div><div> Thanks</div><div><br></div><font color="#888888"><div>
Barry</div></font><div><div></div><div class="h5"><div><br></div><div><br></div><div><div>On Sep 1, 2010, at 10:23 AM, Leo van Kampenhout wrote:</div><br><blockquote type="cite"><br>Using a DA, of equal size and type as the other one. The curious thing is that there are no problems on a single core run. (but that is not what I want ;)<br>
<br><br><br><div class="gmail_quote">2010/9/1 Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div style="word-wrap: break-word;"><div><br></div> Did you get the matrix in the writer code using a DA or some other code that puts the vector from the 2d or 3d problem in the natural ordering on the binary file.<div>
<br></div><font color="#888888"><div> Barry</div></font><div><div></div><div><div><br><div><div>On Sep 1, 2010, at 9:55 AM, Leo van Kampenhout wrote:</div><br><blockquote type="cite">Hi all, <br><br>I'm having trouble with MatLoad on Petsc-Dev. Not sure if this problem also exists in the main PETSc 3.1 since my code won't compile on that. <br>
The problematic calls are as follows: <br><br><span style="font-family: courier new,monospace;"> call PetscViewerBinaryOpen(PETSC_COMM_WORLD,file2,</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> & FILE_MODE_READ,fd,ierr)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> call DAGetMatrix(da,MATMPIAIJ,A,ierr)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> </span><span style="font-family: courier new,monospace;">call MatLoad(A,fd,ierr)</span><br><br>which supposedly will load a matrix A which is stored by the following commands:<br>
<br><span style="font-family: courier new,monospace;"> call PetscViewerBinaryOpen(PETSC_COMM_WORLD,file2,</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> & FILE_MODE_WRITE,fd,ierr)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> call MatView(A,fd,ierr)</span><br><br>by another program (on a single core). The MatLoad() does work perfectly on a single core, however two or more cores make the program crash with the following error (see below)<br>
The programcall is <br><div style="margin-left: 40px;"><span style="font-family: courier new,monospace;"> mpirun -n 2 ./main</span><br></div><br>Is something broken with Petsc or do I need to set some extra parameters? <br>
Thanks in advance, <br><br>Leo van Kampenhout<br><br>PS. I'm using a build from last friday, since when I updated petsc-dev today using Mercurial I couldnt get it to configure or compile :S (which is another subject)<br>
<br><br>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>[0]PETSC ERROR: No support for this operation for this object type!<br>[0]PETSC ERROR: MatLoad is not supported for type: mpiaij!<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: Petsc Development HG revision: 84fec3ca75412473b06251029b58f2a68f1334f3 HG Date: Fri Aug 27 07:06:08 2010 +0200<br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[0]PETSC ERROR: See docs/index.html for manual pages.<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: ./main on a arch-linu named wingtip72 by csg4035 Wed Sep 1 15:52:11 2010<br>[0]PETSC ERROR: Libraries linked from /net/users/csg/csg4035/install/petsc-dev/arch-linux-gnu-c-debug/lib<br>[0]PETSC ERROR: Configure run at Fri Aug 27 08:37:32 2010<br>
[0]PETSC ERROR: Configure options <br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: MatLoad() line 843 in src/mat/interface/matrix.c<br><br><br>
</blockquote></div><br></div></div></div></div></blockquote></div><br>
</blockquote></div><br></div></div></div></blockquote></div><br>
</blockquote></div><br></div></body></html>