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


                size_t mat_out = local_n*length_x*length_y*output_size_z + j;

                elina_double_interval_mul(&tmp1, &tmp2, inf_coeff[mat_out], sup_coeff[mat_out], input_inf[i], input_sup[i]);

                res_inf = res_inf + tmp1;
           }


    for(size_t i = 0; i < expr_size; 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_inf = res_inf + tmp1;
    }
