Nektar::MatrixSubtractionTests Namespace Reference

Functions
	BOOST_AUTO_TEST_CASE (TestLhsFullRhsFullSubtraction)

Function Documentation

BOOST_AUTO_TEST_CASE()

Nektar::MatrixSubtractionTests::BOOST_AUTO_TEST_CASE

(

TestLhsFullRhsFullSubtraction

)

Definition at line 232 of file TestNekMatrixOperations.cpp.

{
    // double lhs_values[] = {2, 4, 6, 8,
    //                        10, 12, 14, 16,
    //                        18, 20, 22, 24,
    //                        26, 28, 30, 32};
 
    double lhs_values[] = {2, 10, 18, 26, 4, 12, 20, 28,
                           6, 14, 22, 30, 8, 16, 24, 32};
 
    std::shared_ptr<NekMatrix<NekDouble, StandardMatrixTag>> lhs1;
    std::shared_ptr<NekMatrix<NekMatrix<NekDouble>, ScaledMatrixTag>> lhs2;
    std::shared_ptr<NekMatrix<NekMatrix<NekDouble>, BlockMatrixTag>> lhs3;
 
    GenerateFullMatrices(lhs_values, 2.0, lhs1, lhs2, lhs3);
    // double rhs_values[] = {4, 8, 12, 16,
    //                        20, 24, 28, 32,
    //                        36, 40, 44, 48,
    //                        52, 56, 60, 64};
    double rhs_values[] = {4,  20, 36, 52, 8,  24, 40, 56,
                           12, 28, 44, 60, 16, 32, 48, 64};
    std::shared_ptr<NekMatrix<NekDouble, StandardMatrixTag>> rhs1;
    std::shared_ptr<NekMatrix<NekMatrix<NekDouble>, ScaledMatrixTag>> rhs2;
    std::shared_ptr<NekMatrix<NekMatrix<NekDouble>, BlockMatrixTag>> rhs3;
    GenerateFullMatrices(rhs_values, 2.0, rhs1, rhs2, rhs3);
 
    double result_values[16];
    std::transform(lhs_values, lhs_values + 16, rhs_values, result_values,
                   std::minus<NekDouble>());
    NekMatrix<NekDouble> result(4, 4, result_values);
    RunAllSubCombinations(*lhs1, *lhs2, *lhs3, *rhs1, *rhs2, *rhs3, result);
}

References Nektar::GenerateFullMatrices(), and Nektar::RunAllSubCombinations().