10 #pragma implementation
17 #include "CLHEP/Matrix/GenMatrix.h"
18 #include "CLHEP/Matrix/SymMatrix.h"
19 #include "CLHEP/Matrix/Matrix.h"
21 #ifdef HEP_DEBUG_INLINE
22 #include "CLHEP/Matrix/GenMatrix.icc"
27 #ifdef HEP_THIS_FUNCTION_IS_NOT_NEEDED
28 static void delete_array(
double *hm)
36 for(
int r=1;r<=hm.
num_row();r++) {
38 for(
int c=1;c<=hm.
num_col();c++) {
48 for(
int c=1;c<=hm.
num_col();c++) {
50 for(
int r=1;r<=hm.
num_row();r++)
62 for(r=1;r<=A.num_row();r++)
64 for(
int i=1;i<=hm.
num_row();i++)
65 A.
fast(r,c)=hm(i,r)*hm(i,c);
67 double max=fabs(A(1,1));
68 for(r=2;r<=A.num_row();r++)
69 if(max<fabs(A(r,r))) max=fabs(A(r,r));
75 std::cerr << es << std::endl;
76 std::cerr <<
"---Exiting to System." << std::endl;
82 for (
int k1=1; k1<=
num_row(); k1++)
83 for (
int k2=1; k2<=
num_col(); k2++)
84 if(o(k1,k2) != (*this)(k1,k2))
return false;
double norm(const HepGenMatrix &m)
virtual int num_row() const =0
const double & fast(int row, int col) const
double norm1(const HepGenMatrix &m)
virtual int num_col() const =0
void delete_m(int size, double *)
static void error(const char *s)
virtual bool operator==(const HepGenMatrix &) const
double norm_infinity(const HepGenMatrix &m)
HepMatrix diagonalize(HepSymMatrix *s)