Correct number of integration points. Dynamic order.

This corrects a bug in calculation of integration points number.
Other changes allow a user to specify the integration order at runtime at
initialization of the local assemblers.