On Wed, Jun 13, 2012 at 3:57 AM, Mark F. Adams <span dir="ltr">&lt;<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word"><div>This seems to be the first call to VecSetBlockSize:</div></div></blockquote><div><br></div><div>If that were true, how could we ever change the block size for any vector?</div><div><br>
</div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Breakpoint 1 at 0x1dbcee: file /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c, line 1523.</div>
<div>(gdb) r</div><div>Starting program: /Users/markadams/Desktop/GAMG/ex10 -f0 binaryoutput-sym -vecload_block_size 3</div><div>Reading symbols for shared libraries .++++++ done</div><div><br></div><div>Breakpoint 1, VecSetBlockSize (v=0x18c1850, bs=3) at /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c:1523</div>
<div>1523      PetscFunctionBegin;</div><div>(gdb) where</div><div>#0  VecSetBlockSize (v=0x18c1850, bs=3) at /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c:1523</div><div>#1  0x00246e6f in VecLoad_Binary (vec=0x18c1850, viewer=0x180d850) at /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c:104</div>
<div>#2  0x00247c85 in VecLoad_Default (newvec=0x18c1850, viewer=0x180d850) at /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c:348</div><div>#3  0x001d7ed8 in VecLoad (newvec=0x18c1850, viewer=0x180d850) at /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c:1111</div>
<div>#4  0x00003a91 in main ()</div><div><div>(gdb) n</div><div>1524      PetscValidHeaderSpecific(v,VEC_CLASSID,1);</div><div>(gdb) </div><div>1525      if (bs == v-&gt;map-&gt;bs) PetscFunctionReturn(0);</div><div>(gdb) </div>
<div>1526      ierr = PetscLayoutSetBlockSize(v-&gt;map,bs);CHKERRQ(ierr);</div><div>(gdb) </div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------------</div><div>[0]PETSC ERROR: Arguments are incompatible!</div>
<div>[0]PETSC ERROR: Cannot change block size 1 to 3!</div><div>[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: Petsc Development HG revision: 73637941c41b51aa3e638c914e6cb05248ece8b0  HG Date: Mon Jun 11 14:55:24 2012 -0400</div>
<div>[0]PETSC ERROR: See docs/changes/index.html for recent updates.</div><div>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.</div><div>[0]PETSC ERROR: See docs/index.html for manual pages.</div><div>
[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[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</div>
<div>[0]PETSC ERROR: Libraries linked from /Users/markadams/Codes/petsc-dev/arch-macosx-gnu-g/lib</div><div>[0]PETSC ERROR: Configure run at Mon Jun 11 17:03:36 2012</div><div>[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</div>
<div>[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: PetscLayoutSetBlockSize() line 460 in /Users/markadams/Codes/petsc-dev/src/vec/vec/impls/mpi/pmap.c</div>
<div>[0]PETSC ERROR: VecSetBlockSize() line 1526 in /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c</div></div><div><br></div><div><br></div><div><div>On Jun 12, 2012, at 3:42 PM, Matthew Knepley wrote:</div>
<br><blockquote type="cite">On Wed, Jun 13, 2012 at 3:45 AM, Mark F. Adams <span dir="ltr">&lt;<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>&gt;</span> wrote:<br><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
&gt;<br>
&gt; I thought there would be a -vec_block_size option, but alas its missing. Can you check in the debugger<br>
&gt; when the block size is getting set? I cannot find it in VecSetFromOptions().<br>
<br>
Its in the stack trace right?<br></blockquote><div><br></div><div>That is the second time, I want the first.</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">


Thanks,<br>
mark<br>
<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: PetscLayoutSetBlockSize() line 460 in /Users/markadams/Codes/petsc-dev/src/vec/vec/impls/mpi/pmap.c<br>
[0]PETSC ERROR: VecSetBlockSize() line 1526 in /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c<br>
[0]PETSC ERROR: VecLoad_Binary() line 104 in /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c<br>
[0]PETSC ERROR: VecLoad_Default() line 348 in /Users/markadams/Codes/petsc-dev/src/vec/vec/utils/vecio.c<br>
[0]PETSC ERROR: VecLoad() line 1111 in /Users/markadams/Codes/petsc-dev/src/vec/vec/interface/vector.c<br>
[0]PETSC ERROR: main() line 132 in src/ksp/ksp/examples/tutorials/ex10.c<br>
<br>
<br>
</blockquote></div><br><br clear="all"><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>

-- Norbert Wiener<br>
</font></span></blockquote></div><br></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>