VtVertex.hh File Reference

#include <iosfwd>
#include <vector>
#include "CMatrix.hh"
#include "VtVector.hh"
#include "VtTrack.hh"
#include "VtMassC_list.hh"
#include "VtVertex.icc"

Include dependency graph for VtVertex.hh:

Go to the source code of this file.

Classes
class	VERTEX::Vertex

Namespaces
	VERTEX

Typedefs
typedef std::vector< Track * >::const_iterator	VERTEX::Track_cit
typedef std::vector< Track * >::iterator	VERTEX::Track_it
typedef std::vector< Track * >	VERTEX::Track_v
typedef std::vector< double >	VERTEX::Vector_d

Functions
void	Dinv1 (int *n, double *a, int *idim, double *ir, int *ifail)
std::ostream &	VERTEX::operator<< (std::ostream &s, const Vertex &a)

Function Documentation

Dinv1()

void Dinv1	(	int *	n,
		double *	a,
		int *	idim,
		double *	ir,
		int *	ifail
	)

{
      static double det = 0;
 
      *ifail = 0;
 
      if (*n < 1 || *n > *idim) {
    *ifail = -1;
    return;
      }
 
      if ( *n == 1)
      {
    if (a[0] == 0.) {
        *ifail = -1;
            return;
    }
    a[0] = 1. / a[0];
    return;
      }
 
      if ( *n == 2)
      {
    det = a[0] * a[3] - a[2] * a[1];
    
    if (det == 0.) { *ifail = -1; return; }
 
    double s = 1. / det;
    double c11 = s * a[3];
 
    a[2] = -s * a[2];
    a[1] = -s * a[1];
    a[3] =  s * a[0];
    a[0] = c11;
    return;
      }
 
 
 
      if ( *n == 3 )
      {
    static double t1, t2, t3, temp, s;
    static double c11, c12, c13, c21, c22, c23, c31, c32, c33;
  
    /*     COMPUTE COFACTORS. */
    c11 = a[4] * a[8] - a[7] * a[5];
    c12 = a[7] * a[2] - a[1] * a[8];
    c13 = a[1] * a[5] - a[4] * a[2];
    c21 = a[5] * a[6] - a[8] * a[3];
    c22 = a[8] * a[0] - a[2] * a[6];
    c23 = a[2] * a[3] - a[5] * a[0];
    c31 = a[3] * a[7] - a[6] * a[4];
    c32 = a[6] * a[1] - a[0] * a[7];
    c33 = a[0] * a[4] - a[3] * a[1];
 
    t1 = fabs(a[0]);
    t2 = fabs(a[1]);
    t3 = fabs(a[2]);
    
    /*     (SET TEMP=PIVOT AND DET=PIVOT*DET.) */
    if(t1 < t2) {
            if (t3 < t2) {
        /*        (PIVOT IS A21) */
        temp = a[1];
        det = c13 * c32 - c12 * c33;
            } else {
        /*     (PIVOT IS A31) */
        temp = a[2];
        det = c23 * c12 - c22 * c13;
            }
    } else {
            if(t3 < t1) {
        /*     (PIVOT IS A11) */
        temp = a[0];
        det = c22 * c33 - c23 * c32;
            } else {
        /*     (PIVOT IS A31) */
        temp = a[2];
        det = c23 * c12 - c22 * c13;
            }
    }
 
    if (det == 0.) {
        *ifail = -1;
            return;
    }
 
    /*     SET ELEMENTS OF INVERSE IN A. */
    s = temp / det;
    a[0] = s * c11;
    a[3] = s * c21;
    a[6] = s * c31;
    a[1] = s * c12;
    a[4] = s * c22;
    a[7] = s * c32;
    a[2] = s * c13;
    a[5] = s * c23;
    a[8] = s * c33;
    return;
      }
 
      /* Initialized data */
 
      static int work[100], jfail;
      for(int i=0; i<*n; ++i) work[i] = 0;
 
      /* Function Body */
      
      /*  N.GT.3 CASES.  FACTORIZE MATRIX AND INVERT. */
      Dfactir1(n, a, idim, work, ifail, &det, &jfail);
      if (*ifail != 0)
      {
    cerr << "Dfactir failed!!" << endl;
    return;
      }
      Dfinv1(n, a, idim, work, ifail);
      return;
} // Dinv1