Hi, I increased the fill value but the problem remains. I'm sending the matrices to petsc-maint at mcs.anl.gov. Thanks alot /H?kan
On Thu, Oct 8, 2009 at 4:03 PM, Hong Zhang <hzhang at mcs.anl.gov> wrote: > > Please try followings: > > 1) only run the second pairs > ? and > 2) use a larger fill value, e.g., fill=5.0 or 10.0 > ? MatPtAP(K2, V2, MAT_INITIAL_MATRIX, fill, &VtKV); let us know what > happens. > > If it still crashes, send us your matrices > K2 and V2 in petsc binary format. > I'll test it. > > Hong > > On Thu, 8 Oct 2009, H??kan Jakobsson wrote: > >> 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 >
