On Wed, Jun 13, 2012 at 3:57 AM, Mark F. Adams <[email protected]>wrote:
This seems to be the first call to VecSetBlockSize:
If that were true, how could we ever change the block size for any vector? Matt
Breakpoint 1 at 0x1dbcee: file /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c, line 1523. (gdb) r Starting program: /Users/markadams/Desktop/GAMG/ex10 -f0 binaryoutput-sym -vecload_block_size 3 Reading symbols for shared libraries .++++++ done
Breakpoint 1, VecSetBlockSize (v=0x18c1850, bs=3) at /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c:1523 1523 PetscFunctionBegin; (gdb) where #0 VecSetBlockSize (v=0x18c1850, bs=3) at /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c:1523 #1 0x00246e6f in VecLoad_Binary (vec=0x18c1850, viewer=0x180d850) at /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c:104 #2 0x00247c85 in VecLoad_Default (newvec=0x18c1850, viewer=0x180d850) at /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c:348 #3 0x001d7ed8 in VecLoad (newvec=0x18c1850, viewer=0x180d850) at /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c:1111 #4 0x00003a91 in main () (gdb) n 1524 PetscValidHeaderSpecific(v,VEC_CLASSID,1); (gdb) 1525 if (bs == v->map->bs) PetscFunctionReturn(0); (gdb) 1526 ierr = PetscLayoutSetBlockSize(v->map,bs);CHKERRQ(ierr); (gdb) [0]PETSC ERROR: --------------------- Error Message ------------------------------------ [0]PETSC ERROR: Arguments are incompatible! [0]PETSC ERROR: Cannot change block size 1 to 3! [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: Petsc Development HG revision: 73637941c41b51aa3e638c914e6cb05248ece8b0 HG Date: Mon Jun 11 14:55:24 2012 -0400 [0]PETSC ERROR: See docs/changes/index.html for recent updates. [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. [0]PETSC ERROR: See docs/index.html for manual pages. [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: /Users/markadams/Desktop/GAMG/ex10 on a arch-maco named madams-macbk-3.local by markadams Tue Jun 12 15:52:56 2012 [0]PETSC ERROR: Libraries linked from /Users/markadams/Codes/petsc-dev/arch-macosx-gnu-g/lib [0]PETSC ERROR: Configure run at Mon Jun 11 17:03:36 2012 [0]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 [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: PetscLayoutSetBlockSize() line 460 in /Users/markadams/Codes/petsc-dev/src/vec/vec/impls/mpi/pmap.c [0]PETSC ERROR: VecSetBlockSize() line 1526 in /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c
On Jun 12, 2012, at 3:42 PM, Matthew Knepley wrote:
On Wed, Jun 13, 2012 at 3:45 AM, Mark F. Adams <[email protected]>wrote:
I thought there would be a -vec_block_size option, but alas its
missing. Can you check in the debugger
when the block size is getting set? I cannot find it in VecSetFromOptions().
Its in the stack trace right?
That is the second time, I want the first.
Matt
Thanks, mark
[0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: PetscLayoutSetBlockSize() line 460 in /Users/markadams/Codes/petsc-dev/src/vec/vec/impls/mpi/pmap.c [0]PETSC ERROR: VecSetBlockSize() line 1526 in /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c [0]PETSC ERROR: VecLoad_Binary() line 104 in /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c [0]PETSC ERROR: VecLoad_Default() line 348 in /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c [0]PETSC ERROR: VecLoad() line 1111 in /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c [0]PETSC ERROR: main() line 132 in src/ksp/ksp/examples/tutorials/ex10.c
-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener