Implement sparse isometry with binary encoding

docs/source/_static/examples/python/circuit.py

@@ -104,7 +104,7 @@
 _, wfn_cas = create("multi_configuration_calculator").run(ham, 1, 1)
 
 # StatePreparation produces a Circuit with a native Q# factory
-state_prep = create("state_prep", "sparse_isometry_gf2x")
+state_prep = create("state_prep", "sparse_isometry")
 circuit = state_prep.run(wfn_cas)
 
 # Inspect the Q# circuit (prune unused qubits for clarity)

docs/source/_static/examples/python/phase_estimation.py

@@ -66,7 +66,7 @@
 qubit_ham = qubit_mapper.run(hamiltonian)
 
 # 6. State preparation
-state_prep = create("state_prep", "sparse_isometry_gf2x")
+state_prep = create("state_prep", "sparse_isometry")
 circuit = state_prep.run(wfn_cas)
 
 # 7. Create and run IQPE with nested algorithm settings

docs/source/_static/examples/python/quickstart.py

@@ -94,7 +94,7 @@
 # start-state-prep-circuit
 
 # Generate state preparation circuit for the sparse state via sparse isometry (GF2 + X)
-state_prep = create("state_prep", "sparse_isometry_gf2x")
+state_prep = create("state_prep", "sparse_isometry")
 sparse_isometry_circuit = state_prep.run(wfn_sparse)
 
 # end-state-prep-circuit

docs/source/_static/examples/python/release_notes_v1_1.py

@@ -199,7 +199,7 @@
 _cas = create("multi_configuration_calculator", "macis_cas")
 _E_cas, _wfn_cas = _cas.run(_hamiltonian, 1, 1)
 
-_state_prep = create("state_prep", "sparse_isometry_gf2x")
+_state_prep = create("state_prep", "sparse_isometry")
 _circuit = _state_prep.run(_wfn_cas)
 
 circuit_executor = create("circuit_executor", "qdk_sparse_state_simulator")

docs/source/_static/examples/python/state_preparation.py

@@ -10,7 +10,7 @@
 from qdk_chemistry.algorithms import create
 
 # Create a StatePreparation instance
-sparse_prep = create("state_prep", "sparse_isometry_gf2x")
+sparse_prep = create("state_prep", "sparse_isometry")
 regular_prep = create("state_prep", "qiskit_regular_isometry")
 # end-cell-create
 ################################################################################
@@ -61,6 +61,6 @@
 from qdk_chemistry.algorithms import registry
 
 print(registry.available("state_prep"))
-# ['sparse_isometry_gf2x', 'qiskit_regular_isometry']
+# ['sparse_isometry', 'qiskit_regular_isometry']
 # end-cell-list-implementations
 ################################################################################

docs/source/conf.py

@@ -59,13 +59,13 @@
     "sphinx.ext.autosummary",  # Create summary tables for modules/classes
     "sphinx.ext.intersphinx",  # Link to other projects' documentation
     "sphinx.ext.viewcode",  # Add links to view source code
+    "sphinx.ext.napoleon",  # Support for Google-style and NumPy-style docstrings
     # Additional extensions
     "sphinx_autodoc_typehints",  # Better support for Python type annotations
     "sphinx_inline_tabs",  # Support for tabbed content in docs
     # C++ documentation
     "breathe",  # Bridge between Sphinx and Doxygen
     # Enable Google-style docstrings parsing
-    "sphinx.ext.napoleon",  # Support for Google-style and NumPy-style docstrings
     "sphinx.ext.todo",  # Support for listing to-dos
     "sphinx.ext.graphviz",  # Support for Graphviz diagrams
     "sphinxcontrib.bibtex",  # Support for bibliographic references

docs/source/user/comprehensive/algorithms/state_preparation.rst

@@ -78,7 +78,7 @@ You can discover available implementations programmatically:
 Sparse Isometry GF2+X
 ~~~~~~~~~~~~~~~~~~~~~
 
-.. rubric:: Factory name: ``"sparse_isometry_gf2x"``
+.. rubric:: Factory name: ``"sparse_isometry"``
 
 This method is an optimized approach that leverages sparsity in the target wavefunction. The GF2+X method, a modification of the original sparse isometry work in :cite:`Malvetti2021`, applies GF(2) Gaussian elimination to the binary matrix representation of the state to determine a reduced space representation of the sparse state. This reduced state is then densely encoded via regular isometry :cite:`Christandl2016` on a smaller number of qubits, and finally scattered to the full qubit space using X and :term:`CNOT` gates. These reductions correspond to efficient gate sequences that simplify the preparation basis. By focusing only on non-zero amplitudes, this approach substantially reduces circuit depth and gate count compared with dense isometry methods. This method is native to QDK/Chemistry and is especially efficient for wavefunctions with sparse amplitude structure.
 

examples/extended_hubbard.ipynb

@@ -502,7 +502,7 @@
     "from qdk.widgets import Circuit\n",
     "\n",
     "# Generate state preparation circuit for the sparse state via GF2+X sparse isometry\n",
-    "state_prep = create(\"state_prep\", \"sparse_isometry_gf2x\")\n",
+    "state_prep = create(\"state_prep\", \"sparse_isometry\")\n",
     "state_prep_circuit = state_prep.run(wfn_trial)\n",
     "\n",
     "# Visualize the sparse isometry circuit\n",

examples/factory_list.ipynb

@@ -1134,28 +1134,28 @@
        "</tr>\n",
        "<tr>\n",
        "<td>state_prep</td>\n",
-       "<td>sparse_isometry_gf2x</td>\n",
+       "<td>sparse_isometry</td>\n",
        "<td>State preparation using sparse isometry with enhanced GF2+X elimination.</td>\n",
        "<td>basis_gates</td>\n",
        "<td>[&#x27;x&#x27;, &#x27;y&#x27;, &#x27;z&#x27;, &#x27;cx&#x27;, &#x27;cz&#x27;, &#x27;id&#x27;, &#x27;h&#x27;, &#x27;s&#x27;, &#x27;sdg&#x27;, &#x27;rz&#x27;]</td>\n",
        "</tr>\n",
        "<tr>\n",
        "<td>state_prep</td>\n",
-       "<td>sparse_isometry_gf2x</td>\n",
+       "<td>sparse_isometry</td>\n",
        "<td>State preparation using sparse isometry with enhanced GF2+X elimination.</td>\n",
        "<td>dense_preparation_method</td>\n",
        "<td>qdk</td>\n",
        "</tr>\n",
        "<tr>\n",
        "<td>state_prep</td>\n",
-       "<td>sparse_isometry_gf2x</td>\n",
+       "<td>sparse_isometry</td>\n",
        "<td>State preparation using sparse isometry with enhanced GF2+X elimination.</td>\n",
        "<td>transpile</td>\n",
        "<td>True</td>\n",
        "</tr>\n",
        "<tr>\n",
        "<td>state_prep</td>\n",
-       "<td>sparse_isometry_gf2x</td>\n",
+       "<td>sparse_isometry</td>\n",
        "<td>State preparation using sparse isometry with enhanced GF2+X elimination.</td>\n",
        "<td>transpile_optimization_level</td>\n",
        "<td>0</td>\n",

examples/interoperability/qiskit/iqpe_no_trotter.py

@@ -228,7 +228,7 @@ def run_iterative_exact_qpe(
     active_hamiltonian, list(top_configurations.keys())
 )
 
-sparse_state_prep = create("state_prep", algorithm_name="sparse_isometry_gf2x")
+sparse_state_prep = create("state_prep", algorithm_name="sparse_isometry")
 state_prep = sparse_state_prep.run(sparse_wavefunction).get_qiskit_circuit()
 state_prep = transpile(
     state_prep,

examples/interoperability/qiskit/iqpe_trotter.py

@@ -94,7 +94,7 @@
     active_hamiltonian, list(top_configurations.keys())
 )
 
-sparse_state_prep = create("state_prep", algorithm_name="sparse_isometry_gf2x")
+sparse_state_prep = create("state_prep", algorithm_name="sparse_isometry")
 state_prep = sparse_state_prep.run(sparse_wavefunction).get_qiskit_circuit()
 state_prep = transpile(
     state_prep,

examples/qpe_stretched_n2.ipynb

@@ -359,7 +359,7 @@
     "from qdk.widgets import Circuit\n",
     "\n",
     "# Generate state preparation circuit for the sparse state via GF2+X sparse isometry\n",
-    "state_prep = create(\"state_prep\", \"sparse_isometry_gf2x\")\n",
+    "state_prep = create(\"state_prep\", \"sparse_isometry\")\n",
     "sparse_isometry_circuit = state_prep.run(wfn_trial)\n",
     "\n",
     "# Visualize the sparse isometry circuit, idle and classical qubits are removed\n",