[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

matrix is singular error II



Hi all,

So I've been fighting with this for awhile now and and I've come seeking
offline guru wisdom. My analysis code generates the following error and
I've been unable to figure out how to correct it or where it comes from.
Any thoughts?

Error in <TDecompLU::DecomposeLUCrout>: matrix is singular
Error in <TDecompLU::DecomposeLUCrout>: matrix is singular
Error in <TDecompLU::InvertLU>: matrix is singular, 3 diag elements <
tolerance of 2.2204e-16
Fatal: b.IsValid() violated at line 176 of `matrix/src/TMatrixT.cxx'
aborting
Generating stack trace...
JApplication.cc:858 didn't sleep full 0.5 seconds!
  0x039c2d8a in TMatrixT<double> operator*<double>(TMatrixT<double> 
const&, TMatrixT<double> const&) + 0x22 from 
/usr/local/root/lib/libMatrix.so
  0x08132021 in DKinFit::FitTwoGammas(float, float) + 0x334d from 
/s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x0812aa4c in DTwoGammaFit_factory::evnt(jana::JEventLoop*, int) + 
0x15c from /s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x0809a391 in 
jana::JFactory<DTwoGammaFit>::Get(std::vector<DTwoGammaFit const*, 
std::allocator<DTwoGammaFit const*> >&) + 0x75 from 
/s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x0809a8a1 in 
_ZN4jana10JEventLoop14GetFromFactoryI12DTwoGammaFitEEPNS_8JFactoryIT_EERSt6vectorIPKS4_SaIS9_EEPKcRN 
+ 0x1d1 from /s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x0809aa88 in jana::JFactory<DTwoGammaFit>* 
jana::JEventLoop::Get<DTwoGammaFit>(std::vector<DTwoGammaFit const*, 
std::allocator<DTwoGammaFit const*> >&, char const*) + 0xdc from 
/s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x0808e524 in FCALTreeProc::evnt(jana::JEventLoop*, int) + 0xf70 from 
/s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x080b5981 in jana::JEventLoop::OneEvent() + 0x28b from 
/s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x080b5b4c in jana::JEventLoop::Loop() + 0x56 from 
/s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x0809f89a in LaunchThread(void*) + 0x36 from 
/s4/leverin/gluex/my_src/bin/Linux_RH-athlon-gcc3.4.4/fcalTree
  0x0026c341 in <unknown> from /lib/tls/libpthread.so.0
  0x02cd16fe in __clone + 0x5e from /lib/tls/libc.so.6
Abort





I included the stack trace this time

thanks,
-Blake