Non-array magnitude are casted to numpy arrays when a numpy function is applied.

fixes #44
Support for mathematical operations on dimensionless (but not simplified quantities)

fixes #45
Updated CHANGES and AUTHORS