#include "petscdt.h" PetscErrorCode PetscDTAltVPullback(PetscInt N, PetscInt M, const PetscReal *L, PetscInt k, const PetscReal *w, PetscReal *Lstarw)
|N||- the dimension of the origin vector space of the linear transformation, M >= 0|
|M||- the dimension of the image vector space of the linear transformation, N >= 0|
|L||- a linear transformation, an [M x N] matrix in row-major format|
|k||- the *signed* degree k of the |k|-form w, -(min(M,N)) <= k <= min(M,N). A negative form degree indicates that the pullback should be conjugated by the Hodge star operator (see note).|
|w||- a |k|-form in the image space, size [M choose |k|]|
|Lstarw||- the pullback of w to a |k|-form in the origin space, size [N choose |k|]: (Lstarw)(v_1,...v_k) = w(L*v_1,...,L*v_k).|
Note: negative form degrees accomodate, e.g., H-div conforming vector fields. An H-div conforming vector field stores its degrees of freedom as (dx, dy, dz), like a 1-form, but its normal trace is integrated on faces, like a 2-form. The correct pullback then is to apply the Hodge star transformation from (M-2)-form to 2-form, pullback as a 2-form, then the inverse Hodge star transformation.