test_matrix

wield.control.algorithms.statespace.dense.test.test_matrix

This is a pytest module needing documentation

pytest-html report

Functions

test_QRH()	This is a pytest needing documentation
test_QRHpivot()	This is a pytest needing documentation

Details

test_QRH()

This is a pytest needing documentation

code

def test_QRH():

    N = 10
    M = np.random.rand(N, N)
    M = np.array(
        [
            [1, 0, 0],
            [1, 1, 0],
            [0, 0, 0],
        ],
        float,
    )

    eye = np.eye(M.shape[0], M.shape[1])

    R, [Q], [QT] = matrix_algorithms.QR(
        mat=M,
        mshadow=None,
        qmul=[eye],
        qAmul=[eye],
        pivoting=False,
        # method   = 'Householder',
        method="Givens",
        Rexact=False,
    )
    R2, [Q], [QT] = matrix_algorithms.QR(
        mat=M,
        mshadow=None,
        qmul=[eye],
        qAmul=[eye],
        pivoting=False,
        # method   = 'Householder',
        method="Givens",
        Rexact=True,
    )

    import tabulate

    print("near", tabulate.tabulate(R))
    print("exact", tabulate.tabulate(R2))
    print(tabulate.tabulate(Q))
    print(tabulate.tabulate(QT))

    numpy.testing.assert_almost_equal(Q @ Q.T, eye)
    numpy.testing.assert_almost_equal(Q @ QT, eye)
    numpy.testing.assert_almost_equal(Q @ R2, M)

pytest information

This code is wrapped in a pytest function using conventions detailed in Pytest Conventions. The full name of this test, as known by the documentation, is:

wield.control.algorithms.statespace.dense.test.test_matrix.test_QRH

The full name is useful when building documentation, to link a reference to this page using :func:`name`, or directly include it with an autofunction directive. The collapse nodes below show every instance of the test run. There may only be one, but if the test was run multiple times through pytest parametrizations, the list can be longer.

test_QRH

output

status: skipped
duration: 0.000s
captured errors:
def test_QRH():

        N = 10
        M = np.random.rand(N, N)
        M = np.array(
            [
                [1, 0, 0],
                [1, 1, 0],
                [0, 0, 0],
            ],
            float,
        )

        eye = np.eye(M.shape[0], M.shape[1])

>       R, [Q], [QT] = matrix_algorithms.QR(
            mat=M,
            mshadow=None,
            qmul=[eye],
            qAmul=[eye],
            pivoting=False,
            # method   = 'Householder',
            method="Givens",
            Rexact=False,
        )
E       TypeError: QR() got an unexpected keyword argument 'qmul'

../../src/wield/control/algorithms/statespace/dense/test/test_matrix.py:44: TypeError

• capture.txt

test_QRHpivot()

This is a pytest needing documentation

code

def test_QRHpivot():

    N = 10
    M = np.random.rand(N, N)
    M = np.array(
        [
            [1, 0, 0],
            [0, 1, 0],
            [0, 1, 0],
        ],
        float,
    )

    eye = np.eye(M.shape[0], M.shape[1])

    R, [Q], [QT], P = matrix_algorithms.QR(
        mat=M,
        mshadow=None,
        qmul=[eye],
        qAmul=[eye],
        pivoting=True,
        # method   = 'Householder',
        method="Givens",
        Rexact=True,
    )

    import tabulate

    print(P)
    print(tabulate.tabulate(R))
    # print(tabulate.tabulate(Q))
    # print(tabulate.tabulate(QT))

    numpy.testing.assert_almost_equal(Q @ Q.T, eye)
    numpy.testing.assert_almost_equal(Q @ QT, eye)
    numpy.testing.assert_almost_equal((Q @ R)[:, P], M)

pytest information

This code is wrapped in a pytest function using conventions detailed in Pytest Conventions. The full name of this test, as known by the documentation, is:

wield.control.algorithms.statespace.dense.test.test_matrix.test_QRHpivot

The full name is useful when building documentation, to link a reference to this page using :func:`name`, or directly include it with an autofunction directive. The collapse nodes below show every instance of the test run. There may only be one, but if the test was run multiple times through pytest parametrizations, the list can be longer.

test_QRHpivot

output

status: skipped
duration: 0.000s
captured errors:
def test_QRHpivot():

        N = 10
        M = np.random.rand(N, N)
        M = np.array(
            [
                [1, 0, 0],
                [0, 1, 0],
                [0, 1, 0],
            ],
            float,
        )

        eye = np.eye(M.shape[0], M.shape[1])

>       R, [Q], [QT], P = matrix_algorithms.QR(
            mat=M,
            mshadow=None,
            qmul=[eye],
            qAmul=[eye],
            pivoting=True,
            # method   = 'Householder',
            method="Givens",
            Rexact=True,
        )
E       TypeError: QR() got an unexpected keyword argument 'qmul'

../../src/wield/control/algorithms/statespace/dense/test/test_matrix.py:92: TypeError

• capture.txt