6 CL_PROVIDE(cl_F_epspos)
18 #include "cl_LF_impl.h"
20 // Bei Floats mit d Bits (incl. Hiddem Bit, also d = ?F_mant_len+1)
21 // ist ?F_epsilon = 2^-d*(1+2^(1-d)), d.h. Mantisse 10...01, Vorzeichen +.
23 static const cl_SF SF_epsilon =
24 make_SF(0,SF_exp_mid-SF_mant_len,bit(SF_mant_len)+1);
26 static const cl_FF FF_epsilon =
27 encode_FF(0,-FF_mant_len,bit(FF_mant_len)+1);
29 static const cl_DF DF_epsilon =
30 #if (cl_word_size==64)
31 encode_DF(0,-DF_mant_len,bit(DF_mant_len)+1);
33 encode_DF(0,-DF_mant_len,bit(DF_mant_len-32),1);
36 inline const cl_LF LF_epsilon (uintC len)
38 var Lfloat erg = allocate_lfloat(len,LF_exp_mid+1-intDsize*(uintL)len,0);
39 var uintD* ptr = &TheLfloat(erg)->data[0];
40 #if CL_DS_BIG_ENDIAN_P
41 *ptr++ = bit(intDsize-1);
42 ptr = clear_loop_up(ptr,len-2);
46 ptr = clear_loop_up(ptr,len-2);
47 *ptr = bit(intDsize-1);
52 const cl_F float_epsilon (cl_float_format_t f)
54 floatformatcase((uintL)f
58 , return LF_epsilon(len);
62 CL_PROVIDE_END(cl_F_epspos)