                     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);
+                    elina_double_interval_mul2(&tmp1, &tmp2, old_inf_coeff, old_sup_coeff, 0, ub);
 
                     inf_cst[n] += tmp1;
                     sup_cst[n] += tmp1;

            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);

            atomicAdd(&inf_cst[n], tmp2);
            atomicAdd(&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);

            atomicAdd(&inf_cst[n], tmp1);
            atomicAdd(&sup_cst[n], tmp1);
