Opt branch

@acroy, thanks for the nice work. I created a pull request, since this facilitates adding line comments to the code. Hence, see below for further comments of further potential comments.
```julia
scale!(copy!(vm[j+1],p),1/hm[j+1,j])
```