From dfeae720d3d0933c498f95a764422a4f914275a8 Mon Sep 17 00:00:00 2001
From: Sebastian Ehlert <28669218+awvwgk@users.noreply.github.com>
Date: Mon, 16 Mar 2020 10:22:46 +0100
Subject: [PATCH] Fix: Parallelisation of numerical Hessian (#146)

- declare strain derivs in omp as private
---
 src/hessian.F90 | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/hessian.F90 b/src/hessian.F90
index fab06dde4..2c9f6ffd8 100644
--- a/src/hessian.F90
+++ b/src/hessian.F90
@@ -235,9 +235,9 @@ subroutine numhess( &
 !! ------------------------------------------------------------------------
       !$ nproc = omp_get_num_threads()
       !$omp parallel default(shared) &
-      !$omp&         firstprivate(mol,calc,et,maxiter,acc,wf0) &
-      !$omp&         private(ia,ic,ii,ja,jc,jj,eel,gr,gl,egap,sccr,sccl,wfx,tmol) &
-      !$omp&         shared (h,dipd,pold,step,step2,t1,t0,w1,w0,xyzsave,env)
+      !$omp firstprivate(mol,calc,et,maxiter,acc,wf0) &
+      !$omp private(ia,ic,ii,ja,jc,jj,eel,gr,gl,egap,sccr,sccl,sr,sl,wfx,tmol) &
+      !$omp shared (h,dipd,pold,step,step2,t1,t0,w1,w0,xyzsave,env)
       !$ call omp_set_num_threads(1)
 #ifdef WITH_MKL
       !$ call mkl_set_num_threads(1)