Awesome, thanks.  I ran into that bug and thought I was doing something wrong.


-Andrew


—
Andrew

On Wed, Dec 3, 2014 at 5:40 PM, Barry Smith <[email protected]> wrote:

>> On Dec 3, 2014, at 5:13 PM, Andrew Spott <[email protected]> wrote:
>> 
>> What is the easiest way to change the block size of a matrix?
>> 
>> I have a matrix that was saved with a block size of 1, and I would like to 
>> increase it upon load to a larger block size.  Is there a simple way of 
>> doing this?
>    Yes, declare it as a BAIJ before calling MatLoad() for example
>     Mat A;
>     MatCreate(comm,&A);
>     MatSetType(A,MATBAIJ);
>     MatLoad(A,viewer);
>     and pass the command line option -matload_block_size 2  
>     In the development version of PETSc you can skip the command line option 
> and instead use 
>     MatSetBlockSize(A,2);   before the call to MatLoad().    There was a bug 
> in the release version of PETSc that prevented this from previously working.
>    Barry
>> 
>> Thanks
>> 
>> -Andrew Spott
>> 

Reply via email to