<html><head><style> body {height: 100%; color:#000000; font-size:12pt; font-family:Times New Roman;}</style></head><body><strongsmith@mcs.anl.gov><br><div>> That code cannot be right. When you copy the App to A it does not have the same nonzero pattern. You should not use a MatCopy() there you should use a MatHeaderMerge() or a MatHeaderReplace() whichever one is right.</div><br><div>Thanks,i was searching for a similar function but did not find it hence used MatCopy. I've fixed this and pushed it.</div><br><div> > Then try running with a big DA to see how long it takes.</div><div> </div><div>I tried MatLoad on Leo's 3D grid (202 X 102 X 102) but it takes forever to load. Hence, I tried it on a relatively smaller grid (100 X 100 X 100) and here are the log summary numbers that i got.</div><br><div><div>Event Count Time (sec) Flops --- Global --- --- Stage --- Total</div><div> Max Ratio Max Ratio Max Ratio Mess Avg len Reduct %T %F %M %L %R %T %F %M %L %R Mflop/s</div><div>------------------------------------------------------------------------------------------------------------------------</div><br><div>--- Event Stage 0: Main Stage</div><br><div>MatAssemblyBegin 4 1.0 2.6508e-0134.6 0.00e+00 0.0 0.0e+00 0.0e+00 6.0e+00 3 0 0 0 5 3 0 0 0 6 0</div><div>MatAssemblyEnd 4 1.0 4.2604e-01 1.0 0.00e+00 0.0 1.2e+01 2.0e+04 3.6e+01 9 0 63 1 28 9 0 63 1 38 0</div><div>MatGetSubMatrice 1 1.0 5.9222e-01 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 5.0e+00 13 0 0 0 4 13 0 0 0 5 0</div><div>MatLoad 1 1.0 3.2231e+00 1.0 0.00e+00 0.0 1.1e+01 4.0e+06 4.9e+01 70 0 58100 38 70 0 58100 52 0</div><div>------------------------------------------------------------------------------------------------------------------------</div></div><br><br><div>Shri</div><div> </div><div> Barry</div><div><br><div><div>On Sep 2, 2010, at 10:17 AM, Leo van Kampenhout wrote:</div><br><blockquote>Hello Shri,<br><br>thanks for the fix. The error is now gone but the loading of the Mat takes now infinitely long.. help!!<br><br>Leo<br><br><br><div class="gmail_quote">2010/9/2 Shri <span dir="ltr"><<a href="mailto:abhyshr@mcs.anl.gov" target="_blank">abhyshr@mcs.anl.gov</a>></span><br><br>><blockquote class="gmail_quote" style="border-left: 1px solid #cccccc; margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div>Pushed. Please let me know if you get any errors.<br><div><font color="#888888">Shri</font><div><br>><div class="h5"><br><br>----- Barry Smith wrote:<br>><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><br>><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><br>><div> }</div><br><br><div> Shri,</div><br><div> Could you please add the support for this? And push to petsc-dev</div><br><div> Thanks</div><br><div> Barry</div><br><br><div><div>On Sep 1, 2010, at 10:23 AM, Leo van Kampenhout wrote:</div><br>><br><blockquote><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><br>><br>><blockquote class="gmail_quote" style="border-left: 1px solid #cccccc; margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div style="word-wrap: break-word;"><br> 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>><br>><br></div><font color="#888888"><div> Barry</div></font><div><div><div><br><div><div>On Sep 1, 2010, at 9:55 AM, Leo van Kampenhout wrote:</div><br><blockquote>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><br>><br>> The problematic calls are as follows: <br><br><span style="font-family: courier new,monospace;"> call PetscViewerBinaryOpen(PETSC_COMM_WORLD,file2,</span><br><br>><span style="font-family: courier new,monospace;"> & FILE_MODE_READ,fd,ierr)</span><br><br>><span style="font-family: courier new,monospace;"> call DAGetMatrix(da,MATMPIAIJ,A,ierr)</span><br><br>> <br>><span style="font-family: courier new,monospace;"> </span><span style="font-family: courier new,monospace;">call MatLoad(A,fd,ierr)</span><br><br>><br>which supposedly will load a matrix A which is stored by the following commands:<br><br>> <br>><br><span style="font-family: courier new,monospace;"> call PetscViewerBinaryOpen(PETSC_COMM_WORLD,file2,</span><br><br>><span style="font-family: courier new,monospace;"> & FILE_MODE_WRITE,fd,ierr)</span><br><br>> <br>><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><br>><br>> <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><br>><br>> <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>> <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><br>><br>> <br>> [0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: Petsc Development HG revision: 84fec3ca75412473b06251029b58f2a68f1334f3 HG Date: Fri Aug 27 07:06:08 2010 +0200<br><br>><br>> <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><br>><br>> <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><br>><br>> <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><br>><br>></blockquote></div><br></div></div></div></div></blockquote></div><br><br>></blockquote></div><br><br></div></div></div></div></blockquote></div><br><br>></blockquote></div><br></div><br></bsmith@mcs.anl.gov></body></html>