Unrolling the dimension loops broke the diagonal_jacobian
transformation since it was looking for index_sums over certain
expressions.

Example: Before this commit the expression

\sum_i a_i (\nabla v)_i + \sum_j b_j (\nabla v)_j

was split into two accumulation terms:

1) a_k with corresponding test function (\nabla v)_k
2) b_l with corresponding test function (\nabla v)_l

Now we split into:

a_k + b_k with corresponding test function (\nabla v)_k

This is possible since we have linearity in the test function.

TODO: Jacobians are not yet working

Very wild, totally WIP.

Various fixes, including diagonal jacobian transformations.

Some work on parameter classes, a bugfix on the splitting algorithm
some simplifications in the generation stage

UFL did not like  my way of replacing some args with 1 as it
destroyed the shape information. It now remains in the expression,
but triggers code generation separately.

and remove old implementation