diff --git a/rvec.c b/rvec.c
index ee3b6a05fc01221ea83ec3eb5096c2ad5421ea2f..5e9ca3c66b9a7e797dad08a61d8366d88e65eef4 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 017d403077f582b14de1100cf9c2ea8d2dff9fda..19957621015e43287a85ece06b3bce0608d0ee95 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 */