[petsc-dev] Merge branch 'rlacroix/fix_bs_mpibaij_to_mpiaij'

Barry Smith bsmith at mcs.anl.gov
Mon Feb 2 11:55:35 CST 2015


  Hong,

   MatConvert is totally broken for bs > 1 and has been for a long time; I've just been living with it being broken. It is on Jed's list to fix 

  Barry


> On Feb 2, 2015, at 11:51 AM, Hong <hzhang at mcs.anl.gov> wrote:
> 
> Barry,
> Your change
> https://bitbucket.org/petsc/petsc/commits/fedbf777941bf3b58323d827ccd2c59330060afe
> 
> causes petsc/src/mat/examples/tests/ex55 (runex55_2, master) to crash.
> 
> In this example, a mpibaij matrix C with bs=2 is created, then converted to aij matrix A with bs=2 (without your change, A has default bs=1).
> 
> Then ex55 converts A (mpiaij with bs=2) to a baij matrix B (line 118) by calling 
> MatConvert_basic() -> crashes at line 154 of pmap.c
> if (!map->range) {
>     ierr = PetscMalloc1(size+1, &map->range);CHKERRQ(ierr);
>   }
> 
> I cannot figure out why this Malloc1() crashes when running with 
> mpiexec -n 2 ./ex55 -verbose 1 (i.e., size=2).
> 
> Any hint/suggestion?
> 
> Hong




More information about the petsc-dev mailing list