From ee14ad272e68d9363202d7f668e0b20302827209 Mon Sep 17 00:00:00 2001
From: hachem <im@hachem.wtf>
Date: Mon, 14 Sep 2026 12:20:52 +0200
Subject: feat: simd + testing + formatting

---
 tester/non_clifford.c | 42 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100644 tester/non_clifford.c

(limited to 'tester/non_clifford.c')

diff --git a/tester/non_clifford.c b/tester/non_clifford.c
new file mode 100644
index 0000000..c361fb5
--- /dev/null
+++ b/tester/non_clifford.c
@@ -0,0 +1,42 @@
+#include "tests.h"
+
+static const double R2 = 0.7071067811865476;
+static const double PI = 3.141592653589793;
+
+void run_non_clifford_tests(void)
+{
+	psi_test_section("Non-Clifford gates");
+
+	struct PsiQuantumCircuit t = psi_new_quantum_circuit(1);
+	psi_apply_x(&t, 0);
+	psi_apply_t(&t, 0);
+	struct PsiComplex t_exp[] = { psi_new_complex(0.0, 0.0), psi_new_complex(R2, R2) };
+	psi_check_circuit("T|1> = e^{i pi/4}|1>", &t, t_exp, 2);
+	psi_free_quantum_circuit(&t);
+
+	struct PsiQuantumCircuit rx = psi_new_quantum_circuit(1);
+	psi_apply_rx(&rx, 0, PI);
+	struct PsiComplex rx_exp[] = { psi_new_complex(0.0, 0.0), psi_new_complex(0.0, -1.0) };
+	psi_check_circuit("Rx(pi)|0> = -i|1>", &rx, rx_exp, 2);
+	psi_free_quantum_circuit(&rx);
+
+	struct PsiQuantumCircuit ry = psi_new_quantum_circuit(1);
+	psi_apply_ry(&ry, 0, PI);
+	struct PsiComplex ry_exp[] = { psi_new_complex(0.0, 0.0), psi_new_complex(1.0, 0.0) };
+	psi_check_circuit("Ry(pi)|0> = |1>", &ry, ry_exp, 2);
+	psi_free_quantum_circuit(&ry);
+
+	struct PsiQuantumCircuit p = psi_new_quantum_circuit(1);
+	psi_apply_x(&p, 0);
+	psi_apply_p(&p, 0, PI / 2.0);
+	struct PsiComplex p_exp[] = { psi_new_complex(0.0, 0.0), psi_new_complex(0.0, 1.0) };
+	psi_check_circuit("P(pi/2)|1> = i|1>", &p, p_exp, 2);
+	psi_free_quantum_circuit(&p);
+
+	struct PsiQuantumCircuit sx = psi_new_quantum_circuit(1);
+	psi_apply_sx(&sx, 0);
+	psi_apply_sx(&sx, 0);
+	struct PsiComplex sx_exp[] = { psi_new_complex(0.0, 0.0), psi_new_complex(1.0, 0.0) };
+	psi_check_circuit("sqrt(X) applied twice = X", &sx, sx_exp, 2);
+	psi_free_quantum_circuit(&sx);
+}
-- 
cgit v1.3