Nektar::NekLinAlgTests Namespace Reference

Functions
	BOOST_AUTO_TEST_CASE (TestGramSchmidtOrthogonalizationBookExample)

Function Documentation

BOOST_AUTO_TEST_CASE()

Nektar::NekLinAlgTests::BOOST_AUTO_TEST_CASE

(

TestGramSchmidtOrthogonalizationBookExample

)

Definition at line 47 of file testNekLinAlgAlgorithms.cpp.

{
    double buf_0[] = {1.0, 0.0, 2.0};
    double buf_1[] = {2.0, 3.0, 0.0};
 
    std::vector<NekVector<double>> x;
    x.push_back(NekVector<double>(3, buf_0));
    x.push_back(NekVector<double>(3, buf_1));
 
    std::vector<NekVector<double>> q = GramSchmidtOrthogonalization(x);
 
    BOOST_CHECK_EQUAL(q.size(), 2);
 
    double epsilon = 1e-1;
    BOOST_CHECK_CLOSE(q[0][0], .4472, epsilon);
    BOOST_CHECK_CLOSE(q[0][1], .0, epsilon);
    BOOST_CHECK_CLOSE(q[0][2], .8942, epsilon);
 
    BOOST_CHECK_CLOSE(q[1][0], .45811, epsilon);
    BOOST_CHECK_CLOSE(q[1][1], .85890, epsilon);
    BOOST_CHECK_CLOSE(q[1][2], -.22890, epsilon);
}

References Nektar::GramSchmidtOrthogonalization(), and Nektar::UnitTests::q().