This routine is part of the liesymm package and is loaded via with(liesymm) .

The routine wsubs() is analogous to powsubs() but is for wedge products.  One or more equations specify replacements that are to be made.

The replacements occur even if only a subset of the $\&\^$ arguments match.

If two or more substitutions are specified, they are completed in the order specified, even if they are given as a set. Simultaneous substitution (as in subs()) has not been implemented.

$\mathrm{with}\left(\mathrm{liesymm}\right)\:$$\mathrm{setup}\left(x\,y\,z\,t\right)\:$

$\mathrm{wsubs}\left(d\left(x\right)=d\left(z\right)\,ad\left(x\right)\&ˆd\left(y\right)+bd\left(x\right)\&ˆd\left(y\right)\right)$

$-ad\left(y\right)\&ˆd\left(z\right)-bd\left(y\right)\&ˆd\left(z\right)$

$\mathrm{getform}\left(ad\left(x\right)\&ˆd\left(y\right)\right)$

$d\left(x\right)\&ˆd\left(y\right)$

$\mathrm{wsubs}\left(=d\left(x\right)\&ˆd\left(t\right)\,\left(d\left(x\right)\&ˆd\left(y\right)\right)\&ˆd\left(z\right)\right)$

$-\mathrm{\&^}\left(d\left(x\right)\,d\left(z\right)\,d\left(t\right)\right)$