Simple_spacer::rod_force (int l, int r, Real dist)
{
Real d = range_ideal_len (l, r);
- Real c = range_stiffness (l, r, d > dist);
+ Real c = range_stiffness (l, r, dist > d);
Real block_stretch = dist - d;
+
+ if (isinf (c)) /* take care of the 0*infinity_f case */
+ return 0;
return c * block_stretch;
}
for (vsize i=0; i < springs_.size (); i++)
inv_hooke += springs_[i].inverse_stretch_strength ();
+ if (inv_hooke == 0.0) /* avoid division by zero. If springs are infinitely stiff */
+ return 0.0; /* anyway, then it makes no difference what the force is */
+
assert (cur_len <= line_len_);
return (line_len_ - cur_len) / inv_hooke + force_;
}
vector<Spring> sorted_springs = springs_;
sort (sorted_springs.begin (), sorted_springs.end (), greater<Spring> ());
+
for (vsize i = 0; i < sorted_springs.size (); i++)
{
Spring sp = sorted_springs[i];
/****************************************************************/
-/*
- TODO: should a add penalty for widely varying spring forces (caused
- by constraints, eg.
-
-
- . =====
- . | |
- .o|o|x ##x
- .
-
- The ## forces the notes apart; we shouldn't allow the O's to touch
- this closely.
-*/
-
struct Rod_description
{
vsize r_;