Cross derivatives implementation in FEniCS

Integrate them by parts and use the Function.dx operator as shown in Expressions depending on spatially varying coefficients and their derivatives in advection-diffusion-PDE

1 Like