PetscErrorCode PotSetHarm(Pot self, PetscScalar a) { Harm *ctx; PetscNew(&ctx); ctx->a = a; PFSet(self, HarmApply, NULL, HarmView, HarmDestroy, ctx); return 0; }