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
