fixup: dilatation compile

Author: streeve

src/CabanaPD_Force.hpp

@@ -224,7 +224,7 @@ class Dilatation<MemorySpace, ModelType<PDModelType, MechanicsType, NoFracture,
             // Get the bond distance, displacement, and stretch.
             double xi, r, s;
             getDistance( x, u, i, j, xi, r, s );
-            theta( i ) += model.dilatation( s, xi, vol( j ), m( i ) );
+            theta( i ) += model.dilatation( i, s, xi, vol( j ), m( i ) );
         };
 
         neighbor.iterate( exec_space{}, dilatation, particles,
@@ -338,8 +338,8 @@ class Dilatation<MemorySpace, ModelType<PDModelType, MechanicsType, Fracture,
                 // broken, because m=0 only occurs when all bonds are broken.
                 // mu is still included to account for individual bond breaking.
                 if ( m( i ) > 0 )
-                    theta( i ) += mu( i, n ) *
-                                  model.dilatation( s, xi, vol( j ), m( i ) );
+                    theta( i ) += mu( i, n ) * model.dilatation(
+                                                   i, s, xi, vol( j ), m( i ) );
             }
         };
 

src/force/CabanaPD_ForceModels_LPS.hpp

@@ -63,8 +63,7 @@ struct ForceModel<LPS, Elastic, NoFracture> : public BaseForceModel
 
     KOKKOS_INLINE_FUNCTION auto dilatation( const int, const double s,
                                             const double xi, const double vol,
-                                            const double m_i,
-                                            const double ) const
+                                            const double m_i ) const
     {
         double theta_i = influenceFunction( xi ) * s * xi * xi * vol;
         return 3.0 * theta_i / m_i;

src/force/CabanaPD_ForceModels_PMB.hpp

@@ -224,8 +224,7 @@ struct ForceModel<PMB, ElasticPerfectlyPlastic, Fracture,
     // Fused density update using plastic dilatation.
     KOKKOS_INLINE_FUNCTION auto dilatation( const int i, const double s,
                                             const double xi, const double vol,
-                                            const double m_i,
-                                            const double theta_i_n ) const
+                                            const double ) const
     {
         double theta_i = coeff * s * xi * vol;