1 files changed, 27 insertions, 0 deletions
diff --git a/src/slalib/dvdv.f b/src/slalib/dvdv.f
new file mode 100644
index 0000000..ab3e411
--- /dev/null
+++ b/src/slalib/dvdv.f
@@ -0,0 +1,27 @@
+      DOUBLE PRECISION FUNCTION sla_DVDV (VA, VB)
+*+
+*     - - - - -
+*      D V D V
+*     - - - - -
+*
+*  Scalar product of two 3-vectors  (double precision)
+*
+*  Given:
+*      VA      dp(3)     first vector
+*      VB      dp(3)     second vector
+*
+*  The result is the scalar product VA.VB (double precision)
+*
+*  P.T.Wallace   Starlink   November 1984
+*
+*  Copyright (C) 1995 Rutherford Appleton Laboratory
+*-
+
+      IMPLICIT NONE
+
+      DOUBLE PRECISION VA(3),VB(3)
+
+
+      sla_DVDV=VA(1)*VB(1)+VA(2)*VB(2)+VA(3)*VB(3)
+
+      END