1. A method of optimizing a Business Process Management (BPM) model, said model being associated with a plurality of key performance indicators (KPI) depending on a control parameter vector comprising a plurality of control parameters, the method comprising: separating, by a processor, said control parameters into a linear control parameter vector (X) and a non-linear control parameter vector (X′), iteratively calculating, by the processor, a set of candidate values for said control parameter vector, each candidate value being determined for a preselected value of a non-linear control parameter belonging to said non-linear control parameter vector (X′) from execution of a mixed integer program (MIP) model associated with said BPM model, said plurality of key performance indicators (KPI) and said preselected value of the non-linear control parameter vector (X′), and adjusting, by the processor, said BPM model based on at least one of said candidate values for said control parameter vector.