@@ -259,14 +259,14 @@ auto SmoothCollisionTemplate<max_vert, PrimitiveA, PrimitiveB>::gradient(
259259}
260260
261261template  <int  max_vert, typename  PrimitiveA, typename  PrimitiveB>
262- Eigen::MatrixXd 
262+ auto 
263263SmoothCollisionTemplate<max_vert, PrimitiveA, PrimitiveB>::hessian(
264264    const  Vector<
265265        double ,
266266        -1 ,
267267        SmoothCollisionTemplate<max_vert, PrimitiveA, PrimitiveB>::max_size>&
268268        positions,
269-     const  ParameterType& params) const 
269+     const  ParameterType& params) const  -> Eigen::MatrixXd 
270270{
271271    const  auto  core_indices = get_core_indices ();
272272
@@ -387,17 +387,15 @@ SmoothCollisionTemplate<max_vert, PrimitiveA, PrimitiveB>::hessian(
387387    //  grad of tangent/normal terms
388388    double  orient = 0 ;
389389    Vector<double , -1 , max_size> gOrient ;
390-     MatrixMax< double , max_size, max_size>  hOrient;
390+     Eigen::MatrixXd  hOrient;
391391    {
392392        Vector<double , -1 , max_size> gA  = Vector<double , -1 , max_size>::Zero (
393393                                         n_dofs ()),
394394                                     gB  = Vector<double , -1 , max_size>::Zero (
395395                                         n_dofs ());
396-         MatrixMax<double , max_size, max_size>
397-             hA =
398-                 MatrixMax<double , max_size, max_size>::Zero (n_dofs (), n_dofs ()),
399-             hB =
400-                 MatrixMax<double , max_size, max_size>::Zero (n_dofs (), n_dofs ());
396+         Eigen::MatrixXd
397+             hA = Eigen::MatrixXd::Zero (n_dofs (), n_dofs ()),
398+             hB = Eigen::MatrixXd::Zero (n_dofs (), n_dofs ());
401399        {
402400            gA (core_indices) =
403401                closest_direction_grad.transpose () * gA_reduced .head (dim);
0 commit comments