Use a composite basis for the brittle fracture simulation

Conceptually the brittle fracture simulation uses two spaces:
one for the deformation and one for the damage fields.
Previously, these were squeezed into a single power basis
with dim+1 coefficients.  This works, but it sometimes requires
manual index manipulations, which are not always obvious.
This patch introduces a proper composite basis that separates
the deformation from the damage.

The multi-indices for the global degrees of freedom remain unchanged!
2 jobs for use-composite-basis-for-brittle-fracture in 10 minutes and 32 seconds (queued for 2 seconds)