[petsc-dev] MetSetBlockSize issue

Mark F. Adams mark.adams at columbia.edu
Mon Apr 23 17:50:45 CDT 2012


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