From 71d868df177cec1ff52a3482e86e5cf92cb29ad9 Mon Sep 17 00:00:00 2001
From: Jeff Hammond <jeff.r.hammond@intel.com>
Date: Tue, 16 Jun 2020 14:22:05 -0700
Subject: [PATCH] replace selci_popcnt32 with popcnt

---
 src/selci/iodiff.F   | 1 -
 src/selci/ointer.F   | 8 --------
 src/selci/popcnt32.F | 9 ---------
 3 files changed, 18 deletions(-)

diff --git a/src/selci/iodiff.F b/src/selci/iodiff.F
index e61d92008b..83f55db789 100644
--- a/src/selci/iodiff.F
+++ b/src/selci/iodiff.F
@@ -82,7 +82,6 @@ c
             enddo
          endif
 c
-         !call selci_vpop32(nk, tmp, out)
          do i = 1,nk
            out(i) = popcnt(tmp(i))
          enddo
diff --git a/src/selci/ointer.F b/src/selci/ointer.F
index 7c2a1d7216..45bbeb91f8 100644
--- a/src/selci/ointer.F
+++ b/src/selci/ointer.F
@@ -21,10 +21,6 @@ c
 c     use ieor and popcnt to compute the occupation difference
 c     between the I and J occupations represented as bit patterns
 c     
-c     strip mine into blocks so can call vpop32 instead of
-c     popcnt(popcnt32 -> 100,000/s, vpop32 -> 540,000/s at vl=32
-c     480,000/s at vl=16, 740,000/s at vl=1000)
-c     
 c     Search backwards since it is most likely to find something
 c     that interacts near to the occupation that generated the
 c     excitation being tested.
@@ -39,7 +35,6 @@ c
          do i = 1,ndo
             iwrk1(i) = ieor(itest,ioconf(ioff+i))
          enddo
-         !call selci_vpop32(ndo,iwrk1,iwrk2)
          do i = 1,ndo
            iwrk2(i) = popcnt(iwrk1(i))
          enddo
@@ -54,7 +49,6 @@ c
             iwrk1(i  ) = ieor(itest1,ioconf(ioff+i  ))
             iwrk1(i+1) = ieor(itest2,ioconf(ioff+i+1))
          enddo
-         !call selci_vpop32(ndo2,iwrk1,iwrk2)
          do i = 1,ndo2
            iwrk2(i) = popcnt(iwrk1(i))
          enddo
@@ -71,7 +65,6 @@ c
             iwrk1(i+1) = ieor(itest2,ioconf(ioff+i+1))
             iwrk1(i+2) = ieor(itest3,ioconf(ioff+i+2))
          enddo
-         !call selci_vpop32(ndo3,iwrk1,iwrk2)
          do i = 1,ndo3
            iwrk2(i) = popcnt(iwrk1(i))
          enddo
@@ -87,7 +80,6 @@ c
             enddo
             k = k + nintpo
          enddo
-         !call selci_vpop32(ndo*nintpo,iwrk1,iwrk2)
          do i = 1,ndo*nintpo
            iwrk2(i) = popcnt(iwrk1(i))
          enddo
diff --git a/src/selci/popcnt32.F b/src/selci/popcnt32.F
index da0944e01d..e69de29bb2 100644
--- a/src/selci/popcnt32.F
+++ b/src/selci/popcnt32.F
@@ -1,9 +0,0 @@
-      subroutine selci_vpop32(n, input, output)
-      integer n
-      integer input(*)
-      integer output(*)
-      integer i
-      do i = 1,n
-        output(i) = popcnt(input(i))
-      enddo
-      end