edits
make opsahl faster; fix bug in np.gaussian