update

src/mesa_pd/kernel/VBondModel/AnisotropicVDWContact.h

@@ -50,9 +50,9 @@ class AnisotropicVDWContact
 {
 public:
    template<typename Accessor>
-   real_t operator()(const size_t p_idx1,
-                     const size_t p_idx2,
-                     Accessor &ac);
+   void operator()(const size_t p_idx1,
+                   const size_t p_idx2,
+                   Accessor &ac);
 
    static constexpr real_t eps_ = 0.03733_r;
    static constexpr real_t A_ = 0.0223_r;
@@ -89,9 +89,9 @@ public:
  */
 template<typename Accessor>
 inline
-real_t AnisotropicVDWContact::operator()(const size_t p_idx1,
-                                         const size_t p_idx2,
-                                         Accessor &ac)
+void AnisotropicVDWContact::operator()(const size_t p_idx1,
+                                       const size_t p_idx2,
+                                       Accessor &ac)
 {
    //===Adaptation of PFC5 vdW contact model implementation====
 
@@ -297,6 +297,7 @@ real_t AnisotropicVDWContact::operator()(const size_t p_idx1,
       real_t W_u = 1_r;
       real_t W_u_L = 0_r;
       real_t W_u_LL = 0_r;
+      WALBERLA_UNUSED(W_u_LL);
 
       // Cubic cutoff function 3T->4T (hardcoded since we do not need to mess w these parameters)
       constexpr auto Q1_ = -80.0_r;
@@ -341,8 +342,10 @@ real_t AnisotropicVDWContact::operator()(const size_t p_idx1,
       real_t GA = 1_r;
       real_t GA_L = 0_r;
       real_t GA_LL = 0_r;
+      WALBERLA_UNUSED(GA_LL);
       real_t GA_G = 0_r;
       real_t GA_GG = 0_r;
+      WALBERLA_UNUSED(GA_GG);
       if (std::abs(sin_gamma) > TOL)
       {
          GA = 1_r + W_ga * (1_r - cos_2gamma);