154359Sroberto/* 254359Sroberto * ntp_fp.h - definitions for NTP fixed/floating-point arithmetic 354359Sroberto */ 454359Sroberto 554359Sroberto#ifndef NTP_FP_H 654359Sroberto#define NTP_FP_H 754359Sroberto 854359Sroberto#include <sys/types.h> 954359Sroberto#include <sys/socket.h> 1054359Sroberto#include <netinet/in.h> 1154359Sroberto 12182007Sroberto#include "ntp_rfc2553.h" 13182007Sroberto 1454359Sroberto#include "ntp_types.h" 1554359Sroberto 1654359Sroberto/* 1754359Sroberto * NTP uses two fixed point formats. The first (l_fp) is the "long" 1854359Sroberto * format and is 64 bits long with the decimal between bits 31 and 32. 1954359Sroberto * This is used for time stamps in the NTP packet header (in network 2054359Sroberto * byte order) and for internal computations of offsets (in local host 2154359Sroberto * byte order). We use the same structure for both signed and unsigned 2254359Sroberto * values, which is a big hack but saves rewriting all the operators 2354359Sroberto * twice. Just to confuse this, we also sometimes just carry the 2454359Sroberto * fractional part in calculations, in both signed and unsigned forms. 2554359Sroberto * Anyway, an l_fp looks like: 2654359Sroberto * 2754359Sroberto * 0 1 2 3 2854359Sroberto * 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2954359Sroberto * +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ 3054359Sroberto * | Integral Part | 3154359Sroberto * +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ 3254359Sroberto * | Fractional Part | 3354359Sroberto * +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ 3454359Sroberto * 3554359Sroberto */ 3654359Srobertotypedef struct { 3754359Sroberto union { 3854359Sroberto u_int32 Xl_ui; 3954359Sroberto int32 Xl_i; 4054359Sroberto } Ul_i; 4154359Sroberto union { 4254359Sroberto u_int32 Xl_uf; 4354359Sroberto int32 Xl_f; 4454359Sroberto } Ul_f; 4554359Sroberto} l_fp; 4654359Sroberto 4754359Sroberto#define l_ui Ul_i.Xl_ui /* unsigned integral part */ 4854359Sroberto#define l_i Ul_i.Xl_i /* signed integral part */ 4954359Sroberto#define l_uf Ul_f.Xl_uf /* unsigned fractional part */ 5054359Sroberto#define l_f Ul_f.Xl_f /* signed fractional part */ 5154359Sroberto 5254359Sroberto/* 5354359Sroberto * Fractional precision (of an l_fp) is actually the number of 5454359Sroberto * bits in a long. 5554359Sroberto */ 5654359Sroberto#define FRACTION_PREC (32) 5754359Sroberto 5854359Sroberto 5954359Sroberto/* 6054359Sroberto * The second fixed point format is 32 bits, with the decimal between 6154359Sroberto * bits 15 and 16. There is a signed version (s_fp) and an unsigned 6254359Sroberto * version (u_fp). This is used to represent synchronizing distance 6354359Sroberto * and synchronizing dispersion in the NTP packet header (again, in 6454359Sroberto * network byte order) and internally to hold both distance and 6554359Sroberto * dispersion values (in local byte order). In network byte order 6654359Sroberto * it looks like: 6754359Sroberto * 6854359Sroberto * 0 1 2 3 6954359Sroberto * 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 7054359Sroberto * +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ 7154359Sroberto * | Integer Part | Fraction Part | 7254359Sroberto * +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ 7354359Sroberto * 7454359Sroberto */ 7554359Srobertotypedef int32 s_fp; 7654359Srobertotypedef u_int32 u_fp; 7754359Sroberto 7854359Sroberto/* 7954359Sroberto * A unit second in fp format. Actually 2**(half_the_bits_in_a_long) 8054359Sroberto */ 8154359Sroberto#define FP_SECOND (0x10000) 8254359Sroberto 8354359Sroberto/* 8454359Sroberto * Byte order conversions 8554359Sroberto */ 8654359Sroberto#define HTONS_FP(x) (htonl(x)) 8754359Sroberto#define HTONL_FP(h, n) do { (n)->l_ui = htonl((h)->l_ui); \ 8854359Sroberto (n)->l_uf = htonl((h)->l_uf); } while (0) 8954359Sroberto#define NTOHS_FP(x) (ntohl(x)) 9054359Sroberto#define NTOHL_FP(n, h) do { (h)->l_ui = ntohl((n)->l_ui); \ 9154359Sroberto (h)->l_uf = ntohl((n)->l_uf); } while (0) 9254359Sroberto#define NTOHL_MFP(ni, nf, hi, hf) \ 9354359Sroberto do { (hi) = ntohl(ni); (hf) = ntohl(nf); } while (0) 9454359Sroberto#define HTONL_MFP(hi, hf, ni, nf) \ 9554359Sroberto do { (ni) = ntohl(hi); (nf) = ntohl(hf); } while (0) 9654359Sroberto 9754359Sroberto/* funny ones. Converts ts fractions to net order ts */ 9854359Sroberto#define HTONL_UF(uf, nts) \ 9954359Sroberto do { (nts)->l_ui = 0; (nts)->l_uf = htonl(uf); } while (0) 10054359Sroberto#define HTONL_F(f, nts) do { (nts)->l_uf = htonl(f); \ 10154359Sroberto if ((f) & 0x80000000) \ 10254359Sroberto (nts)->l_i = -1; \ 10354359Sroberto else \ 10454359Sroberto (nts)->l_i = 0; \ 10554359Sroberto } while (0) 10654359Sroberto 10754359Sroberto/* 10854359Sroberto * Conversions between the two fixed point types 10954359Sroberto */ 11054359Sroberto#define MFPTOFP(x_i, x_f) (((x_i) >= 0x00010000) ? 0x7fffffff : \ 11154359Sroberto (((x_i) <= -0x00010000) ? 0x80000000 : \ 11254359Sroberto (((x_i)<<16) | (((x_f)>>16)&0xffff)))) 11354359Sroberto#define LFPTOFP(v) MFPTOFP((v)->l_i, (v)->l_f) 11454359Sroberto 11554359Sroberto#define UFPTOLFP(x, v) ((v)->l_ui = (u_fp)(x)>>16, (v)->l_uf = (x)<<16) 11654359Sroberto#define FPTOLFP(x, v) (UFPTOLFP((x), (v)), (x) < 0 ? (v)->l_ui -= 0x10000 : 0) 11754359Sroberto 11854359Sroberto#define MAXLFP(v) ((v)->l_ui = 0x7fffffff, (v)->l_uf = 0xffffffff) 11954359Sroberto#define MINLFP(v) ((v)->l_ui = 0x80000000, (v)->l_uf = 0) 12054359Sroberto 12154359Sroberto/* 12254359Sroberto * Primitive operations on long fixed point values. If these are 12354359Sroberto * reminiscent of assembler op codes it's only because some may 12454359Sroberto * be replaced by inline assembler for particular machines someday. 12554359Sroberto * These are the (kind of inefficient) run-anywhere versions. 12654359Sroberto */ 12754359Sroberto#define M_NEG(v_i, v_f) /* v = -v */ \ 12854359Sroberto do { \ 12954359Sroberto if ((v_f) == 0) \ 13054359Sroberto (v_i) = -((s_fp)(v_i)); \ 13154359Sroberto else { \ 13254359Sroberto (v_f) = -((s_fp)(v_f)); \ 13354359Sroberto (v_i) = ~(v_i); \ 13454359Sroberto } \ 13554359Sroberto } while(0) 13654359Sroberto 13754359Sroberto#define M_NEGM(r_i, r_f, a_i, a_f) /* r = -a */ \ 13854359Sroberto do { \ 13954359Sroberto if ((a_f) == 0) { \ 14054359Sroberto (r_f) = 0; \ 14154359Sroberto (r_i) = -(a_i); \ 14254359Sroberto } else { \ 14354359Sroberto (r_f) = -(a_f); \ 14454359Sroberto (r_i) = ~(a_i); \ 14554359Sroberto } \ 14654359Sroberto } while(0) 14754359Sroberto 14854359Sroberto#define M_ADD(r_i, r_f, a_i, a_f) /* r += a */ \ 14954359Sroberto do { \ 15054359Sroberto register u_int32 lo_tmp; \ 15154359Sroberto register u_int32 hi_tmp; \ 15254359Sroberto \ 15354359Sroberto lo_tmp = ((r_f) & 0xffff) + ((a_f) & 0xffff); \ 15454359Sroberto hi_tmp = (((r_f) >> 16) & 0xffff) + (((a_f) >> 16) & 0xffff); \ 15554359Sroberto if (lo_tmp & 0x10000) \ 15654359Sroberto hi_tmp++; \ 15754359Sroberto (r_f) = ((hi_tmp & 0xffff) << 16) | (lo_tmp & 0xffff); \ 15854359Sroberto \ 15954359Sroberto (r_i) += (a_i); \ 16054359Sroberto if (hi_tmp & 0x10000) \ 16154359Sroberto (r_i)++; \ 16254359Sroberto } while (0) 16354359Sroberto 16454359Sroberto#define M_ADD3(r_ovr, r_i, r_f, a_ovr, a_i, a_f) /* r += a, three word */ \ 16554359Sroberto do { \ 16654359Sroberto register u_int32 lo_tmp; \ 16754359Sroberto register u_int32 hi_tmp; \ 16854359Sroberto \ 16954359Sroberto lo_tmp = ((r_f) & 0xffff) + ((a_f) & 0xffff); \ 17054359Sroberto hi_tmp = (((r_f) >> 16) & 0xffff) + (((a_f) >> 16) & 0xffff); \ 17154359Sroberto if (lo_tmp & 0x10000) \ 17254359Sroberto hi_tmp++; \ 17354359Sroberto (r_f) = ((hi_tmp & 0xffff) << 16) | (lo_tmp & 0xffff); \ 17454359Sroberto \ 17554359Sroberto lo_tmp = ((r_i) & 0xffff) + ((a_i) & 0xffff); \ 17654359Sroberto if (hi_tmp & 0x10000) \ 17754359Sroberto lo_tmp++; \ 17854359Sroberto hi_tmp = (((r_i) >> 16) & 0xffff) + (((a_i) >> 16) & 0xffff); \ 17954359Sroberto if (lo_tmp & 0x10000) \ 18054359Sroberto hi_tmp++; \ 18154359Sroberto (r_i) = ((hi_tmp & 0xffff) << 16) | (lo_tmp & 0xffff); \ 18254359Sroberto \ 18354359Sroberto (r_ovr) += (a_ovr); \ 18454359Sroberto if (hi_tmp & 0x10000) \ 18554359Sroberto (r_ovr)++; \ 18654359Sroberto } while (0) 18754359Sroberto 18854359Sroberto#define M_SUB(r_i, r_f, a_i, a_f) /* r -= a */ \ 18954359Sroberto do { \ 19054359Sroberto register u_int32 lo_tmp; \ 19154359Sroberto register u_int32 hi_tmp; \ 19254359Sroberto \ 19354359Sroberto if ((a_f) == 0) { \ 19454359Sroberto (r_i) -= (a_i); \ 19554359Sroberto } else { \ 19654359Sroberto lo_tmp = ((r_f) & 0xffff) + ((-((s_fp)(a_f))) & 0xffff); \ 19754359Sroberto hi_tmp = (((r_f) >> 16) & 0xffff) \ 19854359Sroberto + (((-((s_fp)(a_f))) >> 16) & 0xffff); \ 19954359Sroberto if (lo_tmp & 0x10000) \ 20054359Sroberto hi_tmp++; \ 20154359Sroberto (r_f) = ((hi_tmp & 0xffff) << 16) | (lo_tmp & 0xffff); \ 20254359Sroberto \ 20354359Sroberto (r_i) += ~(a_i); \ 20454359Sroberto if (hi_tmp & 0x10000) \ 20554359Sroberto (r_i)++; \ 20654359Sroberto } \ 20754359Sroberto } while (0) 20854359Sroberto 20954359Sroberto#define M_RSHIFTU(v_i, v_f) /* v >>= 1, v is unsigned */ \ 21054359Sroberto do { \ 21154359Sroberto (v_f) = (u_int32)(v_f) >> 1; \ 21254359Sroberto if ((v_i) & 01) \ 21354359Sroberto (v_f) |= 0x80000000; \ 21454359Sroberto (v_i) = (u_int32)(v_i) >> 1; \ 21554359Sroberto } while (0) 21654359Sroberto 21754359Sroberto#define M_RSHIFT(v_i, v_f) /* v >>= 1, v is signed */ \ 21854359Sroberto do { \ 21954359Sroberto (v_f) = (u_int32)(v_f) >> 1; \ 22054359Sroberto if ((v_i) & 01) \ 22154359Sroberto (v_f) |= 0x80000000; \ 22254359Sroberto if ((v_i) & 0x80000000) \ 22354359Sroberto (v_i) = ((v_i) >> 1) | 0x80000000; \ 22454359Sroberto else \ 22554359Sroberto (v_i) = (v_i) >> 1; \ 22654359Sroberto } while (0) 22754359Sroberto 22854359Sroberto#define M_LSHIFT(v_i, v_f) /* v <<= 1 */ \ 22954359Sroberto do { \ 23054359Sroberto (v_i) <<= 1; \ 23154359Sroberto if ((v_f) & 0x80000000) \ 23254359Sroberto (v_i) |= 0x1; \ 23354359Sroberto (v_f) <<= 1; \ 23454359Sroberto } while (0) 23554359Sroberto 23654359Sroberto#define M_LSHIFT3(v_ovr, v_i, v_f) /* v <<= 1, with overflow */ \ 23754359Sroberto do { \ 23854359Sroberto (v_ovr) <<= 1; \ 23954359Sroberto if ((v_i) & 0x80000000) \ 24054359Sroberto (v_ovr) |= 0x1; \ 24154359Sroberto (v_i) <<= 1; \ 24254359Sroberto if ((v_f) & 0x80000000) \ 24354359Sroberto (v_i) |= 0x1; \ 24454359Sroberto (v_f) <<= 1; \ 24554359Sroberto } while (0) 24654359Sroberto 24754359Sroberto#define M_ADDUF(r_i, r_f, uf) /* r += uf, uf is u_int32 fraction */ \ 24854359Sroberto M_ADD((r_i), (r_f), 0, (uf)) /* let optimizer worry about it */ 24954359Sroberto 25054359Sroberto#define M_SUBUF(r_i, r_f, uf) /* r -= uf, uf is u_int32 fraction */ \ 25154359Sroberto M_SUB((r_i), (r_f), 0, (uf)) /* let optimizer worry about it */ 25254359Sroberto 25354359Sroberto#define M_ADDF(r_i, r_f, f) /* r += f, f is a int32 fraction */ \ 25454359Sroberto do { \ 25554359Sroberto if ((f) > 0) \ 25654359Sroberto M_ADD((r_i), (r_f), 0, (f)); \ 25754359Sroberto else if ((f) < 0) \ 25854359Sroberto M_ADD((r_i), (r_f), (-1), (f));\ 25954359Sroberto } while(0) 26054359Sroberto 26154359Sroberto#define M_ISNEG(v_i, v_f) /* v < 0 */ \ 26254359Sroberto (((v_i) & 0x80000000) != 0) 26354359Sroberto 26454359Sroberto#define M_ISHIS(a_i, a_f, b_i, b_f) /* a >= b unsigned */ \ 26554359Sroberto (((u_int32)(a_i)) > ((u_int32)(b_i)) || \ 26654359Sroberto ((a_i) == (b_i) && ((u_int32)(a_f)) >= ((u_int32)(b_f)))) 26754359Sroberto 26854359Sroberto#define M_ISGEQ(a_i, a_f, b_i, b_f) /* a >= b signed */ \ 26954359Sroberto (((int32)(a_i)) > ((int32)(b_i)) || \ 27054359Sroberto ((a_i) == (b_i) && ((u_int32)(a_f)) >= ((u_int32)(b_f)))) 27154359Sroberto 27254359Sroberto#define M_ISEQU(a_i, a_f, b_i, b_f) /* a == b unsigned */ \ 27354359Sroberto ((a_i) == (b_i) && (a_f) == (b_f)) 27454359Sroberto 27554359Sroberto/* 27654359Sroberto * Operations on the long fp format 27754359Sroberto */ 27854359Sroberto#define L_ADD(r, a) M_ADD((r)->l_ui, (r)->l_uf, (a)->l_ui, (a)->l_uf) 27954359Sroberto#define L_SUB(r, a) M_SUB((r)->l_ui, (r)->l_uf, (a)->l_ui, (a)->l_uf) 28054359Sroberto#define L_NEG(v) M_NEG((v)->l_ui, (v)->l_uf) 28154359Sroberto#define L_ADDUF(r, uf) M_ADDUF((r)->l_ui, (r)->l_uf, (uf)) 28254359Sroberto#define L_SUBUF(r, uf) M_SUBUF((r)->l_ui, (r)->l_uf, (uf)) 28354359Sroberto#define L_ADDF(r, f) M_ADDF((r)->l_ui, (r)->l_uf, (f)) 28454359Sroberto#define L_RSHIFT(v) M_RSHIFT((v)->l_i, (v)->l_uf) 285182007Sroberto#define L_RSHIFTU(v) M_RSHIFTU((v)->l_ui, (v)->l_uf) 28654359Sroberto#define L_LSHIFT(v) M_LSHIFT((v)->l_ui, (v)->l_uf) 28754359Sroberto#define L_CLR(v) ((v)->l_ui = (v)->l_uf = 0) 28854359Sroberto 28954359Sroberto#define L_ISNEG(v) (((v)->l_ui & 0x80000000) != 0) 29054359Sroberto#define L_ISZERO(v) ((v)->l_ui == 0 && (v)->l_uf == 0) 29154359Sroberto#define L_ISHIS(a, b) ((a)->l_ui > (b)->l_ui || \ 29254359Sroberto ((a)->l_ui == (b)->l_ui && (a)->l_uf >= (b)->l_uf)) 29354359Sroberto#define L_ISGEQ(a, b) ((a)->l_i > (b)->l_i || \ 29454359Sroberto ((a)->l_i == (b)->l_i && (a)->l_uf >= (b)->l_uf)) 29554359Sroberto#define L_ISEQU(a, b) M_ISEQU((a)->l_ui, (a)->l_uf, (b)->l_ui, (b)->l_uf) 29654359Sroberto 29754359Sroberto/* 29854359Sroberto * s_fp/double and u_fp/double conversions 29954359Sroberto */ 30054359Sroberto#define FRIC 65536. /* 2^16 as a double */ 30154359Sroberto#define DTOFP(r) ((s_fp)((r) * FRIC)) 30254359Sroberto#define DTOUFP(r) ((u_fp)((r) * FRIC)) 30354359Sroberto#define FPTOD(r) ((double)(r) / FRIC) 30454359Sroberto 30554359Sroberto/* 30654359Sroberto * l_fp/double conversions 30754359Sroberto */ 30854359Sroberto#define FRAC 4294967296. /* 2^32 as a double */ 30954359Sroberto#define M_DTOLFP(d, r_i, r_uf) /* double to l_fp */ \ 31054359Sroberto do { \ 31154359Sroberto register double d_tmp; \ 31254359Sroberto \ 31354359Sroberto d_tmp = (d); \ 31454359Sroberto if (d_tmp < 0) { \ 31554359Sroberto d_tmp = -d_tmp; \ 31654359Sroberto (r_i) = (int32)(d_tmp); \ 31754359Sroberto (r_uf) = (u_int32)(((d_tmp) - (double)(r_i)) * FRAC); \ 31854359Sroberto M_NEG((r_i), (r_uf)); \ 31954359Sroberto } else { \ 32054359Sroberto (r_i) = (int32)(d_tmp); \ 32154359Sroberto (r_uf) = (u_int32)(((d_tmp) - (double)(r_i)) * FRAC); \ 32254359Sroberto } \ 32354359Sroberto } while (0) 32454359Sroberto#define M_LFPTOD(r_i, r_uf, d) /* l_fp to double */ \ 32554359Sroberto do { \ 32654359Sroberto register l_fp l_tmp; \ 32754359Sroberto \ 32854359Sroberto l_tmp.l_i = (r_i); \ 32954359Sroberto l_tmp.l_f = (r_uf); \ 33054359Sroberto if (l_tmp.l_i < 0) { \ 33154359Sroberto M_NEG(l_tmp.l_i, l_tmp.l_uf); \ 33254359Sroberto (d) = -((double)l_tmp.l_i + ((double)l_tmp.l_uf) / FRAC); \ 33354359Sroberto } else { \ 33454359Sroberto (d) = (double)l_tmp.l_i + ((double)l_tmp.l_uf) / FRAC; \ 33554359Sroberto } \ 33654359Sroberto } while (0) 33754359Sroberto#define DTOLFP(d, v) M_DTOLFP((d), (v)->l_ui, (v)->l_uf) 33854359Sroberto#define LFPTOD(v, d) M_LFPTOD((v)->l_ui, (v)->l_uf, (d)) 33954359Sroberto 34054359Sroberto/* 34154359Sroberto * Prototypes 34254359Sroberto */ 343132451Srobertoextern char * dofptoa P((u_fp, int, short, int)); 344132451Srobertoextern char * dolfptoa P((u_long, u_long, int, short, int)); 34554359Sroberto 34654359Srobertoextern int atolfp P((const char *, l_fp *)); 34754359Srobertoextern int buftvtots P((const char *, l_fp *)); 348132451Srobertoextern char * fptoa P((s_fp, short)); 349132451Srobertoextern char * fptoms P((s_fp, short)); 35054359Srobertoextern int hextolfp P((const char *, l_fp *)); 35154359Srobertoextern void gpstolfp P((int, int, unsigned long, l_fp *)); 35254359Srobertoextern int mstolfp P((const char *, l_fp *)); 35354359Srobertoextern char * prettydate P((l_fp *)); 35454359Srobertoextern char * gmprettydate P((l_fp *)); 35554359Srobertoextern char * uglydate P((l_fp *)); 35654359Srobertoextern void mfp_mul P((int32 *, u_int32 *, int32, u_int32, int32, u_int32)); 35754359Sroberto 35854359Srobertoextern void get_systime P((l_fp *)); 35954359Srobertoextern int step_systime P((double)); 36054359Srobertoextern int adj_systime P((double)); 36154359Sroberto 362182007Srobertoextern struct tm * ntp2unix_tm P((u_long ntp, int local)); 363182007Sroberto 36454359Sroberto#define lfptoa(_fpv, _ndec) mfptoa((_fpv)->l_ui, (_fpv)->l_uf, (_ndec)) 36554359Sroberto#define lfptoms(_fpv, _ndec) mfptoms((_fpv)->l_ui, (_fpv)->l_uf, (_ndec)) 36654359Sroberto 367132451Sroberto#define stoa(_sin) socktoa((_sin)) 368132451Sroberto#define stohost(_sin) socktohost((_sin)) 36954359Sroberto 370132451Sroberto#define ntoa(_sin) stoa(_sin) 371132451Sroberto#define ntohost(_sin) stohost(_sin) 372132451Sroberto 37354359Sroberto#define ufptoa(_fpv, _ndec) dofptoa((_fpv), 0, (_ndec), 0) 37454359Sroberto#define ufptoms(_fpv, _ndec) dofptoa((_fpv), 0, (_ndec), 1) 37554359Sroberto#define ulfptoa(_fpv, _ndec) dolfptoa((_fpv)->l_ui, (_fpv)->l_uf, 0, (_ndec), 0) 37654359Sroberto#define ulfptoms(_fpv, _ndec) dolfptoa((_fpv)->l_ui, (_fpv)->l_uf, 0, (_ndec), 1) 37754359Sroberto#define umfptoa(_fpi, _fpf, _ndec) dolfptoa((_fpi), (_fpf), 0, (_ndec), 0) 37854359Sroberto 37954359Sroberto#endif /* NTP_FP_H */ 380