Function assertMatrixIsHermitian

quest/src/core/validation.cpp:2459–2474

// type T can be CompMatr1, CompMatr2, CompMatr, DiagMatr1, DiagMatr2, DiagMatr, FullStateDiagMatr
template <class T> 
void assertMatrixIsHermitian(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 hermiticity check below (do not overwrite  matr.isApproxHermitian)
    if (isNumericalValidationDisabled())
        return;

    // may overwrite matr.isApproxHermitian of heap matrices, otherwise ignores epsilon
    assertThat(util_isHermitian(matr, global_validationEpsilon), report::MATRIX_NOT_HERMITIAN, caller);
}
void validate_matrixIsHermitian(CompMatr1 m, const char* caller) { assertMatrixIsHermitian(m, caller); }
void validate_matrixIsHermitian(CompMatr2 m, const char* caller) { assertMatrixIsHermitian(m, caller); }
void validate_matrixIsHermitian(CompMatr  m, const char* caller) { assertMatrixIsHermitian(m, caller); }