From 2c96a7d724d9b3e14644587f551b29029588df17 Mon Sep 17 00:00:00 2001
From: "Andrew E. Torda" <torda@zbh.uni-hamburg.de>
Date: Mon, 17 Jun 2024 15:04:43 +0200
Subject: [PATCH] rvec_dot() and rvec_sub_pbc() were declared static and not
 listed in rvec.h. This was confusing.

---
 rvec.c | 6 +++---
 rvec.h | 2 ++
 2 files changed, 5 insertions(+), 3 deletions(-)

diff --git a/rvec.c b/rvec.c
index ee3b6a0..5e9ca3c 100755
--- a/rvec.c
+++ b/rvec.c
@@ -20,7 +20,7 @@ void rvec_rand(rvec *u) {
 }
 
 /* dot product of u and v */
-static real rvec_dot(rvec u, rvec v) {
+real rvec_dot(rvec u, rvec v) {
     real r = 0;
     for (uint i = 0; i < NDIM; i++)
         r += u[i] * v[i];
@@ -39,7 +39,7 @@ void rvec_add(rvec *u, rvec v, rvec w) {
 }
 
 /* u = v - w */
-static void rvec_sub(rvec *u, rvec v, rvec w) {
+void rvec_sub(rvec *u, rvec v, rvec w) {
     for (uint i = 0; i < NDIM; i++)
         (*u)[i] = v[i] - w[i];
 }
@@ -58,7 +58,7 @@ void rvec_muladd(rvec *u, rvec v, real r, rvec w) {
 
 /* calculate v - w, taking periodic boundary conditions into
    account */
-static void rvec_sub_pbc(rvec *u, rvec box, rvec v, rvec w) {
+void rvec_sub_pbc(rvec *u, rvec box, rvec v, rvec w) {
     for (uint i = 0; i < NDIM; i++) {
         (*u)[i] = v[i] - w[i];
         while ((*u)[i] < -box[i]/2)
diff --git a/rvec.h b/rvec.h
index 017d403..1995762 100644
--- a/rvec.h
+++ b/rvec.h
@@ -17,5 +17,7 @@ void rvecary_muladd(uint n, rvec *a, rvec *b, real r, rvec *c);
 void rvec_muladd(rvec *u, rvec v, real r, rvec w);
 real rvec_len(rvec u);
 void rvec_rand(rvec *u);
+void rvec_sub_pbc(rvec *u, rvec box, rvec v, rvec w);
+real rvec_dot(rvec u, rvec v);
 #endif /* RVEC_H */
  
-- 
GitLab