 void elina_double_interval_mul_cst_coeff(float_type* const res_inf, float_type* const res_sup, const float_type inf, const float_type sup, const float_type inf_expr, const float_type sup_expr)
 {
     elina_double_interval_mul_expr_coeff(res_inf, res_sup, inf, sup, inf_expr, sup_expr);
-    *res_inf += min_denormal;
+    *res_inf -= min_denormal;
     *res_sup += min_denormal;
 }
 

    elina_double_interval_mul(res_inf, res_sup, inf, sup, inf_expr, sup_expr);
    const float_type maxA = max(fabs(inf_expr), fabs(sup_expr));
    float_type tmp1, tmp2;
    elina_double_interval_mul(&tmp1, &tmp2, inf, sup, maxA*ulp, maxA*ulp);
    *res_inf += tmp1;
    *res_sup += tmp2;
}

                                maxRes = max(fabs(inf_coeff), fabs(sup_coeff));
                                maxMul = max(fabs(tmp1), fabs(tmp2));

                                inf_coeff = inf_coeff + tmp1 - (maxRes + maxMul)*ulp;
                                sup_coeff = sup_coeff + tmp2 + (maxRes + maxMul)*ulp;
                            }
                        }
