problem: problem_logexp_of_spline, type = maximize
logexp_sum(spline_sum(matrix_vars_f1, matrix_data_f1))