added more codes

Author: mhjensen

doc/Programs/OOcodes/codes.do.txt

@@ -2064,3 +2064,9 @@ for label in labels:
 
 
 
+!split
+===== Corresponding Qiskit implementations =====
+
+!bc pycod
+
+!ec

doc/Programs/OOcodes/qiskit/core/circuit.py

@@ -0,0 +1,26 @@
+from qiskit import QuantumCircuit, Aer, execute
+
+class Circuit:
+    def __init__(self, num_qubits):
+        self.num_qubits = num_qubits
+        self.qc = QuantumCircuit(num_qubits)
+    
+    def add_gate(self, gate_func):
+        gate_func(self.qc)
+    
+    def add_measure_all(self):
+        self.qc.measure_all()
+    
+    def run_statevector(self):
+        backend = Aer.get_backend('statevector_simulator')
+        job = execute(self.qc, backend)
+        result = job.result()
+        statevector = result.get_statevector()
+        return statevector
+    
+    def run_measurement(self, shots=1024):
+        backend = Aer.get_backend('qasm_simulator')
+        job = execute(self.qc, backend, shots=shots)
+        result = job.result()
+        counts = result.get_counts()
+        return counts

doc/Programs/OOcodes/qiskit/core/gates.py

@@ -0,0 +1,56 @@
+from qiskit import QuantumCircuit
+import numpy as np
+
+class Gates:
+    """Wrapper class for basic gates"""
+    
+    @staticmethod
+    def I(qc, qubit):
+        qc.id(qubit)
+    
+    @staticmethod
+    def X(qc, qubit):
+        qc.x(qubit)
+    
+    @staticmethod
+    def Y(qc, qubit):
+        qc.y(qubit)
+    
+    @staticmethod
+    def Z(qc, qubit):
+        qc.z(qubit)
+   @staticmethod
+    def H(qc, qubit):
+        qc.h(qubit)
+    
+    @staticmethod
+    def S(qc, qubit):
+        qc.s(qubit)
+    
+    @staticmethod
+    def T(qc, qubit):
+        qc.t(qubit)
+    
+    @staticmethod
+    def RX(qc, qubit, theta):
+        qc.rx(theta, qubit)
+    
+    @staticmethod
+    def RY(qc, qubit, theta):
+        qc.ry(theta, qubit)
+    
+    @staticmethod
+    def RZ(qc, qubit, theta):
+        qc.rz(theta, qubit)
+    
+    @staticmethod
+    def CNOT(qc, control, target):
+        qc.cx(control, target)
+    
+    @staticmethod
+    def CZ(qc, control, target):
+        qc.cz(control, target)
+    
+    @staticmethod
+    def SWAP(qc, qubit1, qubit2):
+        qc.swap(qubit1, qubit2)

doc/Programs/OOcodes/qiskit/examples/bell_test.py

@@ -0,0 +1,23 @@
+from states.bell_states import bell_state
+import matplotlib.pyplot as plt
+from qiskit.visualization import plot_histogram
+
+labels = ["Phi+", "Phi-", "Psi+", "Psi-"]
+shots = 1000
+
+for label in labels:
+    print(f"\nBell state {label}")
+    c = bell_state(label)
+    
+    # Run statevector simulation
+    state = c.run_statevector()
+    print("Statevector:", state)
+    
+    # Run measurements
+    c.add_measure_all()
+    counts = c.run_measurement(shots=shots)
+    print("Measurement counts:", counts)
+    
+    # Plot histogram
+    plot_histogram(counts, title=f"Bell state {label}")
+    plt.show()

doc/Programs/OOcodes/qiskit/readme.md

@@ -0,0 +1,24 @@
+
+# Modular Quantum Simulator (Qiskit Version)
+
+This is a fully modular quantum simulator package implemented with Qiskit.
+
+## Modules
+
+- `core/gates.py` — gate factory (wrapper around Qiskit gates)
+- `core/circuit.py` — modular circuit builder
+- `states/bell_states.py` — Bell state generators
+- `examples/bell_test.py` — Bell state experiment
+
+## Features
+
+- Fully based on Qiskit's QuantumCircuit
+- Supports any number of qubits
+- Native statevector simulation and measurement
+- Measurement histograms
+- Mirrors NumPy version structure for easy comparison
+
+## Dependencies
+
+- Qiskit
+- Matplotlib

doc/Programs/OOcodes/qiskit/states/bell_states.py

@@ -0,0 +1,24 @@
+from core.circuit import Circuit
+from core.gates import Gates
+
+def bell_state(label="Phi+"):
+    c = Circuit(2)
+    
+    # Create standard Bell basis preparation
+    Gates.H(c.qc, 0)
+    Gates.CNOT(c.qc, 0, 1)
+    
+    # Apply additional gates depending on label
+    if label == "Phi+":
+        pass
+    elif label == "Phi-":
+        Gates.Z(c.qc, 0)
+    elif label == "Psi+":
+        Gates.X(c.qc, 1)
+    elif label == "Psi-":
+        Gates.X(c.qc, 1)
+        Gates.Z(c.qc, 0)
+    else:
+        raise ValueError("Unknown Bell state")
+
+    return c