 
                     inf_cst[n] += tmp1;
                     sup_cst[n] += tmp1;
-
-                    //atomicAdd(&inf_cst[n], tmp1);
-                    //atomicAdd(&sup_cst[n], -tmp1);
                 }
             }
         }

    for(size_t i = 0; i < num_in_neurons; i++)
    {
        elina_double_interval_mul(&tmp1, &tmp2, coeffs[n*num_in_neurons + i], coeffs[n*num_in_neurons + i], input_inf[i], input_sup[i]);
        res_inf = res_inf + tmp1;
        res_sup = res_sup + tmp2;
    }

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