[petsc-users] Building a rectangular MPI matrix relating to PetscSections

[Alejandro D Otero]

Hi all, is it posible to build a rectangular matrix that relates two fields
with (posible) different number of DoF per point defining the nonzero
structure from 2 petscsections representing one field each?

As an example, I want to build an operator that given a 2D velocity
computes the corresponding vorticity by means of:

W = Bcurl * V

with V the vector with the velocity components at grid points and W the
vorticity at the same points.

Thanks in advance, regards,

Alejandro
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150713/d663eeb7/attachment.html>