And I now get the following error when I try to run this in parallel. 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)
Mark [1]PETSC ERROR: --------------------- Error Message ------------------------------------ [1]PETSC ERROR: Nonconforming object sizes! [1]PETSC ERROR: Sum of local lengths 29784 does not equal global length 29791, my local length 2482 likely a call to VecSetSizes() or MatSetSizes() is wrong. See http://www.mcs.anl.gov/petsc/documentation/faq.html#split! [1]PETSC ERROR: ------------------------------------------------------------------------ [1]PETSC ERROR: Petsc Development HG revision: 73637941c41b51aa3e638c914e6cb05248ece8b0 HG Date: Mon Jun 11 14:55:24 2012 -0400 [1]PETSC ERROR: See docs/changes/index.html for recent updates. [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting. [1]PETSC ERROR: See docs/index.html for manual pages. [1]PETSC ERROR: ------------------------------------------------------------------------ [1]PETSC ERROR: ./ex10 on a arch-maco named madams-macbk-3.local by markadams Tue Jun 12 17:31:55 2012 [1]PETSC ERROR: Libraries linked from /Users/markadams/Codes/petsc-dev/arch-macosx-gnu-g/lib [1]PETSC ERROR: Configure run at Mon Jun 11 17:03:36 2012 [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 [1]PETSC ERROR: ------------------------------------------------------------------------ [1]PETSC ERROR: PetscSplitOwnership() line 93 in /Users/markadams/Codes/petsc-dev/src/sys/utils/psplit.c [1]PETSC ERROR: PetscLayoutSetUp() line 147 in /Users/markadams/Codes/petsc-dev/src/vec/vec/impls/mpi/pmap.c [1]PETSC ERROR: MatMPIAIJSetPreallocation_MPIAIJ() line 3194 in /Users/markadams/Codes/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c On Jun 12, 2012, at 5:19 PM, Blaise Bourdin wrote: > Hi, > > 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. > > 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 > > If you give me until tonight, I can try to set a prefix one the vecs. > > Thanks a lot for your help. > > Blaise > > On Jun 12, 2012, at 4:06 PM, Mark F. Adams wrote: > >> >> On Jun 12, 2012, at 4:56 PM, Jed Brown wrote: >> >>> On Tue, Jun 12, 2012 at 3:54 PM, Mark F. Adams <mark.adams at columbia.edu> >>> wrote: >>> How is what? The ex10 code or "I set the block size for matrices before >>> SetFromOptions". >>> >>> This how I create matrices with blocks sizes: >>> >>> ierr = MatCreate( wcomm, &Prol ); CHKERRQ(ierr); >>> ierr = >>> MatSetSizes(Prol,nloc*bs,nLocalSelected*col_bs,PETSC_DETERMINE,PETSC_DETERMINE); >>> >>> CHKERRQ(ierr); >>> ierr = MatSetBlockSizes( Prol, bs, col_bs ); CHKERRQ(ierr); >>> >>> Sorry, I thought you were talking about vectors because the error is coming >>> from a vector load. >>> >> >> There is a vector load, at least ex10 loads a vec and its seems to work (so >> I assume Blaise put one in there). >> >>> Where is this input file? I'll grab it sort out the issue. >> >> You can get an archive containing the KSP and the data files at >> https://filestogeaux.lsu.edu/public/download.php?FILE=bourdin/52338Q0at2Y >> >> These are just standard "binaryoutput" files I think, >> >> Mark > > -- > Department of Mathematics and Center for Computation & Technology > Louisiana State University, Baton Rouge, LA 70803, USA > Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 > http://www.math.lsu.edu/~bourdin > > > > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120612/6051ba12/attachment.html>
