Fix #317 (segfault when accessing S matrix property of SR)

Created by: femtobit

This PR changes SR::LastSMatrix to return a reference to optional. This seems to work with pybind and fixes #317 (closed).