CLHEP::HepDiagMatrix Class Reference

#include <DiagMatrix.h>

Inheritance diagram for CLHEP::HepDiagMatrix:

Classes
class	HepDiagMatrix_row
class	HepDiagMatrix_row_const

Public Member Functions
	HepDiagMatrix ()
	HepDiagMatrix (int p)
	HepDiagMatrix (int p, int)
	HepDiagMatrix (int p, HepRandom &r)
	HepDiagMatrix (const HepDiagMatrix &hm1)
virtual	~HepDiagMatrix ()
int	num_row () const
int	num_col () const
double &	operator() (int row, int col)
const double &	operator() (int row, int col) const
double &	fast (int row, int col)
const double &	fast (int row, int col) const
void	assign (const HepMatrix &hm2)
void	assign (const HepSymMatrix &hm2)
void	assign (const HepDiagMatrix &hm2)
HepDiagMatrix &	operator*= (double t)
HepDiagMatrix &	operator/= (double t)
HepDiagMatrix &	operator+= (const HepDiagMatrix &hm2)
HepDiagMatrix &	operator-= (const HepDiagMatrix &hm2)
HepDiagMatrix &	operator= (const HepDiagMatrix &hm2)
HepDiagMatrix	operator- () const
HepDiagMatrix	T () const
HepDiagMatrix	apply (double(*f)(double, int, int)) const
HepSymMatrix	similarity (const HepMatrix &hm1) const
HepSymMatrix	similarityT (const HepMatrix &hm1) const
double	similarity (const HepVector &) const
HepDiagMatrix	sub (int min_row, int max_row) const
HepDiagMatrix	sub (int min_row, int max_row)
void	sub (int row, const HepDiagMatrix &hm1)
HepDiagMatrix	inverse (int &ierr) const
void	invert (int &ierr)
void	invert ()
HepDiagMatrix	inverse () const
double	determinant () const
double	trace () const
HepDiagMatrix_row	operator[] (int)
HepDiagMatrix_row_const	operator[] (int) const
Public Member Functions inherited from CLHEP::HepGenMatrix
virtual	~HepGenMatrix ()
virtual int	num_row () const =0
virtual int	num_col () const =0
virtual const double &	operator() (int row, int col) const =0
virtual double &	operator() (int row, int col)=0
virtual void	invert (int &)=0
HepGenMatrix_row	operator[] (int)
const HepGenMatrix_row_const	operator[] (int) const
virtual bool	operator== (const HepGenMatrix &) const

Protected Member Functions
int	num_size () const
Protected Member Functions inherited from CLHEP::HepGenMatrix
virtual int	num_size () const =0
void	delete_m (int size, double *)
double *	new_m (int size)

Friends
class	HepDiagMatrix_row
class	HepDiagMatrix_row_const
class	HepMatrix
class	HepSymMatrix
HepDiagMatrix	operator* (const HepDiagMatrix &hm1, const HepDiagMatrix &hm2)
HepDiagMatrix	operator+ (const HepDiagMatrix &hm1, const HepDiagMatrix &hm2)
HepDiagMatrix	operator- (const HepDiagMatrix &hm1, const HepDiagMatrix &hm2)
HepMatrix	operator* (const HepDiagMatrix &hm1, const HepMatrix &hm2)
HepMatrix	operator* (const HepMatrix &hm1, const HepDiagMatrix &hm2)
HepVector	operator* (const HepDiagMatrix &hm1, const HepVector &hm2)

Additional Inherited Members
Public Types inherited from CLHEP::HepGenMatrix
enum	{ size_max = 25 }
typedef std::vector< double, Alloc< double, 25 > >::iterator	mIter
typedef std::vector< double, Alloc< double, 25 > >::const_iterator	mcIter
Static Public Member Functions inherited from CLHEP::HepGenMatrix
static void	swap (int &, int &)
static void	swap (std::vector< double, Alloc< double, 25 > > &, std::vector< double, Alloc< double, 25 > > &)
static void	error (const char *s)

Detailed Description

Author

Definition at line 35 of file DiagMatrix.h.

Constructor & Destructor Documentation

HepDiagMatrix() [1/5]

CLHEP::HepDiagMatrix::HepDiagMatrix ( )

inline

HepDiagMatrix() [2/5]

CLHEP::HepDiagMatrix::HepDiagMatrix ( int  p )

explicit

Definition at line 67 of file DiagMatrix.cc.

   : m(p), nrow(p)
{
}

HepDiagMatrix() [3/5]

CLHEP::HepDiagMatrix::HepDiagMatrix	(	int	p,
		int	init
	)

Definition at line 72 of file DiagMatrix.cc.

   : m(p), nrow(p)
{   
   switch(init)
   {
   case 0:
      m.assign(nrow,0);
      break;
 
   case 1:
      {
     HepMatrix::mIter a=m.begin();
     HepMatrix::mIter b=m.begin() + p;
     for( ; a<b; a++) *a = 1.0;
     break;
      }
   default:
      error("DiagMatrix: initialization must be either 0 or 1.");
   }
}

HepDiagMatrix() [4/5]

CLHEP::HepDiagMatrix::HepDiagMatrix	(	int	p,
		HepRandom &	r
	)

Definition at line 93 of file DiagMatrix.cc.

  : m(p), nrow(p)
{
   HepMatrix::mIter a = m.begin();
   HepMatrix::mIter b = m.begin() + num_size();
   for(;a<b;a++) *a = r();
}

HepDiagMatrix() [5/5]

CLHEP::HepDiagMatrix::HepDiagMatrix

(

const HepDiagMatrix &

hm1

)

Definition at line 106 of file DiagMatrix.cc.

   : HepGenMatrix(hm1), m(hm1.nrow), nrow(hm1.nrow)
{
   m = hm1.m;
}

~HepDiagMatrix()

CLHEP::HepDiagMatrix::~HepDiagMatrix ( )

virtual

Definition at line 103 of file DiagMatrix.cc.

                              {
}

Member Function Documentation

apply()

HepDiagMatrix CLHEP::HepDiagMatrix::apply

(

double(*)(double, int, int)

f

)

const

Definition at line 581 of file DiagMatrix.cc.

{
#else
{
  HepDiagMatrix mret(num_row());
#endif
  HepMatrix::mcIter a = m.begin();
  HepMatrix::mIter b = mret.m.begin();
  for(int ir=1;ir<=num_row();ir++) {
    *(b++) = (*f)(*(a++), ir, ir);
  }
  return mret;
}

Referenced by main().

assign() [1/3]

void CLHEP::HepDiagMatrix::assign

(

const HepDiagMatrix &

hm2

)

assign() [2/3]

void CLHEP::HepDiagMatrix::assign

(

const HepMatrix &

hm2

)

Definition at line 598 of file DiagMatrix.cc.

{
   if(hm1.num_row()!=nrow)
   {
      nrow = hm1.num_row();
      m.resize(nrow);
   }
   HepMatrix::mcIter a = hm1.m.begin();
   HepMatrix::mIter b = m.begin();
   for(int r=1;r<=nrow;r++) {
      *(b++) = *a;
      if(r<nrow) a += (nrow+1);
   }
}

Referenced by main().

assign() [3/3]

void CLHEP::HepDiagMatrix::assign

(

const HepSymMatrix &

hm2

)

Definition at line 613 of file DiagMatrix.cc.

{
   if(hm1.num_row()!=nrow)
   {
      nrow = hm1.num_row();
      m.resize(nrow);
   }
   HepMatrix::mcIter a = hm1.m.begin();
   HepMatrix::mIter b = m.begin();
   for(int r=1;r<=nrow;r++) {
      *(b++) = *a;
      if(r<nrow) a += (r+1);
   }
}

determinant()

double CLHEP::HepDiagMatrix::determinant

(

)

const

Definition at line 711 of file DiagMatrix.cc.

                                        {
   double d = 1.0;
   HepMatrix::mcIter end = m.begin() + nrow;
   for (HepMatrix::mcIter p=m.begin(); p < end; p++)
      d *= *p;
   return d;
}

fast() [1/2]

double & CLHEP::HepDiagMatrix::fast	(	int	row,
		int	col
	)

Referenced by main().

fast() [2/2]

const double & CLHEP::HepDiagMatrix::fast	(	int	row,
		int	col
	)		const

inverse() [1/2]

HepDiagMatrix CLHEP::HepDiagMatrix::inverse ( ) const

inline

inverse() [2/2]

HepDiagMatrix CLHEP::HepDiagMatrix::inverse

(

int &

ierr

)

const

invert() [1/2]

void CLHEP::HepDiagMatrix::invert ( )

inline

invert() [2/2]

void CLHEP::HepDiagMatrix::invert ( int &  ierr )

virtual

Implements CLHEP::HepGenMatrix.

Definition at line 695 of file DiagMatrix.cc.

                                    {
  int n = num_row();
  ierr = 1;
  HepMatrix::mIter hmm = m.begin();
  int i;
  for(i=0;i<n;i++) {
    if(*(hmm++)==0) return;
  }
  ierr = 0;
  hmm = m.begin();
  for(i=0;i<n;i++) {
    *hmm = 1.0 / *hmm;
    hmm++;
  }  
}

num_col()

int CLHEP::HepDiagMatrix::num_col ( ) const

inlinevirtual

Implements CLHEP::HepGenMatrix.

Referenced by main(), CLHEP::HepMatrix::operator+=(), operator+=(), CLHEP::HepSymMatrix::operator+=(), CLHEP::HepMatrix::operator-=(), operator-=(), CLHEP::HepSymMatrix::operator-=(), and CLHEP::operator<<().

num_row()

int CLHEP::HepDiagMatrix::num_row ( ) const

inlinevirtual

Implements CLHEP::HepGenMatrix.

Referenced by invert(), main(), CLHEP::HepMatrix::operator+=(), operator+=(), CLHEP::HepSymMatrix::operator+=(), CLHEP::HepMatrix::operator-=(), operator-=(), CLHEP::HepSymMatrix::operator-=(), CLHEP::operator<<(), similarity(), and sub().

num_size()

int CLHEP::HepDiagMatrix::num_size ( ) const

inlineprotectedvirtual

Implements CLHEP::HepGenMatrix.

Referenced by HepDiagMatrix().

operator()() [1/2]

double & CLHEP::HepDiagMatrix::operator() ( int  row, int  col  )

virtual

Implements CLHEP::HepGenMatrix.

operator()() [2/2]

const double & CLHEP::HepDiagMatrix::operator() ( int  row, int  col  ) const

virtual

Implements CLHEP::HepGenMatrix.

operator*=()

HepDiagMatrix & CLHEP::HepDiagMatrix::operator*=

(

double

t

)

Definition at line 521 of file DiagMatrix.cc.

{
  SIMPLE_UOP(*=)
  return (*this);
}

operator+=()

HepDiagMatrix & CLHEP::HepDiagMatrix::operator+=

(

const HepDiagMatrix &

hm2

)

Definition at line 476 of file DiagMatrix.cc.

{
  CHK_DIM_2(num_row(),hm2.num_row(),num_col(),hm2.num_col(),+=);
  SIMPLE_BOP(+=)
  return (*this);
}

operator-()

HepDiagMatrix CLHEP::HepDiagMatrix::operator-

(

)

const

Definition at line 176 of file DiagMatrix.cc.

{
#else
{
   HepDiagMatrix hm2(nrow);
#endif
   HepMatrix::mcIter a=m.begin();
   HepMatrix::mIter b=hm2.m.begin();
   HepMatrix::mcIter e=m.begin()+num_size();
   for(;a<e; a++, b++) (*b) = -(*a);
   return hm2;
}

operator-=()

HepDiagMatrix & CLHEP::HepDiagMatrix::operator-=

(

const HepDiagMatrix &

hm2

)

Definition at line 508 of file DiagMatrix.cc.

{
  CHK_DIM_2(num_row(),hm2.num_row(),num_col(),hm2.num_col(),-=);
  SIMPLE_BOP(-=)
  return (*this);
}

operator/=()

HepDiagMatrix & CLHEP::HepDiagMatrix::operator/=

(

double

t

)

Definition at line 515 of file DiagMatrix.cc.

{
  SIMPLE_UOP(/=)
  return (*this);
}

operator=()

HepDiagMatrix & CLHEP::HepDiagMatrix::operator=

(

const HepDiagMatrix &

hm2

)

Definition at line 547 of file DiagMatrix.cc.

{
   if(hm1.nrow != nrow)
   {
      nrow = hm1.nrow;
      m.resize(nrow);
   }
   m=hm1.m;
   return (*this);
}

operator[]() [1/2]

HepDiagMatrix_row CLHEP::HepDiagMatrix::operator[] ( int  )

inline

operator[]() [2/2]

HepDiagMatrix_row_const CLHEP::HepDiagMatrix::operator[] ( int  ) const

inline

similarity() [1/2]

HepSymMatrix CLHEP::HepDiagMatrix::similarity

(

const HepMatrix &

hm1

)

const

Definition at line 628 of file DiagMatrix.cc.

{
#else
{
  HepSymMatrix mret(hm1.num_row());
#endif
  CHK_DIM_1(num_row(),hm1.num_col(),"similarity");
//  HepMatrix temp = hm1*(*this);
// If hm1*(*this) has correct dimensions, then so will the hm1.T multiplication.
// So there is no need to check dimensions again.
  HepMatrix::mIter mrc = mret.m.begin();
  for(int r=1;r<=mret.num_row();r++) {
    HepMatrix::mcIter mrr = hm1.m.begin()+(r-1)*hm1.num_col();
    HepMatrix::mcIter mc = hm1.m.begin();
    for(int c=1;c<=r;c++) {
      HepMatrix::mcIter mi = m.begin();
      double tmp = 0;
      HepMatrix::mcIter mr = mrr;
      for(int i=0;i<hm1.num_col();i++)
    tmp+=*(mr++) * *(mc++) * *(mi++);
      *(mrc++) = tmp;
    }
  }
  return mret;
}

Referenced by main(), and similarity().

similarity() [2/2]

double CLHEP::HepDiagMatrix::similarity

(

const HepVector &

hm1

)

const

Definition at line 656 of file DiagMatrix.cc.

{
  double mret;
  CHK_DIM_1(num_row(),hm1.num_row(),similarity);
  HepMatrix::mcIter mi = m.begin();
  HepMatrix::mcIter mv = hm1.m.begin();
  mret = *(mv)* *(mv)* *(mi++);
  mv++;
  for(int i=2;i<=hm1.num_row();i++) {
    mret+=*(mv)* *(mv)* *(mi++);
    mv++;
  }
  return mret;
}

similarityT()

HepSymMatrix CLHEP::HepDiagMatrix::similarityT

(

const HepMatrix &

hm1

)

const

Definition at line 671 of file DiagMatrix.cc.

{
#else
{
  HepSymMatrix mret(hm1.num_col());
#endif
  CHK_DIM_1(num_col(),hm1.num_row(),similarityT);
//  Matrix temp = (*this)*hm1;
// If hm1*(*this) has correct dimensions, then so will the hm1.T multiplication.
// So there is no need to check dimensions again.
  for(int r=1;r<=mret.num_row();r++)
    for(int c=1;c<=r;c++)
      {
    HepMatrix::mcIter mi = m.begin();
    double tmp = hm1(1,r)*hm1(1,c)* *(mi++);
    for(int i=2;i<=hm1.num_row();i++)
      tmp+=hm1(i,r)*hm1(i,c)* *(mi++);
    mret.fast(r,c) = tmp;
      }
  return mret;
}

Referenced by main().

sub() [1/3]

HepDiagMatrix CLHEP::HepDiagMatrix::sub	(	int	min_row,
		int	max_row
	)

Definition at line 135 of file DiagMatrix.cc.

{
  HepDiagMatrix mret(max_row-min_row+1);
  if(max_row > num_row())
    error("HepDiagMatrix::sub: Index out of range");
  HepMatrix::mIter a = mret.m.begin();
  HepMatrix::mIter b = m.begin() + min_row - 1;
  HepMatrix::mIter e = mret.m.begin() + mret.num_row();
  for(;a<e;) *(a++) = *(b++);
  return mret;
}

sub() [2/3]

HepDiagMatrix CLHEP::HepDiagMatrix::sub	(	int	min_row,
		int	max_row
	)		const

Definition at line 118 of file DiagMatrix.cc.

{
#else
{
  HepDiagMatrix mret(max_row-min_row+1);
#endif
  if(max_row > num_row())
    error("HepDiagMatrix::sub: Index out of range");
  HepMatrix::mIter a = mret.m.begin();
  HepMatrix::mcIter b = m.begin() + min_row - 1;
  HepMatrix::mIter e = mret.m.begin() + mret.num_row();
  for(;a<e;) *(a++) = *(b++);
  return mret;
}

Referenced by diagmatrix_test(), and main().

sub() [3/3]

void CLHEP::HepDiagMatrix::sub	(	int	row,
		const HepDiagMatrix &	hm1
	)

Definition at line 147 of file DiagMatrix.cc.

{
  if(row <1 || row+hm1.num_row()-1 > num_row() )
    error("HepDiagMatrix::sub: Index out of range");
  HepMatrix::mcIter a = hm1.m.begin();
  HepMatrix::mIter b = m.begin() + row - 1;
  HepMatrix::mcIter e = hm1.m.begin() + hm1.num_row();
  for(;a<e;) *(b++) = *(a++);
}

T()

HepDiagMatrix CLHEP::HepDiagMatrix::T

(

)

const

Referenced by main().

trace()

double CLHEP::HepDiagMatrix::trace

(

)

const

Definition at line 719 of file DiagMatrix.cc.

                                  {
   double d = 0.0;
   HepMatrix::mcIter end = m.begin() + nrow;
   for (HepMatrix::mcIter p=m.begin(); p < end; p++)
      d += *p;
   return d;
}

Friends And Related Function Documentation

HepDiagMatrix_row

friend class HepDiagMatrix_row

friend

Definition at line 170 of file DiagMatrix.h.

HepDiagMatrix_row_const

friend class HepDiagMatrix_row_const

friend

Definition at line 171 of file DiagMatrix.h.

HepMatrix

friend class HepMatrix

friend

Definition at line 172 of file DiagMatrix.h.

HepSymMatrix

friend class HepSymMatrix

friend

Definition at line 173 of file DiagMatrix.h.

operator* [1/4]

HepDiagMatrix operator* ( const HepDiagMatrix &  hm1, const HepDiagMatrix &  hm2  )

friend

Definition at line 413 of file DiagMatrix.cc.

{
#else
{
  HepDiagMatrix mret(hm1.num_row());
#endif
  CHK_DIM_1(hm1.num_col(),hm2.num_row(),*);
  HepMatrix::mIter a = mret.m.begin();
  HepMatrix::mcIter b = hm1.m.begin();
  HepMatrix::mcIter c = hm2.m.begin();
  HepMatrix::mIter e = mret.m.begin() + hm1.num_col();
  for(;a<e;) *(a++) = *(b++) * (*(c++));
  return mret;
}

operator* [2/4]

HepMatrix operator* ( const HepDiagMatrix &  hm1, const HepMatrix &  hm2  )

friend

Definition at line 392 of file DiagMatrix.cc.

{
#else
{
  HepMatrix mret(hm1.num_row(),hm2.num_col());
#endif
  CHK_DIM_1(hm1.num_col(),hm2.num_row(),*);
  HepMatrix::mcIter mit1=hm2.m.begin();
  HepMatrix::mIter mir=mret.m.begin();
  HepMatrix::mcIter mrr = hm1.m.begin();
  for(int irow=1;irow<=hm2.num_row();irow++) {
    for(int icol=1;icol<=hm2.num_col();icol++) {
      *(mir++) = *(mit1++) * (*mrr);
    }
    mrr++;
  }
  return mret;
}

operator* [3/4]

HepVector operator* ( const HepDiagMatrix &  hm1, const HepVector &  hm2  )

friend

Definition at line 430 of file DiagMatrix.cc.

{
#else
{
  HepVector mret(hm1.num_row());
#endif
  CHK_DIM_1(hm1.num_col(),hm2.num_row(),*);
  HepGenMatrix::mIter mir=mret.m.begin();
  HepGenMatrix::mcIter mi1 = hm1.m.begin(), mi2 = hm2.m.begin();
  for(int icol=1;icol<=hm1.num_col();icol++) {
    *(mir++) = *(mi1++) * *(mi2++);
  }
  return mret;
}

operator* [4/4]

HepMatrix operator* ( const HepMatrix &  hm1, const HepDiagMatrix &  hm2  )

friend

Definition at line 372 of file DiagMatrix.cc.

{
#else
  {
    HepMatrix mret(hm1.num_row(),hm2.num_col());
#endif
    CHK_DIM_1(hm1.num_col(),hm2.num_row(),*);
    HepMatrix::mcIter mit1=hm1.m.begin();
    HepMatrix::mIter mir=mret.m.begin();
    for(int irow=1;irow<=hm1.num_row();irow++) {
      HepMatrix::mcIter mcc = hm2.m.begin();
      for(int icol=1;icol<=hm1.num_col();icol++) {
    *(mir++) = *(mit1++) * (*(mcc++));
      }
    }
    return mret;
  }

operator+

HepDiagMatrix operator+ ( const HepDiagMatrix &  hm1, const HepDiagMatrix &  hm2  )

friend

Definition at line 221 of file DiagMatrix.cc.

{
#else
{
  HepDiagMatrix mret(hm1.nrow);
#endif
  CHK_DIM_1(hm1.nrow,hm2.nrow,+);
  SIMPLE_TOP(+)
  return mret;
}

operator-

HepDiagMatrix operator- ( const HepDiagMatrix &  hm1, const HepDiagMatrix &  hm2  )

friend

Definition at line 291 of file DiagMatrix.cc.

{
#else
{
  HepDiagMatrix mret(hm1.nrow);
#endif
  CHK_DIM_1(hm1.num_row(),hm2.num_row(),-);
  SIMPLE_TOP(-)
  return mret;
}

---

The documentation for this class was generated from the following files:

• DiagMatrix.h
• DiagMatrix.cc