Nektar::BandedMatrixMatrixMultiplicationTests Namespace Reference

Functions
	BOOST_AUTO_TEST_CASE (TestMatrixMatrixMultiplication)

Function Documentation

BOOST_AUTO_TEST_CASE()

Nektar::BandedMatrixMatrixMultiplicationTests::BOOST_AUTO_TEST_CASE

(

TestMatrixMatrixMultiplication

)

Definition at line 431 of file TestBandedMatrixOperations.cpp.

{
    // [ 2  10   0  0 ]    [ 30 38  0  0 ]
    // [16   4  12  0 ]    [  0 32 40  0 ]
    // [22  18   6 14 ] *  [  0  0 34 42 ]
    // [ 0  24  20  8 ]    [  0  0  0 36 ]
 
    {
        double lhs_buf[] = {0,  2, 16, 22, 10, 4, 18, 24,
                            12, 6, 20, 0,  14, 8, 0,  0};
        double rhs_buf[] = {0, 30, 38, 32, 40, 34, 42, 36};
 
        NekMatrix<double> lhs1(4, 4, lhs_buf, eBANDED, 2, 1);
        std::shared_ptr<
            NekMatrix<NekMatrix<double, StandardMatrixTag>, ScaledMatrixTag>>
            lhs2;
        std::shared_ptr<NekMatrix<NekMatrix<double>, BlockMatrixTag>> lhs3;
 
        NekMatrix<double> rhs1(4, 4, rhs_buf, eBANDED, 0, 1);
        std::shared_ptr<NekMatrix<NekMatrix<double>, ScaledMatrixTag>> rhs2;
        std::shared_ptr<NekMatrix<NekMatrix<double>, BlockMatrixTag>> rhs3;
 
        GenerateMatrices(lhs1, 2.0, 2, 2, lhs2, lhs3);
        GenerateMatrices(rhs1, 2.0, 2, 2, rhs2, rhs3);
 
        double result_buf[] = {60,  480, 660, 0,    396, 736, 1412, 768,
                               400, 568, 924, 1640, 0,   504, 756,  1128};
        NekMatrix<double> result(4, 4, result_buf);
    }
}

References Nektar::eBANDED, and Nektar::GenerateMatrices().