Skip to content
Snippets Groups Projects
Commit c302c490 authored by Markus Blatt's avatar Markus Blatt
Browse files

right hand side is not scaled by SuperLU. No rescaling needed.

[[Imported from SVN: r1257]]
parent eded540f
Branches
Tags
No related merge requests found
......@@ -153,7 +153,6 @@ namespace Dune
bool first, verbose;
};
template<typename T, typename A, int n, int m>
SuperLU<BCRSMatrix<FieldMatrix<T,n,m>,A> >
::~SuperLU()
......@@ -420,10 +419,6 @@ namespace Dune
&L, &U, work, lwork, &B, &X, &rpg, &rcond, &ferr, &berr,
&memusage, &stat, &info);
if(options.Equil==YES)
// undo scaling of right hand side
std::transform(b, b+mat.M(), C, b, std::divides<T>());
if(verbose) {
dinfo<<"Triangular solve: dgssvx() returns info "<< info<<std::endl;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment