Changed sign on 2D particle ind vel

hjabird · Jun 17, 2019 · c6e8e2f · c6e8e2f
1 parent cfc6b87
commit c6e8e2f
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/P2D.c b/src/P2D.c
@@ -66,7 +66,7 @@ CVTX_EXPORT bsv_V2f cvtx_P2D_S2S_vel(
 	bsv_V2f ret;
 	ret = P2D_vel_inner(self, mes_point, kernel,
 		1.f / fabsf(regularisation_radius));
-	return bsv_V2f_mult(ret, -1.f / (2.f * acosf(-1.f)));
+	return bsv_V2f_mult(ret, 1.f / (2.f * acosf(-1.f)));
 }
 
 CVTX_EXPORT bsv_V2f cvtx_P2D_M2S_vel(
@@ -88,7 +88,7 @@ CVTX_EXPORT bsv_V2f cvtx_P2D_M2S_vel(
 		ry += vel.x[1];
 	}
 	bsv_V2f ret = { (float)rx, (float)ry };
-	return bsv_V2f_mult(ret, -1.f / (2.f * acosf(-1.f)));
+	return bsv_V2f_mult(ret, 1.f / (2.f * acosf(-1.f)));
 }
 
 

diff --git a/src/nbody.cl b/src/nbody.cl
@@ -206,7 +206,7 @@ Definitions for the repeated body of kernels
 /* Fill in g calc here */
 
 "#define CVTX_P2D_VEL_END 											\\\n"
-"	cor = - g;	/*1/2pi term is done by host. */					\\\n"
+"	cor = g;	/*1/2pi term is done by host. */					\\\n"
 "	den = pown(radd, 2);											\\\n"
 "	ret.x = rad.y * (cor * particle_vorts[pidx] / den);				\\\n"
 "	ret.y = -rad.x * (cor * particle_vorts[pidx] / den);			\\\n"