further about PCComputeExplicitOperator()

[Yujie]

when 1 processor is used, the matrix M in PCComputeExplicitOperator(pc,&M)
uses MATSEQDENSE type. Now, I want to use MATSEQAIJ, I change the codes as
follows:
1563 if (size == 1) {
1564 //05/01/08
1565 //ierr = MatSetType(*mat,MATSEQDENSE);CHKERRQ(ierr);
1566 //ierr = MatSeqDenseSetPreallocation(*mat,PETSC_NULL);CHKERRQ(ierr);
1567 ierr = MatSetType(*mat,MATSEQAIJ);CHKERRQ(ierr);
1568 ierr = MatSeqAIJSetPreallocation(*mat,0,PETSC_NULL);CHKERRQ(ierr);
1569
1570 } else {
1571 ierr = MatSetType(*mat,MATMPIAIJ);CHKERRQ(ierr);
1572 ierr =
MatMPIAIJSetPreallocation(*mat,0,PETSC_NULL,0,PETSC_NULL);CHKERRQ(ierr);
1573 }

PCApply is fast when running. However,  MatSetValues() is very very slow
when some arraies need to set. I find that the problem likely lies in
MatSeqXAIJReallocateAIJ(A,A->rmap.n,1,nrow,row,col,rmax,aa,ai,aj,rp,ap,imax,nonew,MatScalar)
in MatSetValues_SeqAIJ() after debugging the codes.
I can't further figure out where is the problem. Because it is
difficult to debug. Could you give me some advice? thanks a lot.

the version of PETSc is 2.3.3-p8.

thanks a lot.

Regards,
Yujie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080501/5d7f3faf/attachment.htm>