From 79a83c2cbee5adca798b8976b2a39ecd6ffd39af Mon Sep 17 00:00:00 2001
From: Franklin Wei <me@fwei.tk>
Date: Sat, 2 Feb 2019 23:02:12 -0500
Subject: initial commit

basic B field integrator
---
 test.cpp | 36 ++++++++++++++++++++++++++++++++++++
 1 file changed, 36 insertions(+)
 create mode 100644 test.cpp

(limited to 'test.cpp')

diff --git a/test.cpp b/test.cpp
new file mode 100644
index 0000000..8d20829
--- /dev/null
+++ b/test.cpp
@@ -0,0 +1,36 @@
+#include <cmath>
+#include <iostream>
+#include "curve.h"
+
+vec3 integrand(vec3 s, vec3 ds)
+{
+    return ds;
+}
+
+vec3 point;
+
+/* dl x r / (|r| ^ 2) */
+vec3 dB(vec3 s, vec3 ds)
+{
+    vec3 r = s - point;
+
+    double r2 = r.magnitudeSquared();
+
+    vec3 rnorm = r / std::sqrt(r2);
+
+    return ds.cross(rnorm) / r2;
+}
+
+int main()
+{
+    LineSegment wire(vec3(0, -100, 0), vec3(0, 100, 0));
+
+    vec3 dx = .001;
+
+    point = .01;
+
+    double I = 1;
+
+    for(int i = 0; i < 100; i++, point += dx)
+        std::cout << point[0] << " " << wire.integrate(dB, 1e-2)[2] << endl;
+}
-- 
cgit v1.1