strain limiting tau -> 0 and eps -> 0.01

examples/hang.ipynb

@@ -44,12 +44,8 @@
     "param = app.session.param()\n",
     "gravity = param.get(\"gravity\")\n",
     "\n",
-    "# set strain-limiting upper bound 1% (0.5%+0.5%)\n",
-    "(\n",
-    "    param.set(\"strain-limit-eps\", 0.005)\n",
-    "    .set(\"strain-limit-tau\", 0.005)\n",
-    "    .set(\"friction\", 0.0)\n",
-    ")\n",
+    "# set strain-limiting upper bound 1% (1%+0%)\n",
+    "param.set(\"strain-limit-eps\", 0.01).set(\"strain-limit-tau\", 0.0)\n",
     "\n",
     "# set the maximal video frame to 200\n",
     "param.set(\"frames\", 200)\n",

examples/trampoline.ipynb

@@ -44,8 +44,8 @@
    "source": [
     "param = (\n",
     "    app.session.param()\n",
-    "    .set(\"strain-limit-eps\", 0.005)\n",
-    "    .set(\"strain-limit-tau\", 0.005)\n",
+    "    .set(\"strain-limit-eps\", 0.01)\n",
+    "    .set(\"strain-limit-tau\", 0.0)\n",
     "    .set(\"area-density\", 1e3)\n",
     "    .set(\"volume-density\", 6e3)\n",
     "    .set(\"area-young-mod\", 10000.0)\n",

src/args.rs

@@ -348,7 +348,7 @@ pub struct Args {
     // This bending stiffness not scaled by the density.
     // If you change the density, you need to adjust this value accordingly.
     // This behavior also change by the material thickness.
-    #[clap(long, default_value_t = 10.0)]
+    #[clap(long, default_value_t = 1.0)]
     pub bend: f32,
 
     // Name: Bend Stiffness for Rods

src/cpp/solver/solver.hpp

@@ -10,9 +10,6 @@
 
 namespace solver {
 
-void apply(const DynCSRMat &A, const FixedCSRMat &B, const Vec<Mat3x3f> &C,
-           float D, const Vec<float> &x, Vec<float> &result);
-
 bool solve(const DynCSRMat &A, const FixedCSRMat &B, const Vec<Mat3x3f> &C,
            Vec<float> b, float tol, unsigned max_iter, Vec<float> x,
            unsigned &iter, float &resid);