Function qrTester

test/qr_dense.cpp:90–137

template<typename T>
void qrTester(const int m, const int n, double eps) {
    try {
        SUPPORTED_TYPE_CHECK(T);
        LAPACK_ENABLED_CHECK();

#if 1
        array in = cpu_randu<T>(dim4(m, n));
#else
        array in = randu(m, n, (dtype)dtype_traits<T>::af_type);
#endif

        //! [ex_qr_unpacked]
        array q, r, tau;
        qr(q, r, tau, in);
        //! [ex_qr_unpacked]

        array qq = matmul(q, q.H());
        array ii = identity(qq.dims(), qq.type());

        ASSERT_NEAR(0, max<double>(abs(real(qq - ii))), eps);
        ASSERT_NEAR(0, max<double>(abs(imag(qq - ii))), eps);

        //! [ex_qr_recon]
        array re = matmul(q, r);
        //! [ex_qr_recon]

        ASSERT_NEAR(0, max<double>(abs(real(re - in))), eps);
        ASSERT_NEAR(0, max<double>(abs(imag(re - in))), eps);

        //! [ex_qr_packed]
        array out = in.copy();
        array tau2;
        qrInPlace(tau2, out);
        //! [ex_qr_packed]

        array r2 = upper(out);

        ASSERT_NEAR(0, max<double>(abs(real(tau - tau2))), eps);
        ASSERT_NEAR(0, max<double>(abs(imag(tau - tau2))), eps);

        ASSERT_NEAR(0, max<double>(abs(real(r2 - r))), eps);
        ASSERT_NEAR(0, max<double>(abs(imag(r2 - r))), eps);

    } catch (exception& ex) {
        cout << ex.what() << endl;
        throw;
    }
}

template<typename T>
double eps();