nmxp_base.c 21.3 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.84 2010-09-15 13:14:33 mtheo Exp $
Matteo Quintiliani's avatar
Matteo Quintiliani committed
11
 *
12 13
 */

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>

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>
37
#endif
Matteo Quintiliani's avatar
Matteo Quintiliani committed
38

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

41
int nmxp_openSocket(char *hostname, int portNum, int (*func_cond)(void))
42
{
43
  /*TODO stefano avoid static*/
44
  static int sleepTime = 1;
45
  int isock = -1;
46 47 48
  struct hostent *hostinfo = NULL;
  struct sockaddr_in psServAddr;
  struct in_addr hostaddr;
49

50
#ifdef HAVE_WINDOWS_H
51
  nmxp_initWinsock();
52
#endif
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);
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);
    }
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);
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];
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
}


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 154 155 156
#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)
{
    fd_set fds;
    int n;
    struct timeval tv;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
157
    static int message_times = 0;
158

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

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

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

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

#ifndef HAVE_BROKEN_SO_RCVTIMEO
    struct timeval timeo;
194
#endif
195

196
#endif
197

198 199 200 201
    if(timeoutsec == 0) {
	timeoutsec = NMXP_HIGHEST_TIMEOUT;
    }

202
    if(timeoutsec > 0) {
203
#ifdef HAVE_WINDOWS_H
204 205 206 207 208 209
	timeos  = timeoutsec * 1000;
	ret = setsockopt(isock, SOL_SOCKET, SO_RCVTIMEO, (char *) &timeos, sizeof(timeos));
	if (ret < 0)
	{
	    perror("setsockopt SO_RCVTIMEO");
	}
210
#else
211 212

#ifndef HAVE_BROKEN_SO_RCVTIMEO
213 214 215 216 217 218 219 220
	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.
221
#endif
222

223
#endif
224
    }
225

226 227
    return ret;
}
228

Matteo Quintiliani's avatar
Matteo Quintiliani committed
229

230 231 232
#define MAXLEN_RECV_ERRNO_STR 200

char *nmxp_strerror(int errno_value) {
233
    char * ret_recv_errno_str;
234
#ifdef HAVE_WINDOWS_H
235
    char *recv_errno_str;
236
#else
237

238
#ifdef HAVE_STRERROR_R
239
    char recv_errno_str[MAXLEN_RECV_ERRNO_STR]="";
240
#else
241
    char *recv_errno_str=NULL; 
242 243 244
#endif

#endif
245
    ret_recv_errno_str= (char *) NMXP_MEM_MALLOC (MAXLEN_RECV_ERRNO_STR * sizeof(char));
246 247 248 249 250 251 252 253
    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);
254
#else
255
    recv_errno_str = strerror(errno_value);
256 257
#endif

258
#endif
259

Matteo Quintiliani's avatar
Matteo Quintiliani committed
260
#ifndef HAVE_STRERROR_R
261
    if(recv_errno_str) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
262
#endif
263
	strncpy(ret_recv_errno_str, recv_errno_str, MAXLEN_RECV_ERRNO_STR);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
264
#ifndef HAVE_STRERROR_R
265
    }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
266
#endif
267 268 269 270 271 272 273 274 275 276 277 278

    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;

279
  nmxp_setsockopt_RCVTIMEO(isock, timeoutsec);
280
  
281
  cc = 1;
282 283
  *recv_errno  = 0;
  recvCount = 0;
284
  while(cc > 0 && *recv_errno == 0  && recvCount < length) {
285

Matteo Quintiliani's avatar
Matteo Quintiliani committed
286
      /* TODO some operating system could not reset errno */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
287
      errno = 0;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
288

289
#ifdef HAVE_BROKEN_SO_RCVTIMEO
290 291 292
      cc = nmxp_recv_select_timeout(isock, buffer_char + recvCount, length - recvCount, timeoutsec);
#else
      cc = recv(isock, buffer_char + recvCount, length - recvCount, 0);
293 294
#endif

295 296 297
#ifdef HAVE_WINDOWS_H
      *recv_errno  = WSAGetLastError();
#else
298 299 300 301 302
      if(cc == -2) {
	  *recv_errno  = EWOULDBLOCK;
      } else {
	  *recv_errno  = errno;
      }
303
#endif
304
      if(cc <= 0) {
305
	  /*
306 307
	  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);
308
	  */
309 310 311
      } else {
	  recvCount += cc;
      }
312
  }
313

314
  nmxp_setsockopt_RCVTIMEO(isock, 0);
315

316
  if (recvCount != length  ||  *recv_errno != 0  ||  cc <= 0) {
317

318
      recv_errno_str = nmxp_strerror(*recv_errno);
319

320
#ifdef HAVE_WINDOWS_H
321
      if(*recv_errno != WSAEWOULDBLOCK  &&  *recv_errno != WSAETIMEDOUT)
322 323 324 325
#else
      if(*recv_errno != EWOULDBLOCK)
#endif
      {
326 327
	  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);
328
      }
329
       NMXP_MEM_FREE(recv_errno_str);
330 331 332 333 334 335
      /* TO IMPROVE 
       * Fixed bug receiving zero byte from recv() 'TCP FIN or EOF received'
       * */
      if(cc == 0  &&  *recv_errno == 0) {
	  *recv_errno = -100;
      }
336

337
#ifdef HAVE_WINDOWS_H
338
      if(recvCount != length || (*recv_errno != WSAEWOULDBLOCK  &&  *recv_errno != WSAETIMEDOUT))
339
#else
340
      if(recvCount != length || *recv_errno != EWOULDBLOCK)
341
#endif
342 343 344
      {
	  return NMXP_SOCKET_ERROR;
      }
345 346 347 348 349 350
  }
  
  return NMXP_SOCKET_OK;
}


351
int nmxp_sendHeader(int isock, NMXP_MSG_CLIENT type, int32_t length)
352
{  
Matteo Quintiliani's avatar
Matteo Quintiliani committed
353
    NMXP_MESSAGE_HEADER msg;
354 355 356 357 358

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

Matteo Quintiliani's avatar
Matteo Quintiliani committed
359
    return nmxp_send_ctrl(isock, &msg, sizeof(NMXP_MESSAGE_HEADER));
360 361 362
}


363
int nmxp_receiveHeader(int isock, NMXP_MSG_SERVER *type, int32_t *length, int timeoutsec, int *recv_errno )
364 365
{  
    int ret ;
366
    NMXP_MESSAGE_HEADER msg={0};
367

368
    ret = nmxp_recv_ctrl(isock, &msg, sizeof(NMXP_MESSAGE_HEADER), timeoutsec, recv_errno);
369 370 371 372

    *type = 0;
    *length = 0;

373
    if((ret == NMXP_SOCKET_OK) && (msg.type != 0)) {
374 375 376 377 378 379 380
	msg.signature = ntohl(msg.signature);
	msg.type      = ntohl(msg.type);
	msg.length    = ntohl(msg.length);

	if (msg.signature != NMX_SIGNATURE)
	{
	    ret = NMXP_SOCKET_ERROR;
381 382
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW,
		    "nmxp_receiveHeader(): signature mismatches. signature = %d, type = %d, length = %d\n",
383 384 385 386 387 388 389 390 391 392 393
		    msg.signature, msg.type, msg.length);
	} else {
	    *type = msg.type;
	    *length = msg.length;
	}
    }

    return ret;
}


394
int nmxp_sendMessage(int isock, NMXP_MSG_CLIENT type, void *buffer, int32_t length) {
395 396 397 398 399 400 401 402 403 404 405
    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;
}


406 407 408 409 410 411 412 413 414 415 416 417 418 419
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;
}

420
int nmxp_receiveMessage(int isock, NMXP_MSG_SERVER *type, void *buffer, int32_t *length, int timeoutsec, int *recv_errno, int buffer_length) {
421 422 423
    int ret;
    *length = 0;

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

Matteo Quintiliani's avatar
Matteo Quintiliani committed
426
    if( ret == NMXP_SOCKET_OK  ) {
427 428 429 430 431 432 433

	if(*length > buffer_length) {
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_ANY, "nmxp_receiveMessage(): size of received messagge is bigger than buffer. (%d > %d). \n",
		    *length, buffer_length);
	    ret = NMXP_SOCKET_ERROR;
	} else if (*length > 0) {
	    ret = nmxp_recv_ctrl(isock, buffer, *length, 0, recv_errno);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
434

435 436
	    if(*type == NMXP_MSG_TERMINATESUBSCRIPTION) {
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Received TerminateSubscritption.\n");
437
		nmxp_display_error_from_server(buffer, *length);
438
	    } else if(*type == NMXP_MSG_ERROR) {
439
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Received ErrorMessage: %s\n", NMXP_LOG_STR(buffer));
Matteo Quintiliani's avatar
Matteo Quintiliani committed
440
	    } else {
441
		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
442 443 444
	    }

	}
Matteo Quintiliani's avatar
Matteo Quintiliani committed
445 446 447
    }

    if(*recv_errno != 0) {
448
#ifdef HAVE_WINDOWS_H
449
	if(*recv_errno == WSAEWOULDBLOCK  ||  *recv_errno == WSAETIMEDOUT) {
450
#else
451
	if(*recv_errno == EWOULDBLOCK) {
452
#endif
453
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_DOD, "Timeout receiving in nmxp_receiveMessage()\n");
Matteo Quintiliani's avatar
Matteo Quintiliani committed
454
	} else {
455 456
	    /* 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"); */
457
	}
458
    }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
459

460 461 462
    return ret;
}

463

464
NMXP_DATA_PROCESS *nmxp_processDecompressedData(char* buffer_data, int length_data, NMXP_CHAN_LIST_NET *channelList, const char *network_code_default, const char *location_code_default)
465 466 467 468 469 470
{
  int32_t   netInt    = 0;
  int32_t   pKey      = 0;
  double    pTime     = 0.0;
  int32_t   pNSamp    = 0;
  int32_t   pSampRate = 0;
471
  int32_t  *pDataPtr  = NULL;
472 473
  int       swap      = 0;
  int       idx;
474 475
  int32_t  *outdata   = NULL;

476

477 478 479
  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];
480
  char location_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
481

482
  char *nmxp_channel_name = NULL;
483
  NMXP_DATA_PROCESS *pd   = NULL;
484

485 486 487 488
  pd= (NMXP_DATA_PROCESS *) NMXP_MEM_MALLOC(sizeof(NMXP_DATA_PROCESS));
  memset(pd,0,sizeof(NMXP_DATA_PROCESS));


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

494 495
  nmxp_data_init(pd);
  outdata = (int32_t *) NMXP_MEM_MALLOC(MAX_OUTDATA*sizeof(int32_t));
496 497
  nmxp_channel_name = nmxp_chan_lookupName(pKey, channelList);

498
  if(nmxp_channel_name != NULL) {
499

500
  memcpy(&pTime, &buffer_data[4], 8);
501
  if ( swap ) { nmxp_data_swap_8b(&pTime); }
502 503 504 505 506 507 508

  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
509 510
  memcpy(outdata , (int32_t *) &buffer_data[20], length_data - 20);
  pDataPtr = outdata;
511 512 513 514 515 516 517

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

518
  if(!nmxp_chan_cpy_sta_chan(nmxp_channel_name, station_code, channel_code, network_code, location_code)) {
519 520
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not in STA.CHAN format: %s\n",
	    NMXP_LOG_STR(nmxp_channel_name));
521
  }
522 523

  pd->key = pKey;
524
  if(network_code[0] != 0) {
525
      strncpy(pd->network, network_code, NMXP_DATA_NETWORK_LENGTH);
526
  } else {
527
      strncpy(pd->network, network_code_default, NMXP_DATA_NETWORK_LENGTH);
528
  }
529
  if(station_code[0] != 0) {
530
      strncpy(pd->station, station_code, NMXP_DATA_STATION_LENGTH);
531
  }
532
  if(channel_code[0] != 0) {
533
      strncpy(pd->channel, channel_code, NMXP_DATA_CHANNEL_LENGTH);
534
  }
535 536 537 538 539
  if(location_code[0] != 0) {
      strncpy(pd->location, location_code, NMXP_DATA_LOCATION_LENGTH);
  } else {
      strncpy(pd->location, location_code_default, NMXP_DATA_LOCATION_LENGTH);
  }
540 541 542 543 544 545 546 547 548 549 550
  pd->packet_type = NMXP_MSG_DECOMPRESSED;
  pd->x0 = -1;
  pd->xn = -1;
  pd->x0n_significant = 0;
  pd->time = pTime;
  pd->nSamp = pNSamp;
  pd->pDataPtr = pDataPtr;
  pd->sampRate = pSampRate;



551 552 553
  /* TODO*/
  /* pd.oldest_seq_no = ;*/
  /* pd.seq_no = ;*/
554

555
  NMXP_MEM_FREE(nmxp_channel_name);
556 557
  } else {
      nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not found for key %d\n", pKey);
558 559
  }

560
  return pd;
561 562 563
}


564
NMXP_DATA_PROCESS *nmxp_processCompressedData(char* buffer_data, int length_data, NMXP_CHAN_LIST_NET *channelList, const char *network_code_default, const char *location_code_default)
565 566 567 568 569
{
    int32_t   pKey      = 0;
    double    pTime     = 0.0;
    int32_t   pNSamp    = 0;
    int32_t   pSampRate = 0;
570
    int32_t  *pDataPtr  = NULL;
571

572 573 574
    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];
575
    char location_code[NMXP_CHAN_MAX_SIZE_STR_PATTERN];
576

577 578 579
    NMXP_DATA_PROCESS *pd = NULL;


580

581
    int32_t nmx_rate_code_to_sample_rate[32] = {
582 583 584 585 586
	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};

587
	int32_t nmx_oldest_sequence_number;
588 589
	char nmx_hdr[25];
	unsigned char nmx_ptype;
590
	int32_t nmx_seconds;
591
	double nmx_seconds_double;
592 593
	int16_t nmx_ticks, nmx_instr_id;
	int32_t nmx_seqno;
594
	unsigned char nmx_sample_rate;
595 596
	int32_t nmx_x0;
	int32_t rate_code, chan_code, this_sample_rate;
597

598
	int32_t comp_bytecount;
599
	unsigned char *indata;
600 601
        int32_t * outdata = NULL;

602 603
	int32_t nout, i, k;
	int32_t prev_xn;
604 605
	const uint32_t high_scale = 4096 * 2048;
	const uint32_t high_scale_p = 4096 * 4096;
606

607 608
	char *nmxp_channel_name = NULL;

609 610 611
        pd= (NMXP_DATA_PROCESS *) NMXP_MEM_MALLOC(sizeof(NMXP_DATA_PROCESS));
        memset(pd,0,sizeof(NMXP_DATA_PROCESS));

612
	/* TOREMOVE int my_order = get_my_wordorder();*/
Matteo Quintiliani's avatar
Matteo Quintiliani committed
613
	int my_host_is_bigendian = nmxp_data_bigendianhost();
614
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "my_host_is_bigendian %d\n", my_host_is_bigendian);
615 616

	memcpy(&nmx_oldest_sequence_number, buffer_data, 4);
617 618 619
	if (my_host_is_bigendian) {
	    nmxp_data_swap_4b (&nmx_oldest_sequence_number);
	}
620
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Oldest sequence number = %d\n", nmx_oldest_sequence_number);
621 622 623 624 625 626

	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.   */
627
	    nmxp_log (NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Filler packet - discarding\n");
628
	    /*m continue;*/
629 630 631 632 633 634 635 636 637 638 639
	    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);

640
	if (my_host_is_bigendian) {
641
	    nmxp_data_swap_4b ((int32_t *)&nmx_seconds);
642 643 644 645
	    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);
646
	}
647 648 649 650 651 652 653

	/* 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;
	}

654 655 656 657 658
	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];

659 660 661
	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);
662

663 664
	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);
665

666 667 668
	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);
669
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "this_sample_rate   = %d\n", this_sample_rate);
670

671 672 673 674 675 676
	pKey = (nmx_instr_id << 16) | ( 1 << 8) | ( chan_code);

	pTime = nmx_seconds_double;

	pSampRate = this_sample_rate;

677 678
	nmxp_data_init(pd);
        outdata = (int32_t *) NMXP_MEM_MALLOC(MAX_OUTDATA*sizeof(int32_t));
679 680 681 682
	nmxp_channel_name = nmxp_chan_lookupName(pKey, channelList);

	if(nmxp_channel_name) {

683
	if(!nmxp_chan_cpy_sta_chan(nmxp_channel_name, station_code, channel_code, network_code, location_code)) {
684 685
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not in STA.CHAN format: %s\n",
		    NMXP_LOG_STR(nmxp_channel_name));
686 687
	}
  
688 689
	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));
690

691 692 693
	comp_bytecount = length_data-21;
	indata = (unsigned char *) buffer_data + 21;

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

696 697 698 699 700 701
	/* 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) {
702
		nmxp_log (NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "comp_bytecount = %d, i+17 = %d\n",
703 704 705 706
			comp_bytecount, i+17);
		exit(1);
	    }
	    if (nout+16 > MAX_OUTDATA)  {
707
		nmxp_log (NMXP_LOG_ERR,  NMXP_LOG_D_PACKETMAN, "Output buffer size too small\n");
708 709 710
		exit(1);
	    }
	    k = nmxp_data_unpack_bundle (outdata+nout,indata+i,&prev_xn);
711 712 713 714 715
	    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__);
716 717 718 719 720 721
	    if (k < 0) break;
	    nout += k;
	    /* prev_xn = outdata[nout-1]; */
	}
	nout--;

722
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Unpacked %d samples.\n", nout);
723 724 725 726

	pDataPtr = outdata;
	pNSamp = nout;

727
        pd->key = pKey;
728
	if(network_code[0] != 0) {
729
	    strncpy(pd->network, network_code, NMXP_DATA_NETWORK_LENGTH);
730
	} else {
731
	    strncpy(pd->network, network_code_default, NMXP_DATA_NETWORK_LENGTH);
732
	}
733
	if(station_code[0] != 0) {
734
	    strncpy(pd->station, station_code, NMXP_DATA_STATION_LENGTH);
735
	}
736
	if(channel_code[0] != 0) {
737
	    strncpy(pd->channel, channel_code, NMXP_DATA_CHANNEL_LENGTH);
738
	}
739 740 741 742 743
	if(location_code[0] != 0) {
	    strncpy(pd->location, location_code, NMXP_DATA_LOCATION_LENGTH);
	} else {
	    strncpy(pd->location, location_code_default, NMXP_DATA_LOCATION_LENGTH);
	}
744 745 746 747 748 749 750 751 752 753 754 755
	pd->packet_type = nmx_ptype;
	pd->x0 = nmx_x0;
	pd->xn = pDataPtr[nout];
	pd->x0n_significant = 1;
	pd->oldest_seq_no = nmx_oldest_sequence_number;
	pd->seq_no = nmx_seqno;
	pd->time = pTime;
	pd->nSamp = pNSamp;
	pd->pDataPtr = pDataPtr;
	pd->sampRate = pSampRate;

	NMXP_MEM_FREE(nmxp_channel_name);
756 757
	} else {
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not found for key %d\n", pKey);
758 759
	}

760 761
	return pd;

762 763
}

764 765 766

unsigned int nmxp_sleep(unsigned int sleep_time) {
#ifdef HAVE_WINDOWS_H
767
    Sleep(sleep_time * 1000);
768
    return 0;
769 770 771 772 773
#else
    return sleep(sleep_time);
#endif
}

Matteo Quintiliani's avatar
Matteo Quintiliani committed
774 775 776 777 778 779 780 781 782
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
}