Corrected small error in the dipole calculation. (#34)

Kostaszki · philipheinisch · commit 310c15f6884e · 2019-04-16T11:36:07.000+02:00
diff --git a/include_testcases/ideal_MHD_far_dipole.h b/include_testcases/ideal_MHD_far_dipole.h
@@ -33,7 +33,7 @@ inline REAL4 b_analytical(uint ix, uint iy, uint iz, REAL time) {
 inline REAL4 b_dipole(REAL4 X, REAL4 M) {
   REAL inorm = rsqrt(X.x*X.x + X.y*X.y + X.z*X.z);
   REAL4 n = inorm * X;
-  REAL4 H = (3 * n * (M * n) - M) * inorm * inorm * inorm;
+  REAL4 H = (3 * n * dot(M, n) - M) * inorm * inorm * inorm;
   return H;
 }
 

Original file line number	Diff line number	Diff line change
`@@ -33,7 +33,7 @@ inline REAL4 b_analytical(uint ix, uint iy, uint iz, REAL time) {`
`33`	`33`	`inline REAL4 b_dipole(REAL4 X, REAL4 M) {`
`34`	`34`	`REAL inorm = rsqrt(X.xX.x + X.yX.y + X.z*X.z);`
`35`	`35`	`REAL4 n = inorm * X;`
`36`		`- REAL4 H = (3 * n * (M * n) - M) * inorm * inorm * inorm;`
	`36`	`+ REAL4 H = (3 * n * dot(M, n) - M) * inorm * inorm * inorm;`
`37`	`37`	`return H;`
`38`	`38`	`}`
`39`	`39`