Optimize forward_pad_doys() in mtan module