6 CL_PROVIDE(cl_F_epspos)
18 #include "cl_LF_impl.h"
22 // Bei Floats mit d Bits (incl. Hiddem Bit, also d = ?F_mant_len+1)
23 // ist ?F_epsilon = 2^-d*(1+2^(1-d)), d.h. Mantisse 10...01, Vorzeichen +.
25 static const cl_SF SF_epsilon =
26 make_SF(0,SF_exp_mid-SF_mant_len,bit(SF_mant_len)+1);
28 static const cl_FF FF_epsilon =
29 encode_FF(0,-FF_mant_len,bit(FF_mant_len)+1);
31 static const cl_DF DF_epsilon =
32 #if (cl_word_size==64)
33 encode_DF(0,-DF_mant_len,bit(DF_mant_len)+1);
35 encode_DF(0,-DF_mant_len,bit(DF_mant_len-32),1);
38 inline const cl_LF LF_epsilon (uintC len)
40 var Lfloat erg = allocate_lfloat(len,LF_exp_mid+1-intDsize*len,0);
41 var uintD* ptr = &TheLfloat(erg)->data[0];
42 #if CL_DS_BIG_ENDIAN_P
43 *ptr++ = bit(intDsize-1);
44 ptr = clear_loop_up(ptr,len-2);
48 ptr = clear_loop_up(ptr,len-2);
49 *ptr = bit(intDsize-1);
54 const cl_F float_epsilon (float_format_t f)
56 floatformatcase((uintC)f
60 , return LF_epsilon(len);
66 CL_PROVIDE_END(cl_F_epspos)