G4ErrorMatrix.hh File Reference

#include <vector>
#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 > ::iterator	G4ErrorMatrixIter
typedef std::vector< G4double > ::const_iterator	G4ErrorMatrixConstIter

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

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

Definition at line 44 of file G4ErrorMatrix.hh.

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

Definition at line 41 of file G4ErrorMatrix.hh.

Function Documentation

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 &
	)

Definition at line 201 of file G4ErrorMatrix.cc.

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

Referenced by gadap(), gadap2(), and gadapf().

{
  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;
}

void givens	(	G4double	a,
		G4double	b,
		G4double *	c,
		G4double *	s
	)

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

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

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

Definition at line 272 of file G4ErrorMatrix.cc.

References CHK_DIM_1.

{
  // 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)
        {
           (*pt) += temp * (*pb);
           pb++;
           pt++;
        }
     }
  }
  return mret;
}

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

Definition at line 265 of file G4ErrorMatrix.cc.

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

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

Definition at line 258 of file G4ErrorMatrix.cc.

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

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

Definition at line 225 of file G4ErrorMatrix.cc.

References CHK_DIM_2, G4ErrorMatrix::num_col(), G4ErrorMatrix::num_row(), and SIMPLE_TOP.

{
  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;
}

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

Definition at line 237 of file G4ErrorMatrix.cc.

References CHK_DIM_2, G4ErrorMatrix::num_col(), G4ErrorMatrix::num_row(), and SIMPLE_TOP.

{
  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;
}

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

Definition at line 251 of file G4ErrorMatrix.cc.

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

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

Definition at line 350 of file G4ErrorMatrix.cc.

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

{
  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;
}

void qr_decomp	(	G4ErrorMatrix *	A,
		G4ErrorMatrix *	hsm
	)

G4ErrorMatrix qr_decomp

(

G4ErrorMatrix *

A

)

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
	)