rvec_dot() and rvec_sub_pbc() were declared static and not listed in

rvec.h. This was confusing.

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)

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 */