OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add additional MILU test.

Browse Source
This commit is contained in:
Markus Blatt 2018-05-17 14:12:48 +02:00
parent 0bae240a42
commit 397fc78405
1 changed files with 15 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
tests/test_milu.cpp
View File

@ -77,6 +77,21 @@ void test_milu0(M& A)
#endif #endif
BOOST_ASSERT(point_difference < 1e-12); BOOST_ASSERT(point_difference < 1e-12);
} }
// Test that (LU)^-1Ae=e
A.mv(e, x1);
bilu_backsolve(ILU, x2, x1);
diff = x2;
diff -= e;
for ( std::size_t i = 0, end = A.N(); i < end; ++i)
{
auto point_difference = diff[i].two_norm();
#ifdef DEBUG
std::cout<<"index "<<i<<" size "<<diff.size()<<" point_difference "<<point_difference<<std::endl;
#endif
BOOST_CHECK_CLOSE(x2[i].two_norm(), e[i].two_norm(), 1e-12);
}
} }
template<class B, class Alloc> template<class B, class Alloc>