ROOT logo
// @(#)root/base:$Id: TMatrixTSymCramerInv.cxx 20882 2007-11-19 11:31:26Z rdm $
// Authors: Fons Rademakers, Eddy Offermann  Oct 2004

/*************************************************************************
 * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers.               *
 * All rights reserved.                                                  *
 *                                                                       *
 * For the licensing terms see $ROOTSYS/LICENSE.                         *
 * For the list of contributors see $ROOTSYS/README/CREDITS.             *
 *************************************************************************/

//////////////////////////////////////////////////////////////////////////
//                                                                      //
// TMatrixTSymCramerInv                                                 //
//                                                                      //
// Encapsulate templates of Cramer Inversion routines.                  //
//                                                                      //
// The 4x4, 5x5 and 6x6 are adapted from routines written by            //
// Mark Fischler and Steven Haywood as part of the CLHEP package        //
//                                                                      //
// Although for sizes <= 6x6 the Cramer Inversion has a gain in speed   //
// compared to factorization schemes (like LU) , one pays a price in    //
// accuracy  .                                                          //
//                                                                      //
// For Example:                                                         //
//  H * H^-1 = U, where H is a 5x5 Hilbert matrix                       //
//                      U is a 5x5 Unity matrix                         //
//                                                                      //
// LU    : |U_jk| < 10e-13 for  j!=k                                    //
// Cramer: |U_jk| < 10e-7  for  j!=k                                    //
//                                                                      //
//  however Cramer algorithm is about 10 (!) times faster               //
//////////////////////////////////////////////////////////////////////////

#include "TMatrixTSymCramerInv.h"

#if !defined(R__ALPHA) && !defined(R__SOLARIS) && !defined(R__ACC) && !defined(R__FBSD)
NamespaceImp(TMatrixTSymCramerInv);
#endif

//______________________________________________________________________________
template<class Element> 
Bool_t TMatrixTSymCramerInv::Inv2x2(TMatrixTSym<Element> &m,Double_t *determ)
{
   if (m.GetNrows() != 2) {
      Error("Inv2x2","matrix should be square 2x2");
      return kFALSE;
   }

   Element *pM = m.GetMatrixArray();

   const Double_t det = pM[0] * pM[3] - pM[1] * pM[1];

   if (determ)
      *determ = det;

   if ( det == 0 ) {
      Error("Inv2x2","matrix is singular");
      return kFALSE;
   }

   const Double_t tmp1 =   pM[3] / det;
   pM[3] = pM[0] / det;
   pM[2] = pM[1] = -pM[1] / det;
   pM[0] = tmp1;

   return kTRUE;
}

//______________________________________________________________________________
template<class Element> 
Bool_t TMatrixTSymCramerInv::Inv3x3(TMatrixTSym<Element> &m,Double_t *determ)
{
   if (m.GetNrows() != 3) {
      Error("Inv3x3","matrix should be square 3x3");
      return kFALSE;
   }

   Element *pM = m.GetMatrixArray();

   const Double_t c00 = pM[4] * pM[8] - pM[5] * pM[5];
   const Double_t c01 = pM[5] * pM[2] - pM[1] * pM[8];
   const Double_t c02 = pM[1] * pM[5] - pM[4] * pM[2];
   const Double_t c11 = pM[8] * pM[0] - pM[2] * pM[2];
   const Double_t c12 = pM[2] * pM[1] - pM[5] * pM[0];
   const Double_t c22 = pM[0] * pM[4] - pM[1] * pM[1];

   const Double_t t0  = TMath::Abs(pM[0]);
   const Double_t t1  = TMath::Abs(pM[1]);
   const Double_t t2  = TMath::Abs(pM[2]);

   Double_t det;
   Double_t tmp;

   if (t0 >= t1) {
      if (t2 >= t0) {
         tmp = pM[2];
         det = c12*c01-c11*c02;
      } else {
         tmp = pM[0];
         det = c11*c22-c12*c12;
      }
   } else if (t2 >= t1) {
      tmp = pM[2];
      det = c12*c01-c11*c02;
   } else {
      tmp = pM[1];
      det = c02*c12-c01*c22;
   }

   if ( det == 0 || tmp == 0) {
      Error("Inv3x3","matrix is singular");
      return kFALSE;
   }

   Double_t s = tmp/det;
   if (determ)
      *determ = 1./s;

   pM[0] = s*c00;
   pM[1] = s*c01;
   pM[2] = s*c02;
   pM[3] = pM[1];
   pM[4] = s*c11;
   pM[5] = s*c12;
   pM[6] = pM[2];
   pM[7] = pM[5];
   pM[8] = s*c22;
  
   return kTRUE;
}

// SFij are indices for a 4x4 symmetric matrix.

#define SF00 0
#define SF01 1
#define SF02 2
#define SF03 3

#define SF10 1
#define SF11 5
#define SF12 6
#define SF13 7

#define SF20 2
#define SF21 6
#define SF22 10
#define SF23 11

#define SF30 3
#define SF31 7
#define SF32 11
#define SF33 15

//______________________________________________________________________________
template<class Element> 
Bool_t TMatrixTSymCramerInv::Inv4x4(TMatrixTSym<Element> &m,Double_t *determ)
{
   if (m.GetNrows() != 4) {
      Error("Inv4x4","matrix should be square 4x4");
      return kFALSE;
   }

   Element *pM = m.GetMatrixArray();

   // Find all NECESSARY 2x2 dets:  (14 of them)

   const Double_t mDet2_12_01 = pM[SF10]*pM[SF21] - pM[SF11]*pM[SF20];
   const Double_t mDet2_12_02 = pM[SF10]*pM[SF22] - pM[SF12]*pM[SF20];
   const Double_t mDet2_12_12 = pM[SF11]*pM[SF22] - pM[SF12]*pM[SF21];
   const Double_t mDet2_13_01 = pM[SF10]*pM[SF31] - pM[SF11]*pM[SF30];
   const Double_t mDet2_13_02 = pM[SF10]*pM[SF32] - pM[SF12]*pM[SF30];
   const Double_t mDet2_13_03 = pM[SF10]*pM[SF33] - pM[SF13]*pM[SF30];
   const Double_t mDet2_13_12 = pM[SF11]*pM[SF32] - pM[SF12]*pM[SF31];
   const Double_t mDet2_13_13 = pM[SF11]*pM[SF33] - pM[SF13]*pM[SF31];
   const Double_t mDet2_23_01 = pM[SF20]*pM[SF31] - pM[SF21]*pM[SF30];
   const Double_t mDet2_23_02 = pM[SF20]*pM[SF32] - pM[SF22]*pM[SF30];
   const Double_t mDet2_23_03 = pM[SF20]*pM[SF33] - pM[SF23]*pM[SF30];
   const Double_t mDet2_23_12 = pM[SF21]*pM[SF32] - pM[SF22]*pM[SF31];
   const Double_t mDet2_23_13 = pM[SF21]*pM[SF33] - pM[SF23]*pM[SF31];
   const Double_t mDet2_23_23 = pM[SF22]*pM[SF33] - pM[SF23]*pM[SF32];

  // SFind all NECESSSFRY 3x3 dets:   (10 of them)
  
   const Double_t mDet3_012_012 = pM[SF00]*mDet2_12_12 - pM[SF01]*mDet2_12_02
                                + pM[SF02]*mDet2_12_01;
   const Double_t mDet3_013_012 = pM[SF00]*mDet2_13_12 - pM[SF01]*mDet2_13_02
                                + pM[SF02]*mDet2_13_01;
   const Double_t mDet3_013_013 = pM[SF00]*mDet2_13_13 - pM[SF01]*mDet2_13_03
                                + pM[SF03]*mDet2_13_01;
   const Double_t mDet3_023_012 = pM[SF00]*mDet2_23_12 - pM[SF01]*mDet2_23_02
                                + pM[SF02]*mDet2_23_01;
   const Double_t mDet3_023_013 = pM[SF00]*mDet2_23_13 - pM[SF01]*mDet2_23_03
                                + pM[SF03]*mDet2_23_01;
   const Double_t mDet3_023_023 = pM[SF00]*mDet2_23_23 - pM[SF02]*mDet2_23_03
                                + pM[SF03]*mDet2_23_02;
   const Double_t mDet3_123_012 = pM[SF10]*mDet2_23_12 - pM[SF11]*mDet2_23_02
                                + pM[SF12]*mDet2_23_01;
   const Double_t mDet3_123_013 = pM[SF10]*mDet2_23_13 - pM[SF11]*mDet2_23_03
                                + pM[SF13]*mDet2_23_01;
   const Double_t mDet3_123_023 = pM[SF10]*mDet2_23_23 - pM[SF12]*mDet2_23_03
                                + pM[SF13]*mDet2_23_02;
   const Double_t mDet3_123_123 = pM[SF11]*mDet2_23_23 - pM[SF12]*mDet2_23_13
                                + pM[SF13]*mDet2_23_12;

   // Find the 4x4 det:

   const Double_t det = pM[SF00]*mDet3_123_123 - pM[SF01]*mDet3_123_023
                      + pM[SF02]*mDet3_123_013 - pM[SF03]*mDet3_123_012;

   if (determ)
      *determ = det;

   if ( det == 0 ) {
      Error("Inv4x4","matrix is singular");
      return kFALSE;
   }

   const Double_t oneOverDet = 1.0/det;
   const Double_t mn1OverDet = - oneOverDet;

   pM[SF00] =  mDet3_123_123 * oneOverDet;
   pM[SF01] =  mDet3_123_023 * mn1OverDet;
   pM[SF02] =  mDet3_123_013 * oneOverDet;
   pM[SF03] =  mDet3_123_012 * mn1OverDet;

   pM[SF11] =  mDet3_023_023 * oneOverDet;
   pM[SF12] =  mDet3_023_013 * mn1OverDet;
   pM[SF13] =  mDet3_023_012 * oneOverDet;

   pM[SF22] =  mDet3_013_013 * oneOverDet;
   pM[SF23] =  mDet3_013_012 * mn1OverDet;

   pM[SF33] =  mDet3_012_012 * oneOverDet;

   for (Int_t irow = 0; irow < 4; irow++) {
      const Int_t rowOff1 = irow*4; 
      for (Int_t icol = 0; icol < irow; icol++) {
         const Int_t rowOff2 = icol*4; 
         pM[rowOff1+icol] = pM[rowOff2+irow];
      }
   } 

   return kTRUE;
}

// Mij are indices for a 5x5 matrix.
    
#define SM00 0
#define SM01 1
#define SM02 2
#define SM03 3
#define SM04 4
  
#define SM10 1
#define SM11 6
#define SM12 7
#define SM13 8
#define SM14 9

#define SM20 2
#define SM21 7
#define SM22 12
#define SM23 13
#define SM24 14

#define SM30 3
#define SM31 8
#define SM32 13
#define SM33 18
#define SM34 19

#define SM40 4
#define SM41 9
#define SM42 14
#define SM43 19
#define SM44 24

//______________________________________________________________________________
template<class Element> 
Bool_t TMatrixTSymCramerInv::Inv5x5(TMatrixTSym<Element> &m,Double_t *determ)
{
   if (m.GetNrows() != 5) {
      Error("Inv5x5","matrix should be square 5x5");
      return kFALSE;
   }

   Element *pM = m.GetMatrixArray();

   // Find all NECESSARY 2x2 dets:  (25 of them)

   const Double_t mDet2_23_01 = pM[SM20]*pM[SM31] - pM[SM21]*pM[SM30];
   const Double_t mDet2_23_02 = pM[SM20]*pM[SM32] - pM[SM22]*pM[SM30];
   const Double_t mDet2_23_03 = pM[SM20]*pM[SM33] - pM[SM23]*pM[SM30];
   const Double_t mDet2_23_12 = pM[SM21]*pM[SM32] - pM[SM22]*pM[SM31];
   const Double_t mDet2_23_13 = pM[SM21]*pM[SM33] - pM[SM23]*pM[SM31];
   const Double_t mDet2_23_23 = pM[SM22]*pM[SM33] - pM[SM23]*pM[SM32];
   const Double_t mDet2_24_01 = pM[SM20]*pM[SM41] - pM[SM21]*pM[SM40];
   const Double_t mDet2_24_02 = pM[SM20]*pM[SM42] - pM[SM22]*pM[SM40];
   const Double_t mDet2_24_03 = pM[SM20]*pM[SM43] - pM[SM23]*pM[SM40];
   const Double_t mDet2_24_04 = pM[SM20]*pM[SM44] - pM[SM24]*pM[SM40];
   const Double_t mDet2_24_12 = pM[SM21]*pM[SM42] - pM[SM22]*pM[SM41];
   const Double_t mDet2_24_13 = pM[SM21]*pM[SM43] - pM[SM23]*pM[SM41];
   const Double_t mDet2_24_14 = pM[SM21]*pM[SM44] - pM[SM24]*pM[SM41];
   const Double_t mDet2_24_23 = pM[SM22]*pM[SM43] - pM[SM23]*pM[SM42];
   const Double_t mDet2_24_24 = pM[SM22]*pM[SM44] - pM[SM24]*pM[SM42];
   const Double_t mDet2_34_01 = pM[SM30]*pM[SM41] - pM[SM31]*pM[SM40];
   const Double_t mDet2_34_02 = pM[SM30]*pM[SM42] - pM[SM32]*pM[SM40];
   const Double_t mDet2_34_03 = pM[SM30]*pM[SM43] - pM[SM33]*pM[SM40];
   const Double_t mDet2_34_04 = pM[SM30]*pM[SM44] - pM[SM34]*pM[SM40];
   const Double_t mDet2_34_12 = pM[SM31]*pM[SM42] - pM[SM32]*pM[SM41];
   const Double_t mDet2_34_13 = pM[SM31]*pM[SM43] - pM[SM33]*pM[SM41];
   const Double_t mDet2_34_14 = pM[SM31]*pM[SM44] - pM[SM34]*pM[SM41];
   const Double_t mDet2_34_23 = pM[SM32]*pM[SM43] - pM[SM33]*pM[SM42];
   const Double_t mDet2_34_24 = pM[SM32]*pM[SM44] - pM[SM34]*pM[SM42];
   const Double_t mDet2_34_34 = pM[SM33]*pM[SM44] - pM[SM34]*pM[SM43];

   // Find all NECESSARY 3x3 dets:   (30 of them)

   const Double_t mDet3_123_012 = pM[SM10]*mDet2_23_12 - pM[SM11]*mDet2_23_02 + pM[SM12]*mDet2_23_01;
   const Double_t mDet3_123_013 = pM[SM10]*mDet2_23_13 - pM[SM11]*mDet2_23_03 + pM[SM13]*mDet2_23_01;
   const Double_t mDet3_123_023 = pM[SM10]*mDet2_23_23 - pM[SM12]*mDet2_23_03 + pM[SM13]*mDet2_23_02;
   const Double_t mDet3_123_123 = pM[SM11]*mDet2_23_23 - pM[SM12]*mDet2_23_13 + pM[SM13]*mDet2_23_12;
   const Double_t mDet3_124_012 = pM[SM10]*mDet2_24_12 - pM[SM11]*mDet2_24_02 + pM[SM12]*mDet2_24_01;
   const Double_t mDet3_124_013 = pM[SM10]*mDet2_24_13 - pM[SM11]*mDet2_24_03 + pM[SM13]*mDet2_24_01;
   const Double_t mDet3_124_014 = pM[SM10]*mDet2_24_14 - pM[SM11]*mDet2_24_04 + pM[SM14]*mDet2_24_01;
   const Double_t mDet3_124_023 = pM[SM10]*mDet2_24_23 - pM[SM12]*mDet2_24_03 + pM[SM13]*mDet2_24_02;
   const Double_t mDet3_124_024 = pM[SM10]*mDet2_24_24 - pM[SM12]*mDet2_24_04 + pM[SM14]*mDet2_24_02;
   const Double_t mDet3_124_123 = pM[SM11]*mDet2_24_23 - pM[SM12]*mDet2_24_13 + pM[SM13]*mDet2_24_12;
   const Double_t mDet3_124_124 = pM[SM11]*mDet2_24_24 - pM[SM12]*mDet2_24_14 + pM[SM14]*mDet2_24_12;
   const Double_t mDet3_134_012 = pM[SM10]*mDet2_34_12 - pM[SM11]*mDet2_34_02 + pM[SM12]*mDet2_34_01;
   const Double_t mDet3_134_013 = pM[SM10]*mDet2_34_13 - pM[SM11]*mDet2_34_03 + pM[SM13]*mDet2_34_01;
   const Double_t mDet3_134_014 = pM[SM10]*mDet2_34_14 - pM[SM11]*mDet2_34_04 + pM[SM14]*mDet2_34_01;
   const Double_t mDet3_134_023 = pM[SM10]*mDet2_34_23 - pM[SM12]*mDet2_34_03 + pM[SM13]*mDet2_34_02;
   const Double_t mDet3_134_024 = pM[SM10]*mDet2_34_24 - pM[SM12]*mDet2_34_04 + pM[SM14]*mDet2_34_02;
   const Double_t mDet3_134_034 = pM[SM10]*mDet2_34_34 - pM[SM13]*mDet2_34_04 + pM[SM14]*mDet2_34_03;
   const Double_t mDet3_134_123 = pM[SM11]*mDet2_34_23 - pM[SM12]*mDet2_34_13 + pM[SM13]*mDet2_34_12;
   const Double_t mDet3_134_124 = pM[SM11]*mDet2_34_24 - pM[SM12]*mDet2_34_14 + pM[SM14]*mDet2_34_12;
   const Double_t mDet3_134_134 = pM[SM11]*mDet2_34_34 - pM[SM13]*mDet2_34_14 + pM[SM14]*mDet2_34_13;
   const Double_t mDet3_234_012 = pM[SM20]*mDet2_34_12 - pM[SM21]*mDet2_34_02 + pM[SM22]*mDet2_34_01;
   const Double_t mDet3_234_013 = pM[SM20]*mDet2_34_13 - pM[SM21]*mDet2_34_03 + pM[SM23]*mDet2_34_01;
   const Double_t mDet3_234_014 = pM[SM20]*mDet2_34_14 - pM[SM21]*mDet2_34_04 + pM[SM24]*mDet2_34_01;
   const Double_t mDet3_234_023 = pM[SM20]*mDet2_34_23 - pM[SM22]*mDet2_34_03 + pM[SM23]*mDet2_34_02;
   const Double_t mDet3_234_024 = pM[SM20]*mDet2_34_24 - pM[SM22]*mDet2_34_04 + pM[SM24]*mDet2_34_02;
   const Double_t mDet3_234_034 = pM[SM20]*mDet2_34_34 - pM[SM23]*mDet2_34_04 + pM[SM24]*mDet2_34_03;
   const Double_t mDet3_234_123 = pM[SM21]*mDet2_34_23 - pM[SM22]*mDet2_34_13 + pM[SM23]*mDet2_34_12;
   const Double_t mDet3_234_124 = pM[SM21]*mDet2_34_24 - pM[SM22]*mDet2_34_14 + pM[SM24]*mDet2_34_12;
   const Double_t mDet3_234_134 = pM[SM21]*mDet2_34_34 - pM[SM23]*mDet2_34_14 + pM[SM24]*mDet2_34_13;
   const Double_t mDet3_234_234 = pM[SM22]*mDet2_34_34 - pM[SM23]*mDet2_34_24 + pM[SM24]*mDet2_34_23;

   // Find all NECESSARY 4x4 dets:   (15 of them)

   const Double_t mDet4_0123_0123 = pM[SM00]*mDet3_123_123 - pM[SM01]*mDet3_123_023
                                  + pM[SM02]*mDet3_123_013 - pM[SM03]*mDet3_123_012;
   const Double_t mDet4_0124_0123 = pM[SM00]*mDet3_124_123 - pM[SM01]*mDet3_124_023
                                  + pM[SM02]*mDet3_124_013 - pM[SM03]*mDet3_124_012;
   const Double_t mDet4_0124_0124 = pM[SM00]*mDet3_124_124 - pM[SM01]*mDet3_124_024
                                  + pM[SM02]*mDet3_124_014 - pM[SM04]*mDet3_124_012;
   const Double_t mDet4_0134_0123 = pM[SM00]*mDet3_134_123 - pM[SM01]*mDet3_134_023
                                  + pM[SM02]*mDet3_134_013 - pM[SM03]*mDet3_134_012;
   const Double_t mDet4_0134_0124 = pM[SM00]*mDet3_134_124 - pM[SM01]*mDet3_134_024
                                  + pM[SM02]*mDet3_134_014 - pM[SM04]*mDet3_134_012;
   const Double_t mDet4_0134_0134 = pM[SM00]*mDet3_134_134 - pM[SM01]*mDet3_134_034
                                  + pM[SM03]*mDet3_134_014 - pM[SM04]*mDet3_134_013;
   const Double_t mDet4_0234_0123 = pM[SM00]*mDet3_234_123 - pM[SM01]*mDet3_234_023
                                  + pM[SM02]*mDet3_234_013 - pM[SM03]*mDet3_234_012;
   const Double_t mDet4_0234_0124 = pM[SM00]*mDet3_234_124 - pM[SM01]*mDet3_234_024
                                  + pM[SM02]*mDet3_234_014 - pM[SM04]*mDet3_234_012;
   const Double_t mDet4_0234_0134 = pM[SM00]*mDet3_234_134 - pM[SM01]*mDet3_234_034
                                  + pM[SM03]*mDet3_234_014 - pM[SM04]*mDet3_234_013;
   const Double_t mDet4_0234_0234 = pM[SM00]*mDet3_234_234 - pM[SM02]*mDet3_234_034
                                  + pM[SM03]*mDet3_234_024 - pM[SM04]*mDet3_234_023;
   const Double_t mDet4_1234_0123 = pM[SM10]*mDet3_234_123 - pM[SM11]*mDet3_234_023
                                  + pM[SM12]*mDet3_234_013 - pM[SM13]*mDet3_234_012;
   const Double_t mDet4_1234_0124 = pM[SM10]*mDet3_234_124 - pM[SM11]*mDet3_234_024
                                  + pM[SM12]*mDet3_234_014 - pM[SM14]*mDet3_234_012;
   const Double_t mDet4_1234_0134 = pM[SM10]*mDet3_234_134 - pM[SM11]*mDet3_234_034
                                  + pM[SM13]*mDet3_234_014 - pM[SM14]*mDet3_234_013;
   const Double_t mDet4_1234_0234 = pM[SM10]*mDet3_234_234 - pM[SM12]*mDet3_234_034
                                  + pM[SM13]*mDet3_234_024 - pM[SM14]*mDet3_234_023;
   const Double_t mDet4_1234_1234 = pM[SM11]*mDet3_234_234 - pM[SM12]*mDet3_234_134
                                  + pM[SM13]*mDet3_234_124 - pM[SM14]*mDet3_234_123;

   // Find the 5x5 det:

   const Double_t det = pM[SM00]*mDet4_1234_1234 - pM[SM01]*mDet4_1234_0234 + pM[SM02]*mDet4_1234_0134
                      - pM[SM03]*mDet4_1234_0124 + pM[SM04]*mDet4_1234_0123;
   if (determ)
      *determ = det;

   if ( det == 0 ) {
      Error("Inv5x5","matrix is singular");
      return kFALSE;
   }

   const Double_t oneOverDet = 1.0/det;
   const Double_t mn1OverDet = - oneOverDet;

   pM[SM00] = mDet4_1234_1234 * oneOverDet;
   pM[SM01] = mDet4_1234_0234 * mn1OverDet;
   pM[SM02] = mDet4_1234_0134 * oneOverDet;
   pM[SM03] = mDet4_1234_0124 * mn1OverDet;
   pM[SM04] = mDet4_1234_0123 * oneOverDet;

   pM[SM11] = mDet4_0234_0234 * oneOverDet;
   pM[SM12] = mDet4_0234_0134 * mn1OverDet;
   pM[SM13] = mDet4_0234_0124 * oneOverDet;
   pM[SM14] = mDet4_0234_0123 * mn1OverDet;

   pM[SM22] = mDet4_0134_0134 * oneOverDet;
   pM[SM23] = mDet4_0134_0124 * mn1OverDet;
   pM[SM24] = mDet4_0134_0123 * oneOverDet;

   pM[SM33] = mDet4_0124_0124 * oneOverDet;
   pM[SM34] = mDet4_0124_0123 * mn1OverDet;

   pM[SM44] = mDet4_0123_0123 * oneOverDet;

   for (Int_t irow = 0; irow < 5; irow++) {
      const Int_t rowOff1 = irow*5; 
      for (Int_t icol = 0; icol < irow; icol++) {
         const Int_t rowOff2 = icol*5; 
         pM[rowOff1+icol] = pM[rowOff2+irow];
      }
   } 

   return kTRUE;
}

// Aij are indices for a 6x6 symmetric matrix.

#define SA00 0
#define SA01 1
#define SA02 2
#define SA03 3
#define SA04 4
#define SA05 5

#define SA10 1
#define SA11 7
#define SA12 8
#define SA13 9
#define SA14 10
#define SA15 11

#define SA20 2
#define SA21 8
#define SA22 14
#define SA23 15
#define SA24 16
#define SA25 17

#define SA30 3
#define SA31 9
#define SA32 15
#define SA33 21
#define SA34 22
#define SA35 23

#define SA40 4
#define SA41 10
#define SA42 16
#define SA43 22
#define SA44 28
#define SA45 29

#define SA50 5
#define SA51 11
#define SA52 17
#define SA53 23
#define SA54 29
#define SA55 35

//______________________________________________________________________________
template<class Element> 
Bool_t TMatrixTSymCramerInv::Inv6x6(TMatrixTSym<Element> &m,Double_t *determ)
{
   if (m.GetNrows() != 6 || m.GetNcols() != 6 || m.GetRowLwb() != m.GetColLwb()) {
      Error("Inv6x6","matrix should be square 6x6");
      return kFALSE;
   }

   Element *pM = m.GetMatrixArray();

   // Find all NECESSSARY 2x2 dets:  (39 of them)

   const Double_t mDet2_34_01 = pM[SA30]*pM[SA41] - pM[SA31]*pM[SA40];
   const Double_t mDet2_34_02 = pM[SA30]*pM[SA42] - pM[SA32]*pM[SA40];
   const Double_t mDet2_34_03 = pM[SA30]*pM[SA43] - pM[SA33]*pM[SA40];
   const Double_t mDet2_34_04 = pM[SA30]*pM[SA44] - pM[SA34]*pM[SA40];
   const Double_t mDet2_34_12 = pM[SA31]*pM[SA42] - pM[SA32]*pM[SA41];
   const Double_t mDet2_34_13 = pM[SA31]*pM[SA43] - pM[SA33]*pM[SA41];
   const Double_t mDet2_34_14 = pM[SA31]*pM[SA44] - pM[SA34]*pM[SA41];
   const Double_t mDet2_34_23 = pM[SA32]*pM[SA43] - pM[SA33]*pM[SA42];
   const Double_t mDet2_34_24 = pM[SA32]*pM[SA44] - pM[SA34]*pM[SA42];
   const Double_t mDet2_34_34 = pM[SA33]*pM[SA44] - pM[SA34]*pM[SA43];
   const Double_t mDet2_35_01 = pM[SA30]*pM[SA51] - pM[SA31]*pM[SA50];
   const Double_t mDet2_35_02 = pM[SA30]*pM[SA52] - pM[SA32]*pM[SA50];
   const Double_t mDet2_35_03 = pM[SA30]*pM[SA53] - pM[SA33]*pM[SA50];
   const Double_t mDet2_35_04 = pM[SA30]*pM[SA54] - pM[SA34]*pM[SA50];
   const Double_t mDet2_35_05 = pM[SA30]*pM[SA55] - pM[SA35]*pM[SA50];
   const Double_t mDet2_35_12 = pM[SA31]*pM[SA52] - pM[SA32]*pM[SA51];
   const Double_t mDet2_35_13 = pM[SA31]*pM[SA53] - pM[SA33]*pM[SA51];
   const Double_t mDet2_35_14 = pM[SA31]*pM[SA54] - pM[SA34]*pM[SA51];
   const Double_t mDet2_35_15 = pM[SA31]*pM[SA55] - pM[SA35]*pM[SA51];
   const Double_t mDet2_35_23 = pM[SA32]*pM[SA53] - pM[SA33]*pM[SA52];
   const Double_t mDet2_35_24 = pM[SA32]*pM[SA54] - pM[SA34]*pM[SA52];
   const Double_t mDet2_35_25 = pM[SA32]*pM[SA55] - pM[SA35]*pM[SA52];
   const Double_t mDet2_35_34 = pM[SA33]*pM[SA54] - pM[SA34]*pM[SA53];
   const Double_t mDet2_35_35 = pM[SA33]*pM[SA55] - pM[SA35]*pM[SA53];
   const Double_t mDet2_45_01 = pM[SA40]*pM[SA51] - pM[SA41]*pM[SA50];
   const Double_t mDet2_45_02 = pM[SA40]*pM[SA52] - pM[SA42]*pM[SA50];
   const Double_t mDet2_45_03 = pM[SA40]*pM[SA53] - pM[SA43]*pM[SA50];
   const Double_t mDet2_45_04 = pM[SA40]*pM[SA54] - pM[SA44]*pM[SA50];
   const Double_t mDet2_45_05 = pM[SA40]*pM[SA55] - pM[SA45]*pM[SA50];
   const Double_t mDet2_45_12 = pM[SA41]*pM[SA52] - pM[SA42]*pM[SA51];
   const Double_t mDet2_45_13 = pM[SA41]*pM[SA53] - pM[SA43]*pM[SA51];
   const Double_t mDet2_45_14 = pM[SA41]*pM[SA54] - pM[SA44]*pM[SA51];
   const Double_t mDet2_45_15 = pM[SA41]*pM[SA55] - pM[SA45]*pM[SA51];
   const Double_t mDet2_45_23 = pM[SA42]*pM[SA53] - pM[SA43]*pM[SA52];
   const Double_t mDet2_45_24 = pM[SA42]*pM[SA54] - pM[SA44]*pM[SA52];
   const Double_t mDet2_45_25 = pM[SA42]*pM[SA55] - pM[SA45]*pM[SA52];
   const Double_t mDet2_45_34 = pM[SA43]*pM[SA54] - pM[SA44]*pM[SA53];
   const Double_t mDet2_45_35 = pM[SA43]*pM[SA55] - pM[SA45]*pM[SA53];
   const Double_t mDet2_45_45 = pM[SA44]*pM[SA55] - pM[SA45]*pM[SA54];

   // Find all NECESSSARY 3x3 dets:  (65 of them)

   const Double_t mDet3_234_012 = pM[SA20]*mDet2_34_12 - pM[SA21]*mDet2_34_02 + pM[SA22]*mDet2_34_01;
   const Double_t mDet3_234_013 = pM[SA20]*mDet2_34_13 - pM[SA21]*mDet2_34_03 + pM[SA23]*mDet2_34_01;
   const Double_t mDet3_234_014 = pM[SA20]*mDet2_34_14 - pM[SA21]*mDet2_34_04 + pM[SA24]*mDet2_34_01;
   const Double_t mDet3_234_023 = pM[SA20]*mDet2_34_23 - pM[SA22]*mDet2_34_03 + pM[SA23]*mDet2_34_02;
   const Double_t mDet3_234_024 = pM[SA20]*mDet2_34_24 - pM[SA22]*mDet2_34_04 + pM[SA24]*mDet2_34_02;
   const Double_t mDet3_234_034 = pM[SA20]*mDet2_34_34 - pM[SA23]*mDet2_34_04 + pM[SA24]*mDet2_34_03;
   const Double_t mDet3_234_123 = pM[SA21]*mDet2_34_23 - pM[SA22]*mDet2_34_13 + pM[SA23]*mDet2_34_12;
   const Double_t mDet3_234_124 = pM[SA21]*mDet2_34_24 - pM[SA22]*mDet2_34_14 + pM[SA24]*mDet2_34_12;
   const Double_t mDet3_234_134 = pM[SA21]*mDet2_34_34 - pM[SA23]*mDet2_34_14 + pM[SA24]*mDet2_34_13;
   const Double_t mDet3_234_234 = pM[SA22]*mDet2_34_34 - pM[SA23]*mDet2_34_24 + pM[SA24]*mDet2_34_23;
   const Double_t mDet3_235_012 = pM[SA20]*mDet2_35_12 - pM[SA21]*mDet2_35_02 + pM[SA22]*mDet2_35_01;
   const Double_t mDet3_235_013 = pM[SA20]*mDet2_35_13 - pM[SA21]*mDet2_35_03 + pM[SA23]*mDet2_35_01;
   const Double_t mDet3_235_014 = pM[SA20]*mDet2_35_14 - pM[SA21]*mDet2_35_04 + pM[SA24]*mDet2_35_01;
   const Double_t mDet3_235_015 = pM[SA20]*mDet2_35_15 - pM[SA21]*mDet2_35_05 + pM[SA25]*mDet2_35_01;
   const Double_t mDet3_235_023 = pM[SA20]*mDet2_35_23 - pM[SA22]*mDet2_35_03 + pM[SA23]*mDet2_35_02;
   const Double_t mDet3_235_024 = pM[SA20]*mDet2_35_24 - pM[SA22]*mDet2_35_04 + pM[SA24]*mDet2_35_02;
   const Double_t mDet3_235_025 = pM[SA20]*mDet2_35_25 - pM[SA22]*mDet2_35_05 + pM[SA25]*mDet2_35_02;
   const Double_t mDet3_235_034 = pM[SA20]*mDet2_35_34 - pM[SA23]*mDet2_35_04 + pM[SA24]*mDet2_35_03;
   const Double_t mDet3_235_035 = pM[SA20]*mDet2_35_35 - pM[SA23]*mDet2_35_05 + pM[SA25]*mDet2_35_03;
   const Double_t mDet3_235_123 = pM[SA21]*mDet2_35_23 - pM[SA22]*mDet2_35_13 + pM[SA23]*mDet2_35_12;
   const Double_t mDet3_235_124 = pM[SA21]*mDet2_35_24 - pM[SA22]*mDet2_35_14 + pM[SA24]*mDet2_35_12;
   const Double_t mDet3_235_125 = pM[SA21]*mDet2_35_25 - pM[SA22]*mDet2_35_15 + pM[SA25]*mDet2_35_12;
   const Double_t mDet3_235_134 = pM[SA21]*mDet2_35_34 - pM[SA23]*mDet2_35_14 + pM[SA24]*mDet2_35_13;
   const Double_t mDet3_235_135 = pM[SA21]*mDet2_35_35 - pM[SA23]*mDet2_35_15 + pM[SA25]*mDet2_35_13;
   const Double_t mDet3_235_234 = pM[SA22]*mDet2_35_34 - pM[SA23]*mDet2_35_24 + pM[SA24]*mDet2_35_23;
   const Double_t mDet3_235_235 = pM[SA22]*mDet2_35_35 - pM[SA23]*mDet2_35_25 + pM[SA25]*mDet2_35_23;
   const Double_t mDet3_245_012 = pM[SA20]*mDet2_45_12 - pM[SA21]*mDet2_45_02 + pM[SA22]*mDet2_45_01;
   const Double_t mDet3_245_013 = pM[SA20]*mDet2_45_13 - pM[SA21]*mDet2_45_03 + pM[SA23]*mDet2_45_01;
   const Double_t mDet3_245_014 = pM[SA20]*mDet2_45_14 - pM[SA21]*mDet2_45_04 + pM[SA24]*mDet2_45_01;
   const Double_t mDet3_245_015 = pM[SA20]*mDet2_45_15 - pM[SA21]*mDet2_45_05 + pM[SA25]*mDet2_45_01;
   const Double_t mDet3_245_023 = pM[SA20]*mDet2_45_23 - pM[SA22]*mDet2_45_03 + pM[SA23]*mDet2_45_02;
   const Double_t mDet3_245_024 = pM[SA20]*mDet2_45_24 - pM[SA22]*mDet2_45_04 + pM[SA24]*mDet2_45_02;
   const Double_t mDet3_245_025 = pM[SA20]*mDet2_45_25 - pM[SA22]*mDet2_45_05 + pM[SA25]*mDet2_45_02;
   const Double_t mDet3_245_034 = pM[SA20]*mDet2_45_34 - pM[SA23]*mDet2_45_04 + pM[SA24]*mDet2_45_03;
   const Double_t mDet3_245_035 = pM[SA20]*mDet2_45_35 - pM[SA23]*mDet2_45_05 + pM[SA25]*mDet2_45_03;
   const Double_t mDet3_245_045 = pM[SA20]*mDet2_45_45 - pM[SA24]*mDet2_45_05 + pM[SA25]*mDet2_45_04;
   const Double_t mDet3_245_123 = pM[SA21]*mDet2_45_23 - pM[SA22]*mDet2_45_13 + pM[SA23]*mDet2_45_12;
   const Double_t mDet3_245_124 = pM[SA21]*mDet2_45_24 - pM[SA22]*mDet2_45_14 + pM[SA24]*mDet2_45_12;
   const Double_t mDet3_245_125 = pM[SA21]*mDet2_45_25 - pM[SA22]*mDet2_45_15 + pM[SA25]*mDet2_45_12;
   const Double_t mDet3_245_134 = pM[SA21]*mDet2_45_34 - pM[SA23]*mDet2_45_14 + pM[SA24]*mDet2_45_13;
   const Double_t mDet3_245_135 = pM[SA21]*mDet2_45_35 - pM[SA23]*mDet2_45_15 + pM[SA25]*mDet2_45_13;
   const Double_t mDet3_245_145 = pM[SA21]*mDet2_45_45 - pM[SA24]*mDet2_45_15 + pM[SA25]*mDet2_45_14;
   const Double_t mDet3_245_234 = pM[SA22]*mDet2_45_34 - pM[SA23]*mDet2_45_24 + pM[SA24]*mDet2_45_23;
   const Double_t mDet3_245_235 = pM[SA22]*mDet2_45_35 - pM[SA23]*mDet2_45_25 + pM[SA25]*mDet2_45_23;
   const Double_t mDet3_245_245 = pM[SA22]*mDet2_45_45 - pM[SA24]*mDet2_45_25 + pM[SA25]*mDet2_45_24;
   const Double_t mDet3_345_012 = pM[SA30]*mDet2_45_12 - pM[SA31]*mDet2_45_02 + pM[SA32]*mDet2_45_01;
   const Double_t mDet3_345_013 = pM[SA30]*mDet2_45_13 - pM[SA31]*mDet2_45_03 + pM[SA33]*mDet2_45_01;
   const Double_t mDet3_345_014 = pM[SA30]*mDet2_45_14 - pM[SA31]*mDet2_45_04 + pM[SA34]*mDet2_45_01;
   const Double_t mDet3_345_015 = pM[SA30]*mDet2_45_15 - pM[SA31]*mDet2_45_05 + pM[SA35]*mDet2_45_01;
   const Double_t mDet3_345_023 = pM[SA30]*mDet2_45_23 - pM[SA32]*mDet2_45_03 + pM[SA33]*mDet2_45_02;
   const Double_t mDet3_345_024 = pM[SA30]*mDet2_45_24 - pM[SA32]*mDet2_45_04 + pM[SA34]*mDet2_45_02;
   const Double_t mDet3_345_025 = pM[SA30]*mDet2_45_25 - pM[SA32]*mDet2_45_05 + pM[SA35]*mDet2_45_02;
   const Double_t mDet3_345_034 = pM[SA30]*mDet2_45_34 - pM[SA33]*mDet2_45_04 + pM[SA34]*mDet2_45_03;
   const Double_t mDet3_345_035 = pM[SA30]*mDet2_45_35 - pM[SA33]*mDet2_45_05 + pM[SA35]*mDet2_45_03;
   const Double_t mDet3_345_045 = pM[SA30]*mDet2_45_45 - pM[SA34]*mDet2_45_05 + pM[SA35]*mDet2_45_04;
   const Double_t mDet3_345_123 = pM[SA31]*mDet2_45_23 - pM[SA32]*mDet2_45_13 + pM[SA33]*mDet2_45_12;
   const Double_t mDet3_345_124 = pM[SA31]*mDet2_45_24 - pM[SA32]*mDet2_45_14 + pM[SA34]*mDet2_45_12;
   const Double_t mDet3_345_125 = pM[SA31]*mDet2_45_25 - pM[SA32]*mDet2_45_15 + pM[SA35]*mDet2_45_12;
   const Double_t mDet3_345_134 = pM[SA31]*mDet2_45_34 - pM[SA33]*mDet2_45_14 + pM[SA34]*mDet2_45_13;
   const Double_t mDet3_345_135 = pM[SA31]*mDet2_45_35 - pM[SA33]*mDet2_45_15 + pM[SA35]*mDet2_45_13;
   const Double_t mDet3_345_145 = pM[SA31]*mDet2_45_45 - pM[SA34]*mDet2_45_15 + pM[SA35]*mDet2_45_14;
   const Double_t mDet3_345_234 = pM[SA32]*mDet2_45_34 - pM[SA33]*mDet2_45_24 + pM[SA34]*mDet2_45_23;
   const Double_t mDet3_345_235 = pM[SA32]*mDet2_45_35 - pM[SA33]*mDet2_45_25 + pM[SA35]*mDet2_45_23;
   const Double_t mDet3_345_245 = pM[SA32]*mDet2_45_45 - pM[SA34]*mDet2_45_25 + pM[SA35]*mDet2_45_24;
   const Double_t mDet3_345_345 = pM[SA33]*mDet2_45_45 - pM[SA34]*mDet2_45_35 + pM[SA35]*mDet2_45_34;

   // Find all NECESSSARY 4x4 dets:  (55 of them)

   const Double_t mDet4_1234_0123 = pM[SA10]*mDet3_234_123 - pM[SA11]*mDet3_234_023
                                  + pM[SA12]*mDet3_234_013 - pM[SA13]*mDet3_234_012;
   const Double_t mDet4_1234_0124 = pM[SA10]*mDet3_234_124 - pM[SA11]*mDet3_234_024
                                  + pM[SA12]*mDet3_234_014 - pM[SA14]*mDet3_234_012;
   const Double_t mDet4_1234_0134 = pM[SA10]*mDet3_234_134 - pM[SA11]*mDet3_234_034
                                  + pM[SA13]*mDet3_234_014 - pM[SA14]*mDet3_234_013;
   const Double_t mDet4_1234_0234 = pM[SA10]*mDet3_234_234 - pM[SA12]*mDet3_234_034
                                  + pM[SA13]*mDet3_234_024 - pM[SA14]*mDet3_234_023;
   const Double_t mDet4_1234_1234 = pM[SA11]*mDet3_234_234 - pM[SA12]*mDet3_234_134
                                  + pM[SA13]*mDet3_234_124 - pM[SA14]*mDet3_234_123;
   const Double_t mDet4_1235_0123 = pM[SA10]*mDet3_235_123 - pM[SA11]*mDet3_235_023
                                  + pM[SA12]*mDet3_235_013 - pM[SA13]*mDet3_235_012;
   const Double_t mDet4_1235_0124 = pM[SA10]*mDet3_235_124 - pM[SA11]*mDet3_235_024
                                  + pM[SA12]*mDet3_235_014 - pM[SA14]*mDet3_235_012;
   const Double_t mDet4_1235_0125 = pM[SA10]*mDet3_235_125 - pM[SA11]*mDet3_235_025
                                  + pM[SA12]*mDet3_235_015 - pM[SA15]*mDet3_235_012;
   const Double_t mDet4_1235_0134 = pM[SA10]*mDet3_235_134 - pM[SA11]*mDet3_235_034
                                  + pM[SA13]*mDet3_235_014 - pM[SA14]*mDet3_235_013;
   const Double_t mDet4_1235_0135 = pM[SA10]*mDet3_235_135 - pM[SA11]*mDet3_235_035
                                  + pM[SA13]*mDet3_235_015 - pM[SA15]*mDet3_235_013;
   const Double_t mDet4_1235_0234 = pM[SA10]*mDet3_235_234 - pM[SA12]*mDet3_235_034
                                  + pM[SA13]*mDet3_235_024 - pM[SA14]*mDet3_235_023;
   const Double_t mDet4_1235_0235 = pM[SA10]*mDet3_235_235 - pM[SA12]*mDet3_235_035
                                  + pM[SA13]*mDet3_235_025 - pM[SA15]*mDet3_235_023;
   const Double_t mDet4_1235_1234 = pM[SA11]*mDet3_235_234 - pM[SA12]*mDet3_235_134
                                  + pM[SA13]*mDet3_235_124 - pM[SA14]*mDet3_235_123;
   const Double_t mDet4_1235_1235 = pM[SA11]*mDet3_235_235 - pM[SA12]*mDet3_235_135
                                  + pM[SA13]*mDet3_235_125 - pM[SA15]*mDet3_235_123;
   const Double_t mDet4_1245_0123 = pM[SA10]*mDet3_245_123 - pM[SA11]*mDet3_245_023
                                  + pM[SA12]*mDet3_245_013 - pM[SA13]*mDet3_245_012;
   const Double_t mDet4_1245_0124 = pM[SA10]*mDet3_245_124 - pM[SA11]*mDet3_245_024
                                  + pM[SA12]*mDet3_245_014 - pM[SA14]*mDet3_245_012;
   const Double_t mDet4_1245_0125 = pM[SA10]*mDet3_245_125 - pM[SA11]*mDet3_245_025
                                  + pM[SA12]*mDet3_245_015 - pM[SA15]*mDet3_245_012;
   const Double_t mDet4_1245_0134 = pM[SA10]*mDet3_245_134 - pM[SA11]*mDet3_245_034
                                  + pM[SA13]*mDet3_245_014 - pM[SA14]*mDet3_245_013;
   const Double_t mDet4_1245_0135 = pM[SA10]*mDet3_245_135 - pM[SA11]*mDet3_245_035
                                  + pM[SA13]*mDet3_245_015 - pM[SA15]*mDet3_245_013;
   const Double_t mDet4_1245_0145 = pM[SA10]*mDet3_245_145 - pM[SA11]*mDet3_245_045
                                  + pM[SA14]*mDet3_245_015 - pM[SA15]*mDet3_245_014;
   const Double_t mDet4_1245_0234 = pM[SA10]*mDet3_245_234 - pM[SA12]*mDet3_245_034
                                  + pM[SA13]*mDet3_245_024 - pM[SA14]*mDet3_245_023;
   const Double_t mDet4_1245_0235 = pM[SA10]*mDet3_245_235 - pM[SA12]*mDet3_245_035
                                  + pM[SA13]*mDet3_245_025 - pM[SA15]*mDet3_245_023;
   const Double_t mDet4_1245_0245 = pM[SA10]*mDet3_245_245 - pM[SA12]*mDet3_245_045
                                  + pM[SA14]*mDet3_245_025 - pM[SA15]*mDet3_245_024;
   const Double_t mDet4_1245_1234 = pM[SA11]*mDet3_245_234 - pM[SA12]*mDet3_245_134
                                  + pM[SA13]*mDet3_245_124 - pM[SA14]*mDet3_245_123;
   const Double_t mDet4_1245_1235 = pM[SA11]*mDet3_245_235 - pM[SA12]*mDet3_245_135
                                  + pM[SA13]*mDet3_245_125 - pM[SA15]*mDet3_245_123;
   const Double_t mDet4_1245_1245 = pM[SA11]*mDet3_245_245 - pM[SA12]*mDet3_245_145
                                  + pM[SA14]*mDet3_245_125 - pM[SA15]*mDet3_245_124;
   const Double_t mDet4_1345_0123 = pM[SA10]*mDet3_345_123 - pM[SA11]*mDet3_345_023
                                  + pM[SA12]*mDet3_345_013 - pM[SA13]*mDet3_345_012;
   const Double_t mDet4_1345_0124 = pM[SA10]*mDet3_345_124 - pM[SA11]*mDet3_345_024
                                  + pM[SA12]*mDet3_345_014 - pM[SA14]*mDet3_345_012;
   const Double_t mDet4_1345_0125 = pM[SA10]*mDet3_345_125 - pM[SA11]*mDet3_345_025
                                  + pM[SA12]*mDet3_345_015 - pM[SA15]*mDet3_345_012;
   const Double_t mDet4_1345_0134 = pM[SA10]*mDet3_345_134 - pM[SA11]*mDet3_345_034
                                  + pM[SA13]*mDet3_345_014 - pM[SA14]*mDet3_345_013;
   const Double_t mDet4_1345_0135 = pM[SA10]*mDet3_345_135 - pM[SA11]*mDet3_345_035
                                  + pM[SA13]*mDet3_345_015 - pM[SA15]*mDet3_345_013;
   const Double_t mDet4_1345_0145 = pM[SA10]*mDet3_345_145 - pM[SA11]*mDet3_345_045
                                  + pM[SA14]*mDet3_345_015 - pM[SA15]*mDet3_345_014;
   const Double_t mDet4_1345_0234 = pM[SA10]*mDet3_345_234 - pM[SA12]*mDet3_345_034
                                  + pM[SA13]*mDet3_345_024 - pM[SA14]*mDet3_345_023;
   const Double_t mDet4_1345_0235 = pM[SA10]*mDet3_345_235 - pM[SA12]*mDet3_345_035
                                  + pM[SA13]*mDet3_345_025 - pM[SA15]*mDet3_345_023;
   const Double_t mDet4_1345_0245 = pM[SA10]*mDet3_345_245 - pM[SA12]*mDet3_345_045
                                  + pM[SA14]*mDet3_345_025 - pM[SA15]*mDet3_345_024;
   const Double_t mDet4_1345_0345 = pM[SA10]*mDet3_345_345 - pM[SA13]*mDet3_345_045
                                  + pM[SA14]*mDet3_345_035 - pM[SA15]*mDet3_345_034;
   const Double_t mDet4_1345_1234 = pM[SA11]*mDet3_345_234 - pM[SA12]*mDet3_345_134
                                  + pM[SA13]*mDet3_345_124 - pM[SA14]*mDet3_345_123;
   const Double_t mDet4_1345_1235 = pM[SA11]*mDet3_345_235 - pM[SA12]*mDet3_345_135
                                  + pM[SA13]*mDet3_345_125 - pM[SA15]*mDet3_345_123;
   const Double_t mDet4_1345_1245 = pM[SA11]*mDet3_345_245 - pM[SA12]*mDet3_345_145
                                  + pM[SA14]*mDet3_345_125 - pM[SA15]*mDet3_345_124;
   const Double_t mDet4_1345_1345 = pM[SA11]*mDet3_345_345 - pM[SA13]*mDet3_345_145
                                  + pM[SA14]*mDet3_345_135 - pM[SA15]*mDet3_345_134;
   const Double_t mDet4_2345_0123 = pM[SA20]*mDet3_345_123 - pM[SA21]*mDet3_345_023
                                  + pM[SA22]*mDet3_345_013 - pM[SA23]*mDet3_345_012;
   const Double_t mDet4_2345_0124 = pM[SA20]*mDet3_345_124 - pM[SA21]*mDet3_345_024
                                  + pM[SA22]*mDet3_345_014 - pM[SA24]*mDet3_345_012;
   const Double_t mDet4_2345_0125 = pM[SA20]*mDet3_345_125 - pM[SA21]*mDet3_345_025
                                  + pM[SA22]*mDet3_345_015 - pM[SA25]*mDet3_345_012;
   const Double_t mDet4_2345_0134 = pM[SA20]*mDet3_345_134 - pM[SA21]*mDet3_345_034
                                  + pM[SA23]*mDet3_345_014 - pM[SA24]*mDet3_345_013;
   const Double_t mDet4_2345_0135 = pM[SA20]*mDet3_345_135 - pM[SA21]*mDet3_345_035
                                  + pM[SA23]*mDet3_345_015 - pM[SA25]*mDet3_345_013;
   const Double_t mDet4_2345_0145 = pM[SA20]*mDet3_345_145 - pM[SA21]*mDet3_345_045
                                  + pM[SA24]*mDet3_345_015 - pM[SA25]*mDet3_345_014;
   const Double_t mDet4_2345_0234 = pM[SA20]*mDet3_345_234 - pM[SA22]*mDet3_345_034
                                  + pM[SA23]*mDet3_345_024 - pM[SA24]*mDet3_345_023;
   const Double_t mDet4_2345_0235 = pM[SA20]*mDet3_345_235 - pM[SA22]*mDet3_345_035
                                  + pM[SA23]*mDet3_345_025 - pM[SA25]*mDet3_345_023;
   const Double_t mDet4_2345_0245 = pM[SA20]*mDet3_345_245 - pM[SA22]*mDet3_345_045
                                  + pM[SA24]*mDet3_345_025 - pM[SA25]*mDet3_345_024;
   const Double_t mDet4_2345_0345 = pM[SA20]*mDet3_345_345 - pM[SA23]*mDet3_345_045
                                  + pM[SA24]*mDet3_345_035 - pM[SA25]*mDet3_345_034;
   const Double_t mDet4_2345_1234 = pM[SA21]*mDet3_345_234 - pM[SA22]*mDet3_345_134
                                  + pM[SA23]*mDet3_345_124 - pM[SA24]*mDet3_345_123;
   const Double_t mDet4_2345_1235 = pM[SA21]*mDet3_345_235 - pM[SA22]*mDet3_345_135
                                  + pM[SA23]*mDet3_345_125 - pM[SA25]*mDet3_345_123;
   const Double_t mDet4_2345_1245 = pM[SA21]*mDet3_345_245 - pM[SA22]*mDet3_345_145
                                  + pM[SA24]*mDet3_345_125 - pM[SA25]*mDet3_345_124;
   const Double_t mDet4_2345_1345 = pM[SA21]*mDet3_345_345 - pM[SA23]*mDet3_345_145
                                  + pM[SA24]*mDet3_345_135 - pM[SA25]*mDet3_345_134;
   const Double_t mDet4_2345_2345 = pM[SA22]*mDet3_345_345 - pM[SA23]*mDet3_345_245
                                  + pM[SA24]*mDet3_345_235 - pM[SA25]*mDet3_345_234;

   // Find all NECESSSARY 5x5 dets:  (19 of them)

   const Double_t mDet5_01234_01234 = pM[SA00]*mDet4_1234_1234 - pM[SA01]*mDet4_1234_0234
                                    + pM[SA02]*mDet4_1234_0134 - pM[SA03]*mDet4_1234_0124 + pM[SA04]*mDet4_1234_0123;
   const Double_t mDet5_01235_01234 = pM[SA00]*mDet4_1235_1234 - pM[SA01]*mDet4_1235_0234
                                    + pM[SA02]*mDet4_1235_0134 - pM[SA03]*mDet4_1235_0124 + pM[SA04]*mDet4_1235_0123;
   const Double_t mDet5_01235_01235 = pM[SA00]*mDet4_1235_1235 - pM[SA01]*mDet4_1235_0235
                                    + pM[SA02]*mDet4_1235_0135 - pM[SA03]*mDet4_1235_0125 + pM[SA05]*mDet4_1235_0123;
   const Double_t mDet5_01245_01234 = pM[SA00]*mDet4_1245_1234 - pM[SA01]*mDet4_1245_0234
                                    + pM[SA02]*mDet4_1245_0134 - pM[SA03]*mDet4_1245_0124 + pM[SA04]*mDet4_1245_0123;
   const Double_t mDet5_01245_01235 = pM[SA00]*mDet4_1245_1235 - pM[SA01]*mDet4_1245_0235
                                    + pM[SA02]*mDet4_1245_0135 - pM[SA03]*mDet4_1245_0125 + pM[SA05]*mDet4_1245_0123;
   const Double_t mDet5_01245_01245 = pM[SA00]*mDet4_1245_1245 - pM[SA01]*mDet4_1245_0245
                                    + pM[SA02]*mDet4_1245_0145 - pM[SA04]*mDet4_1245_0125 + pM[SA05]*mDet4_1245_0124;
   const Double_t mDet5_01345_01234 = pM[SA00]*mDet4_1345_1234 - pM[SA01]*mDet4_1345_0234
                                    + pM[SA02]*mDet4_1345_0134 - pM[SA03]*mDet4_1345_0124 + pM[SA04]*mDet4_1345_0123;
   const Double_t mDet5_01345_01235 = pM[SA00]*mDet4_1345_1235 - pM[SA01]*mDet4_1345_0235
                                    + pM[SA02]*mDet4_1345_0135 - pM[SA03]*mDet4_1345_0125 + pM[SA05]*mDet4_1345_0123;
   const Double_t mDet5_01345_01245 = pM[SA00]*mDet4_1345_1245 - pM[SA01]*mDet4_1345_0245
                                    + pM[SA02]*mDet4_1345_0145 - pM[SA04]*mDet4_1345_0125 + pM[SA05]*mDet4_1345_0124;
   const Double_t mDet5_01345_01345 = pM[SA00]*mDet4_1345_1345 - pM[SA01]*mDet4_1345_0345
                                    + pM[SA03]*mDet4_1345_0145 - pM[SA04]*mDet4_1345_0135 + pM[SA05]*mDet4_1345_0134;
   const Double_t mDet5_02345_01234 = pM[SA00]*mDet4_2345_1234 - pM[SA01]*mDet4_2345_0234
                                    + pM[SA02]*mDet4_2345_0134 - pM[SA03]*mDet4_2345_0124 + pM[SA04]*mDet4_2345_0123;
   const Double_t mDet5_02345_01235 = pM[SA00]*mDet4_2345_1235 - pM[SA01]*mDet4_2345_0235
                                    + pM[SA02]*mDet4_2345_0135 - pM[SA03]*mDet4_2345_0125 + pM[SA05]*mDet4_2345_0123;
   const Double_t mDet5_02345_01245 = pM[SA00]*mDet4_2345_1245 - pM[SA01]*mDet4_2345_0245
                                    + pM[SA02]*mDet4_2345_0145 - pM[SA04]*mDet4_2345_0125 + pM[SA05]*mDet4_2345_0124;
   const Double_t mDet5_02345_01345 = pM[SA00]*mDet4_2345_1345 - pM[SA01]*mDet4_2345_0345
                                    + pM[SA03]*mDet4_2345_0145 - pM[SA04]*mDet4_2345_0135 + pM[SA05]*mDet4_2345_0134;
   const Double_t mDet5_02345_02345 = pM[SA00]*mDet4_2345_2345 - pM[SA02]*mDet4_2345_0345
                                    + pM[SA03]*mDet4_2345_0245 - pM[SA04]*mDet4_2345_0235 + pM[SA05]*mDet4_2345_0234;
   const Double_t mDet5_12345_01234 = pM[SA10]*mDet4_2345_1234 - pM[SA11]*mDet4_2345_0234
                                    + pM[SA12]*mDet4_2345_0134 - pM[SA13]*mDet4_2345_0124 + pM[SA14]*mDet4_2345_0123;
   const Double_t mDet5_12345_01235 = pM[SA10]*mDet4_2345_1235 - pM[SA11]*mDet4_2345_0235
                                    + pM[SA12]*mDet4_2345_0135 - pM[SA13]*mDet4_2345_0125 + pM[SA15]*mDet4_2345_0123;
   const Double_t mDet5_12345_01245 = pM[SA10]*mDet4_2345_1245 - pM[SA11]*mDet4_2345_0245
                                    + pM[SA12]*mDet4_2345_0145 - pM[SA14]*mDet4_2345_0125 + pM[SA15]*mDet4_2345_0124;
   const Double_t mDet5_12345_01345 = pM[SA10]*mDet4_2345_1345 - pM[SA11]*mDet4_2345_0345
                                    + pM[SA13]*mDet4_2345_0145 - pM[SA14]*mDet4_2345_0135 + pM[SA15]*mDet4_2345_0134;
   const Double_t mDet5_12345_02345 = pM[SA10]*mDet4_2345_2345 - pM[SA12]*mDet4_2345_0345
                                    + pM[SA13]*mDet4_2345_0245 - pM[SA14]*mDet4_2345_0235 + pM[SA15]*mDet4_2345_0234;
   const Double_t mDet5_12345_12345 = pM[SA11]*mDet4_2345_2345 - pM[SA12]*mDet4_2345_1345
                                    + pM[SA13]*mDet4_2345_1245 - pM[SA14]*mDet4_2345_1235 + pM[SA15]*mDet4_2345_1234;

   // Find the determinant 

   const Double_t det = pM[SA00]*mDet5_12345_12345 - pM[SA01]*mDet5_12345_02345 + pM[SA02]*mDet5_12345_01345
                      - pM[SA03]*mDet5_12345_01245 + pM[SA04]*mDet5_12345_01235 - pM[SA05]*mDet5_12345_01234;

   if (determ)
      *determ = det;

   if ( det == 0 ) {
      Error("Inv6x6","matrix is singular");
      return kFALSE;
   }

   const Double_t oneOverDet = 1.0/det;
   const Double_t mn1OverDet = - oneOverDet;

   pM[SA00] =  mDet5_12345_12345*oneOverDet;
   pM[SA01] =  mDet5_12345_02345*mn1OverDet;
   pM[SA02] =  mDet5_12345_01345*oneOverDet;
   pM[SA03] =  mDet5_12345_01245*mn1OverDet;
   pM[SA04] =  mDet5_12345_01235*oneOverDet;
   pM[SA05] =  mDet5_12345_01234*mn1OverDet;

   pM[SA11] =  mDet5_02345_02345*oneOverDet;
   pM[SA12] =  mDet5_02345_01345*mn1OverDet;
   pM[SA13] =  mDet5_02345_01245*oneOverDet;
   pM[SA14] =  mDet5_02345_01235*mn1OverDet;
   pM[SA15] =  mDet5_02345_01234*oneOverDet;

   pM[SA22] =  mDet5_01345_01345*oneOverDet;
   pM[SA23] =  mDet5_01345_01245*mn1OverDet;
   pM[SA24] =  mDet5_01345_01235*oneOverDet;
   pM[SA25] =  mDet5_01345_01234*mn1OverDet;

   pM[SA33] =  mDet5_01245_01245*oneOverDet;
   pM[SA34] =  mDet5_01245_01235*mn1OverDet;
   pM[SA35] =  mDet5_01245_01234*oneOverDet;

   pM[SA44] =  mDet5_01235_01235*oneOverDet;
   pM[SA45] =  mDet5_01235_01234*mn1OverDet;

   pM[SA55] =  mDet5_01234_01234*oneOverDet;

   for (Int_t irow = 0; irow < 6; irow++) {
      const Int_t rowOff1 = irow*6; 
      for (Int_t icol = 0; icol < irow; icol++) {
         const Int_t rowOff2 = icol*6; 
         pM[rowOff1+icol] = pM[rowOff2+irow];
      }
   } 

   return kTRUE;
}

#ifndef ROOT_TMatrixFSymfwd
#include "TMatrixFSymfwd.h"
#endif

template Bool_t TMatrixTSymCramerInv::Inv2x2<Float_t>(TMatrixFSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv3x3<Float_t>(TMatrixFSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv4x4<Float_t>(TMatrixFSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv5x5<Float_t>(TMatrixFSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv6x6<Float_t>(TMatrixFSym&,Double_t*);

#ifndef ROOT_TMatrixDSymfwd
#include "TMatrixDSymfwd.h"
#endif

template Bool_t TMatrixTSymCramerInv::Inv2x2<Double_t>(TMatrixDSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv3x3<Double_t>(TMatrixDSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv4x4<Double_t>(TMatrixDSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv5x5<Double_t>(TMatrixDSym&,Double_t*);
template Bool_t TMatrixTSymCramerInv::Inv6x6<Double_t>(TMatrixDSym&,Double_t*);
 TMatrixTSymCramerInv.cxx:1
 TMatrixTSymCramerInv.cxx:2
 TMatrixTSymCramerInv.cxx:3
 TMatrixTSymCramerInv.cxx:4
 TMatrixTSymCramerInv.cxx:5
 TMatrixTSymCramerInv.cxx:6
 TMatrixTSymCramerInv.cxx:7
 TMatrixTSymCramerInv.cxx:8
 TMatrixTSymCramerInv.cxx:9
 TMatrixTSymCramerInv.cxx:10
 TMatrixTSymCramerInv.cxx:11
 TMatrixTSymCramerInv.cxx:12
 TMatrixTSymCramerInv.cxx:13
 TMatrixTSymCramerInv.cxx:14
 TMatrixTSymCramerInv.cxx:15
 TMatrixTSymCramerInv.cxx:16
 TMatrixTSymCramerInv.cxx:17
 TMatrixTSymCramerInv.cxx:18
 TMatrixTSymCramerInv.cxx:19
 TMatrixTSymCramerInv.cxx:20
 TMatrixTSymCramerInv.cxx:21
 TMatrixTSymCramerInv.cxx:22
 TMatrixTSymCramerInv.cxx:23
 TMatrixTSymCramerInv.cxx:24
 TMatrixTSymCramerInv.cxx:25
 TMatrixTSymCramerInv.cxx:26
 TMatrixTSymCramerInv.cxx:27
 TMatrixTSymCramerInv.cxx:28
 TMatrixTSymCramerInv.cxx:29
 TMatrixTSymCramerInv.cxx:30
 TMatrixTSymCramerInv.cxx:31
 TMatrixTSymCramerInv.cxx:32
 TMatrixTSymCramerInv.cxx:33
 TMatrixTSymCramerInv.cxx:34
 TMatrixTSymCramerInv.cxx:35
 TMatrixTSymCramerInv.cxx:36
 TMatrixTSymCramerInv.cxx:37
 TMatrixTSymCramerInv.cxx:38
 TMatrixTSymCramerInv.cxx:39
 TMatrixTSymCramerInv.cxx:40
 TMatrixTSymCramerInv.cxx:41
 TMatrixTSymCramerInv.cxx:42
 TMatrixTSymCramerInv.cxx:43
 TMatrixTSymCramerInv.cxx:44
 TMatrixTSymCramerInv.cxx:45
 TMatrixTSymCramerInv.cxx:46
 TMatrixTSymCramerInv.cxx:47
 TMatrixTSymCramerInv.cxx:48
 TMatrixTSymCramerInv.cxx:49
 TMatrixTSymCramerInv.cxx:50
 TMatrixTSymCramerInv.cxx:51
 TMatrixTSymCramerInv.cxx:52
 TMatrixTSymCramerInv.cxx:53
 TMatrixTSymCramerInv.cxx:54
 TMatrixTSymCramerInv.cxx:55
 TMatrixTSymCramerInv.cxx:56
 TMatrixTSymCramerInv.cxx:57
 TMatrixTSymCramerInv.cxx:58
 TMatrixTSymCramerInv.cxx:59
 TMatrixTSymCramerInv.cxx:60
 TMatrixTSymCramerInv.cxx:61
 TMatrixTSymCramerInv.cxx:62
 TMatrixTSymCramerInv.cxx:63
 TMatrixTSymCramerInv.cxx:64
 TMatrixTSymCramerInv.cxx:65
 TMatrixTSymCramerInv.cxx:66
 TMatrixTSymCramerInv.cxx:67
 TMatrixTSymCramerInv.cxx:68
 TMatrixTSymCramerInv.cxx:69
 TMatrixTSymCramerInv.cxx:70
 TMatrixTSymCramerInv.cxx:71
 TMatrixTSymCramerInv.cxx:72
 TMatrixTSymCramerInv.cxx:73
 TMatrixTSymCramerInv.cxx:74
 TMatrixTSymCramerInv.cxx:75
 TMatrixTSymCramerInv.cxx:76
 TMatrixTSymCramerInv.cxx:77
 TMatrixTSymCramerInv.cxx:78
 TMatrixTSymCramerInv.cxx:79
 TMatrixTSymCramerInv.cxx:80
 TMatrixTSymCramerInv.cxx:81
 TMatrixTSymCramerInv.cxx:82
 TMatrixTSymCramerInv.cxx:83
 TMatrixTSymCramerInv.cxx:84
 TMatrixTSymCramerInv.cxx:85
 TMatrixTSymCramerInv.cxx:86
 TMatrixTSymCramerInv.cxx:87
 TMatrixTSymCramerInv.cxx:88
 TMatrixTSymCramerInv.cxx:89
 TMatrixTSymCramerInv.cxx:90
 TMatrixTSymCramerInv.cxx:91
 TMatrixTSymCramerInv.cxx:92
 TMatrixTSymCramerInv.cxx:93
 TMatrixTSymCramerInv.cxx:94
 TMatrixTSymCramerInv.cxx:95
 TMatrixTSymCramerInv.cxx:96
 TMatrixTSymCramerInv.cxx:97
 TMatrixTSymCramerInv.cxx:98
 TMatrixTSymCramerInv.cxx:99
 TMatrixTSymCramerInv.cxx:100
 TMatrixTSymCramerInv.cxx:101
 TMatrixTSymCramerInv.cxx:102
 TMatrixTSymCramerInv.cxx:103
 TMatrixTSymCramerInv.cxx:104
 TMatrixTSymCramerInv.cxx:105
 TMatrixTSymCramerInv.cxx:106
 TMatrixTSymCramerInv.cxx:107
 TMatrixTSymCramerInv.cxx:108
 TMatrixTSymCramerInv.cxx:109
 TMatrixTSymCramerInv.cxx:110
 TMatrixTSymCramerInv.cxx:111
 TMatrixTSymCramerInv.cxx:112
 TMatrixTSymCramerInv.cxx:113
 TMatrixTSymCramerInv.cxx:114
 TMatrixTSymCramerInv.cxx:115
 TMatrixTSymCramerInv.cxx:116
 TMatrixTSymCramerInv.cxx:117
 TMatrixTSymCramerInv.cxx:118
 TMatrixTSymCramerInv.cxx:119
 TMatrixTSymCramerInv.cxx:120
 TMatrixTSymCramerInv.cxx:121
 TMatrixTSymCramerInv.cxx:122
 TMatrixTSymCramerInv.cxx:123
 TMatrixTSymCramerInv.cxx:124
 TMatrixTSymCramerInv.cxx:125
 TMatrixTSymCramerInv.cxx:126
 TMatrixTSymCramerInv.cxx:127
 TMatrixTSymCramerInv.cxx:128
 TMatrixTSymCramerInv.cxx:129
 TMatrixTSymCramerInv.cxx:130
 TMatrixTSymCramerInv.cxx:131
 TMatrixTSymCramerInv.cxx:132
 TMatrixTSymCramerInv.cxx:133
 TMatrixTSymCramerInv.cxx:134
 TMatrixTSymCramerInv.cxx:135
 TMatrixTSymCramerInv.cxx:136
 TMatrixTSymCramerInv.cxx:137
 TMatrixTSymCramerInv.cxx:138
 TMatrixTSymCramerInv.cxx:139
 TMatrixTSymCramerInv.cxx:140
 TMatrixTSymCramerInv.cxx:141
 TMatrixTSymCramerInv.cxx:142
 TMatrixTSymCramerInv.cxx:143
 TMatrixTSymCramerInv.cxx:144
 TMatrixTSymCramerInv.cxx:145
 TMatrixTSymCramerInv.cxx:146
 TMatrixTSymCramerInv.cxx:147
 TMatrixTSymCramerInv.cxx:148
 TMatrixTSymCramerInv.cxx:149
 TMatrixTSymCramerInv.cxx:150
 TMatrixTSymCramerInv.cxx:151
 TMatrixTSymCramerInv.cxx:152
 TMatrixTSymCramerInv.cxx:153
 TMatrixTSymCramerInv.cxx:154
 TMatrixTSymCramerInv.cxx:155
 TMatrixTSymCramerInv.cxx:156
 TMatrixTSymCramerInv.cxx:157
 TMatrixTSymCramerInv.cxx:158
 TMatrixTSymCramerInv.cxx:159
 TMatrixTSymCramerInv.cxx:160
 TMatrixTSymCramerInv.cxx:161
 TMatrixTSymCramerInv.cxx:162
 TMatrixTSymCramerInv.cxx:163
 TMatrixTSymCramerInv.cxx:164
 TMatrixTSymCramerInv.cxx:165
 TMatrixTSymCramerInv.cxx:166
 TMatrixTSymCramerInv.cxx:167
 TMatrixTSymCramerInv.cxx:168
 TMatrixTSymCramerInv.cxx:169
 TMatrixTSymCramerInv.cxx:170
 TMatrixTSymCramerInv.cxx:171
 TMatrixTSymCramerInv.cxx:172
 TMatrixTSymCramerInv.cxx:173
 TMatrixTSymCramerInv.cxx:174
 TMatrixTSymCramerInv.cxx:175
 TMatrixTSymCramerInv.cxx:176
 TMatrixTSymCramerInv.cxx:177
 TMatrixTSymCramerInv.cxx:178
 TMatrixTSymCramerInv.cxx:179
 TMatrixTSymCramerInv.cxx:180
 TMatrixTSymCramerInv.cxx:181
 TMatrixTSymCramerInv.cxx:182
 TMatrixTSymCramerInv.cxx:183
 TMatrixTSymCramerInv.cxx:184
 TMatrixTSymCramerInv.cxx:185
 TMatrixTSymCramerInv.cxx:186
 TMatrixTSymCramerInv.cxx:187
 TMatrixTSymCramerInv.cxx:188
 TMatrixTSymCramerInv.cxx:189
 TMatrixTSymCramerInv.cxx:190
 TMatrixTSymCramerInv.cxx:191
 TMatrixTSymCramerInv.cxx:192
 TMatrixTSymCramerInv.cxx:193
 TMatrixTSymCramerInv.cxx:194
 TMatrixTSymCramerInv.cxx:195
 TMatrixTSymCramerInv.cxx:196
 TMatrixTSymCramerInv.cxx:197
 TMatrixTSymCramerInv.cxx:198
 TMatrixTSymCramerInv.cxx:199
 TMatrixTSymCramerInv.cxx:200
 TMatrixTSymCramerInv.cxx:201
 TMatrixTSymCramerInv.cxx:202
 TMatrixTSymCramerInv.cxx:203
 TMatrixTSymCramerInv.cxx:204
 TMatrixTSymCramerInv.cxx:205
 TMatrixTSymCramerInv.cxx:206
 TMatrixTSymCramerInv.cxx:207
 TMatrixTSymCramerInv.cxx:208
 TMatrixTSymCramerInv.cxx:209
 TMatrixTSymCramerInv.cxx:210
 TMatrixTSymCramerInv.cxx:211
 TMatrixTSymCramerInv.cxx:212
 TMatrixTSymCramerInv.cxx:213
 TMatrixTSymCramerInv.cxx:214
 TMatrixTSymCramerInv.cxx:215
 TMatrixTSymCramerInv.cxx:216
 TMatrixTSymCramerInv.cxx:217
 TMatrixTSymCramerInv.cxx:218
 TMatrixTSymCramerInv.cxx:219
 TMatrixTSymCramerInv.cxx:220
 TMatrixTSymCramerInv.cxx:221
 TMatrixTSymCramerInv.cxx:222
 TMatrixTSymCramerInv.cxx:223
 TMatrixTSymCramerInv.cxx:224
 TMatrixTSymCramerInv.cxx:225
 TMatrixTSymCramerInv.cxx:226
 TMatrixTSymCramerInv.cxx:227
 TMatrixTSymCramerInv.cxx:228
 TMatrixTSymCramerInv.cxx:229
 TMatrixTSymCramerInv.cxx:230
 TMatrixTSymCramerInv.cxx:231
 TMatrixTSymCramerInv.cxx:232
 TMatrixTSymCramerInv.cxx:233
 TMatrixTSymCramerInv.cxx:234
 TMatrixTSymCramerInv.cxx:235
 TMatrixTSymCramerInv.cxx:236
 TMatrixTSymCramerInv.cxx:237
 TMatrixTSymCramerInv.cxx:238
 TMatrixTSymCramerInv.cxx:239
 TMatrixTSymCramerInv.cxx:240
 TMatrixTSymCramerInv.cxx:241
 TMatrixTSymCramerInv.cxx:242
 TMatrixTSymCramerInv.cxx:243
 TMatrixTSymCramerInv.cxx:244
 TMatrixTSymCramerInv.cxx:245
 TMatrixTSymCramerInv.cxx:246
 TMatrixTSymCramerInv.cxx:247
 TMatrixTSymCramerInv.cxx:248
 TMatrixTSymCramerInv.cxx:249
 TMatrixTSymCramerInv.cxx:250
 TMatrixTSymCramerInv.cxx:251
 TMatrixTSymCramerInv.cxx:252
 TMatrixTSymCramerInv.cxx:253
 TMatrixTSymCramerInv.cxx:254
 TMatrixTSymCramerInv.cxx:255
 TMatrixTSymCramerInv.cxx:256
 TMatrixTSymCramerInv.cxx:257
 TMatrixTSymCramerInv.cxx:258
 TMatrixTSymCramerInv.cxx:259
 TMatrixTSymCramerInv.cxx:260
 TMatrixTSymCramerInv.cxx:261
 TMatrixTSymCramerInv.cxx:262
 TMatrixTSymCramerInv.cxx:263
 TMatrixTSymCramerInv.cxx:264
 TMatrixTSymCramerInv.cxx:265
 TMatrixTSymCramerInv.cxx:266
 TMatrixTSymCramerInv.cxx:267
 TMatrixTSymCramerInv.cxx:268
 TMatrixTSymCramerInv.cxx:269
 TMatrixTSymCramerInv.cxx:270
 TMatrixTSymCramerInv.cxx:271
 TMatrixTSymCramerInv.cxx:272
 TMatrixTSymCramerInv.cxx:273
 TMatrixTSymCramerInv.cxx:274
 TMatrixTSymCramerInv.cxx:275
 TMatrixTSymCramerInv.cxx:276
 TMatrixTSymCramerInv.cxx:277
 TMatrixTSymCramerInv.cxx:278
 TMatrixTSymCramerInv.cxx:279
 TMatrixTSymCramerInv.cxx:280
 TMatrixTSymCramerInv.cxx:281
 TMatrixTSymCramerInv.cxx:282
 TMatrixTSymCramerInv.cxx:283
 TMatrixTSymCramerInv.cxx:284
 TMatrixTSymCramerInv.cxx:285
 TMatrixTSymCramerInv.cxx:286
 TMatrixTSymCramerInv.cxx:287
 TMatrixTSymCramerInv.cxx:288
 TMatrixTSymCramerInv.cxx:289
 TMatrixTSymCramerInv.cxx:290
 TMatrixTSymCramerInv.cxx:291
 TMatrixTSymCramerInv.cxx:292
 TMatrixTSymCramerInv.cxx:293
 TMatrixTSymCramerInv.cxx:294
 TMatrixTSymCramerInv.cxx:295
 TMatrixTSymCramerInv.cxx:296
 TMatrixTSymCramerInv.cxx:297
 TMatrixTSymCramerInv.cxx:298
 TMatrixTSymCramerInv.cxx:299
 TMatrixTSymCramerInv.cxx:300
 TMatrixTSymCramerInv.cxx:301
 TMatrixTSymCramerInv.cxx:302
 TMatrixTSymCramerInv.cxx:303
 TMatrixTSymCramerInv.cxx:304
 TMatrixTSymCramerInv.cxx:305
 TMatrixTSymCramerInv.cxx:306
 TMatrixTSymCramerInv.cxx:307
 TMatrixTSymCramerInv.cxx:308
 TMatrixTSymCramerInv.cxx:309
 TMatrixTSymCramerInv.cxx:310
 TMatrixTSymCramerInv.cxx:311
 TMatrixTSymCramerInv.cxx:312
 TMatrixTSymCramerInv.cxx:313
 TMatrixTSymCramerInv.cxx:314
 TMatrixTSymCramerInv.cxx:315
 TMatrixTSymCramerInv.cxx:316
 TMatrixTSymCramerInv.cxx:317
 TMatrixTSymCramerInv.cxx:318
 TMatrixTSymCramerInv.cxx:319
 TMatrixTSymCramerInv.cxx:320
 TMatrixTSymCramerInv.cxx:321
 TMatrixTSymCramerInv.cxx:322
 TMatrixTSymCramerInv.cxx:323
 TMatrixTSymCramerInv.cxx:324
 TMatrixTSymCramerInv.cxx:325
 TMatrixTSymCramerInv.cxx:326
 TMatrixTSymCramerInv.cxx:327
 TMatrixTSymCramerInv.cxx:328
 TMatrixTSymCramerInv.cxx:329
 TMatrixTSymCramerInv.cxx:330
 TMatrixTSymCramerInv.cxx:331
 TMatrixTSymCramerInv.cxx:332
 TMatrixTSymCramerInv.cxx:333
 TMatrixTSymCramerInv.cxx:334
 TMatrixTSymCramerInv.cxx:335
 TMatrixTSymCramerInv.cxx:336
 TMatrixTSymCramerInv.cxx:337
 TMatrixTSymCramerInv.cxx:338
 TMatrixTSymCramerInv.cxx:339
 TMatrixTSymCramerInv.cxx:340
 TMatrixTSymCramerInv.cxx:341
 TMatrixTSymCramerInv.cxx:342
 TMatrixTSymCramerInv.cxx:343
 TMatrixTSymCramerInv.cxx:344
 TMatrixTSymCramerInv.cxx:345
 TMatrixTSymCramerInv.cxx:346
 TMatrixTSymCramerInv.cxx:347
 TMatrixTSymCramerInv.cxx:348
 TMatrixTSymCramerInv.cxx:349
 TMatrixTSymCramerInv.cxx:350
 TMatrixTSymCramerInv.cxx:351
 TMatrixTSymCramerInv.cxx:352
 TMatrixTSymCramerInv.cxx:353
 TMatrixTSymCramerInv.cxx:354
 TMatrixTSymCramerInv.cxx:355
 TMatrixTSymCramerInv.cxx:356
 TMatrixTSymCramerInv.cxx:357
 TMatrixTSymCramerInv.cxx:358
 TMatrixTSymCramerInv.cxx:359
 TMatrixTSymCramerInv.cxx:360
 TMatrixTSymCramerInv.cxx:361
 TMatrixTSymCramerInv.cxx:362
 TMatrixTSymCramerInv.cxx:363
 TMatrixTSymCramerInv.cxx:364
 TMatrixTSymCramerInv.cxx:365
 TMatrixTSymCramerInv.cxx:366
 TMatrixTSymCramerInv.cxx:367
 TMatrixTSymCramerInv.cxx:368
 TMatrixTSymCramerInv.cxx:369
 TMatrixTSymCramerInv.cxx:370
 TMatrixTSymCramerInv.cxx:371
 TMatrixTSymCramerInv.cxx:372
 TMatrixTSymCramerInv.cxx:373
 TMatrixTSymCramerInv.cxx:374
 TMatrixTSymCramerInv.cxx:375
 TMatrixTSymCramerInv.cxx:376
 TMatrixTSymCramerInv.cxx:377
 TMatrixTSymCramerInv.cxx:378
 TMatrixTSymCramerInv.cxx:379
 TMatrixTSymCramerInv.cxx:380
 TMatrixTSymCramerInv.cxx:381
 TMatrixTSymCramerInv.cxx:382
 TMatrixTSymCramerInv.cxx:383
 TMatrixTSymCramerInv.cxx:384
 TMatrixTSymCramerInv.cxx:385
 TMatrixTSymCramerInv.cxx:386
 TMatrixTSymCramerInv.cxx:387
 TMatrixTSymCramerInv.cxx:388
 TMatrixTSymCramerInv.cxx:389
 TMatrixTSymCramerInv.cxx:390
 TMatrixTSymCramerInv.cxx:391
 TMatrixTSymCramerInv.cxx:392
 TMatrixTSymCramerInv.cxx:393
 TMatrixTSymCramerInv.cxx:394
 TMatrixTSymCramerInv.cxx:395
 TMatrixTSymCramerInv.cxx:396
 TMatrixTSymCramerInv.cxx:397
 TMatrixTSymCramerInv.cxx:398
 TMatrixTSymCramerInv.cxx:399
 TMatrixTSymCramerInv.cxx:400
 TMatrixTSymCramerInv.cxx:401
 TMatrixTSymCramerInv.cxx:402
 TMatrixTSymCramerInv.cxx:403
 TMatrixTSymCramerInv.cxx:404
 TMatrixTSymCramerInv.cxx:405
 TMatrixTSymCramerInv.cxx:406
 TMatrixTSymCramerInv.cxx:407
 TMatrixTSymCramerInv.cxx:408
 TMatrixTSymCramerInv.cxx:409
 TMatrixTSymCramerInv.cxx:410
 TMatrixTSymCramerInv.cxx:411
 TMatrixTSymCramerInv.cxx:412
 TMatrixTSymCramerInv.cxx:413
 TMatrixTSymCramerInv.cxx:414
 TMatrixTSymCramerInv.cxx:415
 TMatrixTSymCramerInv.cxx:416
 TMatrixTSymCramerInv.cxx:417
 TMatrixTSymCramerInv.cxx:418
 TMatrixTSymCramerInv.cxx:419
 TMatrixTSymCramerInv.cxx:420
 TMatrixTSymCramerInv.cxx:421
 TMatrixTSymCramerInv.cxx:422
 TMatrixTSymCramerInv.cxx:423
 TMatrixTSymCramerInv.cxx:424
 TMatrixTSymCramerInv.cxx:425
 TMatrixTSymCramerInv.cxx:426
 TMatrixTSymCramerInv.cxx:427
 TMatrixTSymCramerInv.cxx:428
 TMatrixTSymCramerInv.cxx:429
 TMatrixTSymCramerInv.cxx:430
 TMatrixTSymCramerInv.cxx:431
 TMatrixTSymCramerInv.cxx:432
 TMatrixTSymCramerInv.cxx:433
 TMatrixTSymCramerInv.cxx:434
 TMatrixTSymCramerInv.cxx:435
 TMatrixTSymCramerInv.cxx:436
 TMatrixTSymCramerInv.cxx:437
 TMatrixTSymCramerInv.cxx:438
 TMatrixTSymCramerInv.cxx:439
 TMatrixTSymCramerInv.cxx:440
 TMatrixTSymCramerInv.cxx:441
 TMatrixTSymCramerInv.cxx:442
 TMatrixTSymCramerInv.cxx:443
 TMatrixTSymCramerInv.cxx:444
 TMatrixTSymCramerInv.cxx:445
 TMatrixTSymCramerInv.cxx:446
 TMatrixTSymCramerInv.cxx:447
 TMatrixTSymCramerInv.cxx:448
 TMatrixTSymCramerInv.cxx:449
 TMatrixTSymCramerInv.cxx:450
 TMatrixTSymCramerInv.cxx:451
 TMatrixTSymCramerInv.cxx:452
 TMatrixTSymCramerInv.cxx:453
 TMatrixTSymCramerInv.cxx:454
 TMatrixTSymCramerInv.cxx:455
 TMatrixTSymCramerInv.cxx:456
 TMatrixTSymCramerInv.cxx:457
 TMatrixTSymCramerInv.cxx:458
 TMatrixTSymCramerInv.cxx:459
 TMatrixTSymCramerInv.cxx:460
 TMatrixTSymCramerInv.cxx:461
 TMatrixTSymCramerInv.cxx:462
 TMatrixTSymCramerInv.cxx:463
 TMatrixTSymCramerInv.cxx:464
 TMatrixTSymCramerInv.cxx:465
 TMatrixTSymCramerInv.cxx:466
 TMatrixTSymCramerInv.cxx:467
 TMatrixTSymCramerInv.cxx:468
 TMatrixTSymCramerInv.cxx:469
 TMatrixTSymCramerInv.cxx:470
 TMatrixTSymCramerInv.cxx:471
 TMatrixTSymCramerInv.cxx:472
 TMatrixTSymCramerInv.cxx:473
 TMatrixTSymCramerInv.cxx:474
 TMatrixTSymCramerInv.cxx:475
 TMatrixTSymCramerInv.cxx:476
 TMatrixTSymCramerInv.cxx:477
 TMatrixTSymCramerInv.cxx:478
 TMatrixTSymCramerInv.cxx:479
 TMatrixTSymCramerInv.cxx:480
 TMatrixTSymCramerInv.cxx:481
 TMatrixTSymCramerInv.cxx:482
 TMatrixTSymCramerInv.cxx:483
 TMatrixTSymCramerInv.cxx:484
 TMatrixTSymCramerInv.cxx:485
 TMatrixTSymCramerInv.cxx:486
 TMatrixTSymCramerInv.cxx:487
 TMatrixTSymCramerInv.cxx:488
 TMatrixTSymCramerInv.cxx:489
 TMatrixTSymCramerInv.cxx:490
 TMatrixTSymCramerInv.cxx:491
 TMatrixTSymCramerInv.cxx:492
 TMatrixTSymCramerInv.cxx:493
 TMatrixTSymCramerInv.cxx:494
 TMatrixTSymCramerInv.cxx:495
 TMatrixTSymCramerInv.cxx:496
 TMatrixTSymCramerInv.cxx:497
 TMatrixTSymCramerInv.cxx:498
 TMatrixTSymCramerInv.cxx:499
 TMatrixTSymCramerInv.cxx:500
 TMatrixTSymCramerInv.cxx:501
 TMatrixTSymCramerInv.cxx:502
 TMatrixTSymCramerInv.cxx:503
 TMatrixTSymCramerInv.cxx:504
 TMatrixTSymCramerInv.cxx:505
 TMatrixTSymCramerInv.cxx:506
 TMatrixTSymCramerInv.cxx:507
 TMatrixTSymCramerInv.cxx:508
 TMatrixTSymCramerInv.cxx:509
 TMatrixTSymCramerInv.cxx:510
 TMatrixTSymCramerInv.cxx:511
 TMatrixTSymCramerInv.cxx:512
 TMatrixTSymCramerInv.cxx:513
 TMatrixTSymCramerInv.cxx:514
 TMatrixTSymCramerInv.cxx:515
 TMatrixTSymCramerInv.cxx:516
 TMatrixTSymCramerInv.cxx:517
 TMatrixTSymCramerInv.cxx:518
 TMatrixTSymCramerInv.cxx:519
 TMatrixTSymCramerInv.cxx:520
 TMatrixTSymCramerInv.cxx:521
 TMatrixTSymCramerInv.cxx:522
 TMatrixTSymCramerInv.cxx:523
 TMatrixTSymCramerInv.cxx:524
 TMatrixTSymCramerInv.cxx:525
 TMatrixTSymCramerInv.cxx:526
 TMatrixTSymCramerInv.cxx:527
 TMatrixTSymCramerInv.cxx:528
 TMatrixTSymCramerInv.cxx:529
 TMatrixTSymCramerInv.cxx:530
 TMatrixTSymCramerInv.cxx:531
 TMatrixTSymCramerInv.cxx:532
 TMatrixTSymCramerInv.cxx:533
 TMatrixTSymCramerInv.cxx:534
 TMatrixTSymCramerInv.cxx:535
 TMatrixTSymCramerInv.cxx:536
 TMatrixTSymCramerInv.cxx:537
 TMatrixTSymCramerInv.cxx:538
 TMatrixTSymCramerInv.cxx:539
 TMatrixTSymCramerInv.cxx:540
 TMatrixTSymCramerInv.cxx:541
 TMatrixTSymCramerInv.cxx:542
 TMatrixTSymCramerInv.cxx:543
 TMatrixTSymCramerInv.cxx:544
 TMatrixTSymCramerInv.cxx:545
 TMatrixTSymCramerInv.cxx:546
 TMatrixTSymCramerInv.cxx:547
 TMatrixTSymCramerInv.cxx:548
 TMatrixTSymCramerInv.cxx:549
 TMatrixTSymCramerInv.cxx:550
 TMatrixTSymCramerInv.cxx:551
 TMatrixTSymCramerInv.cxx:552
 TMatrixTSymCramerInv.cxx:553
 TMatrixTSymCramerInv.cxx:554
 TMatrixTSymCramerInv.cxx:555
 TMatrixTSymCramerInv.cxx:556
 TMatrixTSymCramerInv.cxx:557
 TMatrixTSymCramerInv.cxx:558
 TMatrixTSymCramerInv.cxx:559
 TMatrixTSymCramerInv.cxx:560
 TMatrixTSymCramerInv.cxx:561
 TMatrixTSymCramerInv.cxx:562
 TMatrixTSymCramerInv.cxx:563
 TMatrixTSymCramerInv.cxx:564
 TMatrixTSymCramerInv.cxx:565
 TMatrixTSymCramerInv.cxx:566
 TMatrixTSymCramerInv.cxx:567
 TMatrixTSymCramerInv.cxx:568
 TMatrixTSymCramerInv.cxx:569
 TMatrixTSymCramerInv.cxx:570
 TMatrixTSymCramerInv.cxx:571
 TMatrixTSymCramerInv.cxx:572
 TMatrixTSymCramerInv.cxx:573
 TMatrixTSymCramerInv.cxx:574
 TMatrixTSymCramerInv.cxx:575
 TMatrixTSymCramerInv.cxx:576
 TMatrixTSymCramerInv.cxx:577
 TMatrixTSymCramerInv.cxx:578
 TMatrixTSymCramerInv.cxx:579
 TMatrixTSymCramerInv.cxx:580
 TMatrixTSymCramerInv.cxx:581
 TMatrixTSymCramerInv.cxx:582
 TMatrixTSymCramerInv.cxx:583
 TMatrixTSymCramerInv.cxx:584
 TMatrixTSymCramerInv.cxx:585
 TMatrixTSymCramerInv.cxx:586
 TMatrixTSymCramerInv.cxx:587
 TMatrixTSymCramerInv.cxx:588
 TMatrixTSymCramerInv.cxx:589
 TMatrixTSymCramerInv.cxx:590
 TMatrixTSymCramerInv.cxx:591
 TMatrixTSymCramerInv.cxx:592
 TMatrixTSymCramerInv.cxx:593
 TMatrixTSymCramerInv.cxx:594
 TMatrixTSymCramerInv.cxx:595
 TMatrixTSymCramerInv.cxx:596
 TMatrixTSymCramerInv.cxx:597
 TMatrixTSymCramerInv.cxx:598
 TMatrixTSymCramerInv.cxx:599
 TMatrixTSymCramerInv.cxx:600
 TMatrixTSymCramerInv.cxx:601
 TMatrixTSymCramerInv.cxx:602
 TMatrixTSymCramerInv.cxx:603
 TMatrixTSymCramerInv.cxx:604
 TMatrixTSymCramerInv.cxx:605
 TMatrixTSymCramerInv.cxx:606
 TMatrixTSymCramerInv.cxx:607
 TMatrixTSymCramerInv.cxx:608
 TMatrixTSymCramerInv.cxx:609
 TMatrixTSymCramerInv.cxx:610
 TMatrixTSymCramerInv.cxx:611
 TMatrixTSymCramerInv.cxx:612
 TMatrixTSymCramerInv.cxx:613
 TMatrixTSymCramerInv.cxx:614
 TMatrixTSymCramerInv.cxx:615
 TMatrixTSymCramerInv.cxx:616
 TMatrixTSymCramerInv.cxx:617
 TMatrixTSymCramerInv.cxx:618
 TMatrixTSymCramerInv.cxx:619
 TMatrixTSymCramerInv.cxx:620
 TMatrixTSymCramerInv.cxx:621
 TMatrixTSymCramerInv.cxx:622
 TMatrixTSymCramerInv.cxx:623
 TMatrixTSymCramerInv.cxx:624
 TMatrixTSymCramerInv.cxx:625
 TMatrixTSymCramerInv.cxx:626
 TMatrixTSymCramerInv.cxx:627
 TMatrixTSymCramerInv.cxx:628
 TMatrixTSymCramerInv.cxx:629
 TMatrixTSymCramerInv.cxx:630
 TMatrixTSymCramerInv.cxx:631
 TMatrixTSymCramerInv.cxx:632
 TMatrixTSymCramerInv.cxx:633
 TMatrixTSymCramerInv.cxx:634
 TMatrixTSymCramerInv.cxx:635
 TMatrixTSymCramerInv.cxx:636
 TMatrixTSymCramerInv.cxx:637
 TMatrixTSymCramerInv.cxx:638
 TMatrixTSymCramerInv.cxx:639
 TMatrixTSymCramerInv.cxx:640
 TMatrixTSymCramerInv.cxx:641
 TMatrixTSymCramerInv.cxx:642
 TMatrixTSymCramerInv.cxx:643
 TMatrixTSymCramerInv.cxx:644
 TMatrixTSymCramerInv.cxx:645
 TMatrixTSymCramerInv.cxx:646
 TMatrixTSymCramerInv.cxx:647
 TMatrixTSymCramerInv.cxx:648
 TMatrixTSymCramerInv.cxx:649
 TMatrixTSymCramerInv.cxx:650
 TMatrixTSymCramerInv.cxx:651
 TMatrixTSymCramerInv.cxx:652
 TMatrixTSymCramerInv.cxx:653
 TMatrixTSymCramerInv.cxx:654
 TMatrixTSymCramerInv.cxx:655
 TMatrixTSymCramerInv.cxx:656
 TMatrixTSymCramerInv.cxx:657
 TMatrixTSymCramerInv.cxx:658
 TMatrixTSymCramerInv.cxx:659
 TMatrixTSymCramerInv.cxx:660
 TMatrixTSymCramerInv.cxx:661
 TMatrixTSymCramerInv.cxx:662
 TMatrixTSymCramerInv.cxx:663
 TMatrixTSymCramerInv.cxx:664
 TMatrixTSymCramerInv.cxx:665
 TMatrixTSymCramerInv.cxx:666
 TMatrixTSymCramerInv.cxx:667
 TMatrixTSymCramerInv.cxx:668
 TMatrixTSymCramerInv.cxx:669
 TMatrixTSymCramerInv.cxx:670
 TMatrixTSymCramerInv.cxx:671
 TMatrixTSymCramerInv.cxx:672
 TMatrixTSymCramerInv.cxx:673
 TMatrixTSymCramerInv.cxx:674
 TMatrixTSymCramerInv.cxx:675
 TMatrixTSymCramerInv.cxx:676
 TMatrixTSymCramerInv.cxx:677
 TMatrixTSymCramerInv.cxx:678
 TMatrixTSymCramerInv.cxx:679
 TMatrixTSymCramerInv.cxx:680
 TMatrixTSymCramerInv.cxx:681
 TMatrixTSymCramerInv.cxx:682
 TMatrixTSymCramerInv.cxx:683
 TMatrixTSymCramerInv.cxx:684
 TMatrixTSymCramerInv.cxx:685
 TMatrixTSymCramerInv.cxx:686
 TMatrixTSymCramerInv.cxx:687
 TMatrixTSymCramerInv.cxx:688
 TMatrixTSymCramerInv.cxx:689
 TMatrixTSymCramerInv.cxx:690
 TMatrixTSymCramerInv.cxx:691
 TMatrixTSymCramerInv.cxx:692
 TMatrixTSymCramerInv.cxx:693
 TMatrixTSymCramerInv.cxx:694
 TMatrixTSymCramerInv.cxx:695
 TMatrixTSymCramerInv.cxx:696
 TMatrixTSymCramerInv.cxx:697
 TMatrixTSymCramerInv.cxx:698
 TMatrixTSymCramerInv.cxx:699
 TMatrixTSymCramerInv.cxx:700
 TMatrixTSymCramerInv.cxx:701
 TMatrixTSymCramerInv.cxx:702
 TMatrixTSymCramerInv.cxx:703
 TMatrixTSymCramerInv.cxx:704
 TMatrixTSymCramerInv.cxx:705
 TMatrixTSymCramerInv.cxx:706
 TMatrixTSymCramerInv.cxx:707
 TMatrixTSymCramerInv.cxx:708
 TMatrixTSymCramerInv.cxx:709
 TMatrixTSymCramerInv.cxx:710
 TMatrixTSymCramerInv.cxx:711
 TMatrixTSymCramerInv.cxx:712
 TMatrixTSymCramerInv.cxx:713
 TMatrixTSymCramerInv.cxx:714
 TMatrixTSymCramerInv.cxx:715
 TMatrixTSymCramerInv.cxx:716
 TMatrixTSymCramerInv.cxx:717
 TMatrixTSymCramerInv.cxx:718
 TMatrixTSymCramerInv.cxx:719
 TMatrixTSymCramerInv.cxx:720
 TMatrixTSymCramerInv.cxx:721
 TMatrixTSymCramerInv.cxx:722
 TMatrixTSymCramerInv.cxx:723
 TMatrixTSymCramerInv.cxx:724
 TMatrixTSymCramerInv.cxx:725
 TMatrixTSymCramerInv.cxx:726
 TMatrixTSymCramerInv.cxx:727
 TMatrixTSymCramerInv.cxx:728
 TMatrixTSymCramerInv.cxx:729
 TMatrixTSymCramerInv.cxx:730
 TMatrixTSymCramerInv.cxx:731
 TMatrixTSymCramerInv.cxx:732
 TMatrixTSymCramerInv.cxx:733
 TMatrixTSymCramerInv.cxx:734
 TMatrixTSymCramerInv.cxx:735
 TMatrixTSymCramerInv.cxx:736
 TMatrixTSymCramerInv.cxx:737
 TMatrixTSymCramerInv.cxx:738
 TMatrixTSymCramerInv.cxx:739
 TMatrixTSymCramerInv.cxx:740
 TMatrixTSymCramerInv.cxx:741
 TMatrixTSymCramerInv.cxx:742
 TMatrixTSymCramerInv.cxx:743
 TMatrixTSymCramerInv.cxx:744
 TMatrixTSymCramerInv.cxx:745
 TMatrixTSymCramerInv.cxx:746
 TMatrixTSymCramerInv.cxx:747
 TMatrixTSymCramerInv.cxx:748
 TMatrixTSymCramerInv.cxx:749
 TMatrixTSymCramerInv.cxx:750
 TMatrixTSymCramerInv.cxx:751
 TMatrixTSymCramerInv.cxx:752
 TMatrixTSymCramerInv.cxx:753
 TMatrixTSymCramerInv.cxx:754
 TMatrixTSymCramerInv.cxx:755
 TMatrixTSymCramerInv.cxx:756
 TMatrixTSymCramerInv.cxx:757
 TMatrixTSymCramerInv.cxx:758
 TMatrixTSymCramerInv.cxx:759
 TMatrixTSymCramerInv.cxx:760
 TMatrixTSymCramerInv.cxx:761
 TMatrixTSymCramerInv.cxx:762
 TMatrixTSymCramerInv.cxx:763
 TMatrixTSymCramerInv.cxx:764
 TMatrixTSymCramerInv.cxx:765
 TMatrixTSymCramerInv.cxx:766
 TMatrixTSymCramerInv.cxx:767
 TMatrixTSymCramerInv.cxx:768
 TMatrixTSymCramerInv.cxx:769
 TMatrixTSymCramerInv.cxx:770
 TMatrixTSymCramerInv.cxx:771
 TMatrixTSymCramerInv.cxx:772
 TMatrixTSymCramerInv.cxx:773
 TMatrixTSymCramerInv.cxx:774
 TMatrixTSymCramerInv.cxx:775
 TMatrixTSymCramerInv.cxx:776
 TMatrixTSymCramerInv.cxx:777
 TMatrixTSymCramerInv.cxx:778
 TMatrixTSymCramerInv.cxx:779
 TMatrixTSymCramerInv.cxx:780
 TMatrixTSymCramerInv.cxx:781
 TMatrixTSymCramerInv.cxx:782
 TMatrixTSymCramerInv.cxx:783
 TMatrixTSymCramerInv.cxx:784
 TMatrixTSymCramerInv.cxx:785
 TMatrixTSymCramerInv.cxx:786
 TMatrixTSymCramerInv.cxx:787
 TMatrixTSymCramerInv.cxx:788
 TMatrixTSymCramerInv.cxx:789
 TMatrixTSymCramerInv.cxx:790
 TMatrixTSymCramerInv.cxx:791
 TMatrixTSymCramerInv.cxx:792
 TMatrixTSymCramerInv.cxx:793
 TMatrixTSymCramerInv.cxx:794
 TMatrixTSymCramerInv.cxx:795
 TMatrixTSymCramerInv.cxx:796
 TMatrixTSymCramerInv.cxx:797
 TMatrixTSymCramerInv.cxx:798
 TMatrixTSymCramerInv.cxx:799
 TMatrixTSymCramerInv.cxx:800
 TMatrixTSymCramerInv.cxx:801
 TMatrixTSymCramerInv.cxx:802
 TMatrixTSymCramerInv.cxx:803
 TMatrixTSymCramerInv.cxx:804
 TMatrixTSymCramerInv.cxx:805
 TMatrixTSymCramerInv.cxx:806
 TMatrixTSymCramerInv.cxx:807
 TMatrixTSymCramerInv.cxx:808
 TMatrixTSymCramerInv.cxx:809
 TMatrixTSymCramerInv.cxx:810
 TMatrixTSymCramerInv.cxx:811
 TMatrixTSymCramerInv.cxx:812
 TMatrixTSymCramerInv.cxx:813
 TMatrixTSymCramerInv.cxx:814
 TMatrixTSymCramerInv.cxx:815
 TMatrixTSymCramerInv.cxx:816
 TMatrixTSymCramerInv.cxx:817
 TMatrixTSymCramerInv.cxx:818
 TMatrixTSymCramerInv.cxx:819
 TMatrixTSymCramerInv.cxx:820
 TMatrixTSymCramerInv.cxx:821
 TMatrixTSymCramerInv.cxx:822
 TMatrixTSymCramerInv.cxx:823
 TMatrixTSymCramerInv.cxx:824
 TMatrixTSymCramerInv.cxx:825