mat_cumsum_columns_with_vector_ceiling function

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