                                         const float_type aux_coeff = aux_coeffs[filter_index];
                                         elina_double_interval_mul_expr_coeff(&tmp1, &tmp2, prev_inf_coeff, prev_sup_coeff, -aux_coeff, aux_coeff);
 
-                                        maxRes = fmax(fabs(inf_coeff), fabs(sup_coeff));
-                                        maxMul = fmax(fabs(tmp1), fabs(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;


                    elina_double_interval_mul_cst_coeff(&tmp1, &tmp2, prev_inf_coeff, prev_sup_coeff, -aux_cst, aux_cst);

                    maxRes = fmax(fabs(inf_cst), fabs(sup_cst));
                    maxMul = fmax(fabs(tmp1), fabs(tmp2));

                    inf_cst += tmp1 + (maxRes + maxMul)*ulp + min_denormal;
                    sup_cst += tmp2 + (maxRes + maxMul)*ulp + min_denormal;

void elina_double_interval_mul_expr_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(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;
