Let $\mathrm{\π}\:E\to M$be a fiber bundle, with base dimension $n$ and fiber dimension $m$ and let ${\mathrm{\π}}^{\mathrm{\∞}}\:J^{\mathrm{\∞}}\left(E\right) \to M$be the infinite jet bundle of $E$. Let $\(x^i\, u^{\mathrm{\α}}\, u_{i_{}}^{\mathrm{\α}}\, u_{i_{}j}^{\mathrm{\α}}$, ..., $u_{\mathrm{ij} \cdot \cdot \cdot k}^{\mathrm{\α}}\, ....\)$be a local system of jet coordinates and let ${\mathrm{\Θ}}^{\mathrm{\α}} \= {\mathrm{du}}^{\mathrm{\α}}-u_{\mathrm{\ℓ}}^{\mathrm{\α}}{\mathrm{dx}}^{\mathrm{\ℓ}}$. Let ${\mathrm{\Omega }}^{\left(n\,s\right)}\left(J^{\infty }\left(E\right)\right)$ be the space of all differential bi-forms of horizontal degree$n$and vertical degree $s\.$Let $\mathrm{\ω} \in {\mathrm{\Omega }}^{\left(n\,s\right)}\left(J^{\infty }\left(E\right)\right)$and let $E_{\mathrm{\α}}\left(\mathrm{\ω}\right) \in {\mathrm{\Ω}}^{\left(n-1\,s\right)}\left(J^{\infty }\left(E\right)\right)$ be the components of the Euler-Lagrange operator applied to $\mathrm{\ω}$. Then the integration by parts operator $I\: {\mathrm{\Omega }}^{\left(n\,s\right)}\left(J^{\infty }\left(E\right)\right)\to {\mathrm{\Omega }}^{\left(n\,s\right)}\left(J^{\infty }\left(E\right)\right)$ is defined by

The command IntegrationByParts($\mathrm{\omega }$) returns the type$\left(n\, s\right)$bi-form $I\left(\mathrm{\ω}\right)$.

The command IntegrationByParts is part of the DifferentialGeometry:-JetCalculus package. It can be used in the form IntegrationByParts(...) only after executing the commands with(DifferentialGeometry) and with(JetCalculus), but can always be used by executing DifferentialGeometry:-JetCalculus:-IntegrationByParts(...).

with(DifferentialGeometry): with(JetCalculus):

DGsetup([x], [u], E, 3):

PDEtools[declare](a(x), b(x), c(x), quiet):

omega1 := Dx &wedge evalDG(a(x)*Cu[] + b(x)*Cu[1] + c(x)*Cu[1, 1] + d(x)*Cu[1, 1, 1]);

$\mathrm{\ω1}≔\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[1\,1\right]\right]\,\left[\left[\left[1\,2\right]\,a\right]\,\left[\left[1\,3\right]\,b\right]\,\left[\left[1\,4\right]\,c\right]\,\left[\left[1\,5\right]\,d\left(x\right)\right]\right]\right]\right)$

IntegrationByParts(omega1);

$\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[1\,1\right]\right]\,\left[\left[\left[1\,2\right]\,-d_{x,x,x}+c_{x,x}-b_x+a\right]\right]\right]\right)$

omega2 := Dx &wedge evalDG(a(x)*Cu[]&w Cu[1] + b(x)*Cu[] &w Cu[1,1] + c(x)*Cu[1] &w Cu[1,1]);

$\mathrm{\ω2}≔\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[1\,2\right]\right]\,\left[\left[\left[1\,2\,3\right]\,a\right]\,\left[\left[1\,2\,4\right]\,b\right]\,\left[\left[1\,3\,4\right]\,c\right]\right]\right]\right)$

omega3 := IntegrationByParts(omega2);

$\mathrm{\ω3}≔\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[1\,2\right]\right]\,\left[\left[\left[1\,2\,3\right]\,-\frac{c_{x,x}}{2}-b_x+a\right]\,\left[\left[1\,2\,4\right]\,-\frac{3c_x}{2}\right]\,\left[\left[1\,2\,5\right]\,-c\right]\right]\right]\right)$

IntegrationByParts(omega3);

$\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[1\,2\right]\right]\,\left[\left[\left[1\,2\,3\right]\,-\frac{c_{x,x}}{2}-b_x+a\right]\,\left[\left[1\,2\,4\right]\,-\frac{3c_x}{2}\right]\,\left[\left[1\,2\,5\right]\,-c\right]\right]\right]\right)$

DGsetup([x, y], [u, v], E, 3):

PDEtools[declare](a(x, y), b(x, y), c(x, y), d(x, y), e(x, y), f(x, y), quiet):

omega4 := Dx &wedge Dy &wedge evalDG(a(x, y)*Cu[] + b(x, y)*Cv[] + c(x, y)*Cu[1] + d(x, y)*Cu[2] + e(x, y)*Cv[1] + f(x, y)*Cv[2]);

$\mathrm{\ω4}≔\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[2\,1\right]\right]\,\left[\left[\left[1\,2\,3\right]\,a\right]\,\left[\left[1\,2\,4\right]\,b\right]\,\left[\left[1\,2\,5\right]\,c\right]\,\left[\left[1\,2\,6\right]\,d\right]\,\left[\left[1\,2\,7\right]\,e\right]\,\left[\left[1\,2\,8\right]\,f\right]\right]\right]\right)$

IntegrationByParts(omega4);

$\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[2\,1\right]\right]\,\left[\left[\left[1\,2\,3\right]\,-d_y-c_x+a\right]\,\left[\left[1\,2\,4\right]\,-f_y-e_x+b\right]\right]\right]\right)$

omega5 := Dx &wedge Dy &wedge evalDG(a(x, y)*Cu[1] &w Cv[1]);

$\mathrm{\ω5}≔\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[2\,2\right]\right]\,\left[\left[\left[1\,2\,5\,7\right]\,a\right]\right]\right]\right)$

IntegrationByParts(omega5);

$\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[2\,2\right]\right]\,\left[\left[\left[1\,2\,3\,7\right]\,-\frac{a_x}{2}\right]\,\left[\left[1\,2\,3\,12\right]\,-\frac{a}{2}\right]\,\left[\left[1\,2\,4\,5\right]\,\frac{a_x}{2}\right]\,\left[\left[1\,2\,4\,9\right]\,\frac{a}{2}\right]\right]\right]\right)$

eta := evalDG(u[1]*Dx &w Cu[2] &w Cv[1] &w Cu[1, 1]);

$\mathrm{\eta }≔\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[1\,3\right]\right]\,\left[\left[\left[1\,6\,7\,9\right]\,u_1\right]\right]\right]\right)$

omega6 := HorizontalExteriorDerivative(eta);

$\mathrm{\ω6}≔\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[2\,3\right]\right]\,\left[\left[\left[1\,2\,6\,7\,9\right]\,-u_{1,2}\right]\,\left[\left[1\,2\,6\,7\,16\right]\,-u_1\right]\,\left[\left[1\,2\,6\,9\,13\right]\,u_1\right]\,\left[\left[1\,2\,7\,9\,11\right]\,-u_1\right]\right]\right]\right)$

IntegrationByParts(omega6);

$\mathrm{\_DG}\left(\left[\left[''biform''\,E\,\left[2\,3\right]\right]\,\left[\left[\left[1\,2\,3\,4\,5\right]\,0\right]\right]\right]\right)$