// smatrix.cpp

#include "smatrix.h"
#include <CLHEP/Matrix/SymMatrix.h>

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;

}