#include "QuEST.h"
#include "utilities.hpp"
#include "catch.hpp"
#include <random>
#include <algorithm>
#include <bitset>

Data Structures
class	SequenceGenerator< T >

class	SubListGenerator

Macros
#define	DEMAND(cond) if (!(cond)) FAIL( );

#define	macro_copyComplexMatrix(dest, src)
	Copies ComplexMatrix structures into a QMatrix. More...

#define	macro_copyQMatrix(dest, src)

Functions
void	applyReferenceMatrix (QMatrix &state, int ctrls, int numCtrls, int targs, int numTargs, QMatrix op)
	Modifies the density matrix `state` to be the result of left-multiplying the multi-target operator matrix `op`, with the specified control and target qubits (in `ctrls` and `targs` respectively). More...

void	applyReferenceMatrix (QVector &state, int ctrls, int numCtrls, int targs, int numTargs, QMatrix op)
	Modifies the state-vector `state` to be the result of left-multiplying the multi-target operator matrix `op`, with the specified control and target qubits (in `ctrls` and `targs` respectively). More...

void	applyReferenceOp (QMatrix &state, int ctrls, int numCtrls, int targs, int numTargs, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the multi-target operator matrix `op`, with the specified control and target qubits (in `ctrls` and `targs` respectively). More...

void	applyReferenceOp (QMatrix &state, int *ctrls, int numCtrls, int targ1, int targ2, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the two-target operator matrix `op`, with the specified control qubits (in `ctrls`). More...

void	applyReferenceOp (QMatrix &state, int *ctrls, int numCtrls, int target, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the single-target operator matrix `op`, with the specified control qubits (in `ctrls`). More...

void	applyReferenceOp (QMatrix &state, int *targs, int numTargs, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the multi-target operator matrix `op`, with no control qubits. More...

void	applyReferenceOp (QMatrix &state, int ctrl, int *targs, int numTargs, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the multi-target operator matrix `op`, with a single control qubit `ctrl`. More...

void	applyReferenceOp (QMatrix &state, int ctrl, int targ, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the single-control single-target operator matrix `op`. More...

void	applyReferenceOp (QMatrix &state, int ctrl, int targ1, int targ2, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the two-target operator matrix `op`, with a single control qubit `ctrl`. More...

void	applyReferenceOp (QMatrix &state, int targ, QMatrix op)
	Modifies the density matrix `state` to be the result of applying the single-target operator matrix `op`, with no control qubit. More...

void	applyReferenceOp (QVector &state, int ctrls, int numCtrls, int targs, int numTargs, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the multi-target operator matrix `op`, with the specified control and target qubits (in `ctrls` and `targs` respectively). More...

void	applyReferenceOp (QVector &state, int *ctrls, int numCtrls, int targ1, int targ2, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the two-target operator matrix `op`, with the specified control qubits (in `ctrls`). More...

void	applyReferenceOp (QVector &state, int *ctrls, int numCtrls, int target, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the single-target operator matrix `op`, with the specified control qubits (in `ctrls`). More...

void	applyReferenceOp (QVector &state, int *targs, int numTargs, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the multi-target operator matrix `op`, with no contorl qubits. More...

void	applyReferenceOp (QVector &state, int ctrl, int *targs, int numTargs, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the multi-target operator matrix `op`, with a single control qubit (`ctrl`) This updates `state` under. More...

void	applyReferenceOp (QVector &state, int ctrl, int targ, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the single-target operator matrix `op`, with a single control qubit (`ctrl`). More...

void	applyReferenceOp (QVector &state, int ctrl, int targ1, int targ2, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the two-target operator matrix `op`, with a single control qubit (`ctrl`). More...

void	applyReferenceOp (QVector &state, int targ, QMatrix op)
	Modifies the state-vector `state` to be the result of applying the single-target operator matrix `op`, with no contorl qubits. More...

bool	areEqual (QMatrix a, QMatrix b)
	Returns true if the absolute value of the difference between every amplitude in matrices `a` and `b` is less than `REAL_EPS`. More...

bool	areEqual (Qureg qureg, QMatrix matr)
	Performs a hardware-agnostic comparison of density-matrix `qureg` to `matr`, checking whether the difference between the real and imaginary components of every amplitude is smaller than the QuEST_PREC-specific REAL_EPS (defined in QuEST_precision) precision. More...

bool	areEqual (Qureg qureg, QMatrix matr, qreal precision)
	Performs a hardware-agnostic comparison of density-matrix `qureg` to `matr`, checking whether the difference between the real and imaginary components of every amplitude is smaller than `precision`. More...

bool	areEqual (Qureg qureg, QVector vec)
	Performs a hardware-agnostic comparison of state-vector `qureg` to `vec`, checking whether the difference between the real and imaginary components of every amplitude is smaller than the QuEST_PREC-specific REAL_EPS (defined in QuEST_precision) precision. More...

bool	areEqual (Qureg qureg, QVector vec, qreal precision)
	Performs a hardware-agnostic comparison of state-vector `qureg` to `vec`, checking whether the difference between the real and imaginary components of every amplitude is smaller than `precision`. More...

bool	areEqual (Qureg qureg1, Qureg qureg2)
	Performs a hardware-agnostic comparison of the given quregs, checking whether the difference between the real and imaginary components of every amplitude is smaller than the QuEST_PREC-specific REAL_EPS (defined in QuEST_precision) precision. More...

bool	areEqual (Qureg qureg1, Qureg qureg2, qreal precision)
	Performs a hardware-agnostic comparison of the given quregs, checking whether the difference between the real and imaginary components of every amplitude is smaller than `precision`. More...

bool	areEqual (QVector a, QVector b)
	Returns true if the absolute value of the difference between every amplitude in vectors `a` and `b` is less than `REAL_EPS`. More...

bool	areEqual (QVector vec, qreal *reals)
	Returns true if the absolute value of the difference between every element in `vec` (which must be strictly real) and those implied by `reals`, is less than `REAL_EPS`. More...

bool	areEqual (QVector vec, qreal reals, qreal imags)
	Returns true if the absolute value of the difference between every element in `vec` and those implied by `reals` and `imags`, is less than `REAL_EPS`. More...

Catch::Generators::GeneratorWrapper< int * >	bitsets (int numBits)
	Returns a Catch2 generator of every `numBits-length` bit-set, in increasing lexographic order, where left-most (zero index) bit is treated as LEAST significant (opposite typical convention). More...

unsigned int	calcLog2 (long unsigned int res)
	Returns log2 of numbers which must be gauranteed to be 2^n. More...

void	deleteFilesWithPrefixSynch (char *prefix)
	Deletes all files with filename starting with prefix. More...

qcomp	expI (qreal phase)
	Returns the unit-norm complex number exp(i*`phase`). More...

QMatrix	getConjugateTranspose (QMatrix a)
	Returns the conjugate transpose of the complex square matrix `a`. More...

QVector	getDFT (QVector in)
	Returns the discrete fourier transform of vector in. More...

QVector	getDFT (QVector in, int *targs, int numTargs)
	Returns the discrete fourier transform of a sub-partition of the vector in. More...

QMatrix	getExponentialOfDiagonalMatrix (QMatrix a)
	Returns the matrix exponential of a diagonal, square, complex matrix. More...

QMatrix	getExponentialOfPauliMatrix (qreal angle, QMatrix a)
	Returns the matrix exponential of a kronecker product of pauli matrices (or of any involutory matrices), with exponent factor (-i `angle` / 2). More...

QMatrix	getFullOperatorMatrix (int ctrls, int numCtrls, int targs, int numTargs, QMatrix op, int numQubits)
	Takes a 2^`numTargs-by-2^numTargs matrix` `op` and a returns a 2^`numQubits-by-2^numQubits matrix` where `op` is controlled on the given `ctrls` qubits. More...

QMatrix	getIdentityMatrix (size_t dim)
	Returns a dim-by-dim identity matrix. More...

long long int	getIndexOfTargetValues (long long int ref, int *targs, int numTargs, int targVal)

QMatrix	getKetBra (QVector ket, QVector bra)
	Returns the matrix \|`ket><bra\|,` with ith-jth element `ket(i)` conj(`bra(j)`), since \|`ket><bra\| =` sum_i a_i\|i> sum_j b_j* <j\| = sum_{ij} a_i b_j* \|i><j\|. More...

QMatrix	getKroneckerProduct (QMatrix a, QMatrix b)
	Returns the kronecker product of `a` and `b`, where `a` and `b` are square but possibly differently-sized complex matrices. More...

QVector	getKroneckerProduct (QVector b, QVector a)
	Returns b (otimes) a. More...

QVector	getMatrixDiagonal (QMatrix matr)
	Returns the diagonal vector of the given matrix. More...

QMatrix	getMixedDensityMatrix (std::vector< qreal > probs, std::vector< QVector > states)
	Returns a mixed density matrix formed from mixing the given pure states, which are assumed normalised, but not necessarily orthogonal. More...

QVector	getNormalised (QVector vec)
	Returns an L2-normalised copy of `vec`, using Kahan summation for improved accuracy. More...

QMatrix	getPureDensityMatrix (QVector state)
	Returns a density matrix initialised into the given pure state. More...

qcomp	getRandomComplex ()
	Returns a random complex number within the square closing (-1-i) and (1+i), from a distribution uniformly randomising the individual real and imaginary components in their domains. More...

QMatrix	getRandomDensityMatrix (int numQb)
	Returns a random `numQb-by-numQb density` matrix, from an undisclosed distribution, in a very mixed state. More...

int	getRandomInt (int min, int max)
	Returns a random integer between `min` (inclusive) and `max` (exclusive), from the uniform distribution. More...

std::vector< QMatrix >	getRandomKrausMap (int numQb, int numOps)
	Returns a random Kraus map of #`numOps` 2^`numQb-by-2^numQb operators`, from an undisclosed distribution. More...

std::vector< QVector >	getRandomOrthonormalVectors (int numQb, int numStates)
	Returns a list of random orthonormal complex vectors, from an undisclosed distribution. More...

std::vector< qreal >	getRandomProbabilities (int numProbs)
	Returns a list of random real scalars, each in [0, 1], which sum to unity. More...

QMatrix	getRandomPureDensityMatrix (int numQb)
	Returns a random `numQb-by-numQb density` matrix, from an undisclosed distribution, which is pure (corresponds to a random state-vector) More...

QMatrix	getRandomQMatrix (int dim)
	Returns a `dim-by-dim complex` matrix, where the real and imaginary value of each element are independently random, under the standard normal distribution (mean 0, standard deviation 1). More...

QVector	getRandomQVector (int dim)
	Returns a `dim-length` vector with random complex amplitudes in the square joining {-1-i, 1+i}, of an undisclosed distribution. More...

qreal	getRandomReal (qreal min, qreal max)
	Returns a random real between `min` (inclusive) and `max` (exclusive), from the uniform distribution. More...

QVector	getRandomStateVector (int numQb)
	Returns a random `numQb-length` L2-normalised state-vector from an undisclosed distribution. More...

QMatrix	getRandomUnitary (int numQb)
	Returns a uniformly random (under Haar) 2^`numQb-by-2^numQb unitary` matrix. More...

QMatrix	getSwapMatrix (int qb1, int qb2, int numQb)
	Returns the 2^`numQb-by-2^numQb unitary` matrix which swaps qubits `qb1` and `qb2`; the SWAP gate of not-necessarily-adjacent qubits. More...

long long int	getTwosComplement (long long int decimal, int numBits)
	Returns the two's complement signed encoding of the unsigned number decimal, which must be a number between 0 and 2^numBits (exclusive). More...

long long int	getUnsigned (long long int twosComp, int numBits)
	Return the unsigned value of a number, made of `#numBits` bits, which under two's complement, encodes the signed number twosComp. More...

long long int	getValueOfTargets (long long int ind, int *targs, int numTargs)
	Returns the integer value of the targeted sub-register for the given full state index `ind`. More...

QMatrix	getZeroMatrix (size_t dim)
	Returns a dim-by-dim square complex matrix, initialised to all zeroes. More...

QMatrix	operator* (const qcomp &a, const QMatrix &m)

QVector	operator* (const qcomp &a, const QVector &v)

QMatrix	operator* (const QMatrix &m, const qcomp &a)

QVector	operator* (const QMatrix &m, const QVector &v)

QMatrix	operator* (const QMatrix &m1, const QMatrix &m2)

QVector	operator* (const QVector &v, const qcomp &a)

qcomp	operator* (const QVector &v1, const QVector &v2)

void	operator*= (QMatrix &m1, const QMatrix &m2)

void	operator*= (QMatrix &m1, const qreal &a)

void	operator*= (QVector &v1, const qcomp &a)

QMatrix	operator+ (const QMatrix &m1, const QMatrix &m2)

QVector	operator+ (const QVector &v1, const QVector &v2)

void	operator+= (QMatrix &m1, const QMatrix &m2)

void	operator+= (QVector &v1, const QVector &v2)

QMatrix	operator- (const QMatrix &m1, const QMatrix &m2)

QVector	operator- (const QVector &v1, const QVector &v2)

void	operator-= (QMatrix &m1, const QMatrix &m2)

void	operator-= (QVector &v1, const QVector &v2)

QMatrix	operator/ (const QMatrix &m, const qcomp &a)

QVector	operator/ (const QVector &v, const qcomp &a)

void	operator/= (QMatrix &m1, const qreal &a)

void	operator/= (QVector &v1, const qcomp &a)

Catch::Generators::GeneratorWrapper< pauliOpType * >	pauliseqs (int numPaulis)
	Returns a Catch2 generator of every `numPaulis-length` set of Pauli-matrix types (or base-4 integers). More...

Catch::Generators::GeneratorWrapper< int * >	sequences (int base, int numDigits)
	Returns a Catch2 generator of every `numDigits-length` sequence in the given `base`, in increasing lexographic order, where left-most (zero index) bit is treated as LEAST significant (opposite typical convention). More...

long long int	setBit (long long int num, int bitInd, int bitVal)

void	setDiagMatrixOverrides (QMatrix &matr, int numQubitsPerReg, int numRegs, enum bitEncoding encoding, long long int overrideInds, qreal *overridePhases, int numOverrides)
	Modifies the given diagonal matrix such that the diagonal elements which correspond to the coordinates in overrideInds are replaced with exp(i phase), as prescribed by overridePhases. More...

void	setRandomDiagPauliHamil (PauliHamil hamil)
	Populates `hamil` with random coefficients and a random amount number of PAULI_I and PAULI_Z operators. More...

void	setRandomPauliSum (PauliHamil hamil)
	Populates `hamil` with random coefficients and pauli codes. More...

void	setRandomPauliSum (qreal coeffs, pauliOpType codes, int numQubits, int numTerms)
	Populates the `coeffs` array with random qreals in (-5, 5), and populates `codes` with random Pauli codes. More...

void	setSubMatrix (QMatrix &dest, QMatrix sub, size_t r, size_t c)
	Modifies `dest` by overwriting its submatrix (from top-left corner (`r`, `c`) to bottom-right corner (`r` + `dest.size()`, `c` + `dest.size()`) with the complete elements of sub. More...

void	setUniqueFilename (char outFn, char prefix)
	Modifies outFn to be a filename of format prefix_NUM.txt where NUM is a new unique integer so far. More...

Catch::Generators::GeneratorWrapper< int * >	sublists (Catch::Generators::GeneratorWrapper< int > &&gen, int numSamps)

Catch::Generators::GeneratorWrapper< int * >	sublists (Catch::Generators::GeneratorWrapper< int > &&gen, int numSamps, const int *exclude, int numExclude)

Catch::Generators::GeneratorWrapper< int * >	sublists (Catch::Generators::GeneratorWrapper< int > &&gen, int numSamps, int excluded)

Catch::Generators::GeneratorWrapper< int * >	sublists (int *list, int len, int sublen)
	Returns a Catch2 generator of every length-`sublen` sublist of length-`len` `list`, in increasing lexographic order. More...

ComplexMatrix2	toComplexMatrix2 (QMatrix qm)
	Returns a `ComplexMatrix2` copy of QMatix `qm`. More...

ComplexMatrix4	toComplexMatrix4 (QMatrix qm)
	Returns a `ComplexMatrix4` copy of QMatix `qm`. More...

void	toComplexMatrixN (QMatrix qm, ComplexMatrixN cm)
	Initialises `cm` with the values of `qm`. More...

QMatrix	toDiagonalQMatrix (QVector vec)
	Returns a diagonal complex matrix formed by the given vector. More...

QMatrix	toQMatrix (Complex alpha, Complex beta)
	Returns the matrix (where a=`alpha`, b=`beta`) {{a, -conj(b)}, {b, conj(a)}} using the `qcomp` complex type. More...

QMatrix	toQMatrix (ComplexMatrix2 src)
	Returns a copy of the given 2-by-2 matrix. More...

QMatrix	toQMatrix (ComplexMatrix4 src)
	Returns a copy of the given 4-by-4 matrix. More...

QMatrix	toQMatrix (ComplexMatrixN src)
	Returns a copy of the given 2^`N-by-2^N matrix`. More...

QMatrix	toQMatrix (DiagonalOp op)
	Returns a 2^`N-by-2^N complex` diagonal matrix form of the DiagonalOp. More...

QMatrix	toQMatrix (PauliHamil hamil)
	Returns a 2^`N-by-2^N Hermitian` matrix form of the PauliHamil. More...

QMatrix	toQMatrix (qreal coeffs, pauliOpType paulis, int numQubits, int numTerms)
	Returns a 2^`N-by-2^N Hermitian` matrix form of the specified weighted sum of Pauli products. More...

QMatrix	toQMatrix (Qureg qureg)
	Returns an equal-size copy of the given density matrix `qureg`. More...

void	toQureg (Qureg qureg, QMatrix mat)
	Initialises the density matrix `qureg` to have the same amplitudes as `mat`. More...

void	toQureg (Qureg qureg, QVector vec)
	Initialises the state-vector `qureg` to have the same amplitudes as `vec`. More...

QVector	toQVector (DiagonalOp op)
	Returns a vector with the same of the full diagonal operator, populated with `op's` elements. More...

QVector	toQVector (Qureg qureg)
	Returns an equal-size copy of the given state-vector `qureg`. More...

void	updateIndices (int oldEl, int newEl, int list1, int len1, int list2, int len2)

void	writeToFileSynch (char *fn, const string &contents)
	Writes contents to the file with filename fn, which is created and/or overwritten. More...

Variables
static int	fn_unique_suffix_id = 0

Detailed Description

Unoptimised, analytic implementations of matrix operations used by QuEST_unit_tests

Author: Tyson Jones

Definition in file utilities.cpp.

Macro Definition Documentation

◆ DEMAND

#define DEMAND ( cond ) if (!(cond)) FAIL( );

Definition at line 24 of file utilities.cpp.

◆ macro_copyComplexMatrix

#define macro_copyComplexMatrix	(	dest,
		src
	)

Value:

    { \
    for (size_t i=0; i<dest.size(); i++) \
        for (size_t j=0; j<dest.size(); j++) \
            dest[i][j] = qcomp(src.real[i][j], src.imag[i][j]); \
}

Copies ComplexMatrix structures into a QMatrix.

Definition at line 1039 of file utilities.cpp.

◆ macro_copyQMatrix

#define macro_copyQMatrix	(	dest,
		src
	)

Value:

    { \
    for (size_t i=0; i<src.size(); i++) { \
        for (size_t j=0; j<src.size(); j++) { \
            dest.real[i][j] = real(src[i][j]); \
            dest.imag[i][j] = imag(src[i][j]); \
        } \
    } \
}

Definition at line 1013 of file utilities.cpp.

Function Documentation

◆ expI()

qcomp expI ( qreal phase )

Returns the unit-norm complex number exp(i*phase).

This function uses the Euler formula, and avoids problems with calling exp(complex) in a platform agnostic way

Definition at line 417 of file utilities.cpp.

                         {
     return qcomp(cos(phase), sin(phase));
 }

References qcomp.

Referenced by getDFT(), setDiagMatrixOverrides(), and TEST_CASE().

◆ getIndexOfTargetValues()

long long int getIndexOfTargetValues	(	long long int	ref,
		int *	targs,
		int	numTargs,
		int	targVal
	)

Definition at line 689 of file utilities.cpp.

                                                                                                {
     // ref state is the starting index, where the targets can be in any bit state;
     // on the bits of the non-target qubits matter 
     
     for (int t=0; t<numTargs; t++) {
         int bit = (targVal >> t) & 1;
         ref = setBit(ref, targs[t], bit);    
     }
     return ref;
 }

References setBit().

Referenced by getDFT().

◆ operator*() [1/7]

QMatrix operator*	(	const qcomp &	a,
		const QMatrix &	m
	)

Definition at line 92 of file utilities.cpp.

                                                       {
     QMatrix out = m;
     for (size_t r=0; r<m.size(); r++)
         for (size_t c=0; c<m.size(); c++)
             out[r][c] *= a;
     return out;
 }

◆ operator*() [2/7]

QVector operator*	(	const qcomp &	a,
		const QVector &	v
	)

Definition at line 40 of file utilities.cpp.

                                                       {
     QVector out = v;
     for (size_t i=0; i<v.size(); i++)
         out[i] *= a;
     return out;
 }

◆ operator*() [3/7]

QMatrix operator*	(	const QMatrix &	m,
		const qcomp &	a
	)

Definition at line 99 of file utilities.cpp.

                                                       {
     return a * m;
 }

◆ operator*() [4/7]

QVector operator*	(	const QMatrix &	m,
		const QVector &	v
	)

Definition at line 134 of file utilities.cpp.

                                                         {
     DEMAND( m.size() == v.size() );
     QVector prod = QVector(v.size());
     for (size_t r=0; r<v.size(); r++)
         for (size_t c=0; c<v.size(); c++)
             prod[r] += m[r][c] * v[c];
     return prod;
 }

References DEMAND.

◆ operator*() [5/7]

QMatrix operator*	(	const QMatrix &	m1,
		const QMatrix &	m2
	)

Definition at line 109 of file utilities.cpp.

                                                           {
     QMatrix prod = m1; // will be completely overwritten
     for (size_t r=0; r<m1.size(); r++)
         for (size_t c=0; c<m1.size(); c++) {
             prod[r][c] = 0;
             for (size_t k=0; k<m1.size(); k++)
                 prod[r][c] += m1[r][k] * m2[k][c];
         }
     return prod;
 }

◆ operator*() [6/7]

QVector operator*	(	const QVector &	v,
		const qcomp &	a
	)

Definition at line 46 of file utilities.cpp.

                                                       {
     return a * v;
 }

◆ operator*() [7/7]

qcomp operator*	(	const QVector &	v1,
		const QVector &	v2
	)

Definition at line 56 of file utilities.cpp.

                                                         {
     // this is sum_i v1_i conj(v2_i)
     DEMAND( v1.size() == v2.size() );
     qcomp out = 0;
     for (size_t i=0; i<v1.size(); i++)
         out += v1[i] * conj(v2[i]);
     return out;
 }

References DEMAND, and qcomp.

◆ operator*=() [1/3]

void operator*=	(	QMatrix &	m1,
		const QMatrix &	m2
	)

Definition at line 131 of file utilities.cpp.

                                                   {
     m1 = m1 * m2;
 }

◆ operator*=() [2/3]

void operator*=	(	QMatrix &	m1,
		const qreal &	a
	)

Definition at line 125 of file utilities.cpp.

                                                {
     m1 = m1 * a;
 }

◆ operator*=() [3/3]

void operator*=	(	QVector &	v1,
		const qcomp &	a
	)

Definition at line 70 of file utilities.cpp.

                                                {
     v1 = v1 * a;
 }

◆ operator+() [1/2]

QMatrix operator+	(	const QMatrix &	m1,
		const QMatrix &	m2
	)

Definition at line 76 of file utilities.cpp.

                                                           {
     DEMAND( m1.size() == m2.size() );
     QMatrix out = m1;
     for (size_t r=0; r<m1.size(); r++)
         for (size_t c=0; c<m1.size(); c++)
             out[r][c] += m2[r][c];
     return out;
 }

References DEMAND.

◆ operator+() [2/2]

QVector operator+	(	const QVector &	v1,
		const QVector &	v2
	)

Definition at line 26 of file utilities.cpp.

                                                           {
     DEMAND( v1.size() == v2.size() );
     QVector out = v1;
     for (size_t i=0; i<v2.size(); i++)
         out[i] += v2[i];
     return out;
 }

References DEMAND.

◆ operator+=() [1/2]

void operator+=	(	QMatrix &	m1,
		const QMatrix &	m2
	)

Definition at line 119 of file utilities.cpp.

                                                   {
     m1 = m1 + m2;
 }

◆ operator+=() [2/2]

void operator+=	(	QVector &	v1,
		const QVector &	v2
	)

Definition at line 64 of file utilities.cpp.

                                                   { // these violate += returns (fine)
     v1 = v1 + v2;
 }

◆ operator-() [1/2]

QMatrix operator-	(	const QMatrix &	m1,
		const QMatrix &	m2
	)

Definition at line 84 of file utilities.cpp.

                                                           {
     DEMAND( m1.size() == m2.size() );
     QMatrix out = m1;
     for (size_t r=0; r<m1.size(); r++)
         for (size_t c=0; c<m1.size(); c++)
             out[r][c] -= m2[r][c];
     return out;
 }

References DEMAND.

◆ operator-() [2/2]

QVector operator-	(	const QVector &	v1,
		const QVector &	v2
	)

Definition at line 33 of file utilities.cpp.

                                                           {
     DEMAND( v1.size() == v2.size() );
     QVector out = v1;
     for (size_t i=0; i<v2.size(); i++)
         out[i] -= v2[i];
     return out;
 }

References DEMAND.

◆ operator-=() [1/2]

void operator-=	(	QMatrix &	m1,
		const QMatrix &	m2
	)

Definition at line 122 of file utilities.cpp.

                                                   {
     m1 = m1 - m2;
 }

◆ operator-=() [2/2]

void operator-=	(	QVector &	v1,
		const QVector &	v2
	)

Definition at line 67 of file utilities.cpp.

                                                   {
     v1 = v1 - v2;
 }

◆ operator/() [1/2]

QMatrix operator/	(	const QMatrix &	m,
		const qcomp &	a
	)

Definition at line 102 of file utilities.cpp.

                                                       {
     QMatrix out = m;
     for (size_t r=0; r<m.size(); r++)
         for (size_t c=0; c<m.size(); c++)
             out[r][c] /= a;
     return out;
 }

◆ operator/() [2/2]

QVector operator/	(	const QVector &	v,
		const qcomp &	a
	)

Definition at line 49 of file utilities.cpp.

                                                       {
     DEMAND( abs(a) != 0 );
     QVector out = v;
     for (size_t i=0; i<v.size(); i++)
         out[i] /= a;
     return out;
 }

References DEMAND.

◆ operator/=() [1/2]

void operator/=	(	QMatrix &	m1,
		const qreal &	a
	)

Definition at line 128 of file utilities.cpp.

                                                {
     m1 = m1 / a;
 }

◆ operator/=() [2/2]

void operator/=	(	QVector &	v1,
		const qcomp &	a
	)

Definition at line 73 of file utilities.cpp.

                                                {
     v1 = v1 / a;
 }

◆ setBit()

long long int setBit	(	long long int	num,
		int	bitInd,
		int	bitVal
	)

Definition at line 681 of file utilities.cpp.

                                                                 {
     DEMAND( (bitVal == 0 || bitVal == 1) );
     DEMAND( num >= 0 );
     DEMAND( bitInd >= 0 );
     
     return (num & ~(1UL << bitInd)) | (bitVal << bitInd);
 }

References DEMAND.

Referenced by getIndexOfTargetValues().

◆ sublists() [1/3]

Catch::Generators::GeneratorWrapper<int*> sublists	(	Catch::Generators::GeneratorWrapper< int > &&	gen,
		int	numSamps
	)

Definition at line 1510 of file utilities.cpp.

   {
     int exclude[] = {};  
     return Catch::Generators::GeneratorWrapper<int*>(
         std::unique_ptr<Catch::Generators::IGenerator<int*>>(
             new SubListGenerator(std::move(gen), numSamps, exclude, 0)));
 }

◆ sublists() [2/3]

Catch::Generators::GeneratorWrapper<int*> sublists	(	Catch::Generators::GeneratorWrapper< int > &&	gen,
		int	numSamps,
		const int *	exclude,
		int	numExclude
	)

Definition at line 1495 of file utilities.cpp.

   {    
     return Catch::Generators::GeneratorWrapper<int*>(
         std::unique_ptr<Catch::Generators::IGenerator<int*>>(
             new SubListGenerator(std::move(gen), numSamps, exclude, numExclude)));
 }

◆ sublists() [3/3]

Catch::Generators::GeneratorWrapper<int*> sublists	(	Catch::Generators::GeneratorWrapper< int > &&	gen,
		int	numSamps,
		int	excluded
	)

Definition at line 1502 of file utilities.cpp.

   {
     int exclude[] = {excluded};  
     return Catch::Generators::GeneratorWrapper<int*>(
         std::unique_ptr<Catch::Generators::IGenerator<int*>>(
             new SubListGenerator(std::move(gen), numSamps, exclude, 1)));
 }

◆ updateIndices()

void updateIndices	(	int	oldEl,
		int	newEl,
		int *	list1,
		int	len1,
		int *	list2,
		int	len2
	)

Definition at line 287 of file utilities.cpp.

                                                                                      {
     if (list1 != NULL) {
         for (int i=0; i<len1; i++) {
             if (list1[i] == oldEl) {
                 list1[i] = newEl;
                 return;
             }
         }
     }
     for (int i=0; i<len2; i++) {
         if (list2[i] == oldEl) {
             list2[i] = newEl;
             return;
         }
     }
 }

Referenced by getFullOperatorMatrix().

Variable Documentation

◆ fn_unique_suffix_id

int fn_unique_suffix_id = 0

static

Definition at line 1356 of file utilities.cpp.

Referenced by setUniqueFilename().

Data Structures

Macros

Functions

Variables

Detailed Description

Macro Definition Documentation

◆ DEMAND

◆ macro_copyComplexMatrix

◆ macro_copyQMatrix

Function Documentation

◆ expI()

◆ getIndexOfTargetValues()

◆ operator*() [1/7]

◆ operator*() [2/7]

◆ operator*() [3/7]

◆ operator*() [4/7]

◆ operator*() [5/7]

◆ operator*() [6/7]

◆ operator*() [7/7]

◆ operator*=() [1/3]

◆ operator*=() [2/3]

◆ operator*=() [3/3]

◆ operator+() [1/2]

◆ operator+() [2/2]

◆ operator+=() [1/2]

◆ operator+=() [2/2]

◆ operator-() [1/2]

◆ operator-() [2/2]

◆ operator-=() [1/2]

◆ operator-=() [2/2]

◆ operator/() [1/2]

◆ operator/() [2/2]

◆ operator/=() [1/2]

◆ operator/=() [2/2]

◆ setBit()

◆ sublists() [1/3]

◆ sublists() [2/3]

◆ sublists() [3/3]

◆ updateIndices()

Variable Documentation

◆ fn_unique_suffix_id