wrote doc for unit-testing functions

doxyconf

@@ -757,7 +757,7 @@ WARN_LOGFILE           =
 # spaces.
 # Note: If this tag is empty the current directory is searched.
 
-INPUT                  = . README.md ./QuEST/src ./QuEST/src/CPU ./QuEST/src/GPU ./QuEST/include
+INPUT                  = . README.md ./QuEST/src ./QuEST/src/CPU ./QuEST/src/GPU ./QuEST/include ./testing
 
 # This tag can be used to specify the character encoding of the source files
 # that doxygen parses. Internally doxygen uses the UTF-8 encoding. Doxygen uses
@@ -779,7 +779,7 @@ INPUT_ENCODING         = UTF-8
 
 EXTENSION_MAPPING      = cu=c++ cpp=c++ c=c h=c md=markdown
 
-FILE_PATTERNS          = *.cu *.c *.cpp *.h *.md
+FILE_PATTERNS          = *.cu *.c *.cpp *.h *.hpp *.md
 
 # The RECURSIVE tag can be used to specify whether or not subdirectories should
 # be searched for input files as well.

testing/QuEST_test_calculations.cpp

@@ -34,6 +34,7 @@
 using Catch::Matchers::Contains;
 
 
+
 TEST_CASE( "calcDensityInnerProduct", "[calculations]" ) {
 
     QuESTEnv env = createQuESTEnv();

testing/QuEST_test_unitaries.cpp

@@ -9,30 +9,6 @@
  *
  * QuEST's user validation handling is unit tested by redefining exitWithError 
  * (a weak C symbol) to throw a C++ exception, caught by the Catch2 library.
- * 
- * All unitary unit tests follow the master template:
-
-TEST_CASE( "OP", "[unitaries]" ) {
-    
-    PREPARE_TEST( env, quregVec, quregMatr, refVec, refMatr, NUM_QUBITS );
- 
-    SECTION( "state-vector correctness" ) {
-    
-    }
-    
-    SECTION( "density-matrix correctness" ) {
-    
-    }
-    
-    SECTION( "input validation" ) {
-    
-    }
-    CLEANUP_TEST( env, quregVec, quregMatr );
-}
-    
-    CLEANUP_TEST( env, quregVec, quregMatr );
-}
-
  *
  * @author Tyson Jones
  */
@@ -1416,7 +1392,7 @@ TEST_CASE( "multiRotatePauli", "[unitaries]" ) {
         // produces exp(-i param/2 pauliProd), unless pauliProd = I
         QMatrix op;
         if (numRefTargs > 0)
-            op = getExponentialPauliMatrix(param, pauliProd);
+            op = getExponentialOfPauliMatrix(param, pauliProd);
 
         SECTION( "state-vector" ) {
 
@@ -1490,7 +1466,7 @@ TEST_CASE( "multiRotateZ", "[unitaries]" ) {
             getFullOperatorMatrix(NULL, 0, targs, numTargs, zProd, NUM_QUBITS);
 
         // exp( -i param/2 Z . I . Z ...)
-        QMatrix op = getExponentialDiagonalMatrix(expArg);
+        QMatrix op = getExponentialOfDiagonalMatrix(expArg);
 
         // all qubits to specify full operator matrix on reference structures
         int allQubits[NUM_QUBITS];