*! v 1.0.2 11Nov96. STB-34 sbe13
program define xri_slll
        version 4.0
local ll "`1'"
local M "`2'"
local S "`3'"
if "`4'"=="" { local G $S_gamma }
else local G "`4'"

tempname ln2pi small
scalar `ln2pi'=1.837877
scalar `small'=1e-9

quietly {
	if !$S_cv & !$S_lns {
		local s `S' /* sigma parametrization */
	}
	else {
		tempvar s
		if $S_lns { /* log S parametrization */
			if $S_cv { gen double `s' = `M'*exp(`S') }
			else gen double `s' = exp(`S')
		}
		else gen double `s' = `M'*`S'
	}

* s(igma) is scale factor

	tempvar Gz
	cap drop __U
	gen double `Gz' = `G'*($S_mldepn-`M')/`s'
	gen double __U = cond(abs(`Gz')<`small', /*
	 */ (1-0.5*`Gz')*($S_mldepn-`M')/`s', log(1+`Gz')/`G' )
	replace `ll' = -0.5*`ln2pi'-log(`s')-0.5*__U^2-`G'*__U
	if $S_trunc { replace `ll' = `ll'-cond(abs(`G')>`small', /*
	 */ log(1-normprob(log(1-`M'*`G'/`s')/`G')), /*
	 */ log(1-normprob(-`M'/`s')) ) }
}
end