 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_sup += min_denormal;
 }

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

                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_sup = res_sup + tmp2;
            }
        }
