private void SetScaledRate() { lambda = Blas.Amax(ndim, rate) * factor; Blas.Dcopy(ndim, rate, scaledRate); Blas.Dscal(ndim, 1.0 / lambda, scaledRate); }