End result BAP

examples/Decoherence/Decoherence.ino

@@ -0,0 +1,86 @@
+#include <Qbead.h>
+
+Qbead::Qbead bead;
+int rotationState = 0;
+uint32_t stateColor = color(255, 255, 255);
+uint32_t decoherenceColor = color(122, 0, 122);
+const bool toggleAnimationOn = 1;
+Qbead::Coordinates oldCoordinates(0, 0, 1);
+int t = 0;
+bool wasFrozen = false;
+
+void setup()
+{
+    bead.begin();
+    bead.setBrightness(25);
+    Serial.println("testing all pixels discretely");
+    for (int i = 0; i < bead.pixels.numPixels(); i++)
+    {
+        bead.pixels.setPixelColor(i, color(255, 255, 255));
+        bead.pixels.show();
+        delay(5);
+    }
+    Serial.println("testing smooth transition between pixels");
+    for (int phi = 0; phi < 360; phi += 30)
+    {
+        for (int theta = 0; theta < 180; theta += 3)
+        {
+            bead.clear();
+            bead.setBloch_deg(theta, phi, colorWheel_deg(phi));
+            bead.show();
+        }
+    }
+    Serial.println("starting inertial tracking");
+    oldCoordinates = bead.state.getCoordinates();
+    t = millis();
+}
+
+void loop()
+{
+    bead.readIMU(true);
+    stateColor = color(255, 255, 255);
+    Serial.print("rotationState: ");
+    Serial.println(rotationState);
+    if (bead.frozen)
+    {
+        stateColor = color(122, 122, 0);
+        wasFrozen = true;
+    }
+    else
+    {
+        if (wasFrozen)
+        {
+            wasFrozen = false;
+            oldCoordinates = bead.state.getCoordinates();
+        }
+        rotationState = bead.checkMotion();
+        if (rotationState != 0)
+        {
+            bead.frozen = true;
+            bead.T_freeze = millis();
+        }
+        float phi = bead.state.getCoordinates().phi();
+        int dt = millis() - t;
+        float randInt = random(200, 10000);
+        phi += dt / randInt;
+        if (phi > 2 * PI)
+        {
+            phi -= 2 * PI;
+        }
+        Qbead::Coordinates newCoordinates(bead.state.getCoordinates().theta(), phi);
+        bead.state.setCoordinates(newCoordinates);
+        bead.visualState = bead.state.getCoordinates();
+    }
+    t = millis();
+    bead.animateTo(rotationState, 2000);
+    if (bead.T_led + 20 < millis()) {
+        bead.T_led = millis();
+        bead.clear();
+        if (!bead.frozen) {
+            bead.setLed(oldCoordinates, decoherenceColor);
+        }
+        bead.showAxis();
+        bead.setLed(bead.visualState, stateColor, 1);
+        bead.show();
+    }
+}

examples/IMU_reader/IMU_reader.ino

@@ -25,6 +25,6 @@ void setup() {
 void loop() {
   bead.readIMU();
   bead.clear();
-  bead.setBloch_deg_smooth(bead.t_acc, bead.p_acc, color(255, 0, 255));
+  bead.setLed(Qbead::Coordinates(bead.gravityVector), color(255, 0, 255), true);
   bead.show();
 }

examples/PauliGate_detection/PauliGate_detection.ino

@@ -0,0 +1,54 @@
+#include <Qbead.h>
+
+Qbead::Qbead bead;
+int rotationState = 0;
+uint32_t stateColor = color(255, 255, 255);
+const bool toggleAnimationOn = 1;
+
+void setup() {
+  bead.begin();
+  bead.setBrightness(25); // way too bright
+  Serial.println("testing all pixels discretely");
+  for (int i = 0; i < bead.pixels.numPixels(); i++) {
+    bead.pixels.setPixelColor(i, color(255, 255, 255));
+    bead.pixels.show();
+    delay(5);
+  }
+  Serial.println("testing smooth transition between pixels");
+  for (int phi = 0; phi < 360; phi += 30) {
+    for (int theta = 0; theta < 180; theta += 3) {
+      bead.setBloch_deg(theta, phi, colorWheel_deg(phi));
+      bead.show();
+    }
+  }
+  Serial.println("starting inertial tracking");
+}
+
+void loop() {
+  bead.readIMU(false);
+  stateColor = color(255, 255, 255);
+  Serial.print("rotationState: ");
+  Serial.println(rotationState);
+  if (bead.frozen)
+  {
+    stateColor = color(122, 122, 0);
+  }
+  else
+  {
+    rotationState = bead.checkMotion();
+    if (rotationState != 0)
+    {
+      bead.frozen = true;
+      bead.T_freeze = millis();
+    }
+  }
+  bead.animateTo(rotationState, 2000);
+  // The leds can only be updated every 20ms
+  if (bead.T_led + 20 < millis()) {
+    bead.T_led = millis();
+    bead.clear();
+    bead.showAxis();
+    bead.setLed(bead.visualState, stateColor);
+    bead.show();
+  }
+}

generateCoordinates.py

@@ -0,0 +1,46 @@
+# This is a script to convert the pick and place file to a coordinates map
+# This code is very much not optimized, but it does its job.
+
+import math
+
+inputFile = "../pick-and-place.csv"
+outputFile = "../pick-and-place-coordinates.csv"
+
+lines = []
+
+with open(inputFile, 'r') as file:
+    for line in file:
+        if line.strip():  # Check if the line is not empty
+            lines.append(line.strip())
+
+x0 = float(lines[0].split(',')[3])
+y0 = float(lines[0].split(',')[4])
+
+coordinates = []
+for line in lines:
+    parts = line.split(',')
+    if len(parts) >= 2:
+        x = float(parts[3]) - x0
+        y = float(parts[4]) - y0
+        coordinates.append((x, y))
+
+# Convert to spherical coordinates
+spherical_coordinates = []
+for x, y in coordinates:
+    r = (x**2 + y**2)**0.5
+    phi = math.atan2(y, x)  # angle in radians
+    if phi < 0:
+        phi += 2 * math.pi  # Normalize angle to [0, 2π)
+    spherical_coordinates.append((r, phi))
+
+# scale r from 0 to pi
+scaled_coordinates = []
+for r, phi in spherical_coordinates:
+    scaled_r = r / max(r for r, _ in spherical_coordinates) * math.pi
+    scaled_coordinates.append((round(scaled_r, 2), round(phi, 2)))
+
+# Write to output file
+with open(outputFile, 'w') as file:
+    for theta, phi in scaled_coordinates:
+        file.write(f"Coordinates({theta}, {phi}),\n")
+print(f"Coordinates written to {outputFile}")