problem with MatPtAP

[Håkan Jakobsson]

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
>