[petsc-users] Assemble a Jacobian for a DMComposite object.

[(Rebecca) Xuefei YUAN]

Dear all,

I am working on assembling a Jacobian for a DMComposite object.

In FormJacobian(), I am not sure how to use MatSetValuesStencil() in  
this case, so I use MatSetValues() to assemble this Jacobian row by row.

MatSetValues(Mat mat,PetscInt m,const PetscInt idxm[],PetscInt n,const  
PetscInt idxn[],const PetscScalar v[],InsertMode addv)

To specify the row and columns' global indices idxm[] and idxn[], I  
need processor (rank=i) to receive (info.xm*info.ym) from each  
processor (rank<i), totally n(=rank) ints,  and send its  
info.xm*info.ym to each processor (rank<i), totally m (=size-rank-1).

Is there an existing call for that?

Thanks a lot!

Rebecca




-- 
(Rebecca) Xuefei YUAN
Department of Applied Physics and Applied Mathematics
Columbia University
Tel:917-399-8032
www.columbia.edu/~xy2102