From 70afc5a7d871919776a64782e8b93404e6b0defd Mon Sep 17 00:00:00 2001 From: eug-vs Date: Sun, 15 Dec 2024 13:17:43 +0100 Subject: feat!: add raylib rendering --- src/constraint/beam.rs | 46 ---------------------------------------------- 1 file changed, 46 deletions(-) delete mode 100644 src/constraint/beam.rs (limited to 'src/constraint/beam.rs') diff --git a/src/constraint/beam.rs b/src/constraint/beam.rs deleted file mode 100644 index c953920..0000000 --- a/src/constraint/beam.rs +++ /dev/null @@ -1,46 +0,0 @@ -use nalgebra::{DVector, RowDVector}; - -use crate::particle_system::{ParticleSystem, Scalar, N}; - -use super::Constraint; - -pub struct BeamConstraint { - pub particle_ids: [usize; 2], - pub length: Scalar, - - jacobian: RowDVector, -} - -impl ParticleSystem { - pub fn add_beam_constraint(&mut self, particle_ids: [usize; 2]) { - let a = &self.particles[particle_ids[0]]; - let b = &self.particles[particle_ids[1]]; - - self.constraints.push(Box::new(BeamConstraint { - particle_ids, - length: (a.position - b.position).norm(), - jacobian: RowDVector::zeros(self.particles.len() * N), - })); - } -} - -impl Constraint for BeamConstraint { - fn get_particles(&self) -> Vec { - Vec::from(self.particle_ids) - } - - fn c(&self, q: &DVector) -> Scalar { - let a = q.fixed_rows::(self.particle_ids[0] * N); - let b = q.fixed_rows::(self.particle_ids[1] * N); - - (a - b).norm() - self.length - } - - fn set_jacobian(&mut self, jacobian: RowDVector) { - self.jacobian = jacobian - } - - fn jacobian_prev(&self) -> RowDVector { - self.jacobian.clone() - } -} -- cgit v1.2.3