Put back matrix pass by reference & in theta_s
Put back matrix pass by reference & in theta_s

eth-cscs/compression,eth-cscs/compression,eth-cscs/compression,eth-cscs/compression,eth-cscs/compression,eth-cscs/compression