compute_R_inv function

Compute SPDE and prior precision matrices for S