[petsc-dev] MetSetBlockSize issue
Barry Smith
bsmith at mcs.anl.gov
Mon Apr 23 20:09:10 CDT 2012
Yes, look, for example how ex32 is handled at the bottom of the makefile.
Thanks
Barry
On Apr 23, 2012, at 5:50 PM, Mark F. Adams wrote:
>
> On Apr 23, 2012, at 5:59 PM, Barry Smith wrote:
>
>>
>> It was changed a while ago that MatSetBlockSize() couldn't be set after the matrix was full instantiated. I guess that example did not get fixed because it is not listed in the makefile to run in the makeall.
>
> May I add it?
>
>>
>> I have fixed the example to call MatSetBlockSize() at a safe point.
>>
>> Barry
>>
>> On Apr 23, 2012, at 4:29 PM, Mark F. Adams wrote:
>>
>>> ex55.c in ksp is failing with:
>>>
>>> [0]PETSC ERROR: --------------------- Error Message ------------------------------------
>>> [0]PETSC ERROR: Arguments are incompatible!
>>> [0]PETSC ERROR: Cannot change block size 1 to 2!
>>> [0]PETSC ERROR: ------------------------------------------------------------------------
>>>
>>> on this line 57 of ex55.c:
>>>
>>> ierr = MatSetBlockSize(Amat,2); CHKERRQ(ierr);
>>>
>>> Any idea what happened here?
>>>
>>> Mark
>>
>>
>
More information about the petsc-dev
mailing list