                 const size_t mat_in = x_val*input_size_y*input_size_z + y_val*input_size_z + inp_z;
 
                 const size_t filter_index = out_z*filter_size_x*filter_size_y*input_size_z + x_shift*filter_size_y*input_size_z + y_shift*input_size_z + inp_z;
+
                 elina_double_interval_mul(&tmp1, &tmp2, coeffs[filter_index], coeffs[filter_index], input_inf[mat_in], input_sup[mat_in]);
+
                 res_inf = res_inf + tmp1;
                 res_sup = res_sup + tmp2;
             }

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


                    inf_cst[n] += tmp1;
                    sup_cst[n] += tmp1;
                }
            }
        }
