PauliStrSum

Functions for applying, exponentiating or Trotterising a weigthed sum of Pauli tensors. More...

Functions
void	applyTrotterizedPauliStrSumGadget (Qureg qureg, PauliStrSum sum, qreal angle, int order, int reps)
void	multiplyPauliStrSum (Qureg qureg, PauliStrSum sum, Qureg workspace)

Detailed Description

Functions for applying, exponentiating or Trotterising a weigthed sum of Pauli tensors.

Function Documentation

applyTrotterizedPauliStrSumGadget()

void applyTrotterizedPauliStrSumGadget	(	Qureg	qureg,
		PauliStrSum	sum,
		qreal	angle,
		int	order,
		int	reps )

Note

Documentation for this function or struct is under construction!

Warning

This function has not yet been unit tested and may contain bugs. Please use with caution!

Definition at line 893 of file operations.cpp.

                                                                                                       {
    validate_quregFields(qureg, __func__);
    validate_pauliStrSumFields(sum, __func__);
    validate_pauliStrSumTargets(sum, qureg, __func__);
    validate_pauliStrSumIsHermitian(sum, __func__);
    validate_trotterParams(qureg, order, reps, __func__);
 
    /// @todo
    /// the accuracy of Trotterisation is greatly improved by randomisation
    /// or (even sub-optimal) grouping into commuting terms. Should we 
    /// implement these here or into another function?
 
    if (angle == 0)
        return;
 
    for (int r=0; r<reps; r++)
        applyHigherOrderTrotter(qureg, sum, angle/reps, order);
}

multiplyPauliStrSum()

void multiplyPauliStrSum	(	Qureg	qureg,
		PauliStrSum	sum,
		Qureg	workspace )

Note

Documentation for this function or struct is under construction!

Attention

This function's input validation has not yet been tested, so erroneous usage may produce unexpected output. Please use with caution!

Definition at line 836 of file operations.cpp.

                                                                        {
    validate_quregFields(qureg, __func__);
    validate_quregFields(workspace, __func__);
    validate_quregCanBeWorkspace(qureg, workspace, __func__);
    validate_pauliStrSumFields(sum, __func__);
    validate_pauliStrSumTargets(sum, qureg, __func__);
 
    // clone qureg to workspace, set qureg to blank
    localiser_statevec_setQuregToSuperposition(0, workspace, 1, qureg, 0, qureg);
    localiser_statevec_initUniformState(qureg, 0);
 
    // apply each term in-turn, mixing into output qureg, then undo using idempotency
    for (qindex i=0; i<sum.numTerms; i++) {
        localiser_statevec_anyCtrlPauliTensor(workspace, {}, {}, sum.strings[i]);
        localiser_statevec_setQuregToSuperposition(1, qureg, sum.coeffs[i], workspace, 0, workspace);
        localiser_statevec_anyCtrlPauliTensor(workspace, {}, {}, sum.strings[i]);
    }
 
    // workspace -> qureg, and qureg -> sum * qureg
}

Referenced by TEST_CASE().