@@ -444,33 +444,33 @@ void check_smooth_friction_force_jacobian(
444444 // /////////////////////////////////////////////////////////////////////////
445445
446446 // test contact force norm derivative
447- {
448- const Eigen::MatrixXd lagged_positions = X + Ut;
449- Eigen::VectorXd normal_force_jacobian = Eigen::VectorXd::Zero (X.size ());
450- {
451- auto cc = create_smooth_collision (mesh, lagged_positions);
452- SmoothContactPotential<SmoothCollisions<dim>> potential (params);
453- Eigen::VectorXd g = potential.gradient (cc, mesh, lagged_positions);
454- Eigen::SparseMatrix<double > h = potential.hessian (cc, mesh, lagged_positions);
455- normal_force_jacobian = (h * g) / g.norm ();
456- }
457-
458- // finite difference
459- auto F_X = [&](const Eigen::VectorXd& x) {
460- Eigen::MatrixXd fd_X = fd::unflatten (x, dim);
461- const Eigen::MatrixXd fd_lagged_positions = fd_X + Ut;
462-
463- CollisionMesh fd_mesh (fd_X, mesh.edges (), mesh.faces ());
464- auto fd_cc = create_smooth_collision (fd_mesh, fd_lagged_positions);
465-
466- SmoothContactPotential<SmoothCollisions<dim>> potential (params);
467- return potential.gradient (fd_cc, fd_mesh, fd_lagged_positions).norm ();
468- };
469-
470- Eigen::VectorXd fd_normal_force_jacobian;
471- fd::finite_gradient (fd::flatten (X), F_X, fd_normal_force_jacobian);
472- CHECK ((normal_force_jacobian - fd_normal_force_jacobian).norm () <= 1e-7 * std::max (normal_force_jacobian.norm (), 1e-8 ));
473- }
447+ // {
448+ // const Eigen::MatrixXd lagged_positions = X + Ut;
449+ // Eigen::VectorXd normal_force_jacobian = Eigen::VectorXd::Zero(X.size());
450+ // {
451+ // auto cc = create_smooth_collision(mesh, lagged_positions);
452+ // SmoothContactPotential<SmoothCollisions<dim>> potential(params);
453+ // Eigen::VectorXd g = potential.gradient(cc, mesh, lagged_positions);
454+ // Eigen::SparseMatrix<double> h = potential.hessian(cc, mesh, lagged_positions);
455+ // normal_force_jacobian = (h * g) / g.norm();
456+ // }
457+
458+ // // finite difference
459+ // auto F_X = [&](const Eigen::VectorXd& x) {
460+ // Eigen::MatrixXd fd_X = fd::unflatten(x, dim);
461+ // const Eigen::MatrixXd fd_lagged_positions = fd_X + Ut;
462+
463+ // CollisionMesh fd_mesh(fd_X, mesh.edges(), mesh.faces());
464+ // auto fd_cc = create_smooth_collision(fd_mesh, fd_lagged_positions);
465+
466+ // SmoothContactPotential<SmoothCollisions<dim>> potential(params);
467+ // return potential.gradient(fd_cc, fd_mesh, fd_lagged_positions).norm();
468+ // };
469+
470+ // Eigen::VectorXd fd_normal_force_jacobian;
471+ // fd::finite_gradient(fd::flatten(X), F_X, fd_normal_force_jacobian);
472+ // CHECK((normal_force_jacobian - fd_normal_force_jacobian).norm() <= 1e-7 * std::max(normal_force_jacobian.norm(), 1e-8));
473+ // }
474474
475475 // /////////////////////////////////////////////////////////////////////////
476476
0 commit comments