QR_test.hhs

/**
 * @author Jason Reynolds
 * @param
 * @returns
 * 
 * Function unit test for QR, makes sure it decomposes the right answer for both Q and R with hedgehog's built in tolerance level
 */



function QR_test() {

    *import math: QR

    if (QR([[1, -1, 4], [1, 4, -2], [1, 4, 2], [1, -1, 0]]).R.val === [[2, 3, 2], [0, 5, -2], [0, 0, 4], [0, 0, 0]]) {

    }

    else {
        throw 'unit test failed for QR_test.'
    }


    //built in tolerance in hedgehog lab fixes the 0.49999 issue
    if (QR([[1, -1, 4], [1, 4, -2], [1, 4, 2], [1, -1, 0]]).Q.val === [[0.5, -0.5, 0.5, 0.5], [0.5, 0.5, -0.5, 0.5], [0.5, 0.5, 0.5, -0.5], [0.5, -0.5, -0.5, -0.5]]) {

    }

    else {
        throw 'unit test failed for QR_test.'
    }

}