dist = ifelse( ite == 1, llh, (llh-last_llh) / (1+abs(llh)) )
Dist1 = max( (abs(phi-Phi)) / (1+abs(phi)) )
Dist2 = max( (abs(rho-Rho)) / (1+abs(rho)) )
Dist3 = max( (abs(pi-Pi)) / (1+abs(Pi)) )
dist2 = max(Dist1,Dist2,Dist3)
dist = ifelse( ite == 1, llh, (llh-last_llh) / (1+abs(llh)) )
Dist1 = max( (abs(phi-Phi)) / (1+abs(phi)) )
Dist2 = max( (abs(rho-Rho)) / (1+abs(rho)) )
Dist3 = max( (abs(pi-Pi)) / (1+abs(Pi)) )
dist2 = max(Dist1,Dist2,Dist3)