add compute_stl.R