Ext_err_valid.cxx File Reference

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

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 92 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().