updated dynamics examples to use time-evol

in lieu of direct gadgets, for clarity

Author: TysonRayJones

examples/extended/dynamics.c

@@ -1,6 +1,5 @@
 /** @file
- * An example of using QuEST (primarily function
- * applyTrotterizedPauliStrSumGadget()) to perform
+ * An example of using QuEST to perform closed
  * dynamical simulation via Trotterisation of the
  * unitary-time evolution operator.
  * 
@@ -158,7 +157,7 @@ int main() {
     for (int i=0; i<steps; i++) {
 
         // evolve qureg under (approx) exp(-i dt H)
-        applyTrotterizedPauliStrSumGadget(qureg, hamil, -dt, order, reps);
+        applyTrotterizedUnitaryTimeEvolution(qureg, hamil, dt, order, reps);
 
         // calculate and report <O>
         qreal time = dt * (i+1);
@@ -188,7 +187,7 @@ int main() {
 
     // verify results by uninterrupted higher-order simulation to target time
     initPlusState(qureg);
-    applyTrotterizedPauliStrSumGadget(qureg, hamil, -dt*steps, order+2, reps*steps);
+    applyTrotterizedUnitaryTimeEvolution(qureg, hamil, dt*steps, order+2, reps*steps);
     reportScalar("final <O>", calcExpecPauliStrSum(qureg, observ));
 
     // clean up

examples/extended/dynamics.cpp

@@ -1,6 +1,5 @@
 /** @file
- * An example of using QuEST (primarily function
- * applyTrotterizedPauliStrSumGadget()) to perform
+ * An example of using QuEST to perform closed
  * dynamical simulation via Trotterisation of the
  * unitary-time evolution operator.
  * 
@@ -155,7 +154,7 @@ int main() {
     for (int i=0; i<steps; i++) {
 
         // evolve qureg under (approx) exp(-i dt H)
-        applyTrotterizedPauliStrSumGadget(qureg, hamil, -dt, order, reps);
+        applyTrotterizedUnitaryTimeEvolution(qureg, hamil, dt, order, reps);
 
         // calculate and report <O>
         qreal time = dt * (i+1);
@@ -182,7 +181,7 @@ int main() {
 
     // verify results by uninterrupted higher-order simulation to target time
     initPlusState(qureg);
-    applyTrotterizedPauliStrSumGadget(qureg, hamil, -dt*steps, order+2, reps*steps);
+    applyTrotterizedUnitaryTimeEvolution(qureg, hamil, dt*steps, order+2, reps*steps);
     reportScalar("final <O>", calcExpecPauliStrSum(qureg, observ));
 
     // clean up