The support for multidomain (via a separate Implementation template arg) has been only implemented for the combination implicit=true and DiffMethod::numeric. This adapts the pattern for the other specializations.