problem: problem_logexp_sum_of_splines, type = maximize
logexp_sum(spline_sum(matrix_parameters_vars, matrix_data))