[petsc-dev] error reading binary file

[Matthew Knepley]

On Wed, Jun 13, 2012 at 3:57 AM, Mark F. Adams <mark.adams at columbia.edu>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 <mark.adams at columbia.edu>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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120613/956540e0/attachment.html>