[feature/SPO_apply] Small fixes, include <chrono> in the examples

fast_pauli/cpp/examples/03_summed_pauli_op.cpp

@@ -13,6 +13,7 @@
  */
 
 #include <algorithm>
+#include <chrono>
 #include <experimental/mdspan>
 #include <random>
 
@@ -28,10 +29,10 @@ int main()
     // User settings
     //
     size_t const n_operators = 1000;
-    size_t const n_qubits = 10;
+    size_t const n_qubits = 12;
     size_t const weight = 2;
     size_t const n_states = 1000;
-    using fp_type = double;
+    using fp_type = float;
 
     //
     // Setup the summed pauli operator
@@ -63,11 +64,11 @@ int main()
     //
     // Apply the states
     //
-    auto start_seq = std::chrono::high_resolution_clock::now();
-    summed_op.apply_weighted(std::execution::seq, new_states, states, weights);
-    auto end_seq = std::chrono::high_resolution_clock::now();
-    fmt::println("Time taken for sequential execution: {} seconds",
-                 std::chrono::duration_cast<std::chrono::seconds>(end_seq - start_seq).count());
+    // auto start_seq = std::chrono::high_resolution_clock::now();
+    // summed_op.apply_weighted(std::execution::seq, new_states, states, weights);
+    // auto end_seq = std::chrono::high_resolution_clock::now();
+    // fmt::println("Time taken for sequential execution: {} seconds",
+    //              std::chrono::duration_cast<std::chrono::seconds>(end_seq - start_seq).count());
 
     auto start_par = std::chrono::high_resolution_clock::now();
     summed_op.apply_weighted(std::execution::par, new_states, states, weights);

fast_pauli/cpp/include/__pauli_helpers.hpp

@@ -112,8 +112,8 @@ std::vector<PauliString> calcutate_pauli_strings(size_t const n_qubits, size_t c
     size_t n_pauli_strings = nontrivial_paulis.size() * idx_combos.size();
     std::vector<PauliString> result(n_pauli_strings);
 
-    fmt::println("n_qubits = {}  weight = {}  n_nontrivial_paulis = {}  n_combos = {}", n_qubits, weight,
-                 nontrivial_paulis.size(), idx_combos.size());
+    // fmt::println("n_qubits = {}  weight = {}  n_nontrivial_paulis = {}  n_combos = {}", n_qubits, weight,
+    //              nontrivial_paulis.size(), idx_combos.size());
 
     // Iterate through all the nontrivial paulis and all the combinations
     for (size_t i = 0; i < nontrivial_paulis.size(); ++i)
@@ -153,7 +153,7 @@ std::vector<PauliString> calculate_pauli_strings_max_weight(size_t n_qubits, siz
         result.insert(result.end(), ps.begin(), ps.end());
     }
 
-    fmt::println("n_qubits = {}  weight = {}  n_pauli_strings = {}", n_qubits, weight, result.size());
+    // fmt::println("n_qubits = {}  weight = {}  n_pauli_strings = {}", n_qubits, weight, result.size());
     return result;
 }