Nektar::DiagonalMatrixOperationsUnitTests Namespace Reference

Functions
	BOOST_AUTO_TEST_CASE (TestDoubleDiagonalVectorMultiplication)
	BOOST_AUTO_TEST_CASE (TestDoubleScaledDiagonalVectorMultiplication)

Function Documentation

BOOST_AUTO_TEST_CASE() [1/2]

Nektar::DiagonalMatrixOperationsUnitTests::BOOST_AUTO_TEST_CASE

(

TestDoubleDiagonalVectorMultiplication

)

Definition at line 45 of file TestDiagonalMatrixOperations.cpp.

{
    double m_buf[] = {1, 2, 3};
    double v_buf[] = {4, 5, 6};
 
    NekMatrix<double> m(3, 3, m_buf, eDIAGONAL);
    NekVector<double> v(3, v_buf);
 
    double expected_result_buf[] = {4, 10, 18};
    NekVector<double> expected_result(3, expected_result_buf);
 
    NekVector<double> result = m * v;
 
    BOOST_CHECK_EQUAL(expected_result, result);
}

References Nektar::eDIAGONAL.

BOOST_AUTO_TEST_CASE() [2/2]

Nektar::DiagonalMatrixOperationsUnitTests::BOOST_AUTO_TEST_CASE

(

TestDoubleScaledDiagonalVectorMultiplication

)

Definition at line 61 of file TestDiagonalMatrixOperations.cpp.

{
    double m_buf[] = {1, 2, 3};
    double v_buf[] = {4, 5, 6};
 
    std::shared_ptr<NekMatrix<double>> inner(
        new NekMatrix<double>(3, 3, m_buf, eDIAGONAL));
    NekMatrix<NekMatrix<double>, ScaledMatrixTag> m(5.0, inner);
 
    NekVector<double> v(3, v_buf);
 
    double expected_result_buf[] = {20, 50, 90};
    NekVector<double> expected_result(3, expected_result_buf);
 
    NekVector<double> result = m * v;
 
    BOOST_CHECK_EQUAL(expected_result, result);
}

References Nektar::eDIAGONAL.