====== Constrained summations ======

Bert to write this up.

<code>
A[] = whatever()
</code>