6_inv_free

Author: Roushelfy

simulators/6_inv_free/src/NeoHookeanEnergy.cu

@@ -96,7 +96,6 @@ T NeoHookeanEnergy<T, dim>::val()
          NeoHookeanEnergyVal(E,  Mu, Lambda,X,device_IB(i),device_vol(i));
          device_val(i) = E; })
         .wait();
-    display_vec(device_val);
     return devicesum(device_val);
 }
 

simulators/6_inv_free/src/NeoHookean_auto.cu

@@ -1,12 +1,13 @@
 #include "simulator.h"
-
+#include "Eigen/Dense"
+#include <iostream>
 int main()
 {
-	float nu = 0.4, E = 1e5;
-	float Mu = E / (2 * (1 + nu)), Lam = E * nu / ((1 + nu) * (1 - 2 * nu));
-	float rho = 1000,
-		  k = 2e4, initial_stretch = 1, n_seg = 4, h = 0.01, side_len = 1, tol = 0.01, mu = 0.11;
+	double nu = 0.4, E = 1e5;
+	double Mu = E / (2 * (1 + nu)), Lam = E * nu / ((1 + nu) * (1 - 2 * nu));
+	double rho = 1000,
+		   k = 4e4, initial_stretch = 1, n_seg = 10, h = 0.01, side_len = 1, tol = 0.01, mu = 0.11;
 	// printf("Running mass-spring simulator with parameters: rho = %f, k = %f, initial_stretch = %f, n_seg = %d, h = %f, side_len = %f, tol = %f\n", rho, k, initial_stretch, n_seg, h, side_len, tol);
-	InvFreeSimulator<float, 2> simulator(rho, side_len, initial_stretch, k, h, tol, mu, Mu, Lam, n_seg);
+	InvFreeSimulator<double, 2> simulator(rho, side_len, initial_stretch, k, h, tol, mu, Mu, Lam, n_seg);
 	simulator.run();
 }

simulators/6_inv_free/src/main.cpp

@@ -67,7 +67,7 @@ InvFreeSimulator<T, dim>::Impl::Impl(T rho, T side_len, T initial_stretch, T K,
     DBC_limit.push_back(0);
     DBC_limit.push_back(-0.6);
     DBC_v.push_back(0);
-    DBC_v.push_back(-0.5);
+    DBC_v.push_back(-0.3);
     DBC_stiff = 10;
     x.push_back(0);
     x.push_back(side_len * 0.6);
@@ -145,7 +145,7 @@ void InvFreeSimulator<T, dim>::Impl::step_forward()
     update_x_tilde(add_vector<T>(x, v, 1, h));
     frictionenergy.update_mu_lambda(barrierenergy.compute_mu_lambda(mu));
     update_DBC_target();
-    update_DBC_stiff(10);
+    // update_DBC_stiff(10);
     DeviceBuffer<T> x_n = x; // Copy current positions to x_n
     update_v(add_vector<T>(x, x_n, 1 / h, -1 / h));
     int iter = 0;