contrib/ntp/parseutil/testdcf.c

   1 /*
   2  * /src/NTP/REPOSITORY/v4/parseutil/testdcf.c,v 3.18 1996/12/01 16:05:04 kardel Exp
   3  *
   4  * testdcf.c,v 3.18 1996/12/01 16:05:04 kardel Exp
   5  *
   6  * simple DCF77 100/200ms pulse test program (via 50Baud serial line)
   7  *
   8  * Copyright (c) 1993,1994,1995,1996, 1998 by Frank Kardel
   9  * Friedrich-Alexander Universität Erlangen-Nürnberg, Germany
  10  *
  11  * This program is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
  14  *
  15  * This program may not be sold or used for profit without prior
  16  * written consent of the author.
  17  */
  18
  19 #include <unistd.h>
  20 #include <stdio.h>
  21 #include <fcntl.h>
  22 #include <termios.h>
  23 #include <sys/types.h>
  24 #include <sys/time.h>
  25
  26 #include "ntp_stdlib.h"
  27
  28 /*
  29  * state flags
  30  */
  31 #define DCFB_ANNOUNCE           0x0001 /* switch time zone warning (DST switch) */
  32 #define DCFB_DST                0x0002 /* DST in effect */
  33 #define DCFB_LEAP               0x0004 /* LEAP warning (1 hour prior to occurence) */
  34 #define DCFB_ALTERNATE          0x0008 /* alternate antenna used */
  35
  36 struct clocktime                /* clock time broken up from time code */
  37 {
  38         long wday;
  39         long day;
  40         long month;
  41         long year;
  42         long hour;
  43         long minute;
  44         long second;
  45         long usecond;
  46         long utcoffset; /* in minutes */
  47         long flags;             /* current clock status */
  48 };
  49
  50 typedef struct clocktime clocktime_t;
  51
  52 #define TIMES10(_X_) (((_X_) << 3) + ((_X_) << 1))
  53
  54 /*
  55  * parser related return/error codes
  56  */
  57 #define CVT_MASK        0x0000000F /* conversion exit code */
  58 #define   CVT_NONE      0x00000001 /* format not applicable */
  59 #define   CVT_FAIL      0x00000002 /* conversion failed - error code returned */
  60 #define   CVT_OK        0x00000004 /* conversion succeeded */
  61 #define CVT_BADFMT      0x00000010 /* general format error - (unparsable) */
  62
  63 /*
  64  * DCF77 raw time code
  65  *
  66  * From "Zur Zeit", Physikalisch-Technische Bundesanstalt (PTB), Braunschweig
  67  * und Berlin, Maerz 1989
  68  *
  69  * Timecode transmission:
  70  * AM:
  71  *      time marks are send every second except for the second before the
  72  *      next minute mark
  73  *      time marks consist of a reduction of transmitter power to 25%
  74  *      of the nominal level
  75  *      the falling edge is the time indication (on time)
  76  *      time marks of a 100ms duration constitute a logical 0
  77  *      time marks of a 200ms duration constitute a logical 1
  78  * FM:
  79  *      see the spec. (basically a (non-)inverted psuedo random phase shift)
  80  *
  81  * Encoding:
  82  * Second       Contents
  83  * 0  - 10      AM: free, FM: 0
  84  * 11 - 14      free
  85  * 15           R     - alternate antenna
  86  * 16           A1    - expect zone change (1 hour before)
  87  * 17 - 18      Z1,Z2 - time zone
  88  *               0  0 illegal
  89  *               0  1 MEZ  (MET)
  90  *               1  0 MESZ (MED, MET DST)
  91  *               1  1 illegal
  92  * 19           A2    - expect leap insertion/deletion (1 hour before)
  93  * 20           S     - start of time code (1)
  94  * 21 - 24      M1    - BCD (lsb first) Minutes
  95  * 25 - 27      M10   - BCD (lsb first) 10 Minutes
  96  * 28           P1    - Minute Parity (even)
  97  * 29 - 32      H1    - BCD (lsb first) Hours
  98  * 33 - 34      H10   - BCD (lsb first) 10 Hours
  99  * 35           P2    - Hour Parity (even)
 100  * 36 - 39      D1    - BCD (lsb first) Days
 101  * 40 - 41      D10   - BCD (lsb first) 10 Days
 102  * 42 - 44      DW    - BCD (lsb first) day of week (1: Monday -> 7: Sunday)
 103  * 45 - 49      MO    - BCD (lsb first) Month
 104  * 50           MO0   - 10 Months
 105  * 51 - 53      Y1    - BCD (lsb first) Years
 106  * 54 - 57      Y10   - BCD (lsb first) 10 Years
 107  * 58           P3    - Date Parity (even)
 108  * 59                 - usually missing (minute indication), except for leap insertion
 109  */
 110
 111 static struct rawdcfcode
 112 {
 113         char offset;                    /* start bit */
 114 } rawdcfcode[] =
 115 {
 116         {  0 }, { 15 }, { 16 }, { 17 }, { 19 }, { 20 }, { 21 }, { 25 }, { 28 }, { 29 },
 117         { 33 }, { 35 }, { 36 }, { 40 }, { 42 }, { 45 }, { 49 }, { 50 }, { 54 }, { 58 }, { 59 }
 118 };
 119
 120 #define DCF_M   0
 121 #define DCF_R   1
 122 #define DCF_A1  2
 123 #define DCF_Z   3
 124 #define DCF_A2  4
 125 #define DCF_S   5
 126 #define DCF_M1  6
 127 #define DCF_M10 7
 128 #define DCF_P1  8
 129 #define DCF_H1  9
 130 #define DCF_H10 10
 131 #define DCF_P2  11
 132 #define DCF_D1  12
 133 #define DCF_D10 13
 134 #define DCF_DW  14
 135 #define DCF_MO  15
 136 #define DCF_MO0 16
 137 #define DCF_Y1  17
 138 #define DCF_Y10 18
 139 #define DCF_P3  19
 140
 141 static struct partab
 142 {
 143         char offset;                    /* start bit of parity field */
 144 } partab[] =
 145 {
 146         { 21 }, { 29 }, { 36 }, { 59 }
 147 };
 148
 149 #define DCF_P_P1        0
 150 #define DCF_P_P2        1
 151 #define DCF_P_P3        2
 152
 153 #define DCF_Z_MET 0x2
 154 #define DCF_Z_MED 0x1
 155
 156 static unsigned long
 157 ext_bf(
 158         register unsigned char *buf,
 159         register int   idx
 160         )
 161 {
 162         register unsigned long sum = 0;
 163         register int i, first;
 164
 165         first = rawdcfcode[idx].offset;
 166
 167         for (i = rawdcfcode[idx+1].offset - 1; i >= first; i--)
 168         {
 169                 sum <<= 1;
 170                 sum |= (buf[i] != '-');
 171         }
 172         return sum;
 173 }
 174
 175 static unsigned
 176 pcheck(
 177         register unsigned char *buf,
 178         register int   idx
 179         )
 180 {
 181         register int i,last;
 182         register unsigned psum = 1;
 183
 184         last = partab[idx+1].offset;
 185
 186         for (i = partab[idx].offset; i < last; i++)
 187             psum ^= (buf[i] != '-');
 188
 189         return psum;
 190 }
 191
 192 static unsigned long
 193 convert_rawdcf(
 194         register unsigned char   *buffer,
 195         register int              size,
 196         register clocktime_t     *clock_time
 197         )
 198 {
 199         if (size < 57)
 200         {
 201                 printf("%-30s", "*** INCOMPLETE");
 202                 return CVT_NONE;
 203         }
 204
 205         /*
 206          * check Start and Parity bits
 207          */
 208         if ((ext_bf(buffer, DCF_S) == 1) &&
 209             pcheck(buffer, DCF_P_P1) &&
 210             pcheck(buffer, DCF_P_P2) &&
 211             pcheck(buffer, DCF_P_P3))
 212         {
 213                 /*
 214                  * buffer OK
 215                  */
 216
 217                 clock_time->flags  = 0;
 218                 clock_time->usecond= 0;
 219                 clock_time->second = 0;
 220                 clock_time->minute = ext_bf(buffer, DCF_M10);
 221                 clock_time->minute = TIMES10(clock_time->minute) + ext_bf(buffer, DCF_M1);
 222                 clock_time->hour   = ext_bf(buffer, DCF_H10);
 223                 clock_time->hour   = TIMES10(clock_time->hour) + ext_bf(buffer, DCF_H1);
 224                 clock_time->day    = ext_bf(buffer, DCF_D10);
 225                 clock_time->day    = TIMES10(clock_time->day) + ext_bf(buffer, DCF_D1);
 226                 clock_time->month  = ext_bf(buffer, DCF_MO0);
 227                 clock_time->month  = TIMES10(clock_time->month) + ext_bf(buffer, DCF_MO);
 228                 clock_time->year   = ext_bf(buffer, DCF_Y10);
 229                 clock_time->year   = TIMES10(clock_time->year) + ext_bf(buffer, DCF_Y1);
 230                 clock_time->wday   = ext_bf(buffer, DCF_DW);
 231
 232                 switch (ext_bf(buffer, DCF_Z))
 233                 {
 234                     case DCF_Z_MET:
 235                         clock_time->utcoffset = -60;
 236                         break;
 237
 238                     case DCF_Z_MED:
 239                         clock_time->flags     |= DCFB_DST;
 240                         clock_time->utcoffset  = -120;
 241                         break;
 242
 243                     default:
 244                         printf("%-30s", "*** BAD TIME ZONE");
 245                         return CVT_FAIL|CVT_BADFMT;
 246                 }
 247
 248                 if (ext_bf(buffer, DCF_A1))
 249                     clock_time->flags |= DCFB_ANNOUNCE;
 250
 251                 if (ext_bf(buffer, DCF_A2))
 252                     clock_time->flags |= DCFB_LEAP;
 253
 254                 if (ext_bf(buffer, DCF_R))
 255                     clock_time->flags |= DCFB_ALTERNATE;
 256
 257                 return CVT_OK;
 258         }
 259         else
 260         {
 261                 /*
 262                  * bad format - not for us
 263                  */
 264                 printf("%-30s", "*** BAD FORMAT (invalid/parity)");
 265                 return CVT_FAIL|CVT_BADFMT;
 266         }
 267 }
 268
 269 char
 270 type(
 271         unsigned int c
 272         )
 273 {
 274         c ^= 0xFF;
 275         return (c > 0xF);
 276 }
 277
 278 static const char *wday[8] =
 279 {
 280         "??",
 281         "Mo",
 282         "Tu",
 283         "We",
 284         "Th",
 285         "Fr",
 286         "Sa",
 287         "Su"
 288 };
 289
 290 static char pat[] = "-\\|/";
 291
 292 #define LINES (24-2)    /* error lines after which the two headlines are repeated */
 293
 294 int
 295 main(
 296         int argc,
 297         char *argv[]
 298         )
 299 {
 300         if ((argc != 2) && (argc != 3))
 301         {
 302                 fprintf(stderr, "usage: %s [-f|-t|-ft|-tf] <device>\n", argv[0]);
 303                 exit(1);
 304         }
 305         else
 306         {
 307                 unsigned char c;
 308                 char *file;
 309                 int fd;
 310                 int offset = 15;
 311                 int trace = 0;
 312                 int errs = LINES+1;
 313
 314                 /*
 315                  * SIMPLE(!) argument "parser"
 316                  */
 317                 if (argc == 3)
 318                 {
 319                         if (strcmp(argv[1], "-f") == 0)
 320                             offset = 0;
 321                         if (strcmp(argv[1], "-t") == 0)
 322                             trace = 1;
 323                         if ((strcmp(argv[1], "-ft") == 0) ||
 324                             (strcmp(argv[1], "-tf") == 0))
 325                         {
 326                                 offset = 0;
 327                                 trace = 1;
 328                         }
 329                         file = argv[2];
 330                 }
 331                 else
 332                 {
 333                         file = argv[1];
 334                 }
 335
 336                 fd = open(file, O_RDONLY);
 337                 if (fd == -1)
 338                 {
 339                         perror(file);
 340                         exit(1);
 341                 }
 342                 else
 343                 {
 344                         int i;
 345 #ifdef TIOCM_RTS
 346                         int on = TIOCM_RTS;
 347 #endif
 348                         struct timeval t, tt, tlast;
 349                         char buf[61];
 350                         clocktime_t clock_time;
 351                         struct termios term;
 352                         int rtc = CVT_NONE;
 353
 354                         if (tcgetattr(fd,  &term) == -1)
 355                         {
 356                                 perror("tcgetattr");
 357                                 exit(1);
 358                         }
 359
 360                         memset(term.c_cc, 0, sizeof(term.c_cc));
 361                         term.c_cc[VMIN] = 1;
 362 #ifdef NO_PARENB_IGNPAR /* Was: defined(SYS_IRIX4) || defined (SYS_IRIX5) */
 363                         /* somehow doesn't grok PARENB & IGNPAR (mj) */
 364                         term.c_cflag = B50|CS8|CREAD|CLOCAL;
 365 #else
 366                         term.c_cflag = B50|CS8|CREAD|CLOCAL|PARENB;
 367 #endif
 368                         term.c_iflag = IGNPAR;
 369                         term.c_oflag = 0;
 370                         term.c_lflag = 0;
 371
 372                         if (tcsetattr(fd, TCSANOW, &term) == -1)
 373                         {
 374                                 perror("tcsetattr");
 375                                 exit(1);
 376                         }
 377
 378 #ifdef I_POP
 379                         while (ioctl(fd, I_POP, 0) == 0)
 380                             ;
 381 #endif
 382 #if defined(TIOCMBIC) && defined(TIOCM_RTS)
 383                         if (ioctl(fd, TIOCMBIC, (caddr_t)&on) == -1)
 384                         {
 385                                 perror("TIOCM_RTS");
 386                         }
 387 #endif
 388
 389                         printf("  DCF77 monitor - Copyright (C) 1993-1996, Frank Kardel\n\n");
 390
 391                         clock_time.hour = 0;
 392                         clock_time.minute = 0;
 393                         clock_time.day = 0;
 394                         clock_time.wday = 0;
 395                         clock_time.month = 0;
 396                         clock_time.year = 0;
 397                         clock_time.flags = 0;
 398                         buf[60] = '\0';
 399                         for ( i = 0; i < 60; i++)
 400                             buf[i] = '.';
 401
 402                         gettimeofday(&tlast, 0L);
 403                         i = 0;
 404                         while (read(fd, &c, 1) == 1)
 405                         {
 406                                 gettimeofday(&t, 0L);
 407                                 tt = t;
 408                                 t.tv_sec -= tlast.tv_sec;
 409                                 t.tv_usec -= tlast.tv_usec;
 410                                 if (t.tv_usec < 0)
 411                                 {
 412                                         t.tv_usec += 1000000;
 413                                         t.tv_sec  -= 1;
 414                                 }
 415
 416                                 if (errs > LINES)
 417                                 {
 418                                         printf("  %s", &"PTB private....RADMLSMin....PHour..PMDay..DayMonthYear....P\n"[offset]);
 419                                         printf("  %s", &"---------------RADMLS1248124P124812P1248121241248112481248P\n"[offset]);
 420                                         errs = 0;
 421                                 }
 422
 423                                 if (t.tv_sec > 1 ||
 424                                     (t.tv_sec == 1 &&
 425                                      t.tv_usec > 500000))
 426                                 {
 427                                         printf("%c %.*s ", pat[i % (sizeof(pat)-1)], 59 - offset, &buf[offset]);
 428
 429                                         if ((rtc = convert_rawdcf((unsigned char *)buf, i, &clock_time)) != CVT_OK)
 430                                         {
 431                                                 printf("\n");
 432                                                 clock_time.hour = 0;
 433                                                 clock_time.minute = 0;
 434                                                 clock_time.day = 0;
 435                                                 clock_time.wday = 0;
 436                                                 clock_time.month = 0;
 437                                                 clock_time.year = 0;
 438                                                 clock_time.flags = 0;
 439                                                 errs++;
 440                                         }
 441
 442                                         if (((c^0xFF)+1) & (c^0xFF))
 443                                             buf[0] = '?';
 444                                         else
 445                                             buf[0] = type(c) ? '#' : '-';
 446
 447                                         for ( i = 1; i < 60; i++)
 448                                             buf[i] = '.';
 449
 450                                         i = 0;
 451                                 }
 452                                 else
 453                                 {
 454                                         if (((c^0xFF)+1) & (c^0xFF))
 455                                             buf[i] = '?';
 456                                         else
 457                                             buf[i] = type(c) ? '#' : '-';
 458
 459                                         printf("%c %.*s ", pat[i % (sizeof(pat)-1)], 59 - offset, &buf[offset]);
 460                                 }
 461
 462                                 if (rtc == CVT_OK)
 463                                 {
 464                                         printf("%s, %2d:%02d:%02d, %d.%02d.%02d, <%s%s%s%s>",
 465                                                wday[clock_time.wday],
 466                                                (int)clock_time.hour, (int)clock_time.minute, (int)i, (int)clock_time.day, (int)clock_time.month,
 467                                                (int)clock_time.year,
 468                                                (clock_time.flags & DCFB_ALTERNATE) ? "R" : "_",
 469                                                (clock_time.flags & DCFB_ANNOUNCE) ? "A" : "_",
 470                                                (clock_time.flags & DCFB_DST) ? "D" : "_",
 471                                                (clock_time.flags & DCFB_LEAP) ? "L" : "_"
 472                                                );
 473                                         if (trace && (i == 0))
 474                                         {
 475                                                 printf("\n");
 476                                                 errs++;
 477                                         }
 478                                 }
 479
 480                                 printf("\r");
 481
 482                                 if (i < 60)
 483                                 {
 484                                         i++;
 485                                 }
 486
 487                                 tlast = tt;
 488
 489                                 fflush(stdout);
 490                         }
 491                         close(fd);
 492                 }
 493         }
 494         return 0;
 495 }