problem with MatPtAP

[Hong Zhang]

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
>