documented apfCavityOp.h
putting back the gradient by volume function