 
     for(size_t i = 0; i < expr_size; i++)
     {
-        elina_double_interval_mul(&tmp1, &tmp2, inf_coeff[n*expr_size + i], sup_coeff[n*expr_size + i], input_inf[i], input_sup[i]);
+        elina_double_interval_mul2(&tmp1, &tmp2, inf_coeff[n*expr_size + i], sup_coeff[n*expr_size + i], input_inf[i], input_sup[i]);
 
         res_sup = res_sup + tmp2;
     }

                    inf_coeff[a] = 0.0;
                    sup_coeff[a] = 0.0;
                    float_type tmp1, tmp2;
                    elina_double_interval_mul(&tmp1, &tmp2, old_inf_coeff, old_sup_coeff, 0, ub);

                    inf_cst[n] += tmp2;
                    sup_cst[n] += tmp2;

                    inf_coeff[a] = 0.0;
                    sup_coeff[a] = 0.0;
                    float_type tmp1, tmp2;
                    elina_double_interval_mul2(&tmp1, &tmp2, old_inf_coeff, old_sup_coeff, 0, ub);

                    inf_cst[n] += tmp1;
                    sup_cst[n] += tmp1;
