i've been preallocating the full matrix
const PetscScalar zerotol = 1.e-6;
PetscInt M,N,m,n;
ierr = MatGetSize(MatDense,&M,&N);CHKERRQ(ierr);
ierr = MatGetLocalSize(MatDense,&m,&n);CHKERRQ(ierr);
ierr = MatCreateMPIAIJ(PETSC_COMM_WORLD,m,n,M,N,n,PETSC_NULL,
also, MatMatMult_MPIAIJ_MPIDense() and
MatMatMult_MPIDense_MPIAIJ()
seem to work if you don't want/need to convert both matrices.
df
On Sun, 5 Apr 2009, David Fuentes wrote:
> Hi Yujie,
>
>
> as a work around have you tried converting your dense
> matrices to aij format and using MatMatM
Hi Yujie,
as a work around have you tried converting your dense
matrices to aij format and using MatMatMult_MPIAIJ_MPIAIJ()??
df
On Fri, 3 Apr 2009, Yujie wrote:
>
> Dear Barry:
>
> I am trying to debug the codes you have written with ex123.c. After
> commenting the error output
>