testBug7328.cc File Reference

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

Go to the source code of this file.

Macros
#define	DET_ALSO

Functions
int	test_inversion (int N)
void	heapAddresses (double *&hNew, double *&hMalloc, double *&hNew10000, double *&hMalloc80000)
int	checkHeap (double *&hNew, double *&hMalloc, double *&hNew10000, double *&hMalloc80000, double *&xhNew, double *&xhMalloc, double *&xhNew10000, double *&xhMalloc80000)
int	main (int, char **)

Macro Definition Documentation

DET_ALSO

#define DET_ALSO

Function Documentation

checkHeap()

int checkHeap	(	double *&	hNew,
		double *&	hMalloc,
		double *&	hNew10000,
		double *&	hMalloc80000,
		double *&	xhNew,
		double *&	xhMalloc,
		double *&	xhNew10000,
		double *&	xhMalloc80000
	)

Definition at line 99 of file testBug7328.cc.

                                  {
  int ret = 0;
  if (hNew != xhNew) {
    std::cout<< "Leak:\n"
    << "xhNew - hNew = " <<  xhNew - hNew << "\n";
    ret |= 1;
  }
  if (hMalloc != xhMalloc) {
    std::cout<< "Leak:\n"
    << "xhMalloc - hMalloc = " <<  xhMalloc - hMalloc << "\n";
    ret |= 2;
  }
  if (hNew10000 != xhNew10000) {
    std::cout<< "Leak:\n"
    << "xhNew10000 - hNew10000 = " <<  xhNew10000 - hNew10000 << "\n";
    ret |= 4;
  }
  if (hMalloc80000 != xhMalloc80000) {
    std::cout<< "Leak:\n"
    << "xhMalloc80000 - hMalloc80000 = " <<  xhMalloc80000 -hMalloc80000 
    << "\n";
    ret |= 8;
  }
  return ret; 
}

Referenced by main().

heapAddresses()

void heapAddresses	(	double *&	hNew,
		double *&	hMalloc,
		double *&	hNew10000,
		double *&	hMalloc80000
	)

Definition at line 83 of file testBug7328.cc.

                                      {
  hNew = new double;
  hMalloc =  (double*) malloc(sizeof(double));          
  hNew10000 = new double[10000];
  hMalloc80000 =  (double*) malloc(10000*sizeof(double));
//  std::cout << std::hex << hNew << " " << hMalloc<< " " 
//            << hNew10000 << " " << hMalloc80000 << std::endl;
  free (hMalloc80000);
  delete[] hNew10000;
  free (hMalloc);
  delete hNew;
}

Referenced by main().

main()

int main	(	int	,
		char **
	)

Definition at line 132 of file testBug7328.cc.

                       {
  int ret=0;
  int rhp;
  int i,j;
  for ( i = 1; i <= 50; i++) {
    ret = test_inversion(i);
    if (ret) return ret;
  }
  double *hNew, *hMalloc, *hNew10000, *hMalloc80000;
  double *xhNew, *xhMalloc, *xhNew10000, *xhMalloc80000;
  
  for (int count=0; count < 2; ++count) {
 
    int n1 = 400;
    int n2 = 25;
    heapAddresses ( hNew, hMalloc, hNew10000, hMalloc80000 );
    for (i=0; i<n1; i++) {
      for (j=1; j <= n2; j++) {
    ret = test_inversion(j);
    if (ret) return ret;
      }
    }
    heapAddresses ( xhNew, xhMalloc, xhNew10000, xhMalloc80000 );
    rhp = 0;
    if(count != 0) rhp =  checkHeap ( hNew,  hMalloc,  hNew10000,  hMalloc80000,
                      xhNew, xhMalloc, xhNew10000, xhMalloc80000 );
    if (rhp) std::cout << "Above Leak is after " << n1*n2 << " test inversions\n";
    ret |= rhp;
 
    heapAddresses ( hNew, hMalloc, hNew10000, hMalloc80000 );
    for (i=0; i<2; i++) {
      for (j=1; j < 20; j++) {
    rhp = test_inversion(25+2*j);
    if (rhp) return rhp;
      }
    }
    heapAddresses ( xhNew, xhMalloc, xhNew10000, xhMalloc80000 );
    rhp = 0;
    if(count != 0) rhp =  checkHeap ( hNew,  hMalloc,  hNew10000,  hMalloc80000,
                      xhNew, xhMalloc, xhNew10000, xhMalloc80000 );
    if (rhp) std::cout << "Leak after big inversions\n";
    ret |= rhp;
  } 
  return ret;
}

test_inversion()

int test_inversion

(

int

N

)

Definition at line 22 of file testBug7328.cc.

                           {
 
  int i,j;
  CLHEP::HepSymMatrix S(N,0); 
  for(i=1;i<=N;++i) { 
    for(j=1;j<=N;++j) { 
      if(i<=j) { 
    S (i,j) = (10.0*i+j)/10;
      } 
    } 
  }
  CLHEP::HepSymMatrix SS(N,0); 
  SS = S;
  int ierr = 0;
  SS.invert(ierr);
  if (ierr) {
    std::cout<<"SS.invert failed!!!!  N = " << N << 
                " ierr =  "<< ierr <<std::endl; 
    return 2+100*N;
  }
  CLHEP::HepMatrix SI(N,N,0);
  CLHEP::HepMatrix MS(N,N,0);
  CLHEP::HepMatrix MSS(N,N,0); 
  MS = S;
  MSS = SS;
  SI = MSS*MS;
  for(i=1;i<=N;++i) { 
    for(j=1;j<=N;++j) { 
      if(i!=j) { 
    if (fabs(SI(i,j)) > 1.0e-6) {
          std::cout<<"SS.invert incorrect  N = " << N << 
                " error =  "<< fabs(SI(i,j)) <<std::endl; 
      return 3+100*N;
    }
      } 
      if(i==j) { 
    if (fabs(1-SI(i,j)) > 1.0e-6) {
          std::cout<<"SS.invert incorrect  N = " << N << 
                " error =  "<< fabs(1-SI(i,j)) <<std::endl; 
      return 4+100*N;
    }
      } 
    } 
  }
#define DET_ALSO
#ifdef  DET_ALSO
  double detS  =  S.determinant();
//    std::cout<<"Determinant   N = " << N << 
//      "  =  " << detS <<std::endl;   
  double detSS = SS.determinant();
//    std::cout<<"Determinant Inverse  N = " << N << 
//      "  =  " << detSS <<std::endl; 
  if (fabs((detS-1.0/detSS)/detS) > 1.0e-6) {
    std::cout<<"Determinant incorrect  N = " << N << 
                " error =  " << fabs((detS-1.0/detSS)/detS) <<std::endl; 
      return 5+100*N;
  }  
#endif
  return 0;
}

Referenced by main().