[Trilinos-Users] Teuchos::SerialDenseSolver equilibration: unexpected behavior

[Nico Schlömer]

Hi,

I just encountered a rather odd behavior of Teuchos::SerialDenseSolver 
when used with equilibration.
After building a specific 3x3 equation system (code attached), I call

solver.solve();

without equilibration and get the output

================ *snip* ================
WARNING!  SerialDenseSolver<T>::solve: System should be equilibrated!
-0.0416252 0.836261 0.358417
================ *snap* ================

This solution is correct. Using

solver.equilibrateMatrix();
solver.equilibrateRHS();
solver.solve();
solver.unequilibrateLHS();

however, produces

================ *snip* ================
-0.0416252 0.653305 0.323626
================ *snap* ================

I guess I must be making some sort of mistake here but I fail to see 
what the issue could be.
The comments in ./packages/amesos2/src/Amesos2_Lapack_def.hpp suggest 
that there may be more to it.
Any hints?

Cheers,
Nico
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2011-10-19-teuchos-serialdensesolver.tar.bz2
Type: application/x-bzip
Size: 2344 bytes
Desc: not available
Url : https://software.sandia.gov/pipermail/trilinos-users/attachments/20111019/5159d4c6/attachment.bin