New Amplitude-Based Components

demos/simulation/ab_simulation.ipynb

@@ -0,0 +1,141 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "348de144",
+   "metadata": {},
+   "source": [
+    "# Simulating the Amplitude-Based Collisionless QLBM"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "181af34e",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from qiskit_aer import AerSimulator\n",
+    "\n",
+    "from qlbm.components import (\n",
+    "    CQLBM,\n",
+    "    ABGridMeasurement,\n",
+    "    ABInitialConditions,\n",
+    "    EmptyPrimitive,\n",
+    ")\n",
+    "from qlbm.infra import QiskitRunner, SimulationConfig\n",
+    "from qlbm.lattice import ABLattice\n",
+    "from qlbm.tools.utils import create_directory_and_parents\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "05e4915d",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "lattice = ABLattice(\n",
+    "    {\n",
+    "        \"lattice\": {\"dim\": {\"x\": 64, \"y\": 32}, \"velocities\": \"d2q9\"},\n",
+    "        \"geometry\": [\n",
+    "            {\"shape\": \"cuboid\", \"x\": [10, 13], \"y\": [14, 17], \"boundary\": \"bounceback\"},\n",
+    "        ],\n",
+    "    }\n",
+    ")\n",
+    "\n",
+    "\n",
+    "output_dir = \"qlbm-output/ab-d2q9-64x32-1-obstacle-qiskit\"\n",
+    "create_directory_and_parents(output_dir)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "deb472a0",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "cfg = SimulationConfig(\n",
+    "    initial_conditions=ABInitialConditions(lattice),\n",
+    "    algorithm=CQLBM(lattice),\n",
+    "    postprocessing=EmptyPrimitive(lattice),\n",
+    "    measurement=ABGridMeasurement(lattice),\n",
+    "    target_platform=\"QISKIT\",\n",
+    "    compiler_platform=\"QISKIT\",\n",
+    "    optimization_level=0,\n",
+    "    statevector_sampling=True,\n",
+    "    execution_backend=AerSimulator(method=\"statevector\"),\n",
+    "    sampling_backend=AerSimulator(method=\"statevector\"),\n",
+    ")\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "57240678",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "cfg.prepare_for_simulation()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "1da18eb8",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Number of shots to simulate for each timestep when running the circuit\n",
+    "NUM_SHOTS = 2**12\n",
+    "\n",
+    "# Number of timesteps to simulate\n",
+    "NUM_STEPS = 20"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "5c8c7d1e",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "runner = QiskitRunner(\n",
+    "    cfg,\n",
+    "    lattice,\n",
+    ")\n",
+    "\n",
+    "\n",
+    "# Simulate the circuits using both snapshots\n",
+    "runner.run(\n",
+    "    NUM_STEPS,  # Number of time steps\n",
+    "    NUM_SHOTS,  # Number of shots per time step\n",
+    "    output_dir,\n",
+    "    statevector_snapshots=True,\n",
+    ")"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "qlbm-cpu-venv",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.13.9"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

demos/simulation/lqlga_simulation.ipynb

@@ -1,5 +1,13 @@
 {
  "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "4c316e35",
+   "metadata": {},
+   "source": [
+    "# Simulating the Linear Encoding Quantum Lattice Gas Automata"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -10,12 +18,14 @@
     "from qiskit_aer import AerSimulator\n",
     "\n",
     "from qlbm.components import EmptyPrimitive\n",
-    "from qlbm.components.lqlga.initial import LQGLAInitialConditions\n",
-    "from qlbm.components.lqlga.lqlga import LQLGA\n",
-    "from qlbm.components.lqlga.measurement import LQLGAGridVelocityMeasurement\n",
+    "from qlbm.components.lqlga import (\n",
+    "    LQLGA,\n",
+    "    LQGLAInitialConditions,\n",
+    "    LQLGAGridVelocityMeasurement,\n",
+    ")\n",
     "from qlbm.infra import QiskitRunner, SimulationConfig\n",
-    "from qlbm.lattice.lattices.lqlga_lattice import LQLGALattice\n",
-    "from qlbm.tools.utils import create_directory_and_parents\n"
+    "from qlbm.lattice import LQLGALattice\n",
+    "from qlbm.tools.utils import create_directory_and_parents"
    ]
   },
   {
@@ -25,9 +35,6 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "from qlbm.components.lqlga import LQGLAInitialConditions\n",
-    "from qlbm.lattice import LQLGALattice\n",
-    "\n",
     "lattice = LQLGALattice(\n",
     "    {\n",
     "        \"lattice\": {\n",
@@ -69,18 +76,34 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "4d54ef63",
+   "id": "65ec52e8",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "cfg.prepare_for_simulation()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "11c5dbb7",
    "metadata": {},
    "outputs": [],
    "source": [
-    "cfg.prepare_for_simulation()\n",
-    "\n",
     "# Number of shots to simulate for each timestep when running the circuit\n",
     "NUM_SHOTS = 2**10\n",
     "\n",
     "# Number of timesteps to simulate\n",
-    "NUM_STEPS = 20\n",
-    "\n",
+    "NUM_STEPS = 20"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "4d54ef63",
+   "metadata": {},
+   "outputs": [],
+   "source": [
     "# Create a runner object to simulate the circuit\n",
     "runner = QiskitRunner(\n",
     "    cfg,\n",
@@ -121,7 +144,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.13.7"
+   "version": "3.13.9"
   }
  },
  "nbformat": 4,

demos/simulation/collisionless_simulation.ipynb → demos/simulation/ms_simulation.ipynb

@@ -1,71 +1,69 @@
 {
  "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "328661be",
+   "metadata": {},
+   "source": [
+    "# Simulation the Multi-Speed Collisionless QLBM"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "9b1814b3",
+   "id": "be6f3d11",
    "metadata": {},
    "outputs": [],
    "source": [
     "from qiskit_aer import AerSimulator\n",
     "\n",
     "from qlbm.components import (\n",
     "    CQLBM,\n",
-    "    CollisionlessInitialConditions,\n",
     "    EmptyPrimitive,\n",
     "    GridMeasurement,\n",
+    "    MSInitialConditions,\n",
     ")\n",
     "from qlbm.infra import QiskitRunner, SimulationConfig\n",
-    "from qlbm.lattice import CollisionlessLattice\n",
-    "from qlbm.tools.utils import create_directory_and_parents"
+    "from qlbm.lattice import MSLattice\n",
+    "from qlbm.tools.utils import create_directory_and_parents\n"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "e8e77d2f-a778-441c-b436-a61e1f3154cf",
+   "id": "e8d2b671",
    "metadata": {},
    "outputs": [],
    "source": [
-    "# Load example with mixed boundary conditions and create output directory\n",
-    "lattice = CollisionlessLattice(\n",
+    "lattice = MSLattice(\n",
     "    {\n",
-    "        \"lattice\": {\"dim\": {\"x\": 16, \"y\": 16}, \"velocities\": {\"x\": 4, \"y\": 4}},\n",
+    "        \"lattice\": {\n",
+    "            \"dim\": {\"x\": 64, \"y\": 32},\n",
+    "            \"velocities\": {\n",
+    "                \"x\": 4,\n",
+    "                \"y\": 4,\n",
+    "            },\n",
+    "        },\n",
     "        \"geometry\": [\n",
-    "            {\"shape\": \"cuboid\", \"x\": [9, 12], \"y\": [3, 6], \"boundary\": \"specular\"},\n",
-    "            {\"shape\": \"cuboid\", \"x\": [9, 12], \"y\": [9, 12], \"boundary\": \"bounceback\"},\n",
+    "            {\"shape\": \"cuboid\", \"x\": [10, 13], \"y\": [14, 17], \"boundary\": \"bounceback\"},\n",
     "        ],\n",
     "    }\n",
     ")\n",
     "\n",
     "\n",
-    "output_dir = \"qlbm-output/collisionless-2d-16x16-2-obstacle-mixed-qiskit\"\n",
+    "output_dir = \"qlbm-output/ms-d2q9-64x32-1-obstacle-qiskit\"\n",
     "create_directory_and_parents(output_dir)"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "b4df878f",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "# Number of shots to simulate for each timestep when running the circuit\n",
-    "NUM_SHOTS = 2**12\n",
-    "# Number of timesteps to simulate\n",
-    "NUM_STEPS = 20"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "26ce07ff-55f9-4828-a9b3-4e9c6dc6a856",
+   "id": "da32a65c",
    "metadata": {},
    "outputs": [],
    "source": [
-    "# In the simulation configuration function the user can determine the specifications of the run\n",
     "cfg = SimulationConfig(\n",
-    "    initial_conditions=CollisionlessInitialConditions(lattice),\n",
+    "    initial_conditions=MSInitialConditions(lattice),\n",
     "    algorithm=CQLBM(lattice),\n",
     "    postprocessing=EmptyPrimitive(lattice),\n",
     "    measurement=GridMeasurement(lattice),\n",
@@ -75,13 +73,13 @@
     "    statevector_sampling=True,\n",
     "    execution_backend=AerSimulator(method=\"statevector\"),\n",
     "    sampling_backend=AerSimulator(method=\"statevector\"),\n",
-    ")"
+    ")\n"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "b733d245-c344-46db-88bb-864e4cf07474",
+   "id": "92d4f9ab",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -91,16 +89,30 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "f1c02c70-8244-4a86-8357-9f8c80d2d632",
+   "id": "f2ba1f1f",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Number of shots to simulate for each timestep when running the circuit\n",
+    "NUM_SHOTS = 2**12\n",
+    "\n",
+    "# Number of timesteps to simulate\n",
+    "NUM_STEPS = 20"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "7f26c6d9",
    "metadata": {},
    "outputs": [],
    "source": [
-    "# Create a runner object to simulate the circuit\n",
     "runner = QiskitRunner(\n",
     "    cfg,\n",
     "    lattice,\n",
     ")\n",
     "\n",
+    "\n",
     "# Simulate the circuits using both snapshots\n",
     "runner.run(\n",
     "    NUM_STEPS,  # Number of time steps\n",
@@ -113,7 +125,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "a4c15cbe-3766-4ea1-b3cb-f7caec41328c",
+   "id": "a1737979",
    "metadata": {},
    "outputs": [],
    "source": []
@@ -135,7 +147,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.13.7"
+   "version": "3.13.9"
   }
  },
  "nbformat": 4,