bugfix:damping not added to forces
This commit fixes a bug in the numpy implementation of elastica. Damping
forces was not called by the internal forces. Thus for the numpy
implementation damping forces was not added. This commit fixes this
issue, by calling _compute_damping_forces inside the
_compute_internal_forces.
@@ -322,6 +322,7 @@ class CosseratRod(RodBase, _RodSymplecticStepperMixin):
/ self.dilatation # computed in comp_dilatation <- compute_strain <- compute_stress
)
# self._compute_internal_forces()
+ self._compute_damping_forces()
return difference_kernel(cosserat_internal_stress) - self.damping_forces
