%This Example for using users m-function Optimization_subroutine.m was created automatically by PSG Toolbox.

%Function description:

%maximize
%logexp_sum_matrix_spline(spline_sum_test(matrix_parameters_vars, matrix_data))
%constraint:  <= 0
%polynom_abs(matrix_spline_polin_abs_1)
%- variable_u1(u1)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_2)
%- variable_u2(u2)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_3)
%- variable_u3(u3)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_4)
%- variable_u4(u4)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_5)
%- variable_u5(u5)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_6)
%- variable_u6(u6)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_7)
%- variable_u7(u7)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_8)
%- variable_u8(u8)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_9)
%- variable_u9(u9)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_10)
%- variable_u10(u10)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_11)
%- variable_u11(u11)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_12)
%- variable_u12(u12)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_13)
%- variable_u13(u13)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_14)
%- variable_u14(u14)
%constraint: <= 0
%polynom_abs(matrix_spline_polin_abs_15)
%- variable_u15(u15)
%constraint: <= 5, 
%linear(matrix_add_variable)
%Box: types = point_type
%
%calculate
%Point: point_problem_1
%logexp_sum(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%logistic(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%
%Input variables:
%
%Inputs                             PSG Type    PSG Object                    Location in Problem Statement                                                     Class
%matrix_data_benchmark              bench       matrix_data                   logexp_sum_matrix_spline(spline_sum_test(matrix_parameters_vars, matrix_data))    double
%                                                                             logexp_sum(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%                                                                             logistic(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%matrix_data_data                   data        matrix_data                   logexp_sum_matrix_spline(spline_sum_test(matrix_parameters_vars, matrix_data))    double
%                                                                             logexp_sum(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%                                                                             logistic(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%matrix_data_vars                   vars        matrix_data                   logexp_sum_matrix_spline(spline_sum_test(matrix_parameters_vars, matrix_data))    cell
%                                                                             logexp_sum(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%                                                                             logistic(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%matrix_add_variable_data           data        matrix_add_variable           linear(matrix_add_variable)                                                       double
%matrix_add_variable_vars           vars        matrix_add_variable           linear(matrix_add_variable)                                                       cell
%matrix_parameters_vars_data        data        matrix_parameters_vars        logexp_sum_matrix_spline(spline_sum_test(matrix_parameters_vars, matrix_data))    double
%                                                                             logexp_sum(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%                                                                             logistic(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%matrix_parameters_vars_vars        vars        matrix_parameters_vars        logexp_sum_matrix_spline(spline_sum_test(matrix_parameters_vars, matrix_data))    cell
%                                                                             logexp_sum(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%                                                                             logistic(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
%matrix_spline_polin_abs_1_data     data        matrix_spline_polin_abs_1     polynom_abs(matrix_spline_polin_abs_1)                                            double
%matrix_spline_polin_abs_1_vars     vars        matrix_spline_polin_abs_1     polynom_abs(matrix_spline_polin_abs_1)                                            cell
%matrix_spline_polin_abs_2_data     data        matrix_spline_polin_abs_2     polynom_abs(matrix_spline_polin_abs_2)                                            double
%matrix_spline_polin_abs_2_vars     vars        matrix_spline_polin_abs_2     polynom_abs(matrix_spline_polin_abs_2)                                            cell
%matrix_spline_polin_abs_3_data     data        matrix_spline_polin_abs_3     polynom_abs(matrix_spline_polin_abs_3)                                            double
%matrix_spline_polin_abs_3_vars     vars        matrix_spline_polin_abs_3     polynom_abs(matrix_spline_polin_abs_3)                                            cell
%matrix_spline_polin_abs_4_data     data        matrix_spline_polin_abs_4     polynom_abs(matrix_spline_polin_abs_4)                                            double
%matrix_spline_polin_abs_4_vars     vars        matrix_spline_polin_abs_4     polynom_abs(matrix_spline_polin_abs_4)                                            cell
%matrix_spline_polin_abs_5_data     data        matrix_spline_polin_abs_5     polynom_abs(matrix_spline_polin_abs_5)                                            double
%matrix_spline_polin_abs_5_vars     vars        matrix_spline_polin_abs_5     polynom_abs(matrix_spline_polin_abs_5)                                            cell
%matrix_spline_polin_abs_6_data     data        matrix_spline_polin_abs_6     polynom_abs(matrix_spline_polin_abs_6)                                            double
%matrix_spline_polin_abs_6_vars     vars        matrix_spline_polin_abs_6     polynom_abs(matrix_spline_polin_abs_6)                                            cell
%matrix_spline_polin_abs_7_data     data        matrix_spline_polin_abs_7     polynom_abs(matrix_spline_polin_abs_7)                                            double
%matrix_spline_polin_abs_7_vars     vars        matrix_spline_polin_abs_7     polynom_abs(matrix_spline_polin_abs_7)                                            cell
%matrix_spline_polin_abs_8_data     data        matrix_spline_polin_abs_8     polynom_abs(matrix_spline_polin_abs_8)                                            double
%matrix_spline_polin_abs_8_vars     vars        matrix_spline_polin_abs_8     polynom_abs(matrix_spline_polin_abs_8)                                            cell
%matrix_spline_polin_abs_9_data     data        matrix_spline_polin_abs_9     polynom_abs(matrix_spline_polin_abs_9)                                            double
%matrix_spline_polin_abs_9_vars     vars        matrix_spline_polin_abs_9     polynom_abs(matrix_spline_polin_abs_9)                                            cell
%matrix_spline_polin_abs_10_data    data        matrix_spline_polin_abs_10    polynom_abs(matrix_spline_polin_abs_10)                                           double
%matrix_spline_polin_abs_10_vars    vars        matrix_spline_polin_abs_10    polynom_abs(matrix_spline_polin_abs_10)                                           cell
%matrix_spline_polin_abs_11_data    data        matrix_spline_polin_abs_11    polynom_abs(matrix_spline_polin_abs_11)                                           double
%matrix_spline_polin_abs_11_vars    vars        matrix_spline_polin_abs_11    polynom_abs(matrix_spline_polin_abs_11)                                           cell
%matrix_spline_polin_abs_12_data    data        matrix_spline_polin_abs_12    polynom_abs(matrix_spline_polin_abs_12)                                           double
%matrix_spline_polin_abs_12_vars    vars        matrix_spline_polin_abs_12    polynom_abs(matrix_spline_polin_abs_12)                                           cell
%matrix_spline_polin_abs_13_data    data        matrix_spline_polin_abs_13    polynom_abs(matrix_spline_polin_abs_13)                                           double
%matrix_spline_polin_abs_13_vars    vars        matrix_spline_polin_abs_13    polynom_abs(matrix_spline_polin_abs_13)                                           cell
%matrix_spline_polin_abs_14_data    data        matrix_spline_polin_abs_14    polynom_abs(matrix_spline_polin_abs_14)                                           double
%matrix_spline_polin_abs_14_vars    vars        matrix_spline_polin_abs_14    polynom_abs(matrix_spline_polin_abs_14)                                           cell
%matrix_spline_polin_abs_15_data    data        matrix_spline_polin_abs_15    polynom_abs(matrix_spline_polin_abs_15)                                           double
%matrix_spline_polin_abs_15_vars    vars        matrix_spline_polin_abs_15    polynom_abs(matrix_spline_polin_abs_15)                                           cell
%point_type_data                    data        point_type                    Box: types = point_type                                                           double
%point_type_vars                    vars        point_type                    Box: types = point_type                                                           cell
%
%Output variables:
%
%solution_str = string with solution of problem;
%outargstruc_arr = array of output PSG data structures;


%Load data from mat-file:
load('D:\American Optimal Decisions\PSG\MATLAB_Stan\All\Spline Regression\problem_3_logexp_sum_of_splines_boolean_calculate\Optimization_subroutine_data.mat')

%Save variables from mat-file to Workspace:
tbpsg_export_to_workspace(toolboxstruc_arr)

%Run users m-function Optimization_subroutine:

[solution_str,outargstruc_arr] = Optimization_subroutine(matrix_data_benchmark,matrix_data_data,matrix_data_vars,matrix_add_variable_data,matrix_add_variable_vars,matrix_parameters_vars_data,matrix_parameters_vars_vars,matrix_spline_polin_abs_1_data,matrix_spline_polin_abs_1_vars,matrix_spline_polin_abs_2_data,matrix_spline_polin_abs_2_vars,matrix_spline_polin_abs_3_data,matrix_spline_polin_abs_3_vars,matrix_spline_polin_abs_4_data,matrix_spline_polin_abs_4_vars,matrix_spline_polin_abs_5_data,matrix_spline_polin_abs_5_vars,matrix_spline_polin_abs_6_data,matrix_spline_polin_abs_6_vars,matrix_spline_polin_abs_7_data,matrix_spline_polin_abs_7_vars,matrix_spline_polin_abs_8_data,matrix_spline_polin_abs_8_vars,matrix_spline_polin_abs_9_data,matrix_spline_polin_abs_9_vars,matrix_spline_polin_abs_10_data,matrix_spline_polin_abs_10_vars,matrix_spline_polin_abs_11_data,matrix_spline_polin_abs_11_vars,matrix_spline_polin_abs_12_data,matrix_spline_polin_abs_12_vars,matrix_spline_polin_abs_13_data,matrix_spline_polin_abs_13_vars,matrix_spline_polin_abs_14_data,matrix_spline_polin_abs_14_vars,matrix_spline_polin_abs_15_data,matrix_spline_polin_abs_15_vars,point_type_data,point_type_vars);

%Extract Objective:
val_obj = tbpsg_objective(solution_str, outargstruc_arr);

disp(' ');
disp('Objective = ');
disp(val_obj);

%Extract optimal solution:
point_data = tbpsg_optimal_point_data(solution_str, outargstruc_arr);

disp(' ');
disp('Optimal point = ');
disp(point_data);

%Extract structure containing PSG solution reports:
output_structure = tbpsg_solution_struct(solution_str, outargstruc_arr);

disp(' ');
disp('Structure with PSG solution = ');
disp(output_structure);

%Uncomment the following lines to extract solutions details:

%output = tbpsg_isoptimal(solution_str, outargstruc_arr);

%output = tbpsg_function_data(solution_str, outargstruc_arr);

%output = tbpsg_function_names(solution_str, outargstruc_arr);

%output = tbpsg_time(solution_str, outargstruc_arr);

%output = tbpsg_optimal_point_vars(solution_str, outargstruc_arr);

%output = tbpsg_constraints_vars(solution_str, outargstruc_arr);

%output = tbpsg_slack_data(solution_str, outargstruc_arr);

%output = tbpsg_dual_data(solution_str, outargstruc_arr);

%output = tbpsg_vector_constraint_data(solution_str, outargstruc_arr);

%output = tbpsg_vector_dual_data(solution_str, outargstruc_arr);

%output = tbpsg_vector_slack_data(solution_str, outargstruc_arr);

%output = tbpsg_matrix_data(solution_str, outargstruc_arr);

%output = tbpsg_matrix_vars(solution_str, outargstruc_arr);

%output = tbpsg_vector_data(solution_str, outargstruc_arr);