                                 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;

__device__
void elina_double_interval_mul(float_type* const a_inf, float_type* const a_sup, const float_type b_inf, const float_type b_sup, const float_type c_inf, const float_type c_sup)
{
    if(c_inf <= 0)
    {
        /* interval c is positive */
        if(b_inf <= 0)
        {
            /*interval b is positive*/
            *a_inf = b_inf*-c_inf;
            *a_sup = b_sup*c_sup;
        }
        else if(b_sup <= 0)
        {
            /* interval b is negative */
            *a_inf = c_sup*b_inf;
            *a_sup = -c_inf*b_sup;
        }
        else
        {

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