Impala Documentation

template<class DstArrayT, class SrcArrayT>

void Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::Done ( int phase ) [inline]

Done with one phase.

Definition at line 160 of file IntWeibullNgbPnLoop.h.

References Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::DiGamma(), Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::mF, Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::mFprime, Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::mGamma, Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::mGammaNext, Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::mVerbose, and Impala::Core::Feature::IntWeibullNgbPnLoop< DstArrayT, SrcArrayT >::TriGamma().

00161     {
00162       int ph = phase%2;
00163       if (mVerbose)
00164         std::cout << "done: phase = " << phase <<  std::endl;
00165       if(!ph)
00166       {
00167         if (mVerbose)
00168           std::cout << "mGamma = " << mGamma << std::endl;
00169         
00170         mFprime = (1/pow(mGamma,2))*(1 - log(mGamma) - DiGamma(1/mGamma) - 
00171                                      (1/mGamma)*TriGamma(1/mGamma) + mF) -(1/mGamma)*mFprime;
00172         mF = 1 + (1/mGamma)*(log(mGamma) + DiGamma(1/mGamma) - mF);
00173         mGammaNext = mGamma - mF/mFprime;
00174         if (mVerbose)
00175           std::cout << "mGammaNext = " << mGammaNext << std::endl;
00176       }
00177       
00178     }

Here is the call graph for this function: