From 4c813be666c604b5a606f96db65128ca9fae0b3f Mon Sep 17 00:00:00 2001
From: Iain Bethune <iain.bethune@stfc.ac.uk>
Date: Mon, 14 May 2018 08:14:01 +0000
Subject: [PATCH] Bugfix (Vittorio Zecca, #53)

svn-origin-rev: 18422
---
 src/cp_ddapc_forces.F | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/cp_ddapc_forces.F b/src/cp_ddapc_forces.F
index 6b1bef6823..ec939a3957 100644
--- a/src/cp_ddapc_forces.F
+++ b/src/cp_ddapc_forces.F
@@ -563,7 +563,7 @@ CONTAINS
             factor2 = 0.0_dp
             IF (pos1*pos2 /= 0.0_dp) THEN
                pos1i = 1.0_dp/pos1
-               pos2i = 1.0_dp/pos1
+               pos2i = 1.0_dp/pos2
                dpos1 = pos1i*r1
                dpos2 = pos2i*r2
                ptcos = DOT_PRODUCT(r1, r2)