nmxp_base.c 19.8 KB
Newer Older
1
2
/*! \file
 *
Matteo Quintiliani's avatar
Matteo Quintiliani committed
3
 * \brief Base for Nanometrics Protocol Library
4
5
6
7
8
9
 *
 * Author:
 * 	Matteo Quintiliani
 * 	Istituto Nazionale di Geofisica e Vulcanologia - Italy
 *	quintiliani@ingv.it
 *
10
 * $Id: nmxp_base.c,v 1.79 2009-08-17 08:19:46 mtheo Exp $
Matteo Quintiliani's avatar
Matteo Quintiliani committed
11
 *
12
13
 */

Matteo Quintiliani's avatar
Matteo Quintiliani committed
14
#include "config.h"
15
#include "nmxp_base.h"
16
#include "nmxp_memory.h"
Matteo Quintiliani's avatar
Matteo Quintiliani committed
17
#ifdef HAVE_WINDOWS_H
18
#include "nmxp_win.h"
Matteo Quintiliani's avatar
Matteo Quintiliani committed
19
#endif
20
21
22
23
24
25

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <stdarg.h>

Matteo Quintiliani's avatar
Matteo Quintiliani committed
26
27
28
29
#ifdef HAVE_WINDOWS_H
#include "winsock2.h"
#warning You are compiling on Windows MinGW...
#else
30
31
32
33
34
35
36
#include <sys/socket.h>
#include <netinet/in.h>
#include <arpa/inet.h>
#include <netdb.h>
#include <errno.h>
#include <time.h>
#include <unistd.h>
Matteo Quintiliani's avatar
Matteo Quintiliani committed
37
#endif
Matteo Quintiliani's avatar
Matteo Quintiliani committed
38

Matteo Quintiliani's avatar
Matteo Quintiliani committed
39
40
41
#define MAX_OUTDATA 4096


42
int nmxp_openSocket(char *hostname, int portNum, int (*func_cond)(void))
43
44
{
  static int sleepTime = 1;
45
  int isock = -1;
46
47
48
  struct hostent *hostinfo = NULL;
  struct sockaddr_in psServAddr;
  struct in_addr hostaddr;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
49

50
#ifdef HAVE_WINDOWS_H
Matteo Quintiliani's avatar
Matteo Quintiliani committed
51
  nmxp_initWinsock();
52
#endif
Matteo Quintiliani's avatar
Matteo Quintiliani committed
53

54
55
  if (!hostname)
  {
56
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Empty host name?\n");
57
58
59
60
    return -1;
  }

  if ( (hostinfo = gethostbyname(hostname)) == NULL) {
61
62
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Cannot lookup host: %s\n",
	    NMXP_LOG_STR(hostname));
63
64
65
    return -1;
  }

66
  while(!func_cond())
67
68
  {
    isock = socket (AF_INET, SOCK_STREAM, 0);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
69
70
71
72
73
74
75
76

#ifdef HAVE_WINDOWS_H
    if (isock == INVALID_SOCKET) {
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW,"Error at socket()\n");
	    WSACleanup();
	    exit(1);
    }
#else
77
78
    if (isock < 0)
    {
79
      nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Can't open stream socket\n");
80
81
      exit(1);
    }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
82
#endif
83
84
85
86
87
88

    /* Fill in the structure "psServAddr" with the address of server
       that we want to connect with */
    memset (&psServAddr, 0, sizeof(psServAddr));
    psServAddr.sin_family = AF_INET;
    psServAddr.sin_port = htons((unsigned short) portNum);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
89
90
91
92
93
#ifdef HAVE_WINDOWS_H
    unsigned long address;
    memcpy(&address, hostinfo->h_addr, (size_t) hostinfo->h_length);
    psServAddr.sin_addr.s_addr = address;
#else
94
    psServAddr.sin_addr = *(struct in_addr *)hostinfo->h_addr_list[0];
Matteo Quintiliani's avatar
Matteo Quintiliani committed
95
#endif
96
97
98
99

    /* Report action and resolved address */
    memcpy(&hostaddr.s_addr, *hostinfo->h_addr_list,
	   sizeof (hostaddr.s_addr));
100
    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "Attempting to connect to %s port %d\n",
101
	    NMXP_LOG_STR(inet_ntoa(hostaddr)), portNum);
102

103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
    if(connect(isock, (struct sockaddr *)&psServAddr, sizeof(psServAddr)) >= 0) {
	sleepTime = 1;
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "Connection established: socket=%i,IP=%s,port=%d\n",
		isock, NMXP_LOG_STR(inet_ntoa(hostaddr)), portNum);
	return isock;
    } else {
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Connecting to %s port %d. Trying again after %d seconds...\n",
		NMXP_LOG_STR(inet_ntoa(hostaddr)), portNum, sleepTime);
	nmxp_closeSocket(isock);
	isock = -1;

	if(!func_cond()) {
	    nmxp_sleep (sleepTime);
	    sleepTime *= 2;
	    if (sleepTime > NMXP_SLEEPMAX)
		sleepTime = NMXP_SLEEPMAX;
	}
120
    }
121

122
  }
123
  return isock;
124
125
126
}


Matteo Quintiliani's avatar
Matteo Quintiliani committed
127
128
129
130
131
132
133
134
135
136
int nmxp_closeSocket(int isock) {
	int ret;
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "Closed connection.\n");
#ifdef HAVE_WINDOWS_H
	ret = closesocket(isock);
	WSACleanup();
#else
	ret = close(isock);
#endif
	return ret;
137
138
139
140
141
142
143
144
145
146
147
148
149
}


int nmxp_send_ctrl(int isock, void* buffer, int length)
{
  int sendCount = send(isock, (char*) buffer, length, 0);
  
  if (sendCount != length)
    return NMXP_SOCKET_ERROR;

  return NMXP_SOCKET_OK;
}

150
151
152
153
#ifdef HAVE_BROKEN_SO_RCVTIMEO
#warning Managing non-blocking I/O using select()
int nmxp_recv_select_timeout(int s, char *buf, int len, int timeout)
{
Matteo Quintiliani's avatar
Matteo Quintiliani committed
154
#define HIGHEST_TIMEOUT 30
155
156
157
    fd_set fds;
    int n;
    struct timeval tv;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
158
    static int message_times = 0;
159

160
    /* set up the file descriptor set*/
161
162
163
    FD_ZERO(&fds);
    FD_SET(s, &fds);

164
    /* set up the struct timeval for the timeout*/
165
166
167
168
169
170
171
    if(timeout == 0) {
	if(message_times == 0) {
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_ANY, "nmxp_recv_select_timeout(): timeout = %d\n", HIGHEST_TIMEOUT);
	    message_times++;
	}
	timeout = HIGHEST_TIMEOUT;
    }
172
173
174
    tv.tv_sec = timeout;
    tv.tv_usec = 0;

175
    /* wait until timeout or data received*/
Matteo Quintiliani's avatar
Matteo Quintiliani committed
176
    errno = 0;
177
    n = select(s+1, &fds, NULL, NULL, &tv);
178
    if (n == 0) return -2; /* timeout!*/
179
    if(errno == EINTR) return -2; /* timeout! "Interrupted system call" */
180
    if (n == -1) return -1; /* error*/
181
 
182
    /* data must be here, so do a normal recv()*/
183
184
185
    return recv(s, buf, len, 0);
}
#endif
186

187
188
int nmxp_setsockopt_RCVTIMEO(int isock, int timeoutsec) {
    int ret = 0;
189
#ifdef HAVE_WINDOWS_H
190
    int timeos;
191
#else
192
193
194

#ifndef HAVE_BROKEN_SO_RCVTIMEO
    struct timeval timeo;
195
#endif
196

197
#endif
198

199
    if(timeoutsec > 0) {
200
#ifdef HAVE_WINDOWS_H
201
202
203
204
205
206
	timeos  = timeoutsec * 1000;
	ret = setsockopt(isock, SOL_SOCKET, SO_RCVTIMEO, (char *) &timeos, sizeof(timeos));
	if (ret < 0)
	{
	    perror("setsockopt SO_RCVTIMEO");
	}
Matteo Quintiliani's avatar
Matteo Quintiliani committed
207
#else
208
209

#ifndef HAVE_BROKEN_SO_RCVTIMEO
210
211
212
213
214
215
216
217
	timeo.tv_sec  = timeoutsec;
	timeo.tv_usec = 0;
	ret = setsockopt(isock, SOL_SOCKET, SO_RCVTIMEO, &timeo, sizeof(timeo));
	if (ret < 0) {
	    perror("setsockopt SO_RCVTIMEO");
	}
#else
#warning nmxp_setsockopt_RCVTIMEO() do nothing for your system.
Matteo Quintiliani's avatar
Matteo Quintiliani committed
218
#endif
219

220
#endif
221
    }
222

223
224
    return ret;
}
225

Matteo Quintiliani's avatar
Matteo Quintiliani committed
226

227
228
229
230
#define MAXLEN_RECV_ERRNO_STR 200

char *nmxp_strerror(int errno_value) {
    static char ret_recv_errno_str[MAXLEN_RECV_ERRNO_STR];
Matteo Quintiliani's avatar
Matteo Quintiliani committed
231
#ifdef HAVE_WINDOWS_H
232
    char *recv_errno_str;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
233
#else
234

235
#ifdef HAVE_STRERROR_R
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
    char recv_errno_str[MAXLEN_RECV_ERRNO_STR];
#else
    char *recv_errno_str;
#endif

#endif

    ret_recv_errno_str[0] = 0;

#ifdef HAVE_WINDOWS_H
    recv_errno_str = WSAGetLastErrorMessage(errno_value);
#else

#ifdef HAVE_STRERROR_R
    strerror_r(errno_value, recv_errno_str, MAXLEN_RECV_ERRNO_STR);
251
#else
252
    recv_errno_str = strerror(errno_value);
253
254
#endif

Matteo Quintiliani's avatar
Matteo Quintiliani committed
255
#endif
256

257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
    if(recv_errno_str) {
	strncpy(ret_recv_errno_str, recv_errno_str, MAXLEN_RECV_ERRNO_STR);
    }

    return ret_recv_errno_str;
}


int nmxp_recv_ctrl(int isock, void *buffer, int length, int timeoutsec, int *recv_errno )
{
  int recvCount;
  int cc;
  char *buffer_char = buffer;
  char *recv_errno_str = NULL;

272
  nmxp_setsockopt_RCVTIMEO(isock, timeoutsec);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
273
  
Matteo Quintiliani's avatar
Matteo Quintiliani committed
274
  cc = 1;
275
276
  *recv_errno  = 0;
  recvCount = 0;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
277
  while(cc > 0 && *recv_errno == 0  && recvCount < length) {
278

Matteo Quintiliani's avatar
Matteo Quintiliani committed
279
      /* TODO some operating system could not reset errno */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
280
      errno = 0;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
281

282
#ifdef HAVE_BROKEN_SO_RCVTIMEO
283
284
285
      cc = nmxp_recv_select_timeout(isock, buffer_char + recvCount, length - recvCount, timeoutsec);
#else
      cc = recv(isock, buffer_char + recvCount, length - recvCount, 0);
286
287
#endif

288
289
290
#ifdef HAVE_WINDOWS_H
      *recv_errno  = WSAGetLastError();
#else
291
292
293
294
295
      if(cc == -2) {
	  *recv_errno  = EWOULDBLOCK;
      } else {
	  *recv_errno  = errno;
      }
296
#endif
Matteo Quintiliani's avatar
Matteo Quintiliani committed
297
      if(cc <= 0) {
298
	  /*
299
300
	  nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "nmxp_recv_ctrl(): (cc=%d <= 0) errno=%d  recvCount=%d  length=%d\n",
	  cc, *recv_errno, recvCount, length);
301
	  */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
302
303
304
      } else {
	  recvCount += cc;
      }
305
  }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
306

307
  nmxp_setsockopt_RCVTIMEO(isock, 0);
308

Matteo Quintiliani's avatar
Matteo Quintiliani committed
309
  if (recvCount != length  ||  *recv_errno != 0  ||  cc <= 0) {
310

311
      recv_errno_str = nmxp_strerror(*recv_errno);
312

313
#ifdef HAVE_WINDOWS_H
314
      if(*recv_errno != WSAEWOULDBLOCK  &&  *recv_errno != WSAETIMEDOUT)
315
316
317
318
#else
      if(*recv_errno != EWOULDBLOCK)
#endif
      {
319
320
	  nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "nmxp_recv_ctrl(): %s (errno=%d recvCount=%d length=%d cc=%d)\n",
		  NMXP_LOG_STR(recv_errno_str), *recv_errno, recvCount, length, cc);
321
322
323
324
325
326
327
328
      }

      /* TO IMPROVE 
       * Fixed bug receiving zero byte from recv() 'TCP FIN or EOF received'
       * */
      if(cc == 0  &&  *recv_errno == 0) {
	  *recv_errno = -100;
      }
329

Matteo Quintiliani's avatar
Matteo Quintiliani committed
330
#ifdef HAVE_WINDOWS_H
331
      if(recvCount != length || (*recv_errno != WSAEWOULDBLOCK  &&  *recv_errno != WSAETIMEDOUT))
Matteo Quintiliani's avatar
Matteo Quintiliani committed
332
#else
333
      if(recvCount != length || *recv_errno != EWOULDBLOCK)
Matteo Quintiliani's avatar
Matteo Quintiliani committed
334
#endif
335
336
337
      {
	  return NMXP_SOCKET_ERROR;
      }
338
339
340
341
342
343
  }
  
  return NMXP_SOCKET_OK;
}


344
int nmxp_sendHeader(int isock, NMXP_MSG_CLIENT type, int32_t length)
345
{  
Matteo Quintiliani's avatar
Matteo Quintiliani committed
346
    NMXP_MESSAGE_HEADER msg;
347
348
349
350
351

    msg.signature = htonl(NMX_SIGNATURE);
    msg.type      = htonl(type);
    msg.length    = htonl(length);

Matteo Quintiliani's avatar
Matteo Quintiliani committed
352
    return nmxp_send_ctrl(isock, &msg, sizeof(NMXP_MESSAGE_HEADER));
353
354
355
}


Matteo Quintiliani's avatar
Matteo Quintiliani committed
356
int nmxp_receiveHeader(int isock, NMXP_MSG_SERVER *type, int32_t *length, int timeoutsec, int *recv_errno )
357
358
{  
    int ret ;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
359
    NMXP_MESSAGE_HEADER msg;
360

Matteo Quintiliani's avatar
Matteo Quintiliani committed
361
    ret = nmxp_recv_ctrl(isock, &msg, sizeof(NMXP_MESSAGE_HEADER), timeoutsec, recv_errno);
362
363
364
365
366
367
368
369
370
371
372
373

    *type = 0;
    *length = 0;

    if(ret == NMXP_SOCKET_OK) {
	msg.signature = ntohl(msg.signature);
	msg.type      = ntohl(msg.type);
	msg.length    = ntohl(msg.length);

	if (msg.signature != NMX_SIGNATURE)
	{
	    ret = NMXP_SOCKET_ERROR;
374
375
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW,
		    "nmxp_receiveHeader(): signature mismatches. signature = %d, type = %d, length = %d\n",
376
377
378
379
380
381
382
383
384
385
386
		    msg.signature, msg.type, msg.length);
	} else {
	    *type = msg.type;
	    *length = msg.length;
	}
    }

    return ret;
}


387
int nmxp_sendMessage(int isock, NMXP_MSG_CLIENT type, void *buffer, int32_t length) {
388
389
390
391
392
393
394
395
396
397
398
    int ret;
    ret = nmxp_sendHeader(isock, type, length);
    if( ret == NMXP_SOCKET_OK) {
	if(buffer && length > 0) {
	    ret = nmxp_send_ctrl(isock, buffer, length);
	}
    }
    return ret;
}


399
400
401
402
403
404
405
406
407
408
409
410
411
412
int32_t nmxp_display_error_from_server(char *buffer, int32_t length) {
    /* NMXP_MSG_TERMINATESUBSCRIPTION */
    char *str_msg = NULL;
    int32_t reason;
    memcpy(&reason, buffer, sizeof(reason));
    reason = ntohl(reason);
    str_msg = buffer + sizeof(reason);
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_ANY, "%d %s shutdown: %s\n",
	    reason,
	    (reason == 0)? "Normal" : (reason == 1)? "Error" : (reason == 2)? "Timeout" : "Unknown",
	    str_msg);
    return reason;
}

Matteo Quintiliani's avatar
Matteo Quintiliani committed
413
int nmxp_receiveMessage(int isock, NMXP_MSG_SERVER *type, void **buffer, int32_t *length, int timeoutsec, int *recv_errno ) {
414
415
416
417
    int ret;
    *buffer = NULL;
    *length = 0;

Matteo Quintiliani's avatar
Matteo Quintiliani committed
418
    ret = nmxp_receiveHeader(isock, type, length, timeoutsec, recv_errno);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
419

Matteo Quintiliani's avatar
Matteo Quintiliani committed
420
    if( ret == NMXP_SOCKET_OK  ) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
421
	if (*length > 0) {
422
	    *buffer = NMXP_MEM_MALLOC(*length);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
423
	    ret = nmxp_recv_ctrl(isock, *buffer, *length, 0, recv_errno);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
424

425
426
427
428
	    if(*type == NMXP_MSG_TERMINATESUBSCRIPTION) {
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Received TerminateSubscritption.\n");
		nmxp_display_error_from_server(*buffer, *length);
	    } else if(*type == NMXP_MSG_ERROR) {
429
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Received ErrorMessage: %s\n", NMXP_LOG_STR(*buffer));
Matteo Quintiliani's avatar
Matteo Quintiliani committed
430
	    } else {
431
		nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_PACKETMAN, "Received message type: %d  length=%d\n", *type, *length);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
432
433
434
	    }

	}
Matteo Quintiliani's avatar
Matteo Quintiliani committed
435
436
437
    }

    if(*recv_errno != 0) {
438
#ifdef HAVE_WINDOWS_H
439
	if(*recv_errno == WSAEWOULDBLOCK  ||  *recv_errno == WSAETIMEDOUT) {
440
#else
441
	if(*recv_errno == EWOULDBLOCK) {
442
#endif
443
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_DOD, "Timeout receiving in nmxp_receiveMessage()\n");
Matteo Quintiliani's avatar
Matteo Quintiliani committed
444
	} else {
445
446
	    /* Log message is not necessary because managed by nmxp_recv_ctrl() */
	    /* nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Error in nmxp_receiveMessage()\n"); */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
447
	}
448
    }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
449

450
451
452
    return ret;
}

453

454
NMXP_DATA_PROCESS *nmxp_processDecompressedData(char* buffer_data, int length_data, NMXP_CHAN_LIST_NET *channelList, const char *network_code_default)
455
456
457
458
459
460
{
  int32_t   netInt    = 0;
  int32_t   pKey      = 0;
  double    pTime     = 0.0;
  int32_t   pNSamp    = 0;
  int32_t   pSampRate = 0;
461
  int32_t  *pDataPtr  = NULL;
462
463
  int       swap      = 0;
  int       idx;
464
  static int32_t outdata[MAX_OUTDATA];
465

466
467
468
  char station_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
  char channel_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
  char network_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
469

470
  char *nmxp_channel_name = NULL;
471
472
473
474
475
476
477
  static NMXP_DATA_PROCESS pd;

  /* copy the header contents into local fields and swap */
  memcpy(&netInt, &buffer_data[0], 4);
  pKey = ntohl(netInt);
  if ( pKey != netInt ) { swap = 1; }

478
479
480
481
482
483
  nmxp_data_init(&pd);

  nmxp_channel_name = nmxp_chan_lookupName(pKey, channelList);

  if(nmxp_channel_name) {

484
  memcpy(&pTime, &buffer_data[4], 8);
485
  if ( swap ) { nmxp_data_swap_8b(&pTime); }
486
487
488
489
490
491
492

  memcpy(&netInt, &buffer_data[12], 4);
  pNSamp = ntohl(netInt);
  memcpy(&netInt, &buffer_data[16], 4);
  pSampRate = ntohl(netInt);

  /* There should be (length_data - 20) bytes of data as 32-bit ints here */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
493
494
  memcpy(outdata , (int32_t *) &buffer_data[20], length_data - 20);
  pDataPtr = outdata;
495
496
497
498
499
500
501

  /* Swap the data samples to host order */
  for ( idx=0; idx < pNSamp; idx++ ) {
      netInt = ntohl(pDataPtr[idx]);
      pDataPtr[idx] = netInt;
  }

502
  if(!nmxp_chan_cpy_sta_chan(nmxp_channel_name, station_code, channel_code, network_code)) {
503
504
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not in STA.CHAN format: %s\n",
	    NMXP_LOG_STR(nmxp_channel_name));
505
506
507
  }
  
  pd.key = pKey;
508
  if(network_code[0] != 0) {
509
      strncpy(pd.network, network_code, NMXP_DATA_NETWORK_LENGTH);
510
  } else {
511
      strncpy(pd.network, network_code_default, NMXP_DATA_NETWORK_LENGTH);
512
  }
513
  if(station_code[0] != 0) {
514
      strncpy(pd.station, station_code, NMXP_DATA_STATION_LENGTH);
515
  }
516
  if(channel_code[0] != 0) {
517
      strncpy(pd.channel, channel_code, NMXP_DATA_CHANNEL_LENGTH);
518
  }
519
  pd.packet_type = NMXP_MSG_DECOMPRESSED;
520
521
522
  pd.x0 = -1;
  pd.xn = -1;
  pd.x0n_significant = 0;
523
524
525
  /* TODO*/
  /* pd.oldest_seq_no = ;*/
  /* pd.seq_no = ;*/
526
527
528
529
530
  pd.time = pTime;
  pd.nSamp = pNSamp;
  pd.pDataPtr = pDataPtr;
  pd.sampRate = pSampRate;

531
532
  } else {
      nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not found for key %d\n", pKey);
533
534
  }

535
536
537
538
  return &pd;
}


539
NMXP_DATA_PROCESS *nmxp_processCompressedData(char* buffer_data, int length_data, NMXP_CHAN_LIST_NET *channelList, const char *network_code_default)
540
541
542
543
544
{
    int32_t   pKey      = 0;
    double    pTime     = 0.0;
    int32_t   pNSamp    = 0;
    int32_t   pSampRate = 0;
545
    int32_t  *pDataPtr  = NULL;
546

547
548
549
    char station_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
    char channel_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
    char network_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
550
551
552

    static NMXP_DATA_PROCESS pd;

553
    int32_t nmx_rate_code_to_sample_rate[32] = {
554
555
556
557
558
	0,1,2,5,10,20,40,50,
	80,100,125,200,250,500,1000,25,
	120,0,0,0,0,0,0,0,
	0,0,0,0,0,0,0,0};

559
	int32_t nmx_oldest_sequence_number;
560
561
	char nmx_hdr[25];
	unsigned char nmx_ptype;
562
	int32_t nmx_seconds;
563
	double nmx_seconds_double;
564
565
	int16_t nmx_ticks, nmx_instr_id;
	int32_t nmx_seqno;
566
	unsigned char nmx_sample_rate;
567
568
	int32_t nmx_x0;
	int32_t rate_code, chan_code, this_sample_rate;
569

570
	int32_t comp_bytecount;
571
	unsigned char *indata;
572
	static int32_t outdata[MAX_OUTDATA];
573
574
	int32_t nout, i, k;
	int32_t prev_xn;
575
576
	const uint32_t high_scale = 4096 * 2048;
	const uint32_t high_scale_p = 4096 * 4096;
577

578
579
	char *nmxp_channel_name = NULL;

580
	/* TOREMOVE int my_order = get_my_wordorder();*/
Matteo Quintiliani's avatar
Matteo Quintiliani committed
581
	int my_host_is_bigendian = nmxp_data_bigendianhost();
582
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "my_host_is_bigendian %d\n", my_host_is_bigendian);
583
584

	memcpy(&nmx_oldest_sequence_number, buffer_data, 4);
585
586
587
	if (my_host_is_bigendian) {
	    nmxp_data_swap_4b (&nmx_oldest_sequence_number);
	}
588
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Oldest sequence number = %d\n", nmx_oldest_sequence_number);
589
590
591
592
593
594

	memcpy(nmx_hdr, buffer_data+4, 17);
	/* Decode the Nanometrics packet header bundle. */
	memcpy (&nmx_ptype, nmx_hdr+0, 1);
	if ( (nmx_ptype & 0xf) == 9) {
	    /* Filler packet.  Discard entire packet.   */
595
	    nmxp_log (NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Filler packet - discarding\n");
596
	    /*m continue;*/
597
598
599
600
601
602
603
604
605
606
607
	    exit(0);
	}

	nmx_x0 = 0;
	memcpy (&nmx_seconds, nmx_hdr+1, 4);
	memcpy (&nmx_ticks, nmx_hdr+5, 2);
	memcpy (&nmx_instr_id, nmx_hdr+7, 2);
	memcpy (&nmx_seqno, nmx_hdr+9, 4);
	memcpy (&nmx_sample_rate, nmx_hdr+13, 1);
	memcpy (&nmx_x0, nmx_hdr+14, 3);

608
	if (my_host_is_bigendian) {
609
	    nmxp_data_swap_4b ((int32_t *)&nmx_seconds);
610
611
612
613
	    nmxp_data_swap_2b (&nmx_ticks);
	    nmxp_data_swap_2b (&nmx_instr_id);
	    nmxp_data_swap_4b (&nmx_seqno);
	    nmxp_data_swap_4b (&nmx_x0);
614
	}
615
616
617
618
619
620
621

	/* check if nmx_x0 is negative like as signed 3-byte int */
	if( (nmx_x0 & high_scale) ==  high_scale) {
	    /* nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "WARNING: changed nmx_x0, old value = %d\n",  nmx_x0);*/
	    nmx_x0 -= high_scale_p;
	}

622
623
624
625
626
	nmx_seconds_double = (double) nmx_seconds + ( (double) nmx_ticks / 10000.0 );
	rate_code = nmx_sample_rate>>3;
	chan_code = nmx_sample_rate&7;
	this_sample_rate = nmx_rate_code_to_sample_rate[rate_code];

627
628
629
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_ptype          = %d\n", nmx_ptype);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_seconds        = %d\n", nmx_seconds);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_ticks          = %d\n", nmx_ticks);
630

631
632
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_seconds_double = %f\n", nmx_seconds_double);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_x0             = %d\n", nmx_x0);
633

634
635
636
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_instr_id       = %d\n", nmx_instr_id);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_seqno          = %d\n", nmx_seqno);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_sample_rate    = %d\n", nmx_sample_rate);
637
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "this_sample_rate   = %d\n", this_sample_rate);
638

639
640
641
642
643
644
	pKey = (nmx_instr_id << 16) | ( 1 << 8) | ( chan_code);

	pTime = nmx_seconds_double;

	pSampRate = this_sample_rate;

645
646
647
648
649
650
651
	nmxp_data_init(&pd);

	nmxp_channel_name = nmxp_chan_lookupName(pKey, channelList);

	if(nmxp_channel_name) {

	if(!nmxp_chan_cpy_sta_chan(nmxp_channel_name, station_code, channel_code, network_code)) {
652
653
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not in STA.CHAN format: %s\n",
		    NMXP_LOG_STR(nmxp_channel_name));
654
655
	}
  
656
657
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Channel key %d for %s.%s\n",
		pKey, NMXP_LOG_STR(station_code), NMXP_LOG_STR(channel_code));
658

659
660
661
	comp_bytecount = length_data-21;
	indata = (unsigned char *) buffer_data + 21;

662
663
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "comp_bytecount     = %d  (N = %.2f)\n", comp_bytecount, (double) comp_bytecount / 17.0);

664
665
666
667
668
669
	/* Unpack the data bundles, each 17 bytes long. */
	prev_xn = nmx_x0;
	outdata[0] = nmx_x0;
	nout = 1;
	for (i=0; i<comp_bytecount; i+=17) {
	    if (i+17>comp_bytecount) {
670
		nmxp_log (NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "comp_bytecount = %d, i+17 = %d\n",
671
672
673
674
			comp_bytecount, i+17);
		exit(1);
	    }
	    if (nout+16 > MAX_OUTDATA)  {
675
		nmxp_log (NMXP_LOG_ERR,  NMXP_LOG_D_PACKETMAN, "Output buffer size too small\n");
676
677
678
		exit(1);
	    }
	    k = nmxp_data_unpack_bundle (outdata+nout,indata+i,&prev_xn);
679
680
681
682
683
	    if (k < 0) nmxp_log (NMXP_LOG_WARN, NMXP_LOG_D_PACKETMAN, "Null bundle: %s.%s.%s (k=%d) %s %d\n",
		    NMXP_LOG_STR(network_code),
		    NMXP_LOG_STR(station_code),
		    NMXP_LOG_STR(channel_code), k,
		    __FILE__,  __LINE__);
684
685
686
687
688
689
	    if (k < 0) break;
	    nout += k;
	    /* prev_xn = outdata[nout-1]; */
	}
	nout--;

690
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Unpacked %d samples.\n", nout);
691
692
693
694
695
696

	pDataPtr = outdata;

	pNSamp = nout;

	pd.key = pKey;
697
	if(network_code[0] != 0) {
698
	    strncpy(pd.network, network_code, NMXP_DATA_NETWORK_LENGTH);
699
	} else {
700
	    strncpy(pd.network, network_code_default, NMXP_DATA_NETWORK_LENGTH);
701
	}
702
	if(station_code[0] != 0) {
703
	    strncpy(pd.station, station_code, NMXP_DATA_STATION_LENGTH);
704
	}
705
	if(channel_code[0] != 0) {
706
	    strncpy(pd.channel, channel_code, NMXP_DATA_CHANNEL_LENGTH);
707
	}
708
709
	pd.packet_type = nmx_ptype;
	pd.x0 = nmx_x0;
710
	pd.xn = pDataPtr[nout];
711
	pd.x0n_significant = 1;
712
	pd.oldest_seq_no = nmx_oldest_sequence_number;
713
714
715
716
717
718
	pd.seq_no = nmx_seqno;
	pd.time = pTime;
	pd.nSamp = pNSamp;
	pd.pDataPtr = pDataPtr;
	pd.sampRate = pSampRate;

719
720
	} else {
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not found for key %d\n", pKey);
721
722
	}

723
724
725
	return &pd;
}

726
727
728

unsigned int nmxp_sleep(unsigned int sleep_time) {
#ifdef HAVE_WINDOWS_H
729
    Sleep(sleep_time * 1000);
730
    return 0;
731
732
733
734
735
#else
    return sleep(sleep_time);
#endif
}

Matteo Quintiliani's avatar
Matteo Quintiliani committed
736
737
738
739
740
741
742
743
744
unsigned int nmxp_usleep(unsigned int usleep_time) {
#ifdef HAVE_WINDOWS_H
    Sleep((usleep_time+500)/1000);
    return 0;
#else
    return usleep(usleep_time);
#endif
}