                     float_type tmp1, tmp2;
                     elina_double_interval_mul_cst_coeff(&tmp1, &tmp2, mu_inf, mu_sup, old_inf_coeff, old_sup_coeff);
 
-                    inf_cst[n] += tmp1 - min_denormal;
-                    sup_cst[n] += tmp2 + min_denormal;
+                    res_inf_cst += tmp1 - min_denormal;
+                    res_sup_cst += tmp2 + min_denormal;
                 }
                 else if(old_sup_coeff < 0)
                 {

        }
        else
        {
            compute_lb_from_expr_conv_sparse<<<dim3(fp->layers[layerno]->output_size[0], fp->layers[layerno]->output_size[1], fp->layers[layerno]->output_size[2]/num_chunks), 1>>>(lb_array, linf_coeff, lsup_coeff, linf_cst, fp->input_inf, fp->input_sup, num_chunks, chunk_counter, fp->layers[0]->input_size[0], fp->layers[0]->input_size[1], fp->layers[0]->input_size[2], offset_x, offset_y, length_x, length_y, shift_x, shift_y);
            compute_ub_from_expr_conv_sparse<<<dim3(fp->layers[layerno]->output_size[0], fp->layers[layerno]->output_size[1], fp->layers[layerno]->output_size[2]/num_chunks), 1>>>(ub_array, uinf_coeff, usup_coeff, usup_cst, fp->input_inf, fp->input_sup, num_chunks, chunk_counter, fp->layers[0]->input_size[0], fp->layers[0]->input_size[1], fp->layers[0]->input_size[2], offset_x, offset_y, length_x, length_y, shift_x, shift_y);
        }
    }


        lcsts_from_input_poly_sparse<<<dim3(fp->layers[0]->output_size[0], fp->layers[0]->output_size[1], fp->layers[0]->output_size[2]), 1>>>(coeffs, coeffs, csts, csts, linf_cst, lsup_cst, fp->input_lcst, fp->input_ucst, fp->input_inf, fp->input_sup, fp->layers[0]->input_size[0], fp->layers[0]->input_size[1], fp->layers[0]->input_size[2], offset_x, offset_y, length_x, length_y, shift_x, shift_y);
        ucsts_from_input_poly_sparse<<<dim3(fp->layers[0]->output_size[0], fp->layers[0]->output_size[1], fp->layers[0]->output_size[2]), 1>>>(coeffs, coeffs, csts, csts, uinf_cst, usup_cst, fp->input_lcst, fp->input_ucst, fp->input_inf, fp->input_sup, fp->layers[0]->input_size[0], fp->layers[0]->input_size[1], fp->layers[0]->input_size[2], offset_x, offset_y, length_x, length_y, shift_x, shift_y);

        compute_lb_from_expr<<<num_out_neurons_0_layer, num_threads>>>(fp->layers[0]->lb_array, linf_coeff, lsup_coeff, linf_cst, fp->input_inf + num_in_neurons_0_layer, fp->input_sup + num_in_neurons_0_layer, mu);
        compute_ub_from_expr<<<num_out_neurons_0_layer, num_threads>>>(fp->layers[0]->ub_array, uinf_coeff, usup_coeff, usup_cst, fp->input_inf + num_in_neurons_0_layer, fp->input_sup + num_in_neurons_0_layer, mu);

        cudaFree(linf_coeff);
        cudaFree(lsup_coeff);
