2 * /src/NTP/REPOSITORY/ntp4-dev/libparse/clk_rawdcf.c,v 4.18 2006/06/22 18:40:01 kardel RELEASE_20060622_A
4 * clk_rawdcf.c,v 4.18 2006/06/22 18:40:01 kardel RELEASE_20060622_A
6 * Raw DCF77 pulse clock support
8 * Copyright (c) 1995-2015 by Frank Kardel <kardel <AT> ntp.org>
9 * Copyright (c) 1989-1994 by Frank Kardel, Friedrich-Alexander Universitaet Erlangen-Nuernberg, Germany
11 * Redistribution and use in source and binary forms, with or without
12 * modification, are permitted provided that the following conditions
14 * 1. Redistributions of source code must retain the above copyright
15 * notice, this list of conditions and the following disclaimer.
16 * 2. Redistributions in binary form must reproduce the above copyright
17 * notice, this list of conditions and the following disclaimer in the
18 * documentation and/or other materials provided with the distribution.
19 * 3. Neither the name of the author nor the names of its contributors
20 * may be used to endorse or promote products derived from this software
21 * without specific prior written permission.
23 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
24 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
25 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
26 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
27 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
28 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
29 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
30 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
31 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
32 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
41 #if defined(REFCLOCK) && defined(CLOCK_PARSE) && defined(CLOCK_RAWDCF)
44 #include "timevalops.h"
45 #include "ntp_unixtime.h"
46 #include "ntp_calendar.h"
50 # include <sys/parsestreams.h>
54 # include "ntp_stdlib.h"
60 * From "Zur Zeit", Physikalisch-Technische Bundesanstalt (PTB), Braunschweig
61 * und Berlin, Maerz 1989
63 * Timecode transmission:
65 * time marks are send every second except for the second before the
67 * time marks consist of a reduction of transmitter power to 25%
68 * of the nominal level
69 * the falling edge is the time indication (on time)
70 * time marks of a 100ms duration constitute a logical 0
71 * time marks of a 200ms duration constitute a logical 1
73 * see the spec. (basically a (non-)inverted psuedo random phase shift)
77 * 0 - 10 AM: free, FM: 0
79 * 15 R - "call bit" used to signalize irregularities in the control facilities
80 * (until 2003 indicated transmission via alternate antenna)
81 * 16 A1 - expect zone change (1 hour before)
82 * 17 - 18 Z1,Z2 - time zone
85 * 1 0 MESZ (MED, MET DST)
87 * 19 A2 - expect leap insertion/deletion (1 hour before)
88 * 20 S - start of time code (1)
89 * 21 - 24 M1 - BCD (lsb first) Minutes
90 * 25 - 27 M10 - BCD (lsb first) 10 Minutes
91 * 28 P1 - Minute Parity (even)
92 * 29 - 32 H1 - BCD (lsb first) Hours
93 * 33 - 34 H10 - BCD (lsb first) 10 Hours
94 * 35 P2 - Hour Parity (even)
95 * 36 - 39 D1 - BCD (lsb first) Days
96 * 40 - 41 D10 - BCD (lsb first) 10 Days
97 * 42 - 44 DW - BCD (lsb first) day of week (1: Monday -> 7: Sunday)
98 * 45 - 49 MO - BCD (lsb first) Month
100 * 51 - 53 Y1 - BCD (lsb first) Years
101 * 54 - 57 Y10 - BCD (lsb first) 10 Years
102 * 58 P3 - Date Parity (even)
103 * 59 - usually missing (minute indication), except for leap insertion
106 static parse_pps_fnc_t pps_rawdcf;
107 static parse_cvt_fnc_t cvt_rawdcf;
108 static parse_inp_fnc_t inp_rawdcf;
110 typedef struct last_tcode {
111 time_t tcode; /* last converted time code */
112 timestamp_t tminute; /* sample time for minute start */
113 timestamp_t timeout; /* last timeout timestamp */
116 #define BUFFER_MAX 61
118 clockformat_t clock_rawdcf =
120 inp_rawdcf, /* DCF77 input handling */
121 cvt_rawdcf, /* raw dcf input conversion */
122 pps_rawdcf, /* examining PPS information */
123 0, /* no private configuration data */
124 "RAW DCF77 Timecode", /* direct decoding / time synthesis */
126 BUFFER_MAX, /* bit buffer */
130 static struct dcfparam
132 const unsigned char *onebits;
133 const unsigned char *zerobits;
136 (const unsigned char *)"###############RADMLS1248124P124812P1248121241248112481248P??", /* 'ONE' representation */
137 (const unsigned char *)"--------------------s-------p------p----------------------p__" /* 'ZERO' representation */
140 static struct rawdcfcode
142 char offset; /* start bit */
145 { 0 }, { 15 }, { 16 }, { 17 }, { 19 }, { 20 }, { 21 }, { 25 }, { 28 }, { 29 },
146 { 33 }, { 35 }, { 36 }, { 40 }, { 42 }, { 45 }, { 49 }, { 50 }, { 54 }, { 58 }, { 59 }
172 char offset; /* start bit of parity field */
175 { 21 }, { 29 }, { 36 }, { 59 }
182 #define DCF_Z_MET 0x2
183 #define DCF_Z_MED 0x1
189 const unsigned char *zero
195 first = rawdcfcode[idx].offset;
197 for (i = rawdcfcode[idx+1].offset - 1; i >= first; i--)
200 sum |= (buf[i] != zero[i]);
209 const unsigned char *zero
215 last = partab[idx+1].offset;
217 for (i = partab[idx].offset; i < last; i++)
218 psum ^= (buf[i] != zero[i]);
225 unsigned char *buffer,
227 struct dcfparam *dcfprm,
228 clocktime_t *clock_time
231 unsigned char *s = buffer;
232 const unsigned char *b = dcfprm->onebits;
233 const unsigned char *c = dcfprm->zerobits;
236 parseprintf(DD_RAWDCF,("parse: convert_rawdcf: \"%.*s\"\n", size, buffer));
241 msyslog(LOG_ERR, "parse: convert_rawdcf: INCOMPLETE DATA - time code only has %d bits", size);
243 return CVT_FAIL|CVT_BADFMT;
246 for (i = 0; i < size; i++)
248 if ((*s != *b) && (*s != *c))
251 * we only have two types of bytes (ones and zeros)
254 msyslog(LOG_ERR, "parse: convert_rawdcf: BAD DATA - no conversion");
256 return CVT_FAIL|CVT_BADFMT;
264 * check Start and Parity bits
266 if ((ext_bf(buffer, DCF_S, dcfprm->zerobits) == 1) &&
267 pcheck(buffer, DCF_P_P1, dcfprm->zerobits) &&
268 pcheck(buffer, DCF_P_P2, dcfprm->zerobits) &&
269 pcheck(buffer, DCF_P_P3, dcfprm->zerobits))
274 parseprintf(DD_RAWDCF,("parse: convert_rawdcf: parity check passed\n"));
276 clock_time->flags = PARSEB_S_CALLBIT|PARSEB_S_LEAP;
277 clock_time->utctime= 0;
278 clock_time->usecond= 0;
279 clock_time->second = 0;
280 clock_time->minute = ext_bf(buffer, DCF_M10, dcfprm->zerobits);
281 clock_time->minute = TIMES10(clock_time->minute) + ext_bf(buffer, DCF_M1, dcfprm->zerobits);
282 clock_time->hour = ext_bf(buffer, DCF_H10, dcfprm->zerobits);
283 clock_time->hour = TIMES10(clock_time->hour) + ext_bf(buffer, DCF_H1, dcfprm->zerobits);
284 clock_time->day = ext_bf(buffer, DCF_D10, dcfprm->zerobits);
285 clock_time->day = TIMES10(clock_time->day) + ext_bf(buffer, DCF_D1, dcfprm->zerobits);
286 clock_time->month = ext_bf(buffer, DCF_MO0, dcfprm->zerobits);
287 clock_time->month = TIMES10(clock_time->month) + ext_bf(buffer, DCF_MO, dcfprm->zerobits);
288 clock_time->year = ext_bf(buffer, DCF_Y10, dcfprm->zerobits);
289 clock_time->year = TIMES10(clock_time->year) + ext_bf(buffer, DCF_Y1, dcfprm->zerobits);
291 switch (ext_bf(buffer, DCF_Z, dcfprm->zerobits))
294 clock_time->utcoffset = -1*60*60;
298 clock_time->flags |= PARSEB_DST;
299 clock_time->utcoffset = -2*60*60;
303 parseprintf(DD_RAWDCF,("parse: convert_rawdcf: BAD TIME ZONE\n"));
304 return CVT_FAIL|CVT_BADFMT;
307 if (ext_bf(buffer, DCF_A1, dcfprm->zerobits))
308 clock_time->flags |= PARSEB_ANNOUNCE;
310 if (ext_bf(buffer, DCF_A2, dcfprm->zerobits))
311 clock_time->flags |= PARSEB_LEAPADD; /* default: DCF77 data format deficiency */
313 if (ext_bf(buffer, DCF_R, dcfprm->zerobits))
314 clock_time->flags |= PARSEB_CALLBIT;
316 parseprintf(DD_RAWDCF,("parse: convert_rawdcf: TIME CODE OK: %02d:%02d, %02d.%02d.%02d, flags 0x%lx\n",
317 (int)clock_time->hour, (int)clock_time->minute, (int)clock_time->day, (int)clock_time->month,(int) clock_time->year,
318 (u_long)clock_time->flags));
324 * bad format - not for us
327 msyslog(LOG_ERR, "parse: convert_rawdcf: start bit / parity check FAILED for \"%.*s\"", size, buffer);
329 return CVT_FAIL|CVT_BADFMT;
334 * parse_cvt_fnc_t cvt_rawdcf
335 * raw dcf input routine - needs to fix up 50 baud
336 * characters for 1/0 decision
340 unsigned char *buffer,
342 struct format *param,
343 clocktime_t *clock_time,
347 last_tcode_t *t = (last_tcode_t *)local;
348 unsigned char *s = (unsigned char *)buffer;
349 unsigned char *e = s + size;
350 const unsigned char *b = dcfparameter.onebits;
351 const unsigned char *c = dcfparameter.zerobits;
352 u_long rtc = CVT_NONE;
353 unsigned int i, lowmax, highmax, cutoff, span;
355 unsigned char histbuf[BITS];
357 * the input buffer contains characters with runs of consecutive
358 * bits set. These set bits are an indication of the DCF77 pulse
359 * length. We assume that we receive the pulse at 50 Baud. Thus
360 * a 100ms pulse would generate a 4 bit train (20ms per bit and
362 * a 200ms pulse would create all zeroes (and probably a frame error)
365 for (i = 0; i < BITS; i++)
375 unsigned int ch = *s ^ 0xFF;
377 * these lines are left as an excercise to the reader 8-)
379 if (!((ch+1) & ch) || !*s)
387 *s = (unsigned char) i;
394 parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: character check for 0x%x@%d FAILED\n", *s, (int)(s - (unsigned char *)buffer)));
395 *s = (unsigned char)~0;
396 rtc = CVT_FAIL|CVT_BADFMT;
407 cutoff = 4; /* doesn't really matter - it'll fail anyway, but gives error output */
410 parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: average bit count: %d\n", cutoff));
415 parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: histogram:"));
416 for (i = 0; i <= cutoff; i++)
418 lowmax+=histbuf[i] * i;
419 highmax += histbuf[i];
420 parseprintf(DD_RAWDCF,(" %d", histbuf[i]));
422 parseprintf(DD_RAWDCF, (" <M>"));
424 lowmax += highmax / 2;
438 for (; i < BITS; i++)
440 highmax+=histbuf[i] * i;
442 parseprintf(DD_RAWDCF,(" %d", histbuf[i]));
444 parseprintf(DD_RAWDCF,("\n"));
455 span = cutoff = lowmax;
456 for (i = lowmax; i <= highmax; i++)
458 if (histbuf[cutoff] > histbuf[i])
464 if (histbuf[cutoff] == histbuf[i])
470 cutoff = (cutoff + span) / 2;
472 parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: lower maximum %d, higher maximum %d, cutoff %d\n", lowmax, highmax, cutoff));
474 s = (unsigned char *)buffer;
477 if (*s == (unsigned char)~0)
483 *s = (*s >= cutoff) ? *b : *c;
494 rtc = convert_rawdcf(buffer, size, &dcfparameter, clock_time);
499 newtime = parse_to_unixtime(clock_time, &rtc);
500 if ((rtc == CVT_OK) && t)
502 if ((newtime - t->tcode) <= 600) /* require a successful telegram within last 10 minutes */
504 parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: recent timestamp check OK\n"));
505 clock_time->utctime = newtime;
509 parseprintf(DD_RAWDCF,("parse: cvt_rawdcf: recent timestamp check FAIL - ignore timestamp\n"));
521 * parse_pps_fnc_t pps_rawdcf
523 * currently a very stupid version - should be extended to decode
524 * also ones and zeros (which is easy)
534 if (!status) /* negative edge for simpler wiring (Rx->DCD) */
536 parseio->parse_dtime.parse_ptime = *ptime;
537 parseio->parse_dtime.parse_state |= PARSEB_PPS|PARSEB_S_PPS;
550 struct timeval delta;
554 delta.tv_sec = ref->tv.tv_sec - offset - base->tv.tv_sec;
555 delta.tv_usec = ref->tv.tv_usec - base->tv.tv_usec;
556 if (delta.tv_usec < 0)
559 delta.tv_usec += 1000000;
566 L_SUB(&delt, &base->fp);
567 TSTOTV(&delt, &delta);
570 delta_usec = 1000000 * (int32_t)delta.tv_sec + delta.tv_usec;
581 * only synthesize if all of following conditions are met:
582 * - CVT_OK parse_status (we have a time stamp base)
583 * - ABS(ptime - tminute - (parse_index - 1) sec) < 500ms (spaced by 1 sec +- 500ms)
584 * - minute marker is available (confirms minute raster as base)
586 last_tcode_t *t = (last_tcode_t *)parseio->parse_pdata;
587 long delta_usec = -1;
589 if (t != NULL && t->tminute.tv.tv_sec != 0) {
590 delta_usec = calc_usecdiff(ptime, &t->tminute, parseio->parse_index - 1);
592 delta_usec = -delta_usec;
595 parseprintf(DD_RAWDCF,("parse: snt_rawdcf: synth for offset %d seconds - absolute usec error %ld\n",
596 parseio->parse_index - 1, delta_usec));
598 if (((parseio->parse_dtime.parse_status & CVT_MASK) == CVT_OK) &&
599 (delta_usec < 500000 && delta_usec >= 0)) /* only if minute marker is available */
601 parseio->parse_dtime.parse_stime = *ptime;
604 parseio->parse_dtime.parse_time.tv.tv_sec++;
606 parseio->parse_dtime.parse_time.fp.l_ui++;
609 parseprintf(DD_RAWDCF,("parse: snt_rawdcf: time stamp synthesized offset %d seconds\n", parseio->parse_index - 1));
611 return updatetimeinfo(parseio, parseio->parse_lstate);
617 * parse_inp_fnc_t inp_rawdcf
619 * grab DCF77 data from input stream
628 static struct timeval timeout = { 1, 500000 }; /* 1.5 secongs denote second #60 */
630 parseprintf(DD_PARSE, ("inp_rawdcf(0x%p, 0x%x, ...)\n", (void*)parseio, ch));
632 parseio->parse_dtime.parse_stime = *tstamp; /* collect timestamp */
634 if (parse_timedout(parseio, tstamp, &timeout))
636 last_tcode_t *t = (last_tcode_t *)parseio->parse_pdata;
639 parseprintf(DD_RAWDCF, ("inp_rawdcf: time out seen\n"));
640 /* finish collection */
641 (void) parse_end(parseio);
645 /* remember minute start sample time if timeouts occur in minute raster */
646 if (t->timeout.tv.tv_sec != 0)
648 delta_usec = calc_usecdiff(tstamp, &t->timeout, 60);
650 delta_usec = -delta_usec;
657 if (delta_usec < 500000 && delta_usec >= 0)
659 parseprintf(DD_RAWDCF, ("inp_rawdcf: timeout time difference %ld usec - minute marker set\n", delta_usec));
660 /* collect minute markers only if spaced by 60 seconds */
661 t->tminute = *tstamp;
665 parseprintf(DD_RAWDCF, ("inp_rawdcf: timeout time difference %ld usec - minute marker cleared\n", delta_usec));
666 memset((char *)&t->tminute, 0, sizeof(t->tminute));
668 t->timeout = *tstamp;
670 (void) parse_addchar(parseio, ch);
672 /* pass up to higher layers */
673 return PARSE_INP_TIME;
679 rtc = parse_addchar(parseio, ch);
680 if (rtc == PARSE_INP_SKIP)
682 if (snt_rawdcf(parseio, tstamp) == CVT_OK)
683 return PARSE_INP_SYNTH;
689 #else /* not (REFCLOCK && CLOCK_PARSE && CLOCK_RAWDCF) */
691 #endif /* not (REFCLOCK && CLOCK_PARSE && CLOCK_RAWDCF) */
697 * Revision 4.18 2006/06/22 18:40:01 kardel
698 * clean up signedness (gcc 4)
700 * Revision 4.17 2006/01/22 16:01:55 kardel
701 * update version information
703 * Revision 4.16 2006/01/22 15:51:22 kardel
704 * generate reasonable timecode output on invalid input
706 * Revision 4.15 2005/08/06 19:17:06 kardel
709 * Revision 4.14 2005/08/06 17:39:40 kardel
710 * cleanup size handling wrt/ to buffer boundaries
712 * Revision 4.13 2005/04/16 17:32:10 kardel
715 * Revision 4.12 2004/11/14 15:29:41 kardel
716 * support PPSAPI, upgrade Copyright to Berkeley style
718 * Revision 4.9 1999/12/06 13:42:23 kardel
719 * transfer correctly converted time codes always into tcode
721 * Revision 4.8 1999/11/28 09:13:50 kardel
724 * Revision 4.7 1999/04/01 20:07:20 kardel
725 * added checking for minutie increment of timestamps in clk_rawdcf.c
727 * Revision 4.6 1998/06/14 21:09:37 kardel
730 * Revision 4.5 1998/06/13 12:04:16 kardel
731 * fix SYSV clock name clash
733 * Revision 4.4 1998/06/12 15:22:28 kardel
736 * Revision 4.3 1998/06/06 18:33:36 kardel
737 * simplified condidional compile expression
739 * Revision 4.2 1998/05/24 11:04:18 kardel
740 * triggering PPS on negative edge for simpler wiring (Rx->DCD)
742 * Revision 4.1 1998/05/24 09:39:53 kardel
743 * implementation of the new IO handling model
745 * Revision 4.0 1998/04/10 19:45:30 kardel
746 * Start 4.0 release version numbering
748 * from V3 3.24 log info deleted 1998/04/11 kardel