G4ErrorMatrix.hh File Reference

#include <vector>
#include "G4Types.hh"
#include "G4ErrorMatrix.icc"

Go to the source code of this file.

Data Structures
class	G4ErrorMatrix
class	G4ErrorMatrix::G4ErrorMatrix_row
class	G4ErrorMatrix::G4ErrorMatrix_row_const

Typedefs
typedef std::vector< G4double >::const_iterator	G4ErrorMatrixConstIter
typedef std::vector< G4double >::iterator	G4ErrorMatrixIter

Functions
void	back_solve (const G4ErrorMatrix &R, G4ErrorMatrix *b)
void	col_givens (G4ErrorMatrix *A, G4double c, G4double s, G4int k1, G4int k2, G4int row_min=1, G4int row_max=0)
void	col_house (G4ErrorMatrix *a, const G4ErrorMatrix &v, G4double vnormsq, G4int row, G4int col, G4int row_start, G4int col_start)
void	col_house (G4ErrorMatrix *a, const G4ErrorMatrix &v, G4int row, G4int col, G4int row_start, G4int col_start)
G4ErrorMatrix	dsum (const G4ErrorMatrix &, const G4ErrorMatrix &)
void	house_with_update (G4ErrorMatrix *a, G4ErrorMatrix *v, G4int row=1, G4int col=1)
void	house_with_update (G4ErrorMatrix *a, G4int row=1, G4int col=1)
G4ErrorMatrix	operator* (const G4ErrorMatrix &m1, const G4ErrorMatrix &m2)
G4ErrorMatrix	operator* (const G4ErrorMatrix &m1, G4double t)
G4ErrorMatrix	operator* (G4double t, const G4ErrorMatrix &m1)
G4ErrorMatrix	operator+ (const G4ErrorMatrix &m1, const G4ErrorMatrix &m2)
G4ErrorMatrix	operator- (const G4ErrorMatrix &m1, const G4ErrorMatrix &m2)
G4ErrorMatrix	operator/ (const G4ErrorMatrix &m1, G4double t)
std::ostream &	operator<< (std::ostream &s, const G4ErrorMatrix &q)
G4ErrorMatrix	qr_decomp (G4ErrorMatrix *A)
void	qr_decomp (G4ErrorMatrix *A, G4ErrorMatrix *hsm)
G4ErrorMatrix	qr_inverse (const G4ErrorMatrix &A)
G4ErrorMatrix	qr_inverse (G4ErrorMatrix *A)
G4ErrorMatrix	qr_solve (const G4ErrorMatrix &A, const G4ErrorMatrix &b)
G4ErrorMatrix	qr_solve (G4ErrorMatrix *A, const G4ErrorMatrix &b)
void	row_givens (G4ErrorMatrix *A, G4double c, G4double s, G4int k1, G4int k2, G4int col_min=1, G4int col_max=0)
void	row_house (G4ErrorMatrix *a, const G4ErrorMatrix &v, G4double vnormsq, G4int row, G4int col, G4int row_start, G4int col_start)
void	row_house (G4ErrorMatrix *a, const G4ErrorMatrix &v, G4int row, G4int col, G4int row_start, G4int col_start)

Typedef Documentation

G4ErrorMatrixConstIter

typedef std::vector<G4double>::const_iterator G4ErrorMatrixConstIter

Definition at line 44 of file G4ErrorMatrix.hh.

G4ErrorMatrixIter

typedef std::vector<G4double>::iterator G4ErrorMatrixIter

Definition at line 43 of file G4ErrorMatrix.hh.

Function Documentation

back_solve()

void back_solve	(	const G4ErrorMatrix &	R,
		G4ErrorMatrix *	b
	)

col_givens()

void col_givens	(	G4ErrorMatrix *	A,
		G4double	c,
		G4double	s,
		G4int	k1,
		G4int	k2,
		G4int	row_min = 1,
		G4int	row_max = 0
	)

col_house() [1/2]

void col_house	(	G4ErrorMatrix *	a,
		const G4ErrorMatrix &	v,
		G4double	vnormsq,
		G4int	row,
		G4int	col,
		G4int	row_start,
		G4int	col_start
	)

col_house() [2/2]

void col_house	(	G4ErrorMatrix *	a,
		const G4ErrorMatrix &	v,
		G4int	row,
		G4int	col,
		G4int	row_start,
		G4int	col_start
	)

dsum()

G4ErrorMatrix dsum	(	const G4ErrorMatrix &	mat1,
		const G4ErrorMatrix &	mat2
	)

Definition at line 220 of file G4ErrorMatrix.cc.

{
  G4ErrorMatrix mret(mat1.num_row() + mat2.num_row(),
                     mat1.num_col() + mat2.num_col(), 0);
  mret.sub(1, 1, mat1);
  mret.sub(mat1.num_row() + 1, mat1.num_col() + 1, mat2);
  return mret;
}

References G4ErrorMatrix::num_col(), G4ErrorMatrix::num_row(), and G4ErrorMatrix::sub().

house_with_update() [1/2]

void house_with_update	(	G4ErrorMatrix *	a,
		G4ErrorMatrix *	v,
		G4int	row = 1,
		G4int	col = 1
	)

house_with_update() [2/2]

void house_with_update	(	G4ErrorMatrix *	a,
		G4int	row = 1,
		G4int	col = 1
	)

operator*() [1/3]

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

Definition at line 291 of file G4ErrorMatrix.cc.

{
  // initialize matrix to 0.0
  G4ErrorMatrix mret(mat1.nrow, mat2.ncol, 0);
  CHK_DIM_1(mat1.ncol, mat2.nrow, *);
 
  G4int m1cols = mat1.ncol;
  G4int m2cols = mat2.ncol;
 
  for(G4int i = 0; i < mat1.nrow; i++)
  {
    for(G4int j = 0; j < m1cols; j++)
    {
      G4double temp        = mat1.m[i * m1cols + j];
      G4ErrorMatrixIter pt = mret.m.begin() + i * m2cols;
 
      // Loop over k (the column index in matrix mat2)
      G4ErrorMatrixConstIter pb           = mat2.m.begin() + m2cols * j;
      const G4ErrorMatrixConstIter pblast = pb + m2cols;
      while(pb < pblast)  // Loop checking, 06.08.2015, G.Cosmo
      {
        (*pt) += temp * (*pb);
        pb++;
        pt++;
      }
    }
  }
  return mret;
}

operator*() [2/3]

G4ErrorMatrix operator*	(	const G4ErrorMatrix &	m1,
		G4double	t
	)

Definition at line 277 of file G4ErrorMatrix.cc.

{
  G4ErrorMatrix mret(mat1);
  mret *= t;
  return mret;
}

operator*() [3/3]

G4ErrorMatrix operator*	(	G4double	t,
		const G4ErrorMatrix &	m1
	)

Definition at line 284 of file G4ErrorMatrix.cc.

{
  G4ErrorMatrix mret(mat1);
  mret *= t;
  return mret;
}

operator+()

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

Definition at line 245 of file G4ErrorMatrix.cc.

{
  G4ErrorMatrix mret(mat1.nrow, mat1.ncol);
  CHK_DIM_2(mat1.num_row(), mat2.num_row(), mat1.num_col(), mat2.num_col(), +);
  SIMPLE_TOP(+)
  return mret;
}

operator-()

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

Definition at line 257 of file G4ErrorMatrix.cc.

{
  G4ErrorMatrix mret(mat1.num_row(), mat1.num_col());
  CHK_DIM_2(mat1.num_row(), mat2.num_row(), mat1.num_col(), mat2.num_col(), -);
  SIMPLE_TOP(-)
  return mret;
}

operator/()

G4ErrorMatrix operator/	(	const G4ErrorMatrix &	m1,
		G4double	t
	)

Definition at line 270 of file G4ErrorMatrix.cc.

{
  G4ErrorMatrix mret(mat1);
  mret /= t;
  return mret;
}

operator<<()

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

Definition at line 371 of file G4ErrorMatrix.cc.

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

References G4endl, G4ErrorMatrix::num_col(), and G4ErrorMatrix::num_row().

qr_decomp() [1/2]

G4ErrorMatrix qr_decomp

(

G4ErrorMatrix *

A

)

qr_decomp() [2/2]

void qr_decomp	(	G4ErrorMatrix *	A,
		G4ErrorMatrix *	hsm
	)

qr_inverse() [1/2]

G4ErrorMatrix qr_inverse

(

const G4ErrorMatrix &

A

)

qr_inverse() [2/2]

G4ErrorMatrix qr_inverse

(

G4ErrorMatrix *

A

)

qr_solve() [1/2]

G4ErrorMatrix qr_solve	(	const G4ErrorMatrix &	A,
		const G4ErrorMatrix &	b
	)

qr_solve() [2/2]

G4ErrorMatrix qr_solve	(	G4ErrorMatrix *	A,
		const G4ErrorMatrix &	b
	)

row_givens()

void row_givens	(	G4ErrorMatrix *	A,
		G4double	c,
		G4double	s,
		G4int	k1,
		G4int	k2,
		G4int	col_min = 1,
		G4int	col_max = 0
	)

row_house() [1/2]

void row_house	(	G4ErrorMatrix *	a,
		const G4ErrorMatrix &	v,
		G4double	vnormsq,
		G4int	row,
		G4int	col,
		G4int	row_start,
		G4int	col_start
	)

row_house() [2/2]

void row_house	(	G4ErrorMatrix *	a,
		const G4ErrorMatrix &	v,
		G4int	row,
		G4int	col,
		G4int	row_start,
		G4int	col_start
	)