[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