Fixed a bug with neck strengths

egor-demidov · Mar 13, 2024 · 9ea9339 · 9ea9339
1 parent 6d1d9a0
commit 9ea9339
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/include/libgran/sinter_bridge/alt_breaking_sinter_bridge.h b/include/libgran/sinter_bridge/alt_breaking_sinter_bridge.h
@@ -115,7 +115,7 @@ struct alt_breaking_sinter_functor {
         field_value_t f_tot = f_n * n + f_t;
         field_value_t tau_tot = -tau_t + tau_r + tau_o;
 
-        if (f_tot.norm() > bonded_contacts[i * n_part + j])
+        if (f_tot.norm() > neck_strengths[i * n_part + j])
             bonded_contacts[i * n_part + j] = false;
 
         return std::make_pair(f_tot / mass, tau_tot / inertia);