trf++/src/trfutil
diff -N smatrix.cpp
--- smatrix.cpp 6 Dec 2011 03:29:25 -0000 1.3
+++ /dev/null 1 Jan 1970 00:00:00 -0000
@@ -1,31 +0,0 @@
-// smatrix.cpp
-
-#include "trfutil/smatrix.h"
-//#include <CLHEP/Matrix/SymMatrix.h>
-#include "trfclhep/SymMatrix.h"
-using namespace CLHEP;
-
-int invert( smatrix<double>& sma ) {
-
- // copy matrix to HEP object
- int dim = sma.nrow();
- HepSymMatrix hep(dim);
- int i, j;
- for ( i=0; i<dim; ++i )
- for ( j=0; j<=i; ++j )
- hep.fast(i+1,j+1) = sma(i,j);
-
- // invert
- int err;
- hep.invert( err );
-
- // copy back to smatrix
- if ( ! err )
- for ( i=0; i<dim; ++i )
- for ( j=0; j<=i; ++j )
- sma(i,j) = hep.fast(i+1,j+1);
-
- // return error status.
- return err;
-
-}