// cl_LF_to_FF().
// General includes.
#include "cl_sysdep.h"
// Specification.
#include "cl_F.h"
// Implementation.
#include "cl_LF.h"
#include "cl_LF_impl.h"
#include "cl_FF.h"
#include "cl_DS.h"
namespace cln {
const cl_FF cl_LF_to_FF (const cl_LF& x)
{
// x entpacken:
var cl_signean sign;
var sintL exp;
var uintD* ptr;
var uintC len;
LF_decode(x, { return cl_FF_0; }, sign=,exp=,ptr=,len=,);
// intDsize*len-FF_mant_len-1 Bits der Mantisse wegrunden:
// erste k := ceiling(FF_mant_len+2,intDsize) Digits nach mant holen:
#if (intDsize==64)
var uint64 mant = get_max64_Dptr(FF_mant_len+2,ptr);
#else
var uint32 mant = get_max32_Dptr(FF_mant_len+2,ptr);
#endif
ptr = ptr mspop ceiling(FF_mant_len+2,intDsize);
var const int shiftcount = ceiling(FF_mant_len+2,intDsize)*intDsize-(FF_mant_len+1);
if ( ((mant & bit(shiftcount-1)) ==0) // Bit 7 war 0 -> abrunden
|| ( ((mant & (bit(shiftcount-1)-1)) ==0) // war 1, Bits 6..0 >0 -> aufrunden
&& !test_loop_msp(ptr,len-ceiling(FF_mant_len+2,intDsize)) // weitere Bits /=0 -> aufrunden
// round-to-even
&& ((mant & bit(shiftcount)) ==0)
) )
// abrunden
{ mant = mant >> shiftcount; }
else
// aufrunden
{ mant = mant >> shiftcount;
mant = mant+1;
if (mant >= bit(FF_mant_len+1))
// Überlauf durchs Runden
{ mant = mant>>1; exp = exp+1; } // Mantisse rechts schieben
}
return encode_FF(sign,exp,mant);
}
} // namespace cln
syntax highlighted by Code2HTML, v. 0.9.1