Added unit tests requested by Quantum Motion

quest/include/multiplication.h

@@ -746,23 +746,20 @@ extern "C" {
 
 
 /// @notyetdoced
-/// @notyetvalidated
 /// @see
 /// - leftapplyCompMatr1()
 /// - applyQubitProjector()
 void leftapplyQubitProjector(Qureg qureg, int qubit, int outcome);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 /// @see
 /// - leftapplyCompMatr1()
 /// - applyMultiQubitProjector()
 void leftapplyMultiQubitProjector(Qureg qureg, int* qubits, int* outcomes, int numQubits);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 /// @see
 /// - rightapplyCompMatr1()
 /// - applyQubitProjector()
@@ -782,6 +779,25 @@ void rightapplyMultiQubitProjector(Qureg qureg, int* qubits, int* outcomes, int
 }
 #endif
 
+#ifdef __cplusplus
+
+
+/// @notyetdoced
+/// @cppvectoroverload
+/// @see leftapplyMultiQubitProjector()
+void leftapplyMultiQubitProjector(Qureg qureg, std::vector<int> qubits, std::vector<int> outcomes);
+
+
+/// @notyetdoced
+/// @cppvectoroverload
+/// @see rightapplyMultiQubitProjector()
+void rightapplyMultiQubitProjector(Qureg qureg, std::vector<int> qubits, std::vector<int> outcomes);
+
+
+#endif
+
+/** @} */
+
 
 
 /** 

quest/include/operations.h

@@ -2291,32 +2291,26 @@ extern "C" {
 
 
 /// @notyetdoced
-/// @notyetvalidated
 int applyQubitMeasurement(Qureg qureg, int target);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 int applyQubitMeasurementAndGetProb(Qureg qureg, int target, qreal* probability);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 qreal applyForcedQubitMeasurement(Qureg qureg, int target, int outcome);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 qindex applyMultiQubitMeasurement(Qureg qureg, int* qubits, int numQubits);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 qindex applyMultiQubitMeasurementAndGetProb(Qureg qureg, int* qubits, int numQubits, qreal* probability);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 qreal applyForcedMultiQubitMeasurement(Qureg qureg, int* qubits, int* outcomes, int numQubits);
 
 
@@ -2328,16 +2322,18 @@ qreal applyForcedMultiQubitMeasurement(Qureg qureg, int* qubits, int* outcomes,
 #ifdef __cplusplus
 
 
-/// @notyettested
-/// @notyetvalidated
+/// @notyetdoced
+/// @cppvectoroverload
+/// @see applyMultiQubitMeasurement()
+qindex applyMultiQubitMeasurement(Qureg qureg, std::vector<int> qubits);
+
+
 /// @notyetdoced
 /// @cppvectoroverload
 /// @see applyMultiQubitMeasurementAndGetProb()
 qindex applyMultiQubitMeasurementAndGetProb(Qureg qureg, std::vector<int> qubits, qreal* probability);
 
 
-/// @notyettested
-/// @notyetvalidated
 /// @notyetdoced
 /// @cppvectoroverload
 /// @see applyForcedMultiQubitMeasurement()
@@ -2363,12 +2359,10 @@ extern "C" {
 
 
 /// @notyetdoced
-/// @notyetvalidated
 void applyQubitProjector(Qureg qureg, int target, int outcome);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 void applyMultiQubitProjector(Qureg qureg, int* qubits, int* outcomes, int numQubits);
 
 
@@ -2380,8 +2374,6 @@ void applyMultiQubitProjector(Qureg qureg, int* qubits, int* outcomes, int numQu
 #ifdef __cplusplus
 
 
-/// @notyettested
-/// @notyetvalidated
 /// @notyetdoced
 /// @cppvectoroverload
 /// @see applyMultiQubitProjector()
@@ -2407,12 +2399,10 @@ extern "C" {
 
 
 /// @notyetdoced
-/// @notyetvalidated
 void applyQuantumFourierTransform(Qureg qureg, int* targets, int numTargets);
 
 
 /// @notyetdoced
-/// @notyetvalidated
 void applyFullQuantumFourierTransform(Qureg qureg);
 
 
@@ -2424,8 +2414,6 @@ void applyFullQuantumFourierTransform(Qureg qureg);
 #ifdef __cplusplus
 
 
-/// @notyettested
-/// @notyetvalidated
 /// @notyetdoced
 /// @cppvectoroverload
 /// @see applyQuantumFourierTransform()

quest/include/trotterisation.h

@@ -172,12 +172,17 @@ extern "C" {
  *  - applyTrotterizedUnitaryTimeEvolution()
  * 
  * @author Tyson Jones
+ * @author Vasco Ferreira (randomisation)
  */
 void applyTrotterizedPauliStrSumGadget(Qureg qureg, PauliStrSum sum, qreal angle, int order, int reps, bool permutePaulis);
 
 
 /// @notyetdoced
 /// @notyettested
+///
+/// @author Tyson Jones
+/// @author Vasco Ferreira (randomisation)
+///
 /// @see
 ///  - applyTrotterizedPauliStrSumGadget()
 ///  - applyControlledCompMatr1()
@@ -186,6 +191,10 @@ void applyTrotterizedControlledPauliStrSumGadget(Qureg qureg, int control, Pauli
 
 /// @notyetdoced
 /// @notyettested
+///
+/// @author Tyson Jones
+/// @author Vasco Ferreira (randomisation)
+///
 /// @see
 ///  - applyTrotterizedPauliStrSumGadget()
 ///  - applyMultiControlledCompMatr1()
@@ -194,6 +203,10 @@ void applyTrotterizedMultiControlledPauliStrSumGadget(Qureg qureg, int* controls
 
 /// @notyetdoced
 /// @notyettested
+///
+/// @author Tyson Jones
+/// @author Vasco Ferreira (randomisation)
+///
 /// @see
 ///  - applyTrotterizedPauliStrSumGadget()
 ///  - applyMultiStateControlledCompMatr1()
@@ -248,6 +261,7 @@ void applyTrotterizedMultiStateControlledPauliStrSumGadget(Qureg qureg, int* con
  * - if @p reps is not a positive integer.
  * 
  * @author Tyson Jones
+ * @author Vasco Ferreira (randomisation)
  */
 void applyTrotterizedNonUnitaryPauliStrSumGadget(Qureg qureg, PauliStrSum sum, qcomp angle, int order, int reps, bool permutePaulis);
 
@@ -264,6 +278,10 @@ void applyTrotterizedNonUnitaryPauliStrSumGadget(Qureg qureg, PauliStrSum sum, q
 /// @notyetvalidated
 /// @notyetdoced
 /// @cppvectoroverload
+///
+/// @author Tyson Jones
+/// @author Vasco Ferreira (randomisation)
+///
 /// @see applyTrotterizedMultiControlledPauliStrSumGadget()
 void applyTrotterizedMultiControlledPauliStrSumGadget(Qureg qureg, std::vector<int> controls, PauliStrSum sum, qreal angle, int order, int reps, bool permutePaulis);
 
@@ -272,6 +290,10 @@ void applyTrotterizedMultiControlledPauliStrSumGadget(Qureg qureg, std::vector<i
 /// @notyetvalidated
 /// @notyetdoced
 /// @cppvectoroverload
+///
+/// @author Tyson Jones
+/// @author Vasco Ferreira (randomisation)
+///
 /// @see applyTrotterizedMultiStateControlledPauliStrSumGadget()
 void applyTrotterizedMultiStateControlledPauliStrSumGadget(Qureg qureg, std::vector<int> controls, std::vector<int> states, PauliStrSum sum, qreal angle, int order, int reps, bool permutePaulis);
 
@@ -294,9 +316,7 @@ extern "C" {
 #endif
 
 
-/** @notyettested
- * 
- * Unitarily time evolves @p qureg for the duration @p time under the time-independent Hamiltonian @p hamil, 
+/** Unitarily time evolves @p qureg for the duration @p time under the time-independent Hamiltonian @p hamil, 
  * as approximated by symmetrized Trotterisation of the specified @p order and number of cycles @p reps. 
  * 
  * @formulae 
@@ -385,13 +405,12 @@ extern "C" {
  * - if @p reps is not a positive integer.
  * 
  * @author Tyson Jones
+ * @author Vasco Ferreira (randomisation)
  */
 void applyTrotterizedUnitaryTimeEvolution(Qureg qureg, PauliStrSum hamil, qreal time, int order, int reps, bool permutePaulis);
 
 
-/** @notyettested
- * 
- * Simulates imaginary-time evolution of @p qureg for the duration @p tau under the time-independent 
+/** Simulates imaginary-time evolution of @p qureg for the duration @p tau under the time-independent 
  * Hamiltonian @p hamil, as approximated by symmetrized Trotterisation of the specified @p order and
  * number of cycles @p reps. 
  * 
@@ -517,6 +536,7 @@ void applyTrotterizedUnitaryTimeEvolution(Qureg qureg, PauliStrSum hamil, qreal
  * - if @p reps is not a positive integer.
  * 
  * @author Tyson Jones
+ * @author Vasco Ferreira (randomisation)
  */
 void applyTrotterizedImaginaryTimeEvolution(Qureg qureg, PauliStrSum hamil, qreal tau, int order, int reps, bool permutePaulis);
 
@@ -665,6 +685,7 @@ void applyTrotterizedImaginaryTimeEvolution(Qureg qureg, PauliStrSum hamil, qrea
  * - if @p reps is not a positive integer.
  * 
  * @author Tyson Jones
+ * @author Vasco Ferreira (randomisation)
  */
 void applyTrotterizedNoisyTimeEvolution(Qureg qureg, PauliStrSum hamil, qreal* damps, PauliStrSum* jumps, int numJumps, qreal time, int order, int reps, bool permutePaulis);
 

quest/src/api/operations.cpp

@@ -1717,6 +1717,11 @@ qreal applyForcedMultiQubitMeasurement(Qureg qureg, int* qubits, int* outcomes,
 
 } // end de-mangler
 
+qindex applyMultiQubitMeasurement(Qureg qureg, vector<int> qubits) {
+
+    return applyMultiQubitMeasurement(qureg, qubits.data(), qubits.size());
+}
+
 qindex applyMultiQubitMeasurementAndGetProb(Qureg qureg, vector<int> qubits, qreal* probability) {
 
     return applyMultiQubitMeasurementAndGetProb(qureg, qubits.data(), qubits.size(), probability);

quest/src/core/validation.cpp

@@ -3765,7 +3765,7 @@ void validate_measurementOutcomesMatchTargets(int numQubits, int numOutcomes, co
         {"${NUM_QUBITS}",    numQubits},
         {"${NUM_OUTCOMES}",  numOutcomes}};
 
-    assertThat(numQubits == numOutcomes, report::MEASUREMENT_OUTCOMES_MISMATCH_NUM_TARGETS, caller);
+    assertThat(numQubits == numOutcomes, report::MEASUREMENT_OUTCOMES_MISMATCH_NUM_TARGETS, vars, caller);
 }
 
 

tests/unit/initialisations.cpp

@@ -486,8 +486,7 @@ TEST_CASE( "setQuregToWeightedSum", TEST_CATEGORY ) {
 
     SECTION( LABEL_VALIDATION ) {
 
-        // arbitrary existing qureg
-        Qureg qureg = getCachedStatevecs().begin()->second;
+        Qureg qureg = getArbitraryCachedStatevec();
 
         SECTION( "out qureg uninitialised" ) {
 
@@ -702,7 +701,7 @@ TEST_CASE( "setQuregToMixture", TEST_CATEGORY ) {
 
         SECTION( "out qureg is statevector" ) {
 
-            Qureg badQureg = getCachedStatevecs().begin()->second;
+            Qureg badQureg = getArbitraryCachedStatevec();
 
             REQUIRE_THROWS_WITH( 
                 setQuregToMixture(badQureg, nullptr, nullptr, 1), 
@@ -717,7 +716,7 @@ TEST_CASE( "setQuregToMixture", TEST_CATEGORY ) {
 
             // hide a statevector among them
             int badInd = GENERATE_COPY( range(0,numIn) );
-            inQuregs[badInd] = getCachedStatevecs().begin()->second;;
+            inQuregs[badInd] = getArbitraryCachedStatevec();
 
             REQUIRE_THROWS_WITH( 
                 setQuregToMixture(qureg, nullptr, inQuregs.data(), numIn),