             {
                 elina_double_interval_mul_expr_coeff(&tmp1, &tmp2, expr_inf_coeff[a], expr_sup_coeff[a], aux_inf_coeff[c], aux_sup_coeff[c]);
 
-                maxRes = fmax(fabs(res_inf_coeff[b]), fabs(res_sup_coeff[b]));
+                maxRes = fmax(fabs(inf_coeff), fabs(sup_coeff));
                 maxMul = fmax(fabs(tmp1), fabs(tmp2));
 
-                res_inf_coeff[b] = res_inf_coeff[b] + tmp1 + (maxRes + maxMul)*ulp;
-                res_sup_coeff[b] = res_sup_coeff[b] + tmp2 + (maxRes + maxMul)*ulp;
-            }
-        }
+                inf_coeff = inf_coeff + tmp1 + (maxRes + maxMul)*ulp;
+                sup_coeff = sup_coeff + tmp2 + (maxRes + maxMul)*ulp;
+            }
+        }
+
+        res_inf_coeff[b] = inf_coeff;
+        res_sup_coeff[b] = sup_coeff;
     }
 }
 

        const size_t b = n*num_in_neurons_current_layer + j;
        size_t c = i*num_in_neurons_current_layer + j;

        elina_double_interval_mul_expr_coeff(&res_inf_coeff[b], &res_sup_coeff[b], expr_inf_coeff[a], expr_sup_coeff[a], aux_inf_coeff[c], aux_sup_coeff[c]);

        float_type tmp1, tmp2;
        float_type maxRes, maxMul;


    size_t a = n*num_out_neurons_current_layer + i;

    float_type inf_cst;
    float_type sup_cst;

    elina_double_interval_mul_cst_coeff(&inf_cst, &sup_cst, expr_inf_coeff[a], expr_sup_coeff[a], aux_inf_cst[i], aux_sup_cst[i]);

    float_type tmp1, tmp2;
    float_type maxRes, maxMul;
