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

Store columns of ParallelOverlappingILU0::upper in descending order

Browse Source 
This is the order used by ILU in dune-istl
This commit is contained in:
Markus Blatt
2017-06-14 09:09:26 +00:00
parent e3c898030c
commit c2129185ae
1 changed files with 2 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
opm/autodiff/ParallelOverlappingILU0.hpp
View File

@@ -120,20 +120,18 @@ namespace Opm
// relative to ILU // relative to ILU
for (auto i=A.beforeEnd(); i!=rendi; --i, ++ row ) for (auto i=A.beforeEnd(); i!=rendi; --i, ++ row )
{ {
// coliterator is diagonal after the following loop
const auto endij=(*i).end(); // end of row i
const size_type iIndex = i.index(); const size_type iIndex = i.index();
upper.reserveAdditional( (*i).size() ); upper.reserveAdditional( (*i).size() );
// store in reverse row order // store in reverse row order
// eliminate entries left of diagonal; store L factor // eliminate entries left of diagonal; store L factor
for (auto j=(*i).begin(); j != endij; ++j ) for (auto j=(*i).beforeEnd(); j.index()>=iIndex; --j )
{ {
const size_type jIndex = j.index(); const size_type jIndex = j.index();
if( j.index() == iIndex ) if( j.index() == iIndex )
{ {
inv[ row ] = (*j); inv[ row ] = (*j);
break;
} }
else if ( j.index() >= i.index() ) else if ( j.index() >= i.index() )
{ {