%This Example for using users m-function Optimization_subroutine.m was created automatically by PSG Toolbox. %Function description: %minimize %linear(matrix_w) %Constraint: <= parameter_bound, linearize = 1 %prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) %Constraint: == parameter_bound_1 %linearmulti(pmatrix_dab) %Constraint: == parameter_bound_2 %linearmulti(pmatrix_dac) %Constraint: == parameter_bound_3 %linearmulti(pmatrix_dad) %Constraint: == parameter_bound_4 %linearmulti(pmatrix_dae) %Constraint: == parameter_bound_5 %linearmulti(pmatrix_dba) %Constraint: == parameter_bound_6 %linearmulti(pmatrix_dbc) %Constraint: == parameter_bound_7 %linearmulti(pmatrix_dbd) %Constraint: == parameter_bound_8 %linearmulti(pmatrix_dbe) %Constraint: == parameter_bound_9 %linearmulti(pmatrix_dca) %Constraint: == parameter_bound_10 %linearmulti(pmatrix_dcb) %Constraint: == parameter_bound_11 %linearmulti(pmatrix_dcd) %Constraint: == parameter_bound_12 %linearmulti(pmatrix_dce) %Constraint: == parameter_bound_13 %linearmulti(pmatrix_dda) %Constraint: == parameter_bound_14 %linearmulti(pmatrix_ddb) %Constraint: == parameter_bound_15 %linearmulti(pmatrix_ddc) %Constraint: == parameter_bound_16 %linearmulti(pmatrix_dde) %Constraint: == parameter_bound_17 %linearmulti(pmatrix_dea) %Constraint: == parameter_bound_18 %linearmulti(pmatrix_deb) %Constraint: == parameter_bound_19 %linearmulti(pmatrix_dec) %Constraint: == parameter_bound_20 %linearmulti(pmatrix_ded) %Box: >= parameter_bound_21 %Solver: parameter_solver % %Input variables: % %Inputs PSG Type PSG Object Location in Problem Statement Class %matrix_w_data data matrix_w linear(matrix_w) double %matrix_w_vars vars matrix_w linear(matrix_w) cell %pmatrix_dab_data data pmatrix_dab linearmulti(pmatrix_dab) double %pmatrix_dab_vars vars pmatrix_dab linearmulti(pmatrix_dab) cell %pmatrix_dac_data data pmatrix_dac linearmulti(pmatrix_dac) double %pmatrix_dac_vars vars pmatrix_dac linearmulti(pmatrix_dac) cell %pmatrix_dad_data data pmatrix_dad linearmulti(pmatrix_dad) double %pmatrix_dad_vars vars pmatrix_dad linearmulti(pmatrix_dad) cell %pmatrix_dae_data data pmatrix_dae linearmulti(pmatrix_dae) double %pmatrix_dae_vars vars pmatrix_dae linearmulti(pmatrix_dae) cell %pmatrix_dba_data data pmatrix_dba linearmulti(pmatrix_dba) double %pmatrix_dba_vars vars pmatrix_dba linearmulti(pmatrix_dba) cell %pmatrix_dbc_data data pmatrix_dbc linearmulti(pmatrix_dbc) double %pmatrix_dbc_vars vars pmatrix_dbc linearmulti(pmatrix_dbc) cell %pmatrix_dbd_data data pmatrix_dbd linearmulti(pmatrix_dbd) double %pmatrix_dbd_vars vars pmatrix_dbd linearmulti(pmatrix_dbd) cell %pmatrix_dbe_data data pmatrix_dbe linearmulti(pmatrix_dbe) double %pmatrix_dbe_vars vars pmatrix_dbe linearmulti(pmatrix_dbe) cell %pmatrix_dca_data data pmatrix_dca linearmulti(pmatrix_dca) double %pmatrix_dca_vars vars pmatrix_dca linearmulti(pmatrix_dca) cell %pmatrix_dcb_data data pmatrix_dcb linearmulti(pmatrix_dcb) double %pmatrix_dcb_vars vars pmatrix_dcb linearmulti(pmatrix_dcb) cell %pmatrix_dcd_data data pmatrix_dcd linearmulti(pmatrix_dcd) double %pmatrix_dcd_vars vars pmatrix_dcd linearmulti(pmatrix_dcd) cell %pmatrix_dce_data data pmatrix_dce linearmulti(pmatrix_dce) double %pmatrix_dce_vars vars pmatrix_dce linearmulti(pmatrix_dce) cell %pmatrix_dda_data data pmatrix_dda linearmulti(pmatrix_dda) double %pmatrix_dda_vars vars pmatrix_dda linearmulti(pmatrix_dda) cell %pmatrix_ddb_data data pmatrix_ddb linearmulti(pmatrix_ddb) double %pmatrix_ddb_vars vars pmatrix_ddb linearmulti(pmatrix_ddb) cell %pmatrix_ddc_data data pmatrix_ddc linearmulti(pmatrix_ddc) double %pmatrix_ddc_vars vars pmatrix_ddc linearmulti(pmatrix_ddc) cell %pmatrix_dde_data data pmatrix_dde linearmulti(pmatrix_dde) double %pmatrix_dde_vars vars pmatrix_dde linearmulti(pmatrix_dde) cell %pmatrix_dea_data data pmatrix_dea linearmulti(pmatrix_dea) double %pmatrix_dea_vars vars pmatrix_dea linearmulti(pmatrix_dea) cell %pmatrix_deb_data data pmatrix_deb linearmulti(pmatrix_deb) double %pmatrix_deb_vars vars pmatrix_deb linearmulti(pmatrix_deb) cell %pmatrix_dec_data data pmatrix_dec linearmulti(pmatrix_dec) double %pmatrix_dec_vars vars pmatrix_dec linearmulti(pmatrix_dec) cell %pmatrix_ded_data data pmatrix_ded linearmulti(pmatrix_ded) double %pmatrix_ded_vars vars pmatrix_ded linearmulti(pmatrix_ded) cell %pmatrix_xab_data data pmatrix_xab prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xab_vars vars pmatrix_xab prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xac_data data pmatrix_xac prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xac_vars vars pmatrix_xac prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xba_data data pmatrix_xba prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xba_vars vars pmatrix_xba prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xbc_data data pmatrix_xbc prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xbc_vars vars pmatrix_xbc prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xbd_data data pmatrix_xbd prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xbd_vars vars pmatrix_xbd prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xbe_data data pmatrix_xbe prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xbe_vars vars pmatrix_xbe prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xca_data data pmatrix_xca prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xca_vars vars pmatrix_xca prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xcb_data data pmatrix_xcb prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xcb_vars vars pmatrix_xcb prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xcd_data data pmatrix_xcd prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xcd_vars vars pmatrix_xcd prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xdb_data data pmatrix_xdb prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xdb_vars vars pmatrix_xdb prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xdc_data data pmatrix_xdc prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xdc_vars vars pmatrix_xdc prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xde_data data pmatrix_xde prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xde_vars vars pmatrix_xde prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xeb_data data pmatrix_xeb prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xeb_vars vars pmatrix_xeb prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %pmatrix_xed_data data pmatrix_xed prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %pmatrix_xed_vars vars pmatrix_xed prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) cell %parameter_bound_data data parameter_bound Constraint: <= parameter_bound, linearize = 1 double %parameter_alpha_data data parameter_alpha prmulti_pen(parameter_alpha,pmatrix_xab,pmatrix_xac,pmatrix_xba,pmatrix_xbc,pmatrix_xbd,pmatrix_xbe,pmatrix_xca,pmatrix_xcb,pmatrix_xcd,pmatrix_xdb,pmatrix_xdc,pmatrix_xde,pmatrix_xeb,pmatrix_xed) double %parameter_bound_1_data data parameter_bound_1 Constraint: == parameter_bound_1 double %parameter_bound_2_data data parameter_bound_2 Constraint: == parameter_bound_2 double %parameter_bound_3_data data parameter_bound_3 Constraint: == parameter_bound_3 double %parameter_bound_4_data data parameter_bound_4 Constraint: == parameter_bound_4 double %parameter_bound_5_data data parameter_bound_5 Constraint: == parameter_bound_5 double %parameter_bound_6_data data parameter_bound_6 Constraint: == parameter_bound_6 double %parameter_bound_7_data data parameter_bound_7 Constraint: == parameter_bound_7 double %parameter_bound_8_data data parameter_bound_8 Constraint: == parameter_bound_8 double %parameter_bound_9_data data parameter_bound_9 Constraint: == parameter_bound_9 double %parameter_bound_10_data data parameter_bound_10 Constraint: == parameter_bound_10 double %parameter_bound_11_data data parameter_bound_11 Constraint: == parameter_bound_11 double %parameter_bound_12_data data parameter_bound_12 Constraint: == parameter_bound_12 double %parameter_bound_13_data data parameter_bound_13 Constraint: == parameter_bound_13 double %parameter_bound_14_data data parameter_bound_14 Constraint: == parameter_bound_14 double %parameter_bound_15_data data parameter_bound_15 Constraint: == parameter_bound_15 double %parameter_bound_16_data data parameter_bound_16 Constraint: == parameter_bound_16 double %parameter_bound_17_data data parameter_bound_17 Constraint: == parameter_bound_17 double %parameter_bound_18_data data parameter_bound_18 Constraint: == parameter_bound_18 double %parameter_bound_19_data data parameter_bound_19 Constraint: == parameter_bound_19 double %parameter_bound_20_data data parameter_bound_20 Constraint: == parameter_bound_20 double %parameter_bound_21_data data parameter_bound_21 Box: >= parameter_bound_21 double %parameter_solver_data data parameter_solver Solver: parameter_solver char % %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\Stochastic Multicommodity Network Flow Problem\data_problem_p_multicomm_short\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_w_data,matrix_w_vars,pmatrix_dab_data,pmatrix_dab_vars,pmatrix_dac_data,pmatrix_dac_vars,pmatrix_dad_data,pmatrix_dad_vars,pmatrix_dae_data,pmatrix_dae_vars,pmatrix_dba_data,pmatrix_dba_vars,pmatrix_dbc_data,pmatrix_dbc_vars,pmatrix_dbd_data,pmatrix_dbd_vars,pmatrix_dbe_data,pmatrix_dbe_vars,pmatrix_dca_data,pmatrix_dca_vars,pmatrix_dcb_data,pmatrix_dcb_vars,pmatrix_dcd_data,pmatrix_dcd_vars,pmatrix_dce_data,pmatrix_dce_vars,pmatrix_dda_data,pmatrix_dda_vars,pmatrix_ddb_data,pmatrix_ddb_vars,pmatrix_ddc_data,pmatrix_ddc_vars,pmatrix_dde_data,pmatrix_dde_vars,pmatrix_dea_data,pmatrix_dea_vars,pmatrix_deb_data,pmatrix_deb_vars,pmatrix_dec_data,pmatrix_dec_vars,pmatrix_ded_data,pmatrix_ded_vars,pmatrix_xab_data,pmatrix_xab_vars,pmatrix_xac_data,pmatrix_xac_vars,pmatrix_xba_data,pmatrix_xba_vars,pmatrix_xbc_data,pmatrix_xbc_vars,pmatrix_xbd_data,pmatrix_xbd_vars,pmatrix_xbe_data,pmatrix_xbe_vars,pmatrix_xca_data,pmatrix_xca_vars,pmatrix_xcb_data,pmatrix_xcb_vars,pmatrix_xcd_data,pmatrix_xcd_vars,pmatrix_xdb_data,pmatrix_xdb_vars,pmatrix_xdc_data,pmatrix_xdc_vars,pmatrix_xde_data,pmatrix_xde_vars,pmatrix_xeb_data,pmatrix_xeb_vars,pmatrix_xed_data,pmatrix_xed_vars,parameter_bound_data,parameter_alpha_data,parameter_bound_1_data,parameter_bound_2_data,parameter_bound_3_data,parameter_bound_4_data,parameter_bound_5_data,parameter_bound_6_data,parameter_bound_7_data,parameter_bound_8_data,parameter_bound_9_data,parameter_bound_10_data,parameter_bound_11_data,parameter_bound_12_data,parameter_bound_13_data,parameter_bound_14_data,parameter_bound_15_data,parameter_bound_16_data,parameter_bound_17_data,parameter_bound_18_data,parameter_bound_19_data,parameter_bound_20_data,parameter_bound_21_data,parameter_solver_data); %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);