S_star_fc is already several times faster than the R version, so the problem must be with the sampling code