From 3172bbbb5c96105e5b909fb10bb1f67e66c11ce8 Mon Sep 17 00:00:00 2001
From: vijay gopal chilkuri <vijay.gopal.c@gmail.com>
Date: Mon, 8 Mar 2021 22:16:50 +0100
Subject: [PATCH] Fixed a few bugs.

---
 Makefile             |  2 +-
 codelet_factor_een.f | 29 -----------------------------
 el_nuc_el_blas.irp.f |  8 +++-----
 3 files changed, 4 insertions(+), 35 deletions(-)
 delete mode 100644 codelet_factor_een.f

diff --git a/Makefile b/Makefile
index 6ed9503..164b364 100644
--- a/Makefile
+++ b/Makefile
@@ -1,4 +1,4 @@
-IRPF90 = irpf90 --codelet=jastrow_full:1000 #-s nelec:10 -s nnuc:2 -s ncord:5 #-a -d
+IRPF90 = irpf90 --codelet=factor_een_blas:1000 #-s nelec:10 -s nnuc:2 -s ncord:5 #-a -d
 FC     = ifort -xHost -g -mkl=sequential 
 FCFLAGS= -O2 -ffree-line-length-none -I .
 NINJA  = ninja
diff --git a/codelet_factor_een.f b/codelet_factor_een.f
deleted file mode 100644
index cc2b1e6..0000000
--- a/codelet_factor_een.f
+++ /dev/null
@@ -1,29 +0,0 @@
-
-program codelet_factor_een
-  implicit none
-  integer :: i
-  double precision :: ticks_0, ticks_1, cpu_0, cpu_1
-  integer, parameter :: irp_imax = 100000
-
-  
-
-  call provide_factor_een
-
-  double precision :: irp_rdtsc
-
-  call cpu_time(cpu_0)
-  ticks_0 = irp_rdtsc()
-  do i=1,irp_imax
-    call bld_factor_een
-  enddo
-  ticks_1 = irp_rdtsc()
-  call cpu_time(cpu_1)
-  print *, 'factor_een'
-  print *, '-----------'
-  print *, 'Cycles:'
-  print *,  (ticks_1-ticks_0)/dble(irp_imax)
-  print *, 'Seconds:'
-  print *,  (cpu_1-cpu_0)/dble(irp_imax)
-end
-
-  
\ No newline at end of file
diff --git a/el_nuc_el_blas.irp.f b/el_nuc_el_blas.irp.f
index 7640ec7..73f0d6b 100644
--- a/el_nuc_el_blas.irp.f
+++ b/el_nuc_el_blas.irp.f
@@ -28,7 +28,7 @@
  enddo
 
 
- END_PROVIDER
+END_PROVIDER
 
 
  BEGIN_PROVIDER [ double precision, factor_een_blas ]
@@ -41,13 +41,12 @@
 
  integer                        :: i, j, a, p, k, l, lmax, m, n
  double precision               :: accu
- double precision,dimension(:),allocatable       :: cn
- allocate(cn(ncord))
+ double precision               :: cn(ncord)
 
  factor_een_blas = 0.0d0
  factor_een_deriv_e_blas(1:4,1:nelec) = 0.0d0
 
- TOUCH tmp_c dtmp_c
+ PROVIDE tmp_c dtmp_c
 
  do n = 1, dim_cord_vect
 
@@ -88,5 +87,4 @@
    enddo
  enddo
 
- deallocate(cn)
 END_PROVIDER