From 69d402fa7b7e004827c20351b80ffdffe46477f6 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Mon, 5 Dec 2022 18:01:23 -0500
Subject: [PATCH] handle dead code

---
 src/ML-POD/mlpod.cpp | 21 +++++++++------------
 1 file changed, 9 insertions(+), 12 deletions(-)

diff --git a/src/ML-POD/mlpod.cpp b/src/ML-POD/mlpod.cpp
index 2712dd6b65..72c732b701 100644
--- a/src/ML-POD/mlpod.cpp
+++ b/src/ML-POD/mlpod.cpp
@@ -1641,39 +1641,37 @@ void MLPOD::InitSnap()
 
   double rcutmax = 0.0;
   for (int ielem = 0; ielem < ntypes; ielem++)
-  rcutmax = MAX(2.0*elemradius[ielem]*rcutfac,rcutmax);
+    rcutmax = MAX(2.0*elemradius[ielem]*rcutfac,rcutmax);
 
   snapSetup(twojmax, ntypes);
-  //TemplateCopytoDevice(&sna.radelem[1], elemradius, ntypes, backend);
-  //TemplateCopytoDevice(&sna.wjelem[1], elemweight, ntypes, backend);
   for (int i=0; i<ntypes; i++) {
     sna.radelem[1+i] = elemradius[i];
     sna.wjelem[1+i] = elemweight[i];
   }
   podArrayFill(&sna.map[1], (int) 0, ntypes);
 
-  //double cutsq[100];
-  for (int i=0; i<ntypes; i++)
+  for (int i=0; i<ntypes; i++) {
     for (int j=0; j<ntypes; j++) {
       double cut = (elemradius[i] + elemradius[j])*rcutfac;
       sna.rcutsq[j+1 + (i+1)*(ntypes+1)] = cut*cut;
     }
-  //TemplateCopytoDevice(sna.rcutsq, cutsq, (ntypes+1)*(ntypes+1), backend);
+  }
 
+  // bzeroflag is currently always 0
+#if 0
   if (bzeroflag) {
     double www = wself*wself*wself;
-    //double bzero[100];
     for (int j = 0; j <= twojmax; j++)
       if (bnormflag)
-      sna.bzero[j] = www;
+        sna.bzero[j] = www;
       else
-      sna.bzero[j] = www*(j+1);
-    //TemplateCopytoDevice(sna.bzero, bzero, twojmax+1, backend);
+        sna.bzero[j] = www*(j+1);
   }
+#endif
 
   int nelements = ntypes;
   if (!chemflag)
-  nelements = 1;
+    nelements = 1;
 
   sna.nelements = nelements;
   sna.ndoubles = nelements*nelements;   // number of multi-element pairs
@@ -2599,7 +2597,6 @@ void MLPOD::pod3body(double *eatom, double *yij, double *e2ij, double *tmpmem, i
         costhe = xdot/(rij*rik);
         costhe = costhe > 1.0 ? 1.0 : costhe;
         costhe = costhe < -1.0 ? -1.0 : costhe;
-        xdot = costhe*(rij*rik);
         theta = acos(costhe);
 
         for (int p=0; p <nabf1; p++)