Reduce allocations during model fitting

This improves performance only minimally
Merge branch 'sjk/matrix'
Merge branch 'sjk/perf'
Merge branch 'sjk/gamma'

Conflicts:
	src/glmfit.jl