ntp/ntpd/refclock_wwv.c

122  * sync pulse has been acquired. DSYNC is set when the units digit has
158 #define CMPERR		0x1	/* digit or misc bit compare error */
159 #define LOWERR		0x2	/* low bit or digit amplitude or SNR */
195 #define BTHR		1000.	/* digit threshold */
196 #define BSNR		3.	/* digit likelihood threshold (dB) */
197 #define BCMP		3	/* digit compare threshold */
215 #define TCONST		16	/* data bit/digit time constant */
294 #define DECIM9		4	/* BCD digit 0-9 */
295 #define DECIM6		5	/* BCD digit 0-6 */
296 #define DECIM3		6	/* BCD digit 0-3 */
297 #define DECIM2		7	/* BCD digit 0-2 */
379  * BCD coefficients for maximum-likelihood digit decode
459  * The decoding matrix consists of nine row vectors, one for each digit
467 	int digit;		/* current clock digit */
469 	double digprb;		/* max digit probability */
564 	 * Variables used to estimate signal levels and bit/digit
698 	 * Initialize the decoding matrix with the radix for each digit
1736  * digit bits, decodes miscellaneous bits and dances the leap seconds.
1879 		 * not and the minute digit has been found, restart if
1881 		 * not, restart if the unit digit has not been found
1903 	 * given by the argument. Bits not used in the digit are forced
1923 	 * Correlate coefficient vector with each valid digit vector and
2071 		pp->minute = up->decvec[MN].digit + up->decvec[MN +
2072 		    1].digit * 10;
2073 		pp->hour = up->decvec[HR].digit + up->decvec[HR +
2074 		    1].digit * 10;
2075 		pp->day = up->decvec[DA].digit + up->decvec[DA +
2076 		    1].digit * 10 + up->decvec[DA + 2].digit * 100;
2077 		pp->year = up->decvec[YR].digit + up->decvec[YR +
2078 		    1].digit * 10;
2106  * wwv_corr4 - determine maximum-likelihood digit
2108  * This routine correlates the received digit vector with the BCD
2111  * maximum-likelihood digit, while the ratio of this value to the next
2113  * digit is invalid, the likelihood vector is averaged toward a miss.
2128 	int	mldigit;	/* max likelihood digit */
2135 	 * Correlate digit vector with each BCD coefficient vector. If
2136 	 * any BCD digit bit is bad, consider all bits a miss. Until the
2137 	 * minute units digit has been resolved, don't to anything else.
2160 	 * The current maximum-likelihood digit is compared to the last
2161 	 * maximum-likelihood digit. If different, the compare counter
2162 	 * and maximum-likelihood digit are reset.  When the compare
2163 	 * counter reaches the BCMP threshold (3), the digit is assumed
2167 	 * Note that the clock display digit is set before the compare
2176 		if (vp->digit != mldigit) {
2181 				vp->digit = mldigit;
2196 		    up->epomax, vp->radix, vp->digit, mldigit,
2229 	 * the unit minute digit has been found.
2252 	minute = up->decvec[MN].digit + up->decvec[MN + 1].digit *
2253 	    10 + up->decvec[HR].digit * 60 + up->decvec[HR +
2254 	    1].digit * 600;
2255 	day = up->decvec[DA].digit + up->decvec[DA + 1].digit * 10 +
2256 	    up->decvec[DA + 2].digit * 100;
2261 	isleap = up->decvec[YR].digit & 0x3;
2304  * carry - process digit
2307  * the clock digit modulo the radix. It returns the new clock digit or
2308  * zero if a carry occurred. Once synchronized, the clock digit will
2309  * match the maximum-likelihood digit corresponding to that position.
2319 	dp->digit++;
2320 	if (dp->digit == dp->radix)
2321 		dp->digit = 0;
2326 	return (dp->digit);
2479  *	(DATA) without finding the unit seconds digit.
2632 	year = up->decvec[YR].digit + up->decvec[YR + 1].digit * 10 +
2634 	day = up->decvec[DA].digit + up->decvec[DA + 1].digit * 10 +
2635 	    up->decvec[DA + 2].digit * 100;
2636 	hour = up->decvec[HR].digit + up->decvec[HR + 1].digit * 10;
2637 	minute = up->decvec[MN].digit + up->decvec[MN + 1].digit * 10;