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
