Ext_err_valid.cxx File Reference

#include "CLHEP/Matrix/SymMatrix.h"
#include <iostream>

Go to the source code of this file.

Functions
bool	Ext_err_valid (bool msg, HepSymMatrix &error, const int dimension)
bool	Ext_err_valid (bool msg, const HepSymMatrix &error, const int dimension)

Function Documentation

Ext_err_valid() [1/2]

bool Ext_err_valid	(	bool	msg,
		const HepSymMatrix &	error,
		const int	dimension )

Definition at line 96 of file Ext_err_valid.cxx.

                                                                               {
  bool   valid( 1 );
  double trace( 0 );
 
  for ( int i = 1; i <= dimension; i++ )
  {
    double elem( error( i, i ) );
    trace += elem;
    if ( elem < 0.0 )
    {
      valid = 0;
      if ( msg )
      {
        std::cout << "%ERROR detected at Ext_err_valid: error matrix: error(" << i << "," << i
                  << ")= " << elem << " < 0.0." << std::endl;
      }
    }
    else if ( elem > Large )
    {
      valid = 0;
      if ( msg )
      {
        std::cout << "%ERROR detected at Ext_err_valid: error matrix: error(" << i << "," << i
                  << ")= " << elem << " > " << Large << std::endl;
      }
    }
  }
  if ( !trace ) valid = 0;
  return ( valid );
}

Ext_err_valid() [2/2]

bool Ext_err_valid	(	bool	msg,
		HepSymMatrix &	error,
		const int	dimension )

Definition at line 58 of file Ext_err_valid.cxx.

                                                                         {
  bool   valid( 1 );
  double trace( 0 );
 
  for ( int i = 1; i <= dimension; i++ )
  {
    double elem( error( i, i ) );
    trace += elem;
    if ( elem < 0.0 )
    {
      valid = 0;
      if ( msg )
      {
        std::cout << "%ERROR detected at Ext_err_valid: error(" << i << "," << i
                  << ") = " << elem << " < 0.0. "
                  << "Force to 0.0." << std::endl;
      }
      error( i, i ) = 0.0;
    }
    else if ( elem > Large )
    {
      valid = 0;
      if ( msg )
      {
        std::cout << "%ERROR detected at Ext_err_valid: error(" << i << "," << i
                  << ") = " << elem << " > " << Large << ". Force to " << Large << std::endl;
      }
      error( i, i ) = Large;
    }
  }
  if ( !trace ) valid = 0;
  return ( valid );
}

Referenced by Ext_errmx::valid().