CLHEP::HepVector Class Reference

#include <Vector.h>

Inheritance diagram for CLHEP::HepVector:

Public Member Functions
	HepVector ()
	HepVector (int p)
	HepVector (int p, int)
	HepVector (int p, HepRandom &r)
	HepVector (const HepVector &v)
	HepVector (const HepMatrix &m)
virtual	~HepVector ()
const double &	operator() (int row) const
double &	operator() (int row)
const double &	operator[] (int row) const
double &	operator[] (int row)
virtual const double &	operator() (int row, int col) const
virtual double &	operator() (int row, int col)
HepVector &	operator*= (double t)
HepVector &	operator/= (double t)
HepVector &	operator+= (const HepMatrix &v2)
HepVector &	operator+= (const HepVector &v2)
HepVector &	operator-= (const HepMatrix &v2)
HepVector &	operator-= (const HepVector &v2)
HepVector &	operator= (const HepVector &hm2)
HepVector &	operator= (const HepMatrix &)
HepVector &	operator= (const Hep3Vector &)
HepVector	operator- () const
HepVector	apply (double(*f)(double, int)) const
HepVector	sub (int min_row, int max_row) const
HepVector	sub (int min_row, int max_row)
void	sub (int row, const HepVector &v1)
double	normsq () const
double	norm () const
virtual int	num_row () const
virtual int	num_col () const
HepMatrix	T () 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
virtual 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
class	HepSymMatrix
class	HepMatrix
void	swap (HepVector &v1, HepVector &v2)
double	dot (const HepVector &v1, const HepVector &v2)
HepVector	operator+ (const HepVector &v1, const HepVector &v2)
HepVector	operator- (const HepVector &v1, const HepVector &v2)
HepVector	operator* (const HepSymMatrix &hm1, const HepVector &hm2)
HepVector	operator* (const HepDiagMatrix &hm1, const HepVector &hm2)
HepMatrix	operator* (const HepVector &hm1, const HepMatrix &hm2)
HepVector	operator* (const HepMatrix &hm1, const HepVector &hm2)
HepVector	solve (const HepMatrix &a, const HepVector &v)
void	tridiagonal (HepSymMatrix *a, HepMatrix *hsm)
void	row_house (HepMatrix *, const HepMatrix &, double, int, int, int, int)
void	row_house (HepMatrix *, const HepVector &, double, int, int)
void	back_solve (const HepMatrix &R, HepVector *b)
void	col_house (HepMatrix *, const HepMatrix &, double, int, int, int, int)
HepVector	house (const HepSymMatrix &a, int row, int col)
HepVector	house (const HepMatrix &a, int row, int col)
void	house_with_update (HepMatrix *a, int row, int col)
HepSymMatrix	vT_times_v (const HepVector &v)
HepVector	qr_solve (HepMatrix *, const HepVector &)

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 Vector.h.

Constructor & Destructor Documentation

HepVector() [1/6]

CLHEP::HepVector::HepVector ( )

inline

HepVector() [2/6]

CLHEP::HepVector::HepVector ( int  p )

explicit

Definition at line 53 of file Vector.cc.

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

HepVector() [3/6]

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

Definition at line 58 of file Vector.cc.

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

HepVector() [4/6]

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

Definition at line 79 of file Vector.cc.

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

HepVector() [5/6]

CLHEP::HepVector::HepVector

(

const HepVector &

v

)

Definition at line 94 of file Vector.cc.

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

HepVector() [6/6]

CLHEP::HepVector::HepVector

(

const HepMatrix &

m

)

Definition at line 105 of file Vector.cc.

   : m(hm1.nrow), nrow(hm1.nrow)
{
   if (hm1.num_col() != 1)
      error("Vector::Vector(Matrix) : Matrix is not Nx1");
   
   m = hm1.m;
}

~HepVector()

CLHEP::HepVector::~HepVector ( )

virtual

Definition at line 91 of file Vector.cc.

                      {
}

Member Function Documentation

apply()

HepVector CLHEP::HepVector::apply

(

double(*)(double, int)

f

)

const

Definition at line 554 of file Vector.cc.

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

Referenced by main().

norm()

double CLHEP::HepVector::norm ( ) const

inline

Referenced by main().

normsq()

double CLHEP::HepVector::normsq ( ) const

inline

Referenced by main().

num_col()

int CLHEP::HepVector::num_col ( ) const

virtual

Implements CLHEP::HepGenMatrix.

Definition at line 118 of file Vector.cc.

{ return 1; }

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

num_row()

int CLHEP::HepVector::num_row ( ) const

virtual

Implements CLHEP::HepGenMatrix.

Definition at line 116 of file Vector.cc.

{return nrow;} 

Referenced by CLHEP::RandMultiGauss::fire(), CLHEP::RandMultiGauss::fireArray(), main(), operator+=(), CLHEP::HepMatrix::operator+=(), CLHEP::operator-(), operator-=(), CLHEP::HepMatrix::operator-=(), CLHEP::operator<<(), CLHEP::RandMultiGauss::RandMultiGauss(), CLHEP::HepDiagMatrix::similarity(), CLHEP::HepSymMatrix::similarity(), and sub().

num_size()

int CLHEP::HepVector::num_size ( ) const

protectedvirtual

Implements CLHEP::HepGenMatrix.

Definition at line 117 of file Vector.cc.

{return nrow;} 

operator()() [1/4]

double & CLHEP::HepVector::operator() ( int  row )

inline

operator()() [2/4]

const double & CLHEP::HepVector::operator() ( int  row ) const

inline

operator()() [3/4]

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

virtual

Implements CLHEP::HepGenMatrix.

Definition at line 128 of file Vector.cc.

{
#endif
 
  return *(m.begin()+(row-1));
}

operator()() [4/4]

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

virtual

Implements CLHEP::HepGenMatrix.

Definition at line 141 of file Vector.cc.

{
#endif
 
  return *(m.begin()+(row-1));
}

operator*=()

HepVector & CLHEP::HepVector::operator*=

(

double

t

)

Definition at line 449 of file Vector.cc.

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

operator+=() [1/2]

HepVector & CLHEP::HepVector::operator+=

(

const HepMatrix &

v2

)

Definition at line 408 of file Vector.cc.

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

operator+=() [2/2]

HepVector & CLHEP::HepVector::operator+=

(

const HepVector &

v2

)

Definition at line 415 of file Vector.cc.

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

operator-()

HepVector CLHEP::HepVector::operator-

(

)

const

Definition at line 212 of file Vector.cc.

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

operator-=() [1/2]

HepVector & CLHEP::HepVector::operator-=

(

const HepMatrix &

v2

)

Definition at line 429 of file Vector.cc.

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

operator-=() [2/2]

HepVector & CLHEP::HepVector::operator-=

(

const HepVector &

v2

)

Definition at line 436 of file Vector.cc.

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

operator/=()

HepVector & CLHEP::HepVector::operator/=

(

double

t

)

Definition at line 443 of file Vector.cc.

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

operator=() [1/3]

HepVector & CLHEP::HepVector::operator=

(

const Hep3Vector &

v

)

Definition at line 493 of file Vector.cc.

{
   if(nrow != 3)
   {
      nrow = 3;
      m.resize(nrow);
   }
   m[0] = v.x();
   m[1] = v.y();
   m[2] = v.z();
   return (*this);
}

operator=() [2/3]

HepVector & CLHEP::HepVector::operator=

(

const HepMatrix &

hm1

)

Definition at line 479 of file Vector.cc.

{
   if (hm1.num_col() != 1)
      error("Vector::operator=(Matrix) : Matrix is not Nx1");
   
   if(hm1.nrow != nrow)
   {
      nrow = hm1.nrow;
      m.resize(nrow);
   }
   m = hm1.m;
   return (*this);
}

operator=() [3/3]

HepVector & CLHEP::HepVector::operator=

(

const HepVector &

hm2

)

Definition at line 468 of file Vector.cc.

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

operator[]() [1/2]

double & CLHEP::HepVector::operator[] ( int  row )

inline

operator[]() [2/2]

const double & CLHEP::HepVector::operator[] ( int  row ) const

inline

sub() [1/3]

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

Definition at line 167 of file Vector.cc.

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

sub() [2/3]

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

Definition at line 150 of file Vector.cc.

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

Referenced by main(), and vector_test().

sub() [3/3]

void CLHEP::HepVector::sub	(	int	row,
		const HepVector &	v1
	)

Definition at line 179 of file Vector.cc.

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

T()

HepMatrix CLHEP::HepVector::T

(

)

const

Definition at line 530 of file Vector.cc.

{
#else
{
  HepMatrix mret(1,num_row());
#endif
  mret.m = m;
  return mret;
}

Referenced by main(), and symmatrix_test().

Friends And Related Function Documentation

back_solve

void back_solve ( const HepMatrix &  R, HepVector *  b  )

friend

Definition at line 63 of file MatrixLinear.cc.

{
   (*b)(b->num_row()) /= R(b->num_row(),b->num_row());
   int n = R.num_col();
   int nb = b->num_row();
   HepMatrix::mIter br = b->m.begin() + b->num_row() - 2;
   HepMatrix::mcIter Rrr = R.m.begin() + (nb-2) * (n+1);
   for (int r=b->num_row()-1;r>=1;--r) {
      HepMatrix::mIter bc = br+1;
      HepMatrix::mcIter Rrc = Rrr+1;
      for (int c=r+1;c<=b->num_row();c++) {
     (*br)-=(*(Rrc++))*(*(bc++));
      }
      (*br)/=(*Rrr);
      if(r>1) {
     br--;
     Rrr -= n+1;
      }
   }
}

col_house

void col_house ( HepMatrix *  a, const HepMatrix &  v, double  vnormsq, int  row, int  col, int  row_start, int  col_start  )

friend

Definition at line 154 of file MatrixLinear.cc.

                                                          {
   double beta=-2/vnormsq;
 
   // This is a fast way of calculating w=beta*A.sub(row,n,col,n)*v
 
   HepVector w(a->num_col()-col+1,0);
/* not tested */
   HepMatrix::mIter wptr = w.m.begin();
   int n = a->num_col();
   int nv = v.num_col();
   HepMatrix::mIter acrb = a->m.begin() + (col-1) * n + (row-1);
   int c;
   for (c=col;c<=a->num_col();c++) {
      HepMatrix::mcIter vp = v.m.begin() + (row_start-1) * nv + (col_start-1);
      HepMatrix::mcIter acr = acrb;
      for (int r=row;r<=a->num_row();r++) {
     (*wptr)+=(*(acr++))*(*vp);
     vp += nv;
      }
      wptr++;
      if(c<a->num_col()) acrb += n;
   }
   w*=beta;
  
   // Fast way of calculating A.sub=A.sub+w*v.T()
 
   HepMatrix::mIter arcb = a->m.begin() + (row-1) * n + col-1;
   wptr = w.m.begin();
   for (int r=row; r<=a->num_row();r++) {
      HepMatrix::mIter arc = arcb;
      HepMatrix::mcIter vp = v.m.begin() + (row_start-1) * nv + col_start;
      for (c=col;c<=a->num_col();c++) {
     (*(arc++))+=(*vp)*(*wptr);
     vp += nv;
      }
      wptr++;
      if(r<a->num_row()) arcb += n;
   }
}

dot

double dot ( const HepVector &  v1, const HepVector &  v2  )

friend

Definition at line 542 of file Vector.cc.

{
  if(v1.num_row()!=v2.num_row())
     HepGenMatrix::error("v1 and v2 need to be the same size in dot(HepVector, HepVector)");
  double d= 0;
  HepGenMatrix::mcIter a = v1.m.begin();
  HepGenMatrix::mcIter b = v2.m.begin();
  HepGenMatrix::mcIter e = a + v1.num_size();
  for(;a<e;) d += (*(a++)) * (*(b++));
  return d;
}

HepDiagMatrix

friend class HepDiagMatrix

friend

Definition at line 132 of file Vector.h.

HepMatrix

friend class HepMatrix

friend

Definition at line 134 of file Vector.h.

HepSymMatrix

friend class HepSymMatrix

friend

Definition at line 133 of file Vector.h.

house [1/2]

HepVector house ( const HepMatrix &  a, int  row = 1, int  col = 1  )

friend

Definition at line 371 of file MatrixLinear.cc.

{
   HepVector v(a.num_row()-row+1);
/* not tested */
   int n = a.num_col();
   HepMatrix::mcIter aic = a.m.begin() + (row - 1) * n + (col - 1) ;
   HepMatrix::mIter vp = v.m.begin();
   for (int i=row;i<=a.num_row();i++) {
      (*(vp++))=(*aic);
      aic += n;
   }
   v(1)+=sign(a(row,col))*v.norm();
   return v;
}

house [2/2]

HepVector house ( const HepSymMatrix &  a, int  row = 1, int  col = 1  )

friend

Definition at line 353 of file MatrixLinear.cc.

{
   HepVector v(a.num_row()-row+1);
/* not tested */
   HepMatrix::mIter vp = v.m.begin();
   HepMatrix::mcIter aci = a.m.begin() + col * (col - 1) / 2 + row - 1;
   int i;
   for (i=row;i<=col;i++) {
      (*(vp++))=(*(aci++));
   }
   for (;i<=a.num_row();i++) {
      (*(vp++))=(*aci);
      aci += i;
   }
   v(1)+=sign(a(row,col))*v.norm();
   return v;
}

house_with_update

void house_with_update ( HepMatrix *  a, int  row = 1, int  col = 1  )

friend

Definition at line 396 of file MatrixLinear.cc.

{
   HepVector v(a->num_row()-row+1);
/* not tested */
   HepMatrix::mIter vp = v.m.begin();
   int n = a->num_col();
   HepMatrix::mIter arc = a->m.begin() + (row-1) * n + (col-1);
   int r;
   for (r=row;r<=a->num_row();r++) {
      (*(vp++))=(*arc);
      if(r<a->num_row()) arc += n;
   }
   double normsq=v.normsq();
   double norm=sqrt(normsq);
   normsq-=v(1)*v(1);
   v(1)+=sign((*a)(row,col))*norm;
   normsq+=v(1)*v(1);
   (*a)(row,col)=-sign((*a)(row,col))*norm;
   if (row<a->num_row()) {
      arc = a->m.begin() + row * n + (col-1);
      for (r=row+1;r<=a->num_row();r++) {
     (*arc)=0;
     if(r<a->num_row()) arc += n;
      }
      row_house(a,v,normsq,row,col+1);
   }
}

operator* [1/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* [2/4]

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

friend

Definition at line 354 of file Vector.cc.

{
#else
{
  HepVector mret(hm1.num_row());
#endif
  CHK_DIM_1(hm1.num_col(),hm2.num_row(),*);
  HepGenMatrix::mcIter hm1p,hm2p,vp;
  HepGenMatrix::mIter m3p;
  double temp;
  m3p=mret.m.begin();
  for(hm1p=hm1.m.begin();hm1p<hm1.m.begin()+hm1.num_row()*hm1.num_col();hm1p=hm2p)
    {
      temp=0;
      vp=hm2.m.begin();
      hm2p=hm1p;
      while(hm2p<hm1p+hm1.num_col())
    temp+=(*(hm2p++))*(*(vp++));
      *(m3p++)=temp;
    }
  return mret;
}

operator* [3/4]

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

friend

Definition at line 507 of file SymMatrix.cc.

{
#else
{
  HepVector mret(hm1.num_row());
#endif
  CHK_DIM_1(hm1.num_col(),hm2.num_row(),*);
  HepMatrix::mcIter sp,snp,vpt;
  double temp;
  int step,stept;
  HepMatrix::mIter vrp=mret.m.begin();
  for(step=1,snp=hm1.m.begin();step<=hm1.num_row();++step)
    {
      sp=snp;
      vpt=hm2.m.begin();
      snp+=step;
      temp=0;
      while(sp<snp)
    temp+=*(sp++)*(*(vpt++));
      if(step<hm1.num_row()) sp+=step-1;
      for(stept=step+1;stept<=hm1.num_row();stept++)
    { 
      temp+=*sp*(*(vpt++));
      if(stept<hm1.num_row()) sp+=stept;
    }
      *(vrp++)=temp;
    }   // for(step
  return mret;
}

operator* [4/4]

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

friend

Definition at line 379 of file Vector.cc.

{
#else
{
  HepMatrix mret(hm1.num_row(),hm2.num_col());
#endif
  CHK_DIM_1(1,hm2.num_row(),*);
  HepGenMatrix::mcIter hm1p;
  HepMatrix::mcIter hm2p;
  HepMatrix::mIter mrp=mret.m.begin();
  for(hm1p=hm1.m.begin();hm1p<hm1.m.begin()+hm1.num_row();hm1p++)
    for(hm2p=hm2.m.begin();hm2p<hm2.m.begin()+hm2.num_col();hm2p++)
      *(mrp++)=*hm1p*(*hm2p);
  return mret;
}

operator+

HepVector operator+ ( const HepVector &  v1, const HepVector &  v2  )

friend

Definition at line 255 of file Vector.cc.

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

operator-

HepVector operator- ( const HepVector &  v1, const HepVector &  v2  )

friend

Definition at line 299 of file Vector.cc.

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

qr_solve

HepVector qr_solve ( HepMatrix *  A, const HepVector &  b  )

friend

Definition at line 710 of file MatrixLinear.cc.

{
   HepMatrix Q=qr_decomp(A);
   // Quick way to to Q.T*b.
   HepVector b2(Q.num_col(),0);
   HepMatrix::mIter b2r = b2.m.begin();
   HepMatrix::mIter Qr = Q.m.begin();
   int n = Q.num_col();
   for (int r=1;r<=b2.num_row();r++) {
      HepMatrix::mcIter bc = b.m.begin();
      HepMatrix::mIter Qcr = Qr;
      for (int c=1;c<=b.num_row();c++) {
     *b2r += (*Qcr) * (*(bc++));
     if(c<b.num_row()) Qcr += n;
      }
      b2r++;
      Qr++;
   }
   back_solve(*A,&b2);
   return b2;
}

row_house [1/2]

void row_house ( HepMatrix *  a, const HepMatrix &  v, double  vnormsq, int  row, int  col, int  row_start, int  col_start  )

friend

Definition at line 652 of file MatrixLinear.cc.

                                                           {
   double beta=-2/vnormsq;
 
   // This is a fast way of calculating w=beta*A.sub(row,n,col,n).T()*v
   HepVector w(a->num_col()-col+1,0);
   int na = a->num_col();
   int nv = v.num_col();
   HepMatrix::mIter wptr = w.m.begin();
   HepMatrix::mIter arcb = a->m.begin() + (row-1) * na + (col-1);
   HepMatrix::mcIter vpcb = v.m.begin() + (row_start-1) * nv + (col_start-1);
   int c;
   for (c=col;c<=a->num_col();c++) {
      HepMatrix::mIter arc = arcb;
      HepMatrix::mcIter vpc = vpcb;
      for (int r=row;r<=a->num_row();r++) {
     (*wptr)+=(*arc)*(*vpc);
         if(r<a->num_row()) {
       arc += na;
       vpc += nv;
     }
      }
      wptr++;
      arcb++;
   }
   w*=beta;
 
   arcb = a->m.begin() + (row-1) * na + (col-1);
   HepMatrix::mcIter vpc = v.m.begin() + (row_start-1) * nv + (col_start-1);
   for (int r=row; r<=a->num_row();r++) {
      HepMatrix::mIter arc = arcb;
      HepMatrix::mIter wptr2 = w.m.begin();
      for (c=col;c<=a->num_col();c++) {
     (*(arc++))+=(*vpc)*(*(wptr2++));
      }
      if(r<a->num_row()) {
    arcb += na;
    vpc += nv;
      }
   }
}

row_house [2/2]

void row_house ( HepMatrix *  a, const HepVector &  v, double  vnormsq, int  row = 1, int  col = 1  )

friend

Definition at line 613 of file MatrixLinear.cc.

                             {
   double beta=-2/vnormsq;
 
   // This is a fast way of calculating w=beta*A.sub(row,n,col,n).T()*v
 
   HepVector w(a->num_col()-col+1,0);
/* not tested */
   int na = a->num_col();
   HepMatrix::mIter wptr = w.m.begin();
   HepMatrix::mIter arcb = a->m.begin() + (row-1) * na + (col-1);
   int c;
   for (c=col;c<=a->num_col();c++) {
      HepMatrix::mcIter vp = v.m.begin();
      HepMatrix::mIter arc = arcb;
      for (int r=row;r<=a->num_row();r++) {
     (*wptr)+=(*arc)*(*(vp++));
     if(r<a->num_row()) arc += na;
      }
      wptr++;
      arcb++;
   }
   w*=beta;
  
   // Fast way of calculating A.sub=A.sub+v*w.T()
 
   arcb = a->m.begin() + (row-1) * na + (col-1);
   HepMatrix::mcIter vp = v.m.begin();
   for (int r=row; r<=a->num_row();r++) {
      HepMatrix::mIter wptr2 = w.m.begin();
      HepMatrix::mIter arc = arcb;
      for (c=col;c<=a->num_col();c++) {
     (*(arc++))+=(*vp)*(*(wptr2++));
      }
      vp++;
      if(r<a->num_row()) arcb += na;
   }
}

solve

HepVector solve ( const HepMatrix &  a, const HepVector &  v  )

friend

Definition at line 575 of file Vector.cc.

{
#else
{
  HepVector vret(v);
#endif
  static CLHEP_THREAD_LOCAL int max_array = 20;
  static CLHEP_THREAD_LOCAL int *ir = new int [max_array+1];
 
  if(a.ncol != a.nrow)
     HepGenMatrix::error("Matrix::solve Matrix is not NxN");
  if(a.ncol != v.nrow)
     HepGenMatrix::error("Matrix::solve Vector has wrong number of rows");
 
  int n = a.ncol;
  if (n > max_array) {
    delete [] ir;
    max_array = n;
    ir = new int [max_array+1];
  }
  double det;
  HepMatrix mt(a);
  int i = mt.dfact_matrix(det, ir);
  if (i!=0) {
    for (i=1;i<=n;i++) vret(i) = 0;
    return vret;
  }
  double s21, s22;
  int nxch = ir[n];
  if (nxch!=0) {
    for (int hmm=1;hmm<=nxch;hmm++) {
      int ij = ir[hmm];
      i = ij >> 12;
      int j = ij%4096;
      double te = vret(i);
      vret(i) = vret(j);
      vret(j) = te;
    }
  }
  vret(1) = mt(1,1) * vret(1);
  if (n!=1) {
    for (i=2;i<=n;i++) {
      s21 = -vret(i);
      for (int j=1;j<i;j++) {
    s21 += mt(i,j) * vret(j); 
      }
      vret(i) = -mt(i,i)*s21;
    }
    for (i=1;i<n;i++) {
      int nmi = n-i;
      s22 = -vret(nmi);
      for (int j=1;j<=i;j++) {
    s22 += mt(nmi,n-j+1) * vret(n-j+1);
      }
      vret(nmi) = -s22;
    }
  }
  return vret;
}

swap

void swap ( HepVector &  v1, HepVector &  v2  )

friend

tridiagonal

void tridiagonal ( HepSymMatrix *  a, HepMatrix *  hsm  )

friend

Definition at line 777 of file MatrixLinear.cc.

{
   int nh = hsm->num_col();
   for (int k=1;k<=a->num_col()-2;k++) {
      
      // If this row is already in tridiagonal form, we can skip the
      // transformation.
 
      double scale=0;
      HepMatrix::mIter ajk = a->m.begin() + k * (k+5) / 2;
      int j;
      for (j=k+2;j<=a->num_row();j++) {
     scale+=fabs(*ajk);
     if(j<a->num_row()) ajk += j;
      }
      if (scale ==0) {
     HepMatrix::mIter hsmjkp = hsm->m.begin() + k * (nh+1) - 1;
     for (j=k+1;j<=hsm->num_row();j++) {
        *hsmjkp = 0;
        if(j<hsm->num_row()) hsmjkp += nh;
     }
      } else {
     house_with_update2(a,hsm,k+1,k);
     double normsq=0;
     HepMatrix::mIter hsmrptrkp = hsm->m.begin() + k * (nh+1) - 1;
     int rptr;
     for (rptr=k+1;rptr<=hsm->num_row();rptr++) {
        normsq+=(*hsmrptrkp)*(*hsmrptrkp);
        if(rptr<hsm->num_row()) hsmrptrkp += nh;
     }
     HepVector p(a->num_row()-k,0);
     rptr=k+1;
     HepMatrix::mIter pr = p.m.begin();
     int r;
     for (r=1;r<=p.num_row();r++) {
        HepMatrix::mIter hsmcptrkp = hsm->m.begin() + k * (nh+1) - 1;
        int cptr;
        for (cptr=k+1;cptr<=rptr;cptr++) {
           (*pr)+=a->fast(rptr,cptr)*(*hsmcptrkp);
           if(cptr<a->num_col()) hsmcptrkp += nh;
        }
        for (;cptr<=a->num_col();cptr++) {
           (*pr)+=a->fast(cptr,rptr)*(*hsmcptrkp);
           if(cptr<a->num_col()) hsmcptrkp += nh;
        }
        (*pr)*=2/normsq;
        rptr++;
        pr++;
     }
     double pdotv=0;
     pr = p.m.begin();
     hsmrptrkp = hsm->m.begin() + k * (nh+1) - 1;
     for (r=1;r<=p.num_row();r++) {
        pdotv+=(*(pr++))*(*hsmrptrkp);
        if(r<p.num_row()) hsmrptrkp += nh;
     }
     pr = p.m.begin();
     hsmrptrkp = hsm->m.begin() + k * (nh+1) - 1;
     for (r=1;r<=p.num_row();r++) {
        (*(pr++))-=pdotv*(*hsmrptrkp)/normsq;
        if(r<p.num_row()) hsmrptrkp += nh;
     }
     rptr=k+1;
     pr = p.m.begin();
     hsmrptrkp = hsm->m.begin() + k * (nh+1) - 1;
     for (r=1;r<=p.num_row();r++) {
        int cptr=k+1;
        HepMatrix::mIter mpc = p.m.begin();
        HepMatrix::mIter hsmcptrkp = hsm->m.begin() + k * (nh+1) - 1;
        for (int c=1;c<=r;c++) {
           a->fast(rptr,cptr)-= (*hsmrptrkp)*(*(mpc++))+(*pr)*(*hsmcptrkp);
           cptr++;
           if(c<r) hsmcptrkp += nh;
        }
        pr++;
        rptr++;
        if(r<p.num_row()) hsmrptrkp += nh;
     }
      }
   }
}

vT_times_v

HepSymMatrix vT_times_v ( const HepVector &  v )

friend

Definition at line 539 of file SymMatrix.cc.

{
#else
{
  HepSymMatrix mret(v.num_row());
#endif
  HepMatrix::mIter mr=mret.m.begin();
  HepMatrix::mcIter vt1,vt2;
  for(vt1=v.m.begin();vt1<v.m.begin()+v.num_row();vt1++)
    for(vt2=v.m.begin();vt2<=vt1;vt2++)
      *(mr++)=(*vt1)*(*vt2);
  return mret;
}

---

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

• Vector.h
• Vector.cc