#223 Inversion of FieldMatrix is buggy

Metadata

Property	Value
Reported by	Oliver Sander (oliver.sander@tu-dresden.de)
Reported at	Dec 18, 2006 09:26
Type	Bug Report
Version	Git (pre2.4) [autotools]
Operating System	Unspecified / All
Last edited by	Oliver Sander (oliver.sander@tu-dresden.de)
Last edited at	Dec 21, 2006 10:54
Closed by	Markus Blatt (markus@dr-blatt.de)
Closed at	Dec 21, 2006 13:35
Closed in version	1.0
Resolution	Fixed
Comment	Fixed in rev 4824

Description

Consider the following test program:

#include "config.h"

#include <dune/common/fmatrix.hh>

int main() { Dune::FieldMatrix<double,3,3> A(0);

A[0][1] = 1;
A[1][0] = 1;
A[2][2] = 1;

std::cout << "A:" << std::endl << A << std::endl;

A.invert();

std::cout << "A inverse:" << std::endl << A << std::endl;

}

And here is the corresponding output:

A: 0 1 0 1 0 0 0 0 1

A inverse: 1 0 0 0 1 0 0 0 1

Not quite correct :-)