anonymous_namespace{selfContact.cpp} Namespace Reference

Functions
util::Point	cappedRepulsive (const util::Point &yji, double volj, double Kn, double natural_R)

Function Documentation

cappedRepulsive()

util::Point anonymous_namespace{selfContact.cpp}::cappedRepulsive	(	const util::Point &	yji,
		double	volj,
		double	Kn,
		double	natural_R
	)

Definition at line 17 of file selfContact.cpp.

                                              {
  const double Rji = yji.length();
  if (!(Rji > 0.) || !(natural_R > 0.))
    return util::Point();
 
  double gap = Rji - natural_R;
  if (gap > 0.)
    gap = 0.;
  else {
    // Cap penetration so a bond that breaks while already deep inside
    // natural_R cannot inject a discontinuous Kn*(R-natural_R) kick.
    const double gap_cap = -0.25 * natural_R;
    if (gap < gap_cap)
      gap = gap_cap;
  }
 
  const double scalar_f = Kn * volj * gap / Rji;
  return scalar_f * yji;
}

References util::Point::length().

Here is the call graph for this function: