mat_cumsum_columns_with_scalar_ceiling function

Return matrix with column-wise cumulative sum replacing cells larger than ceil with NA. No checks for arguments is done.