xref: /freebsd/contrib/ntp/libparse/clk_rawdcf.c (revision f5f40dd63bc7acbb5312b26ac1ea1103c12352a6)
1c0b746e5SOllivier Robert /*
2ea906c41SOllivier Robert  * /src/NTP/REPOSITORY/ntp4-dev/libparse/clk_rawdcf.c,v 4.18 2006/06/22 18:40:01 kardel RELEASE_20060622_A
3c0b746e5SOllivier Robert  *
4ea906c41SOllivier Robert  * clk_rawdcf.c,v 4.18 2006/06/22 18:40:01 kardel RELEASE_20060622_A
5c0b746e5SOllivier Robert  *
6c0b746e5SOllivier Robert  * Raw DCF77 pulse clock support
7c0b746e5SOllivier Robert  *
8276da39aSCy Schubert  * Copyright (c) 1995-2015 by Frank Kardel <kardel <AT> ntp.org>
9a25439b6SCy Schubert  * Copyright (c) 1989-1994 by Frank Kardel, Friedrich-Alexander Universitaet Erlangen-Nuernberg, Germany
10c0b746e5SOllivier Robert  *
11ea906c41SOllivier Robert  * Redistribution and use in source and binary forms, with or without
12ea906c41SOllivier Robert  * modification, are permitted provided that the following conditions
13ea906c41SOllivier Robert  * are met:
14ea906c41SOllivier Robert  * 1. Redistributions of source code must retain the above copyright
15ea906c41SOllivier Robert  *    notice, this list of conditions and the following disclaimer.
16ea906c41SOllivier Robert  * 2. Redistributions in binary form must reproduce the above copyright
17ea906c41SOllivier Robert  *    notice, this list of conditions and the following disclaimer in the
18ea906c41SOllivier Robert  *    documentation and/or other materials provided with the distribution.
19ea906c41SOllivier Robert  * 3. Neither the name of the author nor the names of its contributors
20ea906c41SOllivier Robert  *    may be used to endorse or promote products derived from this software
21ea906c41SOllivier Robert  *    without specific prior written permission.
22ea906c41SOllivier Robert  *
23ea906c41SOllivier Robert  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
24ea906c41SOllivier Robert  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
25ea906c41SOllivier Robert  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
26ea906c41SOllivier Robert  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
27ea906c41SOllivier Robert  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
28ea906c41SOllivier Robert  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
29ea906c41SOllivier Robert  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
30ea906c41SOllivier Robert  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
31ea906c41SOllivier Robert  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
32ea906c41SOllivier Robert  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
33ea906c41SOllivier Robert  * SUCH DAMAGE.
34c0b746e5SOllivier Robert  *
35c0b746e5SOllivier Robert  */
36c0b746e5SOllivier Robert 
37c0b746e5SOllivier Robert #ifdef HAVE_CONFIG_H
38c0b746e5SOllivier Robert # include <config.h>
39c0b746e5SOllivier Robert #endif
40c0b746e5SOllivier Robert 
41c0b746e5SOllivier Robert #if defined(REFCLOCK) && defined(CLOCK_PARSE) && defined(CLOCK_RAWDCF)
42c0b746e5SOllivier Robert 
43c0b746e5SOllivier Robert #include "ntp_fp.h"
44276da39aSCy Schubert #include "timevalops.h"
45c0b746e5SOllivier Robert #include "ntp_unixtime.h"
46c0b746e5SOllivier Robert #include "ntp_calendar.h"
47c0b746e5SOllivier Robert 
48c0b746e5SOllivier Robert #include "parse.h"
49c0b746e5SOllivier Robert #ifdef PARSESTREAM
50c0b746e5SOllivier Robert # include <sys/parsestreams.h>
51c0b746e5SOllivier Robert #endif
52c0b746e5SOllivier Robert 
53c0b746e5SOllivier Robert #ifndef PARSEKERNEL
54c0b746e5SOllivier Robert # include "ntp_stdlib.h"
55c0b746e5SOllivier Robert #endif
56c0b746e5SOllivier Robert 
57c0b746e5SOllivier Robert /*
58c0b746e5SOllivier Robert  * DCF77 raw time code
59c0b746e5SOllivier Robert  *
60c0b746e5SOllivier Robert  * From "Zur Zeit", Physikalisch-Technische Bundesanstalt (PTB), Braunschweig
61c0b746e5SOllivier Robert  * und Berlin, Maerz 1989
62c0b746e5SOllivier Robert  *
63c0b746e5SOllivier Robert  * Timecode transmission:
64c0b746e5SOllivier Robert  * AM:
65c0b746e5SOllivier Robert  *	time marks are send every second except for the second before the
66c0b746e5SOllivier Robert  *	next minute mark
67c0b746e5SOllivier Robert  *	time marks consist of a reduction of transmitter power to 25%
68c0b746e5SOllivier Robert  *	of the nominal level
69c0b746e5SOllivier Robert  *	the falling edge is the time indication (on time)
70c0b746e5SOllivier Robert  *	time marks of a 100ms duration constitute a logical 0
71c0b746e5SOllivier Robert  *	time marks of a 200ms duration constitute a logical 1
72c0b746e5SOllivier Robert  * FM:
73c0b746e5SOllivier Robert  *	see the spec. (basically a (non-)inverted psuedo random phase shift)
74c0b746e5SOllivier Robert  *
75c0b746e5SOllivier Robert  * Encoding:
76c0b746e5SOllivier Robert  * Second	Contents
77c0b746e5SOllivier Robert  * 0  - 10	AM: free, FM: 0
78c0b746e5SOllivier Robert  * 11 - 14	free
79a25439b6SCy Schubert  * 15		R     - "call bit" used to signalize irregularities in the control facilities
80a25439b6SCy Schubert  *		        (until 2003 indicated transmission via alternate antenna)
81c0b746e5SOllivier Robert  * 16		A1    - expect zone change (1 hour before)
82c0b746e5SOllivier Robert  * 17 - 18	Z1,Z2 - time zone
83c0b746e5SOllivier Robert  *		 0  0 illegal
84c0b746e5SOllivier Robert  *		 0  1 MEZ  (MET)
85c0b746e5SOllivier Robert  *		 1  0 MESZ (MED, MET DST)
86c0b746e5SOllivier Robert  *		 1  1 illegal
87c0b746e5SOllivier Robert  * 19		A2    - expect leap insertion/deletion (1 hour before)
88c0b746e5SOllivier Robert  * 20		S     - start of time code (1)
89c0b746e5SOllivier Robert  * 21 - 24	M1    - BCD (lsb first) Minutes
90c0b746e5SOllivier Robert  * 25 - 27	M10   - BCD (lsb first) 10 Minutes
91c0b746e5SOllivier Robert  * 28		P1    - Minute Parity (even)
92c0b746e5SOllivier Robert  * 29 - 32	H1    - BCD (lsb first) Hours
93c0b746e5SOllivier Robert  * 33 - 34      H10   - BCD (lsb first) 10 Hours
94c0b746e5SOllivier Robert  * 35		P2    - Hour Parity (even)
95c0b746e5SOllivier Robert  * 36 - 39	D1    - BCD (lsb first) Days
96c0b746e5SOllivier Robert  * 40 - 41	D10   - BCD (lsb first) 10 Days
97c0b746e5SOllivier Robert  * 42 - 44	DW    - BCD (lsb first) day of week (1: Monday -> 7: Sunday)
98c0b746e5SOllivier Robert  * 45 - 49	MO    - BCD (lsb first) Month
99c0b746e5SOllivier Robert  * 50           MO0   - 10 Months
100c0b746e5SOllivier Robert  * 51 - 53	Y1    - BCD (lsb first) Years
101c0b746e5SOllivier Robert  * 54 - 57	Y10   - BCD (lsb first) 10 Years
102c0b746e5SOllivier Robert  * 58 		P3    - Date Parity (even)
103c0b746e5SOllivier Robert  * 59		      - usually missing (minute indication), except for leap insertion
104c0b746e5SOllivier Robert  */
105c0b746e5SOllivier Robert 
106a25439b6SCy Schubert static parse_pps_fnc_t pps_rawdcf;
107a25439b6SCy Schubert static parse_cvt_fnc_t cvt_rawdcf;
108a25439b6SCy Schubert static parse_inp_fnc_t inp_rawdcf;
109c0b746e5SOllivier Robert 
110a151a66cSOllivier Robert typedef struct last_tcode {
111a151a66cSOllivier Robert 	time_t      tcode;	/* last converted time code */
112276da39aSCy Schubert         timestamp_t tminute;	/* sample time for minute start */
113276da39aSCy Schubert         timestamp_t timeout;	/* last timeout timestamp */
114a151a66cSOllivier Robert } last_tcode_t;
115a151a66cSOllivier Robert 
116ea906c41SOllivier Robert #define BUFFER_MAX	61
117ea906c41SOllivier Robert 
118c0b746e5SOllivier Robert clockformat_t clock_rawdcf =
119c0b746e5SOllivier Robert {
120c0b746e5SOllivier Robert   inp_rawdcf,			/* DCF77 input handling */
121c0b746e5SOllivier Robert   cvt_rawdcf,			/* raw dcf input conversion */
122c0b746e5SOllivier Robert   pps_rawdcf,			/* examining PPS information */
123c0b746e5SOllivier Robert   0,				/* no private configuration data */
124c0b746e5SOllivier Robert   "RAW DCF77 Timecode",		/* direct decoding / time synthesis */
125c0b746e5SOllivier Robert 
126ea906c41SOllivier Robert   BUFFER_MAX,			/* bit buffer */
127a151a66cSOllivier Robert   sizeof(last_tcode_t)
128c0b746e5SOllivier Robert };
129c0b746e5SOllivier Robert 
130c0b746e5SOllivier Robert static struct dcfparam
131c0b746e5SOllivier Robert {
1322b15cb3dSCy Schubert 	const unsigned char *onebits;
1332b15cb3dSCy Schubert 	const unsigned char *zerobits;
134c0b746e5SOllivier Robert } dcfparameter =
135c0b746e5SOllivier Robert {
1362b15cb3dSCy Schubert 	(const unsigned char *)"###############RADMLS1248124P124812P1248121241248112481248P??", /* 'ONE' representation */
1372b15cb3dSCy Schubert 	(const unsigned char *)"--------------------s-------p------p----------------------p__"  /* 'ZERO' representation */
138c0b746e5SOllivier Robert };
139c0b746e5SOllivier Robert 
140c0b746e5SOllivier Robert static struct rawdcfcode
141c0b746e5SOllivier Robert {
142c0b746e5SOllivier Robert 	char offset;			/* start bit */
143c0b746e5SOllivier Robert } rawdcfcode[] =
144c0b746e5SOllivier Robert {
145c0b746e5SOllivier Robert 	{  0 }, { 15 }, { 16 }, { 17 }, { 19 }, { 20 }, { 21 }, { 25 }, { 28 }, { 29 },
146c0b746e5SOllivier Robert 	{ 33 }, { 35 }, { 36 }, { 40 }, { 42 }, { 45 }, { 49 }, { 50 }, { 54 }, { 58 }, { 59 }
147c0b746e5SOllivier Robert };
148c0b746e5SOllivier Robert 
149c0b746e5SOllivier Robert #define DCF_M	0
150c0b746e5SOllivier Robert #define DCF_R	1
151c0b746e5SOllivier Robert #define DCF_A1	2
152c0b746e5SOllivier Robert #define DCF_Z	3
153c0b746e5SOllivier Robert #define DCF_A2	4
154c0b746e5SOllivier Robert #define DCF_S	5
155c0b746e5SOllivier Robert #define DCF_M1	6
156c0b746e5SOllivier Robert #define DCF_M10	7
157c0b746e5SOllivier Robert #define DCF_P1	8
158c0b746e5SOllivier Robert #define DCF_H1	9
159c0b746e5SOllivier Robert #define DCF_H10	10
160c0b746e5SOllivier Robert #define DCF_P2	11
161c0b746e5SOllivier Robert #define DCF_D1	12
162c0b746e5SOllivier Robert #define DCF_D10	13
163c0b746e5SOllivier Robert #define DCF_DW	14
164c0b746e5SOllivier Robert #define DCF_MO	15
165c0b746e5SOllivier Robert #define DCF_MO0	16
166c0b746e5SOllivier Robert #define DCF_Y1	17
167c0b746e5SOllivier Robert #define DCF_Y10	18
168c0b746e5SOllivier Robert #define DCF_P3	19
169c0b746e5SOllivier Robert 
170c0b746e5SOllivier Robert static struct partab
171c0b746e5SOllivier Robert {
172c0b746e5SOllivier Robert 	char offset;			/* start bit of parity field */
173c0b746e5SOllivier Robert } partab[] =
174c0b746e5SOllivier Robert {
175c0b746e5SOllivier Robert 	{ 21 }, { 29 }, { 36 }, { 59 }
176c0b746e5SOllivier Robert };
177c0b746e5SOllivier Robert 
178c0b746e5SOllivier Robert #define DCF_P_P1	0
179c0b746e5SOllivier Robert #define DCF_P_P2	1
180c0b746e5SOllivier Robert #define DCF_P_P3	2
181c0b746e5SOllivier Robert 
182c0b746e5SOllivier Robert #define DCF_Z_MET 0x2
183c0b746e5SOllivier Robert #define DCF_Z_MED 0x1
184c0b746e5SOllivier Robert 
185c0b746e5SOllivier Robert static u_long
186c0b746e5SOllivier Robert ext_bf(
187ea906c41SOllivier Robert 	unsigned char *buf,
188ea906c41SOllivier Robert 	int   idx,
1892b15cb3dSCy Schubert 	const unsigned char *zero
190c0b746e5SOllivier Robert 	)
191c0b746e5SOllivier Robert {
192ea906c41SOllivier Robert 	u_long sum = 0;
193ea906c41SOllivier Robert 	int i, first;
194c0b746e5SOllivier Robert 
195c0b746e5SOllivier Robert 	first = rawdcfcode[idx].offset;
196c0b746e5SOllivier Robert 
197c0b746e5SOllivier Robert 	for (i = rawdcfcode[idx+1].offset - 1; i >= first; i--)
198c0b746e5SOllivier Robert 	{
199c0b746e5SOllivier Robert 		sum <<= 1;
200c0b746e5SOllivier Robert 		sum |= (buf[i] != zero[i]);
201c0b746e5SOllivier Robert 	}
202c0b746e5SOllivier Robert 	return sum;
203c0b746e5SOllivier Robert }
204c0b746e5SOllivier Robert 
205c0b746e5SOllivier Robert static unsigned
206c0b746e5SOllivier Robert pcheck(
207c0b746e5SOllivier Robert        unsigned char *buf,
208c0b746e5SOllivier Robert        int   idx,
2092b15cb3dSCy Schubert        const unsigned char *zero
210c0b746e5SOllivier Robert        )
211c0b746e5SOllivier Robert {
212c0b746e5SOllivier Robert 	int i,last;
213c0b746e5SOllivier Robert 	unsigned psum = 1;
214c0b746e5SOllivier Robert 
215c0b746e5SOllivier Robert 	last = partab[idx+1].offset;
216c0b746e5SOllivier Robert 
217c0b746e5SOllivier Robert 	for (i = partab[idx].offset; i < last; i++)
218c0b746e5SOllivier Robert 	    psum ^= (buf[i] != zero[i]);
219c0b746e5SOllivier Robert 
220c0b746e5SOllivier Robert 	return psum;
221c0b746e5SOllivier Robert }
222c0b746e5SOllivier Robert 
2232d4e511cSCy Schubert static int/*BOOL*/
2242d4e511cSCy Schubert zeller_expand(
2252d4e511cSCy Schubert 	clocktime_t     *clock_time,
2262d4e511cSCy Schubert 	unsigned int	wd
2272d4e511cSCy Schubert 	)
2282d4e511cSCy Schubert {
2292d4e511cSCy Schubert         unsigned int  y = (unsigned int)clock_time->year;
2302d4e511cSCy Schubert         unsigned int  m = (unsigned int)clock_time->month - 1u;
2312d4e511cSCy Schubert         unsigned int  d = (unsigned int)clock_time->day - 1u;
2322d4e511cSCy Schubert 	unsigned int  c;
2332d4e511cSCy Schubert 
2342d4e511cSCy Schubert 	/* Check basic constraints first. */
2352d4e511cSCy Schubert         if ((y >= 100u) || (m >= 12u) || (d >= 31u) || (--wd >= 7u))
2362d4e511cSCy Schubert 		return FALSE;
2372d4e511cSCy Schubert 
2382d4e511cSCy Schubert 	/* Get weekday of date in 1st century by a variation on Zeller's
2392d4e511cSCy Schubert 	 * congruence. All operands are non-negative, and the month
2402d4e511cSCy Schubert 	 * formula is adjusted to use a divider of 32, so we can do a
2412d4e511cSCy Schubert 	 * shift instead of a 'true' division:
2422d4e511cSCy Schubert 	 */
2432d4e511cSCy Schubert 	if ((m += 10u) >= 12u)		/* shift base to 0000-03-01 */
2442d4e511cSCy Schubert 		m -= 12u;
2452d4e511cSCy Schubert 	else if (--y >= 100u)
2462d4e511cSCy Schubert 		y += 100;
2472d4e511cSCy Schubert 	d += y + (y >> 2) + 2u;		/* year-related share */
2482d4e511cSCy Schubert 	d += (m * 83u + 16u) >> 5;	/* month-related share */
2492d4e511cSCy Schubert 
2502d4e511cSCy Schubert 	/* The next step combines the exact division by modular inverse
2512d4e511cSCy Schubert 	 * with the (mod 7) step in such way that no true division and
2522d4e511cSCy Schubert 	 * only one multiplication is needed. The multiplier is
2532d4e511cSCy Schubert 	 *      M <- ceil((3*8)/7 * 2**29)
2542d4e511cSCy Schubert 	 * and combines multiplication by invmod(5, 7) -> 3 and modulus
2552d4e511cSCy Schubert 	 * by 7 transformation to (mod 8) in one step.
2562d4e511cSCy Schubert 	 *   Note that 252 == 0 (mod 7) and that 'd' is less than 185,
2572d4e511cSCy Schubert 	 * so the number to invert and reduce is strictly positive. In
2582d4e511cSCy Schubert 	 * the end, 'c' is number of centuries since start of a great
2592d4e511cSCy Schubert 	 * cycle and must be in [0..3] or we had bad input.
2602d4e511cSCy Schubert 	 */
2612d4e511cSCy Schubert 	c = (((252u + wd - d) * 0x6db6db6eU) >> 29) & 7u;
2622d4e511cSCy Schubert 	if (c >= 4)
2632d4e511cSCy Schubert 		return FALSE;
2642d4e511cSCy Schubert 	/* undo calendar base shift now */
2652d4e511cSCy Schubert 	if ((m > 9u) && (++y >= 100u)) {
2662d4e511cSCy Schubert 		y -= 100u;
2672d4e511cSCy Schubert 		c = (c + 1u) & 3u;
2682d4e511cSCy Schubert 	}
2692d4e511cSCy Schubert 	/* combine year with centuries & map to [1970..2369] */
2702d4e511cSCy Schubert 	y += (c * 100u);
2712d4e511cSCy Schubert 	clock_time->year = (int)y + ((y < 370u) ? 2000 : 1600);
2722d4e511cSCy Schubert 	return TRUE;
2732d4e511cSCy Schubert }
2742d4e511cSCy Schubert 
275c0b746e5SOllivier Robert static u_long
276c0b746e5SOllivier Robert convert_rawdcf(
277c0b746e5SOllivier Robert 	       unsigned char   *buffer,
278c0b746e5SOllivier Robert 	       int              size,
279c0b746e5SOllivier Robert 	       struct dcfparam *dcfprm,
280c0b746e5SOllivier Robert 	       clocktime_t     *clock_time
281c0b746e5SOllivier Robert 	       )
282c0b746e5SOllivier Robert {
283ea906c41SOllivier Robert 	unsigned char *s = buffer;
2842b15cb3dSCy Schubert 	const unsigned char *b = dcfprm->onebits;
2852b15cb3dSCy Schubert 	const unsigned char *c = dcfprm->zerobits;
286ea906c41SOllivier Robert 	int i;
287c0b746e5SOllivier Robert 
288276da39aSCy Schubert 	parseprintf(DD_RAWDCF,("parse: convert_rawdcf: \"%.*s\"\n", size, buffer));
289c0b746e5SOllivier Robert 
290c0b746e5SOllivier Robert 	if (size < 57)
291c0b746e5SOllivier Robert 	{
292c0b746e5SOllivier Robert #ifndef PARSEKERNEL
2932b15cb3dSCy Schubert 		msyslog(LOG_ERR, "parse: convert_rawdcf: INCOMPLETE DATA - time code only has %d bits", size);
294c0b746e5SOllivier Robert #endif
295276da39aSCy Schubert 		return CVT_FAIL|CVT_BADFMT;
296c0b746e5SOllivier Robert 	}
297c0b746e5SOllivier Robert 
298ea906c41SOllivier Robert 	for (i = 0; i < size; i++)
299c0b746e5SOllivier Robert 	{
300c0b746e5SOllivier Robert 		if ((*s != *b) && (*s != *c))
301c0b746e5SOllivier Robert 		{
302c0b746e5SOllivier Robert 			/*
303c0b746e5SOllivier Robert 			 * we only have two types of bytes (ones and zeros)
304c0b746e5SOllivier Robert 			 */
305c0b746e5SOllivier Robert #ifndef PARSEKERNEL
306ea906c41SOllivier Robert 			msyslog(LOG_ERR, "parse: convert_rawdcf: BAD DATA - no conversion");
307c0b746e5SOllivier Robert #endif
308276da39aSCy Schubert 			return CVT_FAIL|CVT_BADFMT;
309c0b746e5SOllivier Robert 		}
310ea906c41SOllivier Robert 		if (*b) b++;
311ea906c41SOllivier Robert 		if (*c) c++;
312c0b746e5SOllivier Robert 		s++;
313c0b746e5SOllivier Robert 	}
314c0b746e5SOllivier Robert 
315c0b746e5SOllivier Robert 	/*
316c0b746e5SOllivier Robert 	 * check Start and Parity bits
317c0b746e5SOllivier Robert 	 */
318c0b746e5SOllivier Robert 	if ((ext_bf(buffer, DCF_S, dcfprm->zerobits) == 1) &&
319c0b746e5SOllivier Robert 	    pcheck(buffer, DCF_P_P1, dcfprm->zerobits) &&
320c0b746e5SOllivier Robert 	    pcheck(buffer, DCF_P_P2, dcfprm->zerobits) &&
321c0b746e5SOllivier Robert 	    pcheck(buffer, DCF_P_P3, dcfprm->zerobits))
322c0b746e5SOllivier Robert 	{
323c0b746e5SOllivier Robert 		/*
324c0b746e5SOllivier Robert 		 * buffer OK
325c0b746e5SOllivier Robert 		 */
326c0b746e5SOllivier Robert 		parseprintf(DD_RAWDCF,("parse: convert_rawdcf: parity check passed\n"));
327c0b746e5SOllivier Robert 
328276da39aSCy Schubert 		clock_time->flags  = PARSEB_S_CALLBIT|PARSEB_S_LEAP;
329c0b746e5SOllivier Robert 		clock_time->utctime= 0;
330c0b746e5SOllivier Robert 		clock_time->usecond= 0;
331c0b746e5SOllivier Robert 		clock_time->second = 0;
332c0b746e5SOllivier Robert 		clock_time->minute = ext_bf(buffer, DCF_M10, dcfprm->zerobits);
333c0b746e5SOllivier Robert 		clock_time->minute = TIMES10(clock_time->minute) + ext_bf(buffer, DCF_M1, dcfprm->zerobits);
334c0b746e5SOllivier Robert 		clock_time->hour   = ext_bf(buffer, DCF_H10, dcfprm->zerobits);
335c0b746e5SOllivier Robert 		clock_time->hour   = TIMES10(clock_time->hour) + ext_bf(buffer, DCF_H1, dcfprm->zerobits);
336c0b746e5SOllivier Robert 		clock_time->day    = ext_bf(buffer, DCF_D10, dcfprm->zerobits);
337c0b746e5SOllivier Robert 		clock_time->day    = TIMES10(clock_time->day) + ext_bf(buffer, DCF_D1, dcfprm->zerobits);
338c0b746e5SOllivier Robert 		clock_time->month  = ext_bf(buffer, DCF_MO0, dcfprm->zerobits);
339c0b746e5SOllivier Robert 		clock_time->month  = TIMES10(clock_time->month) + ext_bf(buffer, DCF_MO, dcfprm->zerobits);
340c0b746e5SOllivier Robert 		clock_time->year   = ext_bf(buffer, DCF_Y10, dcfprm->zerobits);
341c0b746e5SOllivier Robert 		clock_time->year   = TIMES10(clock_time->year) + ext_bf(buffer, DCF_Y1, dcfprm->zerobits);
342c0b746e5SOllivier Robert 
3432d4e511cSCy Schubert 		if (!zeller_expand(clock_time, ext_bf(buffer, DCF_DW, dcfprm->zerobits)))
3442d4e511cSCy Schubert 		    return CVT_FAIL|CVT_BADFMT;
3452d4e511cSCy Schubert 
346c0b746e5SOllivier Robert 		switch (ext_bf(buffer, DCF_Z, dcfprm->zerobits))
347c0b746e5SOllivier Robert 		{
348c0b746e5SOllivier Robert 		    case DCF_Z_MET:
349c0b746e5SOllivier Robert 			clock_time->utcoffset = -1*60*60;
350c0b746e5SOllivier Robert 			break;
351c0b746e5SOllivier Robert 
352c0b746e5SOllivier Robert 		    case DCF_Z_MED:
353c0b746e5SOllivier Robert 			clock_time->flags     |= PARSEB_DST;
354c0b746e5SOllivier Robert 			clock_time->utcoffset  = -2*60*60;
355c0b746e5SOllivier Robert 			break;
356c0b746e5SOllivier Robert 
357c0b746e5SOllivier Robert 		    default:
358c0b746e5SOllivier Robert 			parseprintf(DD_RAWDCF,("parse: convert_rawdcf: BAD TIME ZONE\n"));
359c0b746e5SOllivier Robert 			return CVT_FAIL|CVT_BADFMT;
360c0b746e5SOllivier Robert 		}
361c0b746e5SOllivier Robert 
362c0b746e5SOllivier Robert 		if (ext_bf(buffer, DCF_A1, dcfprm->zerobits))
363c0b746e5SOllivier Robert 		    clock_time->flags |= PARSEB_ANNOUNCE;
364c0b746e5SOllivier Robert 
365c0b746e5SOllivier Robert 		if (ext_bf(buffer, DCF_A2, dcfprm->zerobits))
366c0b746e5SOllivier Robert 		    clock_time->flags |= PARSEB_LEAPADD; /* default: DCF77 data format deficiency */
367c0b746e5SOllivier Robert 
368c0b746e5SOllivier Robert 		if (ext_bf(buffer, DCF_R, dcfprm->zerobits))
369a25439b6SCy Schubert 		    clock_time->flags |= PARSEB_CALLBIT;
370c0b746e5SOllivier Robert 
371276da39aSCy Schubert 		parseprintf(DD_RAWDCF,("parse: convert_rawdcf: TIME CODE OK: %02d:%02d, %02d.%02d.%02d, flags 0x%lx\n",
372c0b746e5SOllivier Robert 				       (int)clock_time->hour, (int)clock_time->minute, (int)clock_time->day, (int)clock_time->month,(int) clock_time->year,
373c0b746e5SOllivier Robert 				       (u_long)clock_time->flags));
374c0b746e5SOllivier Robert 		return CVT_OK;
375c0b746e5SOllivier Robert 	}
376c0b746e5SOllivier Robert 	else
377c0b746e5SOllivier Robert 	{
378c0b746e5SOllivier Robert 		/*
379c0b746e5SOllivier Robert 		 * bad format - not for us
380c0b746e5SOllivier Robert 		 */
381c0b746e5SOllivier Robert #ifndef PARSEKERNEL
382276da39aSCy Schubert 		msyslog(LOG_ERR, "parse: convert_rawdcf: start bit / parity check FAILED for \"%.*s\"", size, buffer);
383c0b746e5SOllivier Robert #endif
384c0b746e5SOllivier Robert 		return CVT_FAIL|CVT_BADFMT;
385c0b746e5SOllivier Robert 	}
386c0b746e5SOllivier Robert }
387c0b746e5SOllivier Robert 
388c0b746e5SOllivier Robert /*
389a25439b6SCy Schubert  * parse_cvt_fnc_t cvt_rawdcf
390c0b746e5SOllivier Robert  * raw dcf input routine - needs to fix up 50 baud
391c0b746e5SOllivier Robert  * characters for 1/0 decision
392c0b746e5SOllivier Robert  */
393c0b746e5SOllivier Robert static u_long
394c0b746e5SOllivier Robert cvt_rawdcf(
395c0b746e5SOllivier Robert 	   unsigned char   *buffer,
396c0b746e5SOllivier Robert 	   int              size,
397c0b746e5SOllivier Robert 	   struct format   *param,
398c0b746e5SOllivier Robert 	   clocktime_t     *clock_time,
399c0b746e5SOllivier Robert 	   void            *local
400c0b746e5SOllivier Robert 	   )
401c0b746e5SOllivier Robert {
402a151a66cSOllivier Robert 	last_tcode_t  *t = (last_tcode_t *)local;
403ea906c41SOllivier Robert 	unsigned char *s = (unsigned char *)buffer;
404ea906c41SOllivier Robert 	unsigned char *e = s + size;
4052b15cb3dSCy Schubert 	const unsigned char *b = dcfparameter.onebits;
4062b15cb3dSCy Schubert 	const unsigned char *c = dcfparameter.zerobits;
407a151a66cSOllivier Robert 	u_long       rtc = CVT_NONE;
408ea906c41SOllivier Robert 	unsigned int i, lowmax, highmax, cutoff, span;
409c0b746e5SOllivier Robert #define BITS 9
410c0b746e5SOllivier Robert 	unsigned char     histbuf[BITS];
411c0b746e5SOllivier Robert 	/*
412c0b746e5SOllivier Robert 	 * the input buffer contains characters with runs of consecutive
413c0b746e5SOllivier Robert 	 * bits set. These set bits are an indication of the DCF77 pulse
414c0b746e5SOllivier Robert 	 * length. We assume that we receive the pulse at 50 Baud. Thus
415c0b746e5SOllivier Robert 	 * a 100ms pulse would generate a 4 bit train (20ms per bit and
416c0b746e5SOllivier Robert 	 * start bit)
417c0b746e5SOllivier Robert 	 * a 200ms pulse would create all zeroes (and probably a frame error)
418c0b746e5SOllivier Robert 	 */
419c0b746e5SOllivier Robert 
420c0b746e5SOllivier Robert 	for (i = 0; i < BITS; i++)
421c0b746e5SOllivier Robert 	{
422c0b746e5SOllivier Robert 		histbuf[i] = 0;
423c0b746e5SOllivier Robert 	}
424c0b746e5SOllivier Robert 
425c0b746e5SOllivier Robert 	cutoff = 0;
426c0b746e5SOllivier Robert 	lowmax = 0;
427c0b746e5SOllivier Robert 
428c0b746e5SOllivier Robert 	while (s < e)
429c0b746e5SOllivier Robert 	{
430ea906c41SOllivier Robert 		unsigned int ch = *s ^ 0xFF;
431c0b746e5SOllivier Robert 		/*
432c0b746e5SOllivier Robert 		 * these lines are left as an excercise to the reader 8-)
433c0b746e5SOllivier Robert 		 */
434c0b746e5SOllivier Robert 		if (!((ch+1) & ch) || !*s)
435c0b746e5SOllivier Robert 		{
436c0b746e5SOllivier Robert 
437c0b746e5SOllivier Robert 			for (i = 0; ch; i++)
438c0b746e5SOllivier Robert 			{
439c0b746e5SOllivier Robert 				ch >>= 1;
440c0b746e5SOllivier Robert 			}
441c0b746e5SOllivier Robert 
442a25439b6SCy Schubert 			*s = (unsigned char) i;
443c0b746e5SOllivier Robert 			histbuf[i]++;
444c0b746e5SOllivier Robert 			cutoff += i;
445c0b746e5SOllivier Robert 			lowmax++;
446c0b746e5SOllivier Robert 		}
447c0b746e5SOllivier Robert 		else
448c0b746e5SOllivier Robert 		{
449c0b746e5SOllivier Robert 			parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: character check for 0x%x@%d FAILED\n", *s, (int)(s - (unsigned char *)buffer)));
450c0b746e5SOllivier Robert 			*s = (unsigned char)~0;
451c0b746e5SOllivier Robert 			rtc = CVT_FAIL|CVT_BADFMT;
452c0b746e5SOllivier Robert 		}
453c0b746e5SOllivier Robert 		s++;
454c0b746e5SOllivier Robert 	}
455c0b746e5SOllivier Robert 
456c0b746e5SOllivier Robert 	if (lowmax)
457c0b746e5SOllivier Robert 	{
458c0b746e5SOllivier Robert 		cutoff /= lowmax;
459c0b746e5SOllivier Robert 	}
460c0b746e5SOllivier Robert 	else
461c0b746e5SOllivier Robert 	{
462c0b746e5SOllivier Robert 		cutoff = 4;	/* doesn't really matter - it'll fail anyway, but gives error output */
463c0b746e5SOllivier Robert 	}
464c0b746e5SOllivier Robert 
465c0b746e5SOllivier Robert 	parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: average bit count: %d\n", cutoff));
466c0b746e5SOllivier Robert 
467c0b746e5SOllivier Robert 	lowmax = 0;
468c0b746e5SOllivier Robert 	highmax = 0;
469c0b746e5SOllivier Robert 
470c0b746e5SOllivier Robert 	parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: histogram:"));
471c0b746e5SOllivier Robert 	for (i = 0; i <= cutoff; i++)
472c0b746e5SOllivier Robert 	{
473c0b746e5SOllivier Robert 		lowmax+=histbuf[i] * i;
474c0b746e5SOllivier Robert 		highmax += histbuf[i];
475c0b746e5SOllivier Robert 		parseprintf(DD_RAWDCF,(" %d", histbuf[i]));
476c0b746e5SOllivier Robert 	}
477c0b746e5SOllivier Robert 	parseprintf(DD_RAWDCF, (" <M>"));
478c0b746e5SOllivier Robert 
479c0b746e5SOllivier Robert 	lowmax += highmax / 2;
480c0b746e5SOllivier Robert 
481c0b746e5SOllivier Robert 	if (highmax)
482c0b746e5SOllivier Robert 	{
483c0b746e5SOllivier Robert 		lowmax /= highmax;
484c0b746e5SOllivier Robert 	}
485c0b746e5SOllivier Robert 	else
486c0b746e5SOllivier Robert 	{
487c0b746e5SOllivier Robert 		lowmax = 0;
488c0b746e5SOllivier Robert 	}
489c0b746e5SOllivier Robert 
490c0b746e5SOllivier Robert 	highmax = 0;
491c0b746e5SOllivier Robert 	cutoff = 0;
492c0b746e5SOllivier Robert 
493c0b746e5SOllivier Robert 	for (; i < BITS; i++)
494c0b746e5SOllivier Robert 	{
495c0b746e5SOllivier Robert 		highmax+=histbuf[i] * i;
496c0b746e5SOllivier Robert 		cutoff +=histbuf[i];
497c0b746e5SOllivier Robert 		parseprintf(DD_RAWDCF,(" %d", histbuf[i]));
498c0b746e5SOllivier Robert 	}
499c0b746e5SOllivier Robert 	parseprintf(DD_RAWDCF,("\n"));
500c0b746e5SOllivier Robert 
501c0b746e5SOllivier Robert 	if (cutoff)
502c0b746e5SOllivier Robert 	{
503c0b746e5SOllivier Robert 		highmax /= cutoff;
504c0b746e5SOllivier Robert 	}
505c0b746e5SOllivier Robert 	else
506c0b746e5SOllivier Robert 	{
507c0b746e5SOllivier Robert 		highmax = BITS-1;
508c0b746e5SOllivier Robert 	}
509c0b746e5SOllivier Robert 
510c0b746e5SOllivier Robert 	span = cutoff = lowmax;
511c0b746e5SOllivier Robert 	for (i = lowmax; i <= highmax; i++)
512c0b746e5SOllivier Robert 	{
513c0b746e5SOllivier Robert 		if (histbuf[cutoff] > histbuf[i])
514c0b746e5SOllivier Robert 		{
515c0b746e5SOllivier Robert 			cutoff = i;
516c0b746e5SOllivier Robert 			span = i;
517c0b746e5SOllivier Robert 		}
518c0b746e5SOllivier Robert 		else
519c0b746e5SOllivier Robert 		    if (histbuf[cutoff] == histbuf[i])
520c0b746e5SOllivier Robert 		    {
521c0b746e5SOllivier Robert 			    span = i;
522c0b746e5SOllivier Robert 		    }
523c0b746e5SOllivier Robert 	}
524c0b746e5SOllivier Robert 
525c0b746e5SOllivier Robert 	cutoff = (cutoff + span) / 2;
526c0b746e5SOllivier Robert 
527c0b746e5SOllivier Robert 	parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: lower maximum %d, higher maximum %d, cutoff %d\n", lowmax, highmax, cutoff));
528c0b746e5SOllivier Robert 
529c0b746e5SOllivier Robert 	s = (unsigned char *)buffer;
530ea906c41SOllivier Robert 	while (s < e)
531c0b746e5SOllivier Robert 	{
532c0b746e5SOllivier Robert 		if (*s == (unsigned char)~0)
533c0b746e5SOllivier Robert 		{
534c0b746e5SOllivier Robert 			*s = '?';
535c0b746e5SOllivier Robert 		}
536c0b746e5SOllivier Robert 		else
537c0b746e5SOllivier Robert 		{
538c0b746e5SOllivier Robert 			*s = (*s >= cutoff) ? *b : *c;
539c0b746e5SOllivier Robert 		}
540c0b746e5SOllivier Robert 		s++;
541ea906c41SOllivier Robert 		if (*b) b++;
542ea906c41SOllivier Robert 		if (*c) c++;
543c0b746e5SOllivier Robert 	}
544c0b746e5SOllivier Robert 
545276da39aSCy Schubert 	*s = '\0';
546276da39aSCy Schubert 
547a151a66cSOllivier Robert         if (rtc == CVT_NONE)
548a151a66cSOllivier Robert         {
549a151a66cSOllivier Robert 	       rtc = convert_rawdcf(buffer, size, &dcfparameter, clock_time);
550a151a66cSOllivier Robert 	       if (rtc == CVT_OK)
551a151a66cSOllivier Robert 	       {
552a151a66cSOllivier Robert 			time_t newtime;
553a151a66cSOllivier Robert 
554a151a66cSOllivier Robert 			newtime = parse_to_unixtime(clock_time, &rtc);
555a151a66cSOllivier Robert 			if ((rtc == CVT_OK) && t)
556a151a66cSOllivier Robert 			{
557276da39aSCy Schubert 				if ((newtime - t->tcode) <= 600) /* require a successful telegram within last 10 minutes */
558a151a66cSOllivier Robert 				{
559276da39aSCy Schubert 				        parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: recent timestamp check OK\n"));
560a151a66cSOllivier Robert 					clock_time->utctime = newtime;
561a151a66cSOllivier Robert 				}
562a151a66cSOllivier Robert 				else
563a151a66cSOllivier Robert 				{
564276da39aSCy Schubert 					parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: recent timestamp check FAIL - ignore timestamp\n"));
565276da39aSCy Schubert 					rtc = CVT_SKIP;
566a151a66cSOllivier Robert 				}
567a151a66cSOllivier Robert 				t->tcode            = newtime;
568a151a66cSOllivier Robert 			}
569a151a66cSOllivier Robert 	       }
570a151a66cSOllivier Robert         }
571a151a66cSOllivier Robert 
572a151a66cSOllivier Robert     	return rtc;
573c0b746e5SOllivier Robert }
574c0b746e5SOllivier Robert 
575c0b746e5SOllivier Robert /*
576a25439b6SCy Schubert  * parse_pps_fnc_t pps_rawdcf
577c0b746e5SOllivier Robert  *
578c0b746e5SOllivier Robert  * currently a very stupid version - should be extended to decode
579c0b746e5SOllivier Robert  * also ones and zeros (which is easy)
580c0b746e5SOllivier Robert  */
581c0b746e5SOllivier Robert /*ARGSUSED*/
582c0b746e5SOllivier Robert static u_long
583c0b746e5SOllivier Robert pps_rawdcf(
584ea906c41SOllivier Robert 	parse_t *parseio,
585ea906c41SOllivier Robert 	int status,
586ea906c41SOllivier Robert 	timestamp_t *ptime
587c0b746e5SOllivier Robert 	)
588c0b746e5SOllivier Robert {
589c0b746e5SOllivier Robert 	if (!status)		/* negative edge for simpler wiring (Rx->DCD) */
590c0b746e5SOllivier Robert 	{
591c0b746e5SOllivier Robert 		parseio->parse_dtime.parse_ptime  = *ptime;
592c0b746e5SOllivier Robert 		parseio->parse_dtime.parse_state |= PARSEB_PPS|PARSEB_S_PPS;
593c0b746e5SOllivier Robert 	}
594c0b746e5SOllivier Robert 
595c0b746e5SOllivier Robert 	return CVT_NONE;
596c0b746e5SOllivier Robert }
597c0b746e5SOllivier Robert 
598276da39aSCy Schubert static long
599276da39aSCy Schubert calc_usecdiff(
600276da39aSCy Schubert 	timestamp_t *ref,
601276da39aSCy Schubert 	timestamp_t *base,
602276da39aSCy Schubert 	long         offset
603276da39aSCy Schubert 	)
604276da39aSCy Schubert {
605276da39aSCy Schubert 	struct timeval delta;
606276da39aSCy Schubert 	long delta_usec = 0;
607276da39aSCy Schubert 
608276da39aSCy Schubert #ifdef PARSEKERNEL
609276da39aSCy Schubert 	delta.tv_sec = ref->tv.tv_sec - offset - base->tv.tv_sec;
610276da39aSCy Schubert 	delta.tv_usec = ref->tv.tv_usec - base->tv.tv_usec;
611276da39aSCy Schubert 	if (delta.tv_usec < 0)
612276da39aSCy Schubert 	{
613276da39aSCy Schubert 		delta.tv_sec  -= 1;
614276da39aSCy Schubert 		delta.tv_usec += 1000000;
615276da39aSCy Schubert 	}
616276da39aSCy Schubert #else
617276da39aSCy Schubert 	l_fp delt;
618276da39aSCy Schubert 
619276da39aSCy Schubert 	delt = ref->fp;
620276da39aSCy Schubert 	delt.l_i -= offset;
621276da39aSCy Schubert 	L_SUB(&delt, &base->fp);
622276da39aSCy Schubert 	TSTOTV(&delt, &delta);
623276da39aSCy Schubert #endif
624276da39aSCy Schubert 
625276da39aSCy Schubert 	delta_usec = 1000000 * (int32_t)delta.tv_sec + delta.tv_usec;
626276da39aSCy Schubert 	return delta_usec;
627276da39aSCy Schubert }
628276da39aSCy Schubert 
629c0b746e5SOllivier Robert static u_long
630c0b746e5SOllivier Robert snt_rawdcf(
631ea906c41SOllivier Robert 	parse_t *parseio,
632ea906c41SOllivier Robert 	timestamp_t *ptime
633c0b746e5SOllivier Robert 	)
634c0b746e5SOllivier Robert {
635276da39aSCy Schubert 	/*
636276da39aSCy Schubert 	 * only synthesize if all of following conditions are met:
637276da39aSCy Schubert 	 * - CVT_OK parse_status (we have a time stamp base)
638276da39aSCy Schubert 	 * - ABS(ptime - tminute - (parse_index - 1) sec) < 500ms (spaced by 1 sec +- 500ms)
639276da39aSCy Schubert 	 * - minute marker is available (confirms minute raster as base)
640276da39aSCy Schubert 	 */
641276da39aSCy Schubert 	last_tcode_t  *t = (last_tcode_t *)parseio->parse_pdata;
642276da39aSCy Schubert 	long delta_usec = -1;
643276da39aSCy Schubert 
644276da39aSCy Schubert 	if (t != NULL && t->tminute.tv.tv_sec != 0) {
645276da39aSCy Schubert 		delta_usec = calc_usecdiff(ptime, &t->tminute, parseio->parse_index - 1);
646276da39aSCy Schubert 		if (delta_usec < 0)
647276da39aSCy Schubert 			delta_usec = -delta_usec;
648276da39aSCy Schubert 	}
649276da39aSCy Schubert 
650276da39aSCy Schubert 	parseprintf(DD_RAWDCF,("parse: snt_rawdcf: synth for offset %d seconds - absolute usec error %ld\n",
651276da39aSCy Schubert 			       parseio->parse_index - 1, delta_usec));
652276da39aSCy Schubert 
653276da39aSCy Schubert 	if (((parseio->parse_dtime.parse_status & CVT_MASK) == CVT_OK) &&
654276da39aSCy Schubert 	    (delta_usec < 500000 && delta_usec >= 0)) /* only if minute marker is available */
655c0b746e5SOllivier Robert 	{
656c0b746e5SOllivier Robert 		parseio->parse_dtime.parse_stime = *ptime;
657c0b746e5SOllivier Robert 
658c0b746e5SOllivier Robert #ifdef PARSEKERNEL
659c0b746e5SOllivier Robert 		parseio->parse_dtime.parse_time.tv.tv_sec++;
660c0b746e5SOllivier Robert #else
661c0b746e5SOllivier Robert 		parseio->parse_dtime.parse_time.fp.l_ui++;
662c0b746e5SOllivier Robert #endif
663c0b746e5SOllivier Robert 
664c0b746e5SOllivier Robert 		parseprintf(DD_RAWDCF,("parse: snt_rawdcf: time stamp synthesized offset %d seconds\n", parseio->parse_index - 1));
665c0b746e5SOllivier Robert 
666c0b746e5SOllivier Robert 		return updatetimeinfo(parseio, parseio->parse_lstate);
667c0b746e5SOllivier Robert 	}
668c0b746e5SOllivier Robert 	return CVT_NONE;
669c0b746e5SOllivier Robert }
670c0b746e5SOllivier Robert 
671c0b746e5SOllivier Robert /*
672a25439b6SCy Schubert  * parse_inp_fnc_t inp_rawdcf
673c0b746e5SOllivier Robert  *
674ea906c41SOllivier Robert  * grab DCF77 data from input stream
675c0b746e5SOllivier Robert  */
676c0b746e5SOllivier Robert static u_long
677c0b746e5SOllivier Robert inp_rawdcf(
678c0b746e5SOllivier Robert 	  parse_t      *parseio,
679a25439b6SCy Schubert 	  char         ch,
680c0b746e5SOllivier Robert 	  timestamp_t  *tstamp
681c0b746e5SOllivier Robert 	  )
682c0b746e5SOllivier Robert {
683c0b746e5SOllivier Robert 	static struct timeval timeout = { 1, 500000 }; /* 1.5 secongs denote second #60 */
684c0b746e5SOllivier Robert 
6853311ff84SXin LI 	parseprintf(DD_PARSE, ("inp_rawdcf(0x%p, 0x%x, ...)\n", (void*)parseio, ch));
686c0b746e5SOllivier Robert 
687c0b746e5SOllivier Robert 	parseio->parse_dtime.parse_stime = *tstamp; /* collect timestamp */
688c0b746e5SOllivier Robert 
689c0b746e5SOllivier Robert 	if (parse_timedout(parseio, tstamp, &timeout))
690c0b746e5SOllivier Robert 	{
691276da39aSCy Schubert 		last_tcode_t *t = (last_tcode_t *)parseio->parse_pdata;
692276da39aSCy Schubert 		long delta_usec;
693c0b746e5SOllivier Robert 
694276da39aSCy Schubert 		parseprintf(DD_RAWDCF, ("inp_rawdcf: time out seen\n"));
695276da39aSCy Schubert 		/* finish collection */
696c0b746e5SOllivier Robert 		(void) parse_end(parseio);
697276da39aSCy Schubert 
698276da39aSCy Schubert 		if (t != NULL)
699276da39aSCy Schubert 		{
700276da39aSCy Schubert 			/* remember minute start sample time if timeouts occur in minute raster */
701276da39aSCy Schubert 			if (t->timeout.tv.tv_sec != 0)
702276da39aSCy Schubert 			{
703276da39aSCy Schubert 				delta_usec = calc_usecdiff(tstamp, &t->timeout, 60);
704276da39aSCy Schubert 				if (delta_usec < 0)
705276da39aSCy Schubert 					delta_usec = -delta_usec;
706276da39aSCy Schubert 			}
707276da39aSCy Schubert 			else
708276da39aSCy Schubert 			{
709276da39aSCy Schubert 				delta_usec = -1;
710276da39aSCy Schubert 			}
711276da39aSCy Schubert 
712276da39aSCy Schubert 			if (delta_usec < 500000 && delta_usec >= 0)
713276da39aSCy Schubert 			{
714276da39aSCy Schubert 				parseprintf(DD_RAWDCF, ("inp_rawdcf: timeout time difference %ld usec - minute marker set\n", delta_usec));
715276da39aSCy Schubert 				/* collect minute markers only if spaced by 60 seconds */
716276da39aSCy Schubert 				t->tminute = *tstamp;
717276da39aSCy Schubert 			}
718276da39aSCy Schubert 			else
719276da39aSCy Schubert 			{
720276da39aSCy Schubert 				parseprintf(DD_RAWDCF, ("inp_rawdcf: timeout time difference %ld usec - minute marker cleared\n", delta_usec));
721276da39aSCy Schubert 				memset((char *)&t->tminute, 0, sizeof(t->tminute));
722276da39aSCy Schubert 			}
723276da39aSCy Schubert 			t->timeout = *tstamp;
724276da39aSCy Schubert 		}
725c0b746e5SOllivier Robert 		(void) parse_addchar(parseio, ch);
726276da39aSCy Schubert 
727276da39aSCy Schubert 		/* pass up to higher layers */
728c0b746e5SOllivier Robert 		return PARSE_INP_TIME;
729c0b746e5SOllivier Robert 	}
730c0b746e5SOllivier Robert 	else
731c0b746e5SOllivier Robert 	{
732c0b746e5SOllivier Robert 		unsigned int rtc;
733c0b746e5SOllivier Robert 
734c0b746e5SOllivier Robert 		rtc = parse_addchar(parseio, ch);
735c0b746e5SOllivier Robert 		if (rtc == PARSE_INP_SKIP)
736c0b746e5SOllivier Robert 		{
737c0b746e5SOllivier Robert 			if (snt_rawdcf(parseio, tstamp) == CVT_OK)
738c0b746e5SOllivier Robert 				return PARSE_INP_SYNTH;
739c0b746e5SOllivier Robert 		}
740c0b746e5SOllivier Robert 		return rtc;
741c0b746e5SOllivier Robert 	}
742c0b746e5SOllivier Robert }
743c0b746e5SOllivier Robert 
744c0b746e5SOllivier Robert #else /* not (REFCLOCK && CLOCK_PARSE && CLOCK_RAWDCF) */
745*f5f40dd6SCy Schubert NONEMPTY_TRANSLATION_UNIT
746c0b746e5SOllivier Robert #endif /* not (REFCLOCK && CLOCK_PARSE && CLOCK_RAWDCF) */
747c0b746e5SOllivier Robert 
748c0b746e5SOllivier Robert /*
749c0b746e5SOllivier Robert  * History:
750c0b746e5SOllivier Robert  *
751c0b746e5SOllivier Robert  * clk_rawdcf.c,v
752ea906c41SOllivier Robert  * Revision 4.18  2006/06/22 18:40:01  kardel
753ea906c41SOllivier Robert  * clean up signedness (gcc 4)
754ea906c41SOllivier Robert  *
755ea906c41SOllivier Robert  * Revision 4.17  2006/01/22 16:01:55  kardel
756ea906c41SOllivier Robert  * update version information
757ea906c41SOllivier Robert  *
758ea906c41SOllivier Robert  * Revision 4.16  2006/01/22 15:51:22  kardel
759ea906c41SOllivier Robert  * generate reasonable timecode output on invalid input
760ea906c41SOllivier Robert  *
761ea906c41SOllivier Robert  * Revision 4.15  2005/08/06 19:17:06  kardel
762ea906c41SOllivier Robert  * clean log output
763ea906c41SOllivier Robert  *
764ea906c41SOllivier Robert  * Revision 4.14  2005/08/06 17:39:40  kardel
765ea906c41SOllivier Robert  * cleanup size handling wrt/ to buffer boundaries
766ea906c41SOllivier Robert  *
767ea906c41SOllivier Robert  * Revision 4.13  2005/04/16 17:32:10  kardel
768ea906c41SOllivier Robert  * update copyright
769ea906c41SOllivier Robert  *
770ea906c41SOllivier Robert  * Revision 4.12  2004/11/14 15:29:41  kardel
771ea906c41SOllivier Robert  * support PPSAPI, upgrade Copyright to Berkeley style
772ea906c41SOllivier Robert  *
773a151a66cSOllivier Robert  * Revision 4.9  1999/12/06 13:42:23  kardel
774a151a66cSOllivier Robert  * transfer correctly converted time codes always into tcode
775a151a66cSOllivier Robert  *
776a151a66cSOllivier Robert  * Revision 4.8  1999/11/28 09:13:50  kardel
777a151a66cSOllivier Robert  * RECON_4_0_98F
778a151a66cSOllivier Robert  *
779a151a66cSOllivier Robert  * Revision 4.7  1999/04/01 20:07:20  kardel
780a151a66cSOllivier Robert  * added checking for minutie increment of timestamps in clk_rawdcf.c
781a151a66cSOllivier Robert  *
782c0b746e5SOllivier Robert  * Revision 4.6  1998/06/14 21:09:37  kardel
783c0b746e5SOllivier Robert  * Sun acc cleanup
784c0b746e5SOllivier Robert  *
785c0b746e5SOllivier Robert  * Revision 4.5  1998/06/13 12:04:16  kardel
786c0b746e5SOllivier Robert  * fix SYSV clock name clash
787c0b746e5SOllivier Robert  *
788c0b746e5SOllivier Robert  * Revision 4.4  1998/06/12 15:22:28  kardel
789c0b746e5SOllivier Robert  * fix prototypes
790c0b746e5SOllivier Robert  *
791c0b746e5SOllivier Robert  * Revision 4.3  1998/06/06 18:33:36  kardel
792c0b746e5SOllivier Robert  * simplified condidional compile expression
793c0b746e5SOllivier Robert  *
794c0b746e5SOllivier Robert  * Revision 4.2  1998/05/24 11:04:18  kardel
795c0b746e5SOllivier Robert  * triggering PPS on negative edge for simpler wiring (Rx->DCD)
796c0b746e5SOllivier Robert  *
797c0b746e5SOllivier Robert  * Revision 4.1  1998/05/24 09:39:53  kardel
798c0b746e5SOllivier Robert  * implementation of the new IO handling model
799c0b746e5SOllivier Robert  *
800c0b746e5SOllivier Robert  * Revision 4.0  1998/04/10 19:45:30  kardel
801c0b746e5SOllivier Robert  * Start 4.0 release version numbering
802c0b746e5SOllivier Robert  *
803c0b746e5SOllivier Robert  * from V3 3.24 log info deleted 1998/04/11 kardel
804c0b746e5SOllivier Robert  *
805c0b746e5SOllivier Robert  */
806