Function assertMatrixIsUnitary

quest/src/core/validation.cpp:2419–2434

// type T can be CompMatr1, CompMatr2, CompMatr, DiagMatr1, DiagMatr2, DiagMatr, FullStateDiagMatr
template <class T> 
void assertMatrixIsUnitary(T matr, const char* caller) {

    // validate both stack and heap matrices have been correctly initialised
    validate_matrixFields(matr, caller);

    // validate heap matrices have ever written to their GPU memories (if exists)
    if constexpr (util_isHeapMatrixType<T>())
        validate_matrixIsSynced(matr, caller);

    // avoid superfluous expensive unitarity check below (do not overwrite .isApproxUnitary)
    if (isNumericalValidationDisabled())
        return;

    // may overwrite matr.isApproxUnitary of heap matrices, otherwise ignores epsilon
    assertThat(util_isUnitary(matr, global_validationEpsilon), report::MATRIX_NOT_UNITARY, caller);
}
void validate_matrixIsUnitary(CompMatr1 m, const char* caller) { assertMatrixIsUnitary(m, caller); }
void validate_matrixIsUnitary(CompMatr2 m, const char* caller) { assertMatrixIsUnitary(m, caller); }
void validate_matrixIsUnitary(CompMatr  m, const char* caller) { assertMatrixIsUnitary(m, caller); }