problem with MatPtAP

[Håkan Jakobsson]

Hi,

I have run into some problems with using MatPtAP. I have two similar
matrix/projection matrix pairs (K1,V1) and (K2,V2). Dimensions are

K1 - 10227x10227, V1 - 10227x255
K2 - 10146x10146, V2 - 10146x336

For the first pair MatPtAP returns ok, but for the second one

ierr=MatPtAP(K, V, MAT_INITIAL_MATRIX, 1.0, &VtKV); CHKERRQ(ierr);

fails with

[0]PETSC ERROR: MatPtAPNumeric_SeqAIJ_SeqAIJ() line 245 in
src/mat/impls/aij/seq/matptap.c
[0]PETSC ERROR: MatPtAPNumeric_SeqAIJ() line 36 in
src/mat/impls/aij/seq/matptap.c
[0]PETSC ERROR: MatPtAPNumeric() line 7279 in src/mat/interface/matrix.c
[0]PETSC ERROR: MatPtAP_Basic() line 15 in src/mat/utils/ptap.c
[0]PETSC ERROR: MatPtAP() line 7223 in src/mat/interface/matrix.c

The same also applies when using MatPtAPNumeric. For

Mat VtKV;
ierr=MatPtAPSymbolic(K, V, 1.0, &VtKV); CHKERRQ(ierr);
ierr=MatPtAPNumeric(K, V, VtKV); CHKERRQ(ierr);

MatPtAPSymbolic returns ok, but MatPtAPNumeric fails with error message

[0]PETSC ERROR: MatPtAPNumeric_SeqAIJ_SeqAIJ() line 245 in
src/mat/impls/aij/seq/matptap.c
[0]PETSC ERROR: MatPtAPNumeric_SeqAIJ() line 36 in
src/mat/impls/aij/seq/matptap.c
[0]PETSC ERROR: MatPtAPNumeric() line 7279 in src/mat/interface/matrix.c

When splitting the operation into a MatMatMultTranspose and a
MatMatMult I don't get any errors.

Any ideas on what might be the going on here? PETSc version is 3.0.0-p8.

/Håkan