updated notebooks after change

Author: TomerGoldfriend

applications/cfd/qls_for_hybrid_solvers/banded_be.py

@@ -223,3 +223,47 @@ def block_encode_banded_sym(
             ),
         ),
     )
+
+
+def get_banded_diags_be(mat_raw_scr):
+    """
+    Get relevant block-encoding properties for `block_encode_banded` block encoding,
+
+    Parameters
+    ----------
+    mat_raw_scr : scipy.sparse.spmatrix
+        Square sparse matrix of shape (N, N), real or complex, to be block-encoded.
+
+    Returns
+    -------
+    data_size : int
+       Size of the data variable.
+    block_size : int
+        Size of the block variable.
+    be_scaling_factor : float
+        The scaling factor of the block-encoding unitary
+    BlockEncodedState : QStruct
+        QSVT-compatible QStruct holding the quantum variables, with fields:
+          - data  : QNum[data_size]
+          - block : QNum[block_size]
+    be_qfunc : qfunc
+        Quantum function that implements the block encoding. Signature:
+        be_qfunc(be: BlockEncodedState) → None
+    """
+    raw_size = mat_raw_scr.shape[0]
+    data_size = max(1, (raw_size - 1).bit_length())
+    offsets, diags, diags_maxima, prepare_norm = get_be_banded_data(mat_raw_scr)
+    block_size = int(np.ceil(np.log2(len(offsets)))) + 1
+    be_scaling_factor = prepare_norm
+
+    @qfunc
+    def be_qfunc(block: QNum, data: QNum):
+        block_encode_banded(
+            offsets=offsets,
+            diags=diags,
+            prep_diag=diags_maxima,
+            block=block,
+            data=data,
+        )
+
+    return data_size, block_size, be_scaling_factor, be_qfunc

applications/cfd/qls_for_hybrid_solvers/banded_be.synthesis_options.json

@@ -6,28 +6,28 @@
   "preferences": {
     "custom_hardware_settings": {
       "basis_gates": [
-        "sx",
-        "p",
-        "tdg",
+        "z",
+        "ry",
+        "h",
         "x",
-        "cz",
-        "u2",
+        "u1",
         "id",
         "sxdg",
-        "y",
-        "cy",
-        "s",
-        "ry",
-        "cx",
-        "t",
-        "r",
-        "u",
+        "sx",
         "rz",
-        "u1",
-        "z",
         "rx",
         "sdg",
-        "h"
+        "tdg",
+        "cy",
+        "cz",
+        "cx",
+        "u",
+        "y",
+        "s",
+        "r",
+        "t",
+        "u2",
+        "p"
       ],
       "is_symmetric_connectivity": true
     },
@@ -36,7 +36,7 @@
     "optimization_level": 1,
     "output_format": ["qasm"],
     "pretty_qasm": true,
-    "random_seed": 2454004842,
+    "random_seed": 2459810836,
     "synthesize_all_separately": false,
     "timeout_seconds": 300,
     "transpilation_option": "auto optimize"

applications/cfd/qls_for_hybrid_solvers/banded_sym_be.synthesis_options.json

@@ -6,28 +6,28 @@
   "preferences": {
     "custom_hardware_settings": {
       "basis_gates": [
-        "sx",
-        "p",
-        "tdg",
+        "z",
+        "ry",
+        "h",
         "x",
-        "cz",
-        "u2",
+        "u1",
         "id",
         "sxdg",
-        "y",
-        "cy",
-        "s",
-        "ry",
-        "cx",
-        "t",
-        "r",
-        "u",
+        "sx",
         "rz",
-        "u1",
-        "z",
         "rx",
         "sdg",
-        "h"
+        "tdg",
+        "cy",
+        "cz",
+        "cx",
+        "u",
+        "y",
+        "s",
+        "r",
+        "t",
+        "u2",
+        "p"
       ],
       "is_symmetric_connectivity": true
     },
@@ -36,7 +36,7 @@
     "optimization_level": 1,
     "output_format": ["qasm"],
     "pretty_qasm": true,
-    "random_seed": 3750507641,
+    "random_seed": 1840897986,
     "synthesize_all_separately": false,
     "timeout_seconds": 300,
     "transpilation_option": "auto optimize"