We introduce and study twelve multivariable theta functions defined by
pfaffians with elliptic function entries. We show that, when the crossing
parameter is a cubic root of unity, the domain wall partition function for the
eight-vertex-solid-on-solid model can be written as a sum of two of these
pfaffians. As a limit case, we express the domain wall partition function for
the three-colour model as a sum of two Hankel determinants. We also show that
certain solutions of the TQ-equation for the supersymmetric eight-vertex model
can be expressed in terms of elliptic pfaffians.Comment: 34 page