---

SprMatrix.cc File Reference

#include <cfloat>
#include <cmath>
#include <cstdlib>
#include "SprMatrix.hh"
#include "SprSymMatrix.hh"
#include "SprVector.hh"

Go to the source code of this file.

Defines
#define	CHK_DIM_1(c1, r2, fun)
#define	CHK_DIM_2(r1, r2, c1, c2, fun)
#define	SIMPLE_BOP(OPER)
#define	SIMPLE_TOP(OPER)
#define	SIMPLE_UOP(OPER)
Functions
SprMatrix	dsum (const SprMatrix &m1, const SprMatrix &m2)
SprMatrix	operator * (const SprMatrix &m1, const SprMatrix &m2)
SprMatrix	operator * (double t, const SprMatrix &m1)
SprMatrix	operator * (const SprMatrix &m1, double t)
SprMatrix	operator+ (const SprMatrix &m1, const SprMatrix &m2)
SprMatrix	operator- (const SprMatrix &m1, const SprMatrix &m2)
SprMatrix	operator/ (const SprMatrix &m1, double t)
std::ostream &	operator<< (std::ostream &s, const SprMatrix &q)

---

Define Documentation

#define CHK_DIM_1	(	c1,
		r2,
		fun		)

Value:

if (c1!=r2) { \
     SprGenMatrix::error("Range error in Matrix function " #fun "(2)."); \
   }

Definition at line 74 of file SprMatrix.cc.

Referenced by operator *(), operator+(), and operator-().

#define CHK_DIM_2	(	r1,
		r2,
		c1,
		c2,
		fun		)

Value:

if (r1!=r2 || c1!=c2)  { \
     SprGenMatrix::error("Range error in Matrix function " #fun "(1)."); \
   }

Definition at line 69 of file SprMatrix.cc.

Referenced by operator+(), and operator-().

#define SIMPLE_BOP

(

OPER

)

Value:

register SprMatrix::mIter a=m.begin();                      \
   register SprMatrix::mcIter b=m2.m.begin();                  \
   register SprMatrix::mIter e=m.end();                        \
   for(;a!=e; a++, b++) (*a) OPER (*b);

Definition at line 54 of file SprMatrix.cc.

#define SIMPLE_TOP

(

OPER

)

Value:

register SprMatrix::mcIter a=m1.m.begin();       \
   register SprMatrix::mcIter b=m2.m.begin();       \
   register SprMatrix::mIter t=mret.m.begin();      \
   register SprMatrix::mcIter e=m1.m.end();         \
   for(;a!=e; a++, b++, t++) (*t) = (*a) OPER (*b);

Definition at line 60 of file SprMatrix.cc.

Referenced by operator+(), and operator-().

#define SIMPLE_UOP

(

OPER

)

Value:

register mIter a=m.begin();                      \
   register mIter e=m.end();                        \
   for(;a!=e; a++) (*a) OPER t;

Definition at line 49 of file SprMatrix.cc.

---

Function Documentation

SprMatrix dsum	(	const SprMatrix &	m1,
		const SprMatrix &	m2
	)

Definition at line 209 of file SprMatrix.cc.

Referenced by hcaltb::HcalTBTDCUnpacker::reconstructWC().

{
  SprMatrix mret(m1.num_row() + m2.num_row(), m1.num_col() + m2.num_col(),
                 0);
  mret.sub(1,1,m1);
  mret.sub(m1.num_row()+1,m1.num_col()+1,m2);
  return mret;
}

SprMatrix operator *	(	const SprMatrix &	m1,
		const SprMatrix &	m2
	)

Definition at line 282 of file SprMatrix.cc.

References CHK_DIM_1, i, j, and pyDBSRunClass::temp.

{
  // initialize matrix to 0.0
  SprMatrix mret(m1.nrow,m2.ncol,0);
  CHK_DIM_1(m1.ncol,m2.nrow,*);

  int m1cols = m1.ncol;
  int m2cols = m2.ncol;

  for (int i=0; i<m1.nrow; i++)
  {
     for (int j=0; j<m1cols; j++) 
     {
        register double temp = m1.m[i*m1cols+j];
        register SprMatrix::mIter pt = mret.m.begin() + i*m2cols;
        
        // Loop over k (the column index in matrix m2)
        register SprMatrix::mcIter pb = m2.m.begin() + m2cols*j;
        const SprMatrix::mcIter pblast = pb + m2cols;
        while (pb < pblast)
        {
           (*pt) += temp * (*pb);
           pb++;
           pt++;
        }
     }
  }

  return mret;
}

SprMatrix operator *	(	double	t,
		const SprMatrix &	m1
	)

Definition at line 275 of file SprMatrix.cc.

{
  SprMatrix mret(m1);
  mret *= t;
  return mret;
}

SprMatrix operator *	(	const SprMatrix &	m1,
		double	t
	)

Definition at line 268 of file SprMatrix.cc.

{
  SprMatrix mret(m1);
  mret *= t;
  return mret;
}

SprMatrix operator+	(	const SprMatrix &	m1,
		const SprMatrix &	m2
	)

Definition at line 234 of file SprMatrix.cc.

References CHK_DIM_2, and SIMPLE_TOP.

{
  SprMatrix mret(m1.nrow, m1.ncol);
  CHK_DIM_2(m1.num_row(),m2.num_row(), m1.num_col(),m2.num_col(),+);
  SIMPLE_TOP(+)
  return mret;
}

SprMatrix operator-	(	const SprMatrix &	m1,
		const SprMatrix &	m2
	)

Definition at line 246 of file SprMatrix.cc.

References CHK_DIM_2, and SIMPLE_TOP.

{
  SprMatrix mret(m1.num_row(), m1.num_col());
  CHK_DIM_2(m1.num_row(),m2.num_row(),
                         m1.num_col(),m2.num_col(),-);
  SIMPLE_TOP(-)
  return mret;
}

SprMatrix operator/	(	const SprMatrix &	m1,
		double	t
	)

Definition at line 260 of file SprMatrix.cc.

{
  SprMatrix mret(m1);
  mret /= t;
  return mret;
}

std::ostream& operator<<	(	std::ostream &	s,
		const SprMatrix &	q
	)

Definition at line 361 of file SprMatrix.cc.

References lat::endl(), icol, and width.

{
  s << "\n";
/* Fixed format needs 3 extra characters for field, while scientific needs 7 */
  int width;
  if(s.flags() & std::ios::fixed)
    width = s.precision()+3;
  else
    width = s.precision()+7;
  for(int irow = 1; irow<= q.num_row(); irow++)
    {
      for(int icol = 1; icol <= q.num_col(); icol++)
        {
          s.width(width);
          s << q(irow,icol) << " ";
        }
      s << std::endl;
    }
  return s;
}

---