<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">And I now get the following error when I try to run this in parallel.<div><br></div><div>Note, Blaise, I do not have coordinates, so perhaps you could check that you still get divergence w/o coordinates.  (I don't need anything else form you)</div><div><br></div><div>Mark</div><div><br></div><div><div>[1]PETSC ERROR: --------------------- Error Message ------------------------------------</div><div>[1]PETSC ERROR: Nonconforming object sizes!</div><div>[1]PETSC ERROR: Sum of local lengths 29784 does not equal global length 29791, my local length 2482</div><div>  likely a call to VecSetSizes() or MatSetSizes() is wrong.</div><div>See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#split!">http://www.mcs.anl.gov/petsc/documentation/faq.html#split!</a></div><div>[1]PETSC ERROR: ------------------------------------------------------------------------</div><div>[1]PETSC ERROR: Petsc Development HG revision: 73637941c41b51aa3e638c914e6cb05248ece8b0  HG Date: Mon Jun 11 14:55:24 2012 -0400</div><div>[1]PETSC ERROR: See docs/changes/index.html for recent updates.</div><div>[1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.</div><div>[1]PETSC ERROR: See docs/index.html for manual pages.</div><div>[1]PETSC ERROR: ------------------------------------------------------------------------</div><div>[1]PETSC ERROR: ./ex10 on a arch-maco named madams-macbk-3.local by markadams Tue Jun 12 17:31:55 2012</div><div>[1]PETSC ERROR: Libraries linked from /Users/markadams/Codes/petsc-dev/arch-macosx-gnu-g/lib</div><div>[1]PETSC ERROR: Configure run at Mon Jun 11 17:03:36 2012</div><div>[1]PETSC ERROR: Configure options CXX=/sw/lib/gcc4.4/bin/g++-4 CC=/sw/lib/gcc4.4/bin/gcc-4 CFLAGS=-malign-double CXXFLAGS=-malign-double FC=/sw/lib/gcc4.4/bin/gfortran --download-hypre=1 --download-parmetis=1 --download-metis=1 --download-triangle=1 --with-x=0 --with-debugging=1 --download-mpich=1 PETSC_ARCH=arch-macosx-gnu-g</div><div>[1]PETSC ERROR: ------------------------------------------------------------------------</div><div>[1]PETSC ERROR: PetscSplitOwnership() line 93 in /Users/markadams/Codes/petsc-dev/src/sys/utils/psplit.c</div><div>[1]PETSC ERROR: PetscLayoutSetUp() line 147 in /Users/markadams/Codes/petsc-dev/src/vec/vec/impls/mpi/pmap.c</div><div>[1]PETSC ERROR: MatMPIAIJSetPreallocation_MPIAIJ() line 3194 in /Users/markadams/Codes/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c</div></div><div><br></div><div><br></div><div><br><div><div>On Jun 12, 2012, at 5:19 PM, Blaise Bourdin wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi,<div><br></div><div>Sorry, I should have mentioned that the .h5 and xml files are not for public eyes. I am uploading the binary files created with -U_ksp_view_binary, but my internet connection sucks.</div><div><br></div><div>Jed, I put the file on your account on schur. can you access them? The datafiles.tgz archive contains the  the result files in hdf5 format. ksp_view_binary.tgz contains only the binary output of ksp_view_binary</div><div><br></div><div>If you give me until tonight, I can try to set a prefix one the vecs.</div><div><br></div><div>Thanks a lot for your help.</div><div><br></div><div>Blaise</div><div><br><div><div>On Jun 12, 2012, at 4:06 PM, Mark F. Adams wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On Jun 12, 2012, at 4:56 PM, Jed Brown wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div class="gmail_quote">On Tue, Jun 12, 2012 at 3:54 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>How is what?  The ex10 code or "I set the block size for matrices before SetFromOptions".</div><div><br></div><div>This how I create matrices with blocks sizes:</div><div><br></div><div><div>  ierr = MatCreate( wcomm, &Prol ); CHKERRQ(ierr);</div>
<div>  ierr = MatSetSizes(Prol,nloc*bs,nLocalSelected*col_bs,PETSC_DETERMINE,PETSC_DETERMINE); </div><div>  CHKERRQ(ierr);</div><div>  ierr = MatSetBlockSizes( Prol, bs, col_bs ); CHKERRQ(ierr);</div></div></blockquote></div>
<br><div>Sorry, I thought you were talking about vectors because the error is coming from a vector load.</div><div><br></div></blockquote><div><br></div><div>There is a vector load, at least ex10 loads a vec and its seems to work (so I assume Blaise put one in there).</div><br><blockquote type="cite"><div>Where is this input file? I'll grab it sort out the issue.</div>
</blockquote></div><br><div><div>You can get an archive containing the KSP and the data files at <a href="https://filestogeaux.lsu.edu/public/download.php?FILE=bourdin/52338Q0at2Y">https://filestogeaux.lsu.edu/public/download.php?FILE=bourdin/52338Q0at2Y</a></div></div><div><br></div><div>These are just standard "binaryoutput" files I think,</div><div><br></div><div>Mark</div></div></blockquote></div><br><div apple-content-edited="true">
<span class="Apple-style-span" style="font-size: 12px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">-- </div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Department of Mathematics and Center for Computation & Technology</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Louisiana State University, Baton Rouge, LA 70803, USA</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Tel. +1 (225) 578 1612, Fax  +1 (225) 578 4276 <a href="http://www.math.lsu.edu/~bourdin">http://www.math.lsu.edu/~bourdin</a></div><div><br class="khtml-block-placeholder"></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></div></span></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br></div></div></blockquote></div><br></div></body></html>