[Insight-developers] vnl_svd on Borland compiler

Aljaz Noe noe@grasp.cis.upenn.edu
Tue, 1 Oct 2002 18:57:39 -0400


Hi!

Consider the following piece of code to compute the determinant of a
matrix:

  vnl_matrix<double> m(2,2,0.0);
  m[0][0]=2.0; m[1][0]=1.0;
  m[0][1]=1.0; m[1][1]=2.0;
  double det=vnl_svd<double>(m).determinant_magnitude();
  cout<<det<<"\n";

The result should be det=3. But if the same code is compiled on Borland
bcc32 compiler the result is det=2. Comments???

This also seems to be the reason why many of the FEM tests fail on
Borland compiler.

Aljaz