Merge pull request #44 from mwcraig/mul-etc-take-2

Override basic arithmetic operators
Change linear gain to use automatic error propagation and better handle units
Add bounds checking to FITS slices

Closes #41