nmxp.c 34.4 KB
Newer Older
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
1
2
/*! \file
 *
Matteo Quintiliani's avatar
Matteo Quintiliani committed
3
 * \brief Nanometrics Protocol Library
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
4
5
6
7
8
9
 *
 * Author:
 * 	Matteo Quintiliani
 * 	Istituto Nazionale di Geofisica e Vulcanologia - Italy
 *	quintiliani@ingv.it
 *
10
 * $Id: nmxp.c,v 1.87 2008-04-08 07:57:23 mtheo Exp $
Matteo Quintiliani's avatar
Matteo Quintiliani committed
11
 *
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
12
13
14
 */

#include "nmxp.h"
15
#include "nmxp_memory.h"
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
16

Matteo Quintiliani's avatar
Matteo Quintiliani committed
17
18
#include "config.h"

Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
19
#include <stdio.h>
20
#include <stdlib.h>
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
21
22
#include <string.h>

23
24
25
26
#ifdef HAVE_GETTIMEOFDAY
#include <sys/time.h>
#endif

Matteo Quintiliani's avatar
Matteo Quintiliani committed
27
28
29
30
#ifdef HAVE_ARPA_INET_H
#include <arpa/inet.h>
#endif

31
32
33
34
#ifdef HAVE_WINDOWS_H
#include "winsock2.h"
#endif

35
int nmxp_sendConnect(int isock) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
36
    return nmxp_sendMessage(isock, NMXP_MSG_CONNECT, NULL, 0);
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
37
38
}

Matteo Quintiliani's avatar
Matteo Quintiliani committed
39
40
int nmxp_sendTerminateSubscription(int isock, NMXP_SHUTDOWN_REASON reason, char *message) {
    return nmxp_sendMessage(isock, NMXP_MSG_TERMINATESUBSCRIPTION, message, ((message)? strlen(message)-1 : 0));
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
41
42
}

Matteo Quintiliani's avatar
Matteo Quintiliani committed
43
int nmxp_receiveChannelList(int isock, NMXP_CHAN_LIST **pchannelList) {
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
44
45
    int ret;
    int i;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
46
    int recv_errno;
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
47

Matteo Quintiliani's avatar
Matteo Quintiliani committed
48
    NMXP_MSG_SERVER type;
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
49
    void *buffer;
50
    int32_t length;
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
51
52
53

    *pchannelList = NULL;

Matteo Quintiliani's avatar
Matteo Quintiliani committed
54
    ret = nmxp_receiveMessage(isock, &type, &buffer, &length, 0, &recv_errno);
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
55

Matteo Quintiliani's avatar
Matteo Quintiliani committed
56
    if(type != NMXP_MSG_CHANNELLIST) {
57
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Type %d is not NMXP_MSG_CHANNELLIST!\n", type);
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
58
59
60
61
62
    } else {

	*pchannelList = buffer;
	(*pchannelList)->number = ntohl((*pchannelList)->number);

63
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CHANNEL, "number of channels %d\n", (*pchannelList)->number);
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
64
	
65
	/* TODO check*/
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
66
67
68

	for(i=0; i < (*pchannelList)->number; i++) {
	    (*pchannelList)->channel[i].key = ntohl((*pchannelList)->channel[i].key);
69
70
71
	    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CHANNEL, "%12d %s\n",
		    (*pchannelList)->channel[i].key,
		    NMXP_LOG_STR((*pchannelList)->channel[i].name));
Matteo Quintiliani's avatar
Start    
Matteo Quintiliani committed
72
73
74
75
76
77
78
79
	}

    }

    return ret;
}


80
int nmxp_sendAddTimeSeriesChannel_raw(int isock, NMXP_CHAN_LIST_NET *channelList, int32_t shortTermCompletion, int32_t out_format, NMXP_BUFFER_FLAG buffer_flag) {
81
    int ret;
82
    int32_t buffer_length = 16 + (4 * channelList->number); 
83
    char *buffer = NMXP_MEM_MALLOC(buffer_length);
84
    int32_t app, i, disp;
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109

    disp=0;

    app = htonl(channelList->number);
    memcpy(&buffer[disp], &app, 4);
    disp+=4;

    for(i=0; i < channelList->number; i++) {
	app = htonl(channelList->channel[i].key);
	memcpy(&buffer[disp], &app, 4);
	disp+=4;
    }
    
    app = htonl(shortTermCompletion);
    memcpy(&buffer[disp], &app, 4);
    disp+=4;

    app = htonl(out_format);
    memcpy(&buffer[disp], &app, 4);
    disp+=4;

    app = htonl(buffer_flag);
    memcpy(&buffer[disp], &app, 4);
    disp+=4;

Matteo Quintiliani's avatar
Matteo Quintiliani committed
110
    ret = nmxp_sendMessage(isock, NMXP_MSG_ADDTIMESERIESCHANNELS, buffer, buffer_length);
111
112

    if(buffer) {
113
	NMXP_MEM_FREE(buffer);
114
	buffer = NULL;
115
116
117
118
    }
    return ret;
}

119
#define MAX_LEN_S_CHANNELS 4096
120
int nmxp_sendAddTimeSeriesChannel(int isock, NMXP_CHAN_LIST_NET *channelList, int32_t shortTermCompletion, int32_t out_format, NMXP_BUFFER_FLAG buffer_flag, int n_channel, int n_usec, int flag_restart) {
121
122
123
124
    static int i = 0;
    static int first_time = 1;
    static struct timeval last_tp_now;

125
    char s_channels[MAX_LEN_S_CHANNELS];
126
127
128
129
130
    int j;
    int ret = 0;
    NMXP_CHAN_LIST_NET split_channelList;
    long diff_usec;
    struct timeval tp_now;
131
132
133
134
135
136
137
    double estimated_time = 0.0;

    if(n_usec == 0  &&  n_channel == 0) {
	n_channel = channelList->number;
    }

    estimated_time = (double) channelList->number * ( ((double) n_usec / 1000000.0) / (double) n_channel);
138

139
    if(flag_restart) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
140
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW,
141
142
143
144
145
		"Estimated time for channel requests: %d * (%d/%d) = %.3f sec.\n",
		channelList->number, n_usec / 1000, n_channel,
		estimated_time);


146
147
148
149
	first_time = 1;
	i = 0;
    }

Matteo Quintiliani's avatar
Matteo Quintiliani committed
150
    /* Check if requests could be satisfied within NMXP_MAX_MSCHAN_MSEC */
151
152
153
154
155
156
157
158
159
    if(estimated_time > ((double) NMXP_MAX_MSCHAN_MSEC / 1000.0)) {
	n_usec = ( (double) NMXP_MAX_MSCHAN_MSEC * 1000.0 ) * ( (double) n_channel / (double) channelList->number);
	estimated_time = (double) channelList->number * ( ((double) n_usec / 1000000.0) / (double) n_channel);
	if(flag_restart) {
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_ANY, "Estimated time exceeds. New values %d/%d and estimated time %.3f sec.\n",
		    n_usec/1000, n_channel, estimated_time);
	}
    }

160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
#ifdef HAVE_GETTIMEOFDAY
    gettimeofday(&tp_now, NULL);
#else
    TODO gettimeofday not found
#endif

    if(i <  channelList->number) {
	    if(first_time) {
		    diff_usec = n_usec + 1;
		    first_time = 0;
		    last_tp_now.tv_sec = 0;
		    last_tp_now.tv_usec = 0;
	    } else {
		    diff_usec = (tp_now.tv_sec - last_tp_now.tv_sec) * 1000000;
		    diff_usec += (tp_now.tv_usec - last_tp_now.tv_usec);
	    }
	    if(diff_usec >= n_usec) {
		    /* while(ret == 0  &&  i <  channelList->number) { */
		    split_channelList.number = 0;
		    while(split_channelList.number < n_channel  &&  i < channelList->number) {
			    split_channelList.channel[split_channelList.number].key = channelList->channel[i].key;
			    /* Not necessary, but it could help for debugging */
			    strcpy(split_channelList.channel[split_channelList.number].name, channelList->channel[i].name);
			    split_channelList.number++;
			    i++;
		    }
		    if(split_channelList.number > 0) {
187
			snprintf(s_channels, MAX_LEN_S_CHANNELS, "%.0f/%d chan %d of %d:",
188
				(double)diff_usec/1000.0, split_channelList.number, i, channelList->number);
189
			    for(j=0; j < split_channelList.number; j++) {
190
191
				strncat(s_channels, " ", MAX_LEN_S_CHANNELS);
				strncat(s_channels, NMXP_LOG_STR(split_channelList.channel[j].name), MAX_LEN_S_CHANNELS);
192
			    }
193
			    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "%s\n", s_channels);
194
195
196
197
198
199
200
201
202
203
204
			    ret = nmxp_sendAddTimeSeriesChannel_raw(isock, &split_channelList, shortTermCompletion, out_format, buffer_flag);
		    }
		    /* } */
		    last_tp_now.tv_sec = tp_now.tv_sec;
		    last_tp_now.tv_usec = tp_now.tv_usec;
	    }
    }

    return ret;
}

205

Matteo Quintiliani's avatar
Matteo Quintiliani committed
206
NMXP_DATA_PROCESS *nmxp_receiveData(int isock, NMXP_CHAN_LIST_NET *channelList, const char *network_code, int timeoutsec, int *recv_errno ) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
207
    NMXP_MSG_SERVER type;
208
    void *buffer = NULL;
209
    int32_t length;
210
    NMXP_DATA_PROCESS *pd = NULL;
211

Matteo Quintiliani's avatar
Matteo Quintiliani committed
212
    if(nmxp_receiveMessage(isock, &type, &buffer, &length, timeoutsec, recv_errno) == NMXP_SOCKET_OK) {
213
	if(type == NMXP_MSG_COMPRESSED) {
214
	    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Type %d is NMXP_MSG_COMPRESSED!\n", type);
215
	    pd = nmxp_processCompressedData(buffer, length, channelList, network_code);
216
	} else if(type == NMXP_MSG_DECOMPRESSED) {
217
	    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Type %d is NMXP_MSG_DECOMPRESSED!\n", type);
218
	    pd = nmxp_processDecompressedData(buffer, length, channelList, network_code);
219
	} else {
220
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Type %d is not NMXP_MSG_COMPRESSED or NMXP_MSG_DECOMPRESSED!\n", type);
221
	}
Matteo Quintiliani's avatar
Matteo Quintiliani committed
222
223
    }

224
    return pd;
225
226
}

227

228
int nmxp_sendConnectRequest(int isock, char *naqs_username, char *naqs_password, int32_t connection_time) {
229
    int ret;
230
    int i;
231
    char crc32buf[100];
232
233
    char *pcrc32buf = crc32buf;
    int crc32buf_length = 0;
234
    NMXP_CONNECT_REQUEST connectRequest;
235
    int naqs_username_length, naqs_password_length;
236
    int32_t protocol_version = 0;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
237
    char *pp = NULL;
238

239
240
    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CRC, "%s - %s\n",
	    NMXP_LOG_STR(naqs_username), NMXP_LOG_STR(naqs_password));
241

242
243
244
    naqs_username_length = (naqs_username)? strlen(naqs_username) : 0;
    naqs_password_length = (naqs_password)? strlen(naqs_password) : 0;

245
246
247
248
249
250
251
252
    for(i=0; i < 100; i++) {
	crc32buf[i] = 0;
    }

    for(i=0; i<12; i++) {
	connectRequest.username[i] = 0;
    }
    if(naqs_username_length != 0) {
253
	strcpy(connectRequest.username, naqs_username);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
254
    }
255
256
257
258

    connectRequest.version = protocol_version;
    connectRequest.connection_time = connection_time;

259
260
261
    if(!nmxp_data_bigendianhost()) {
	nmxp_data_swap_4b ((int32_t *)&connection_time);
    }
262

263
    if(naqs_username_length == 0  &&  naqs_password_length == 0 ) {
264
265
	/* sprintf(crc32buf, "%d%d", protocol_version, connection_time); */
	/* TODO */
266
    } else if(naqs_username_length != 0  &&  naqs_password_length != 0 ) {
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
	/* sprintf(crc32buf, "%s%d%d%s", naqs_username, protocol_version,
		connection_time, naqs_password); */

	memcpy(pcrc32buf, naqs_username, naqs_username_length);
	crc32buf_length += naqs_username_length;
	pcrc32buf = crc32buf + crc32buf_length;

	memcpy(pcrc32buf, &(protocol_version), sizeof(protocol_version));
	crc32buf_length += sizeof(protocol_version);
	pcrc32buf = crc32buf + crc32buf_length;

	memcpy(pcrc32buf, &(connection_time), sizeof(connection_time));
	crc32buf_length += sizeof(connection_time);
	pcrc32buf = crc32buf + crc32buf_length;

	memcpy(pcrc32buf, naqs_password, naqs_password_length);
	crc32buf_length += naqs_password_length;
	pcrc32buf = crc32buf + crc32buf_length;

286
    } else if(naqs_username_length != 0 ) {
287
288
	/* sprintf(crc32buf, "%s%d%d", naqs_username, protocol_version, connection_time); */
	/* TODO */
289
    } else if(naqs_password_length != 0 ) {
290
291
	/* sprintf(crc32buf, "%d%d%s", protocol_version, connection_time, naqs_password); */
	/* TODO */
292
    }
293
294
295
    connectRequest.version = htonl(connectRequest.version);
    connectRequest.connection_time = htonl(connectRequest.connection_time);
    connectRequest.crc32 = htonl(crc32(0L, crc32buf, crc32buf_length));
296
297
298
299

    ret = nmxp_sendMessage(isock, NMXP_MSG_CONNECTREQUEST, &connectRequest, sizeof(NMXP_CONNECT_REQUEST));

    if(ret == NMXP_SOCKET_OK) {
300
301
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CRC, "Send a ConnectRequest crc32buf length %d, crc32 = %d\n",
		crc32buf_length, connectRequest.crc32);
302
303
304
	for(i=0; i < crc32buf_length; i++) {
	    nmxp_log(NMXP_LOG_NORM_NO, NMXP_LOG_D_CRC, "%d ", crc32buf[i]);
	}
Matteo Quintiliani's avatar
Matteo Quintiliani committed
305
	pp = (char *) &connectRequest.crc32;
306
307
308
309
	for(i=0; i < sizeof(connectRequest.crc32); i++) {
	    nmxp_log(NMXP_LOG_NORM_NO, NMXP_LOG_D_CRC, "%d ", pp[i]);
	}
	nmxp_log(NMXP_LOG_NORM_NO, NMXP_LOG_D_CRC, "\n");
310
    } else {
311
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CRC, "Send a ConnectRequest.\n");
312
313
314
315
316
317
    }

    return ret;
}


318
int nmxp_readConnectionTime(int isock, int32_t *connection_time) {
319
    int ret;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
320
321
    int recv_errno;
    ret = nmxp_recv_ctrl(isock, connection_time, sizeof(int32_t), 0, &recv_errno);
322
    *connection_time = ntohl(*connection_time);
323
    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "Read connection time from socket %d.\n", *connection_time);
324
    if(ret != NMXP_SOCKET_OK) {
325
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Read connection time from socket.\n");
326
327
328
329
330
    }
    return ret;
}


331
int nmxp_waitReady(int isock) {
332
333
    int times = 0;
    int rc = NMXP_SOCKET_OK;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
334
335
336
    int32_t signature;
    int32_t type = 0;
    int32_t length;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
337
    int recv_errno;
338
339

    while(rc == NMXP_SOCKET_OK  &&  type != NMXP_MSG_READY) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
340
	rc = nmxp_recv_ctrl(isock, &signature, sizeof(signature), 0, &recv_errno);
341
342
343
	if(rc != NMXP_SOCKET_OK) return rc;
	signature = ntohl(signature);
	if(signature == 0) {
344
	    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "signature is equal to zero. receive again.\n");
Matteo Quintiliani's avatar
Matteo Quintiliani committed
345
	    rc = nmxp_recv_ctrl(isock, &signature, sizeof(signature), 0, &recv_errno);
346
347
348
	    signature = ntohl(signature);
	}
	if(signature != NMX_SIGNATURE) {
349
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "signature is not valid. signature = %d\n", signature);
350
351
352
353
	    if(signature == 200) {
		    int32_t err_length;
		    int32_t err_reason;
		    char err_buff[200];
Matteo Quintiliani's avatar
Matteo Quintiliani committed
354
		    rc = nmxp_recv_ctrl(isock, &err_length, sizeof(err_length), 0, &recv_errno);
355
		    err_length = ntohl(err_length);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
356
		    rc = nmxp_recv_ctrl(isock, &err_reason, sizeof(err_reason), 0, &recv_errno);
357
358
		    err_reason = ntohl(err_reason);
		    if(err_length > 4) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
359
			    rc = nmxp_recv_ctrl(isock, err_buff, err_length-4, 0, &recv_errno);
360
361
			    err_buff[err_length] = 0;
		    }
362
		    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "TerminateMessage from Server: %s (%d).\n",
363
			    NMXP_LOG_STR(err_buff), err_reason);
364
	    }
365
366
367
	    return NMXP_SOCKET_ERROR;
	}

Matteo Quintiliani's avatar
Matteo Quintiliani committed
368
	rc = nmxp_recv_ctrl(isock, &type, sizeof(type), 0, &recv_errno);
369
370
371
	if(rc != NMXP_SOCKET_OK) return rc;
	type = ntohl(type);
	if(type != NMXP_MSG_READY) {
372
	    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "type is not READY. type = %d\n", type);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
373
	    rc = nmxp_recv_ctrl(isock, &length, sizeof(length), 0, &recv_errno);
374
375
376
	    if(rc != NMXP_SOCKET_OK) return rc;
	    length = ntohl(length);
	    if(length > 0) {
377
378
		if(type == NMXP_MSG_TERMINATESUBSCRIPTION) {
		    char *str_msg = NULL;
379
		    char *buf_app = (char *) NMXP_MEM_MALLOC(sizeof(char) * length);
380
381
382
383
384
385
386
387
388
389
		    int32_t reason;
		    memcpy(&reason, buf_app, sizeof(reason));
		    reason = ntohl(reason);
		    str_msg = buf_app + sizeof(reason);
		    rc = nmxp_recv_ctrl(isock, buf_app, length, 0, &recv_errno);
		    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);
		    if(buf_app) {
390
			NMXP_MEM_FREE(buf_app);
391
			buf_app = NULL;
392
		    }
393
394
395
		    /* Close the socket*/
		    nmxp_closeSocket(isock);
		    exit(-1);
396
		} else if(length == 4) {
397
		    int32_t app;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
398
		    rc = nmxp_recv_ctrl(isock, &app, length, 0, &recv_errno);
399
400
		    if(rc != NMXP_SOCKET_OK) return rc;
		    app = ntohl(app);
401
		    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "value = %d\n", app);
402
		} else {
403
		    char *buf_app = (char *) NMXP_MEM_MALLOC(sizeof(char) * length);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
404
		    rc = nmxp_recv_ctrl(isock, buf_app, length, 0, &recv_errno);
405
		    if(buf_app) {
406
			NMXP_MEM_FREE(buf_app);
407
			buf_app = NULL;
408
409
410
411
		    }
		}
	    }
	} else {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
412
	    rc = nmxp_recv_ctrl(isock, &length, sizeof(length), 0, &recv_errno);
413
414
415
	    if(rc != NMXP_SOCKET_OK) return rc;
	    length = ntohl(length);
	    if(length != 0) {
416
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "length is not equal to zero. length = %d\n", length);
417
418
419
420
421
422
		return NMXP_SOCKET_ERROR;
	    }
	}

	times++;
	if(times > 10) {
423
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "waiting_ready_message. times > 10\n");
424
425
426
427
428
429
430
431
432
	    rc = NMXP_SOCKET_ERROR;
	}

    }

    return rc;
}


433
int nmxp_sendDataRequest(int isock, int32_t key, int32_t start_time, int32_t end_time) {
434
435
436
437
438
439
440
441
442
443
    int ret;
    NMXP_DATA_REQUEST dataRequest;

    dataRequest.chan_key = htonl(key);
    dataRequest.start_time = htonl(start_time);
    dataRequest.end_time = htonl(end_time);

    ret = nmxp_sendMessage(isock, NMXP_MSG_DATAREQUEST, &dataRequest, sizeof(dataRequest));

    if(ret != NMXP_SOCKET_OK) {
444
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Send a Request message\n");
445
446
447
448
449
450
    }

    return ret;
}


451
452
453
NMXP_CHAN_LIST *nmxp_getAvailableChannelList(char * hostname, int portnum, NMXP_DATATYPE datatype) {
    int naqssock;
    NMXP_CHAN_LIST *channelList = NULL, *channelList_subset = NULL;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
454
    /* int i; */
455

456
    /* 1. Open a socket*/
457
458
459
460
    naqssock = nmxp_openSocket(hostname, portnum);

    if(naqssock != NMXP_SOCKET_ERROR) {

461
	/* 2. Send a Connect*/
462
463
	if(nmxp_sendConnect(naqssock) == NMXP_SOCKET_OK) {

464
	    /* 3. Receive ChannelList*/
465
466
467
	     if(nmxp_receiveChannelList(naqssock, &channelList) == NMXP_SOCKET_OK) {

		 channelList_subset = nmxp_chan_getType(channelList, datatype);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
468
		 nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CHANNEL, "%d / %d are DataType channel.\n", channelList_subset->number, channelList->number);
469

470
		 /* nmxp_chan_sortByKey(channelList_subset);*/
471
472
		 nmxp_chan_sortByName(channelList_subset);

Matteo Quintiliani's avatar
Matteo Quintiliani committed
473
		 /*
474
		 for(i=0; i < channelList_subset->number; i++) {
475
476
477
		     nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CHANNEL, "%12d %s\n",
			     channelList_subset->channel[i].key,
			     NMXP_LOG_STR(channelList_subset->channel[i].name));
478
		 }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
479
		 */
480

481
		 /* 4. Send a Request Pending (optional)*/
482

483
		 /* 5. Send AddChannels*/
484

485
		 /* 6. Repeat until finished: receive and handle packets*/
486

487
		 /* 7. Send Terminate Subscription*/
488
489
490
		 nmxp_sendTerminateSubscription(naqssock, NMXP_SHUTDOWN_NORMAL, "Good Bye!");

	     } else {
491
		 nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CHANNEL, "Error on receiveChannelList()\n");
492
493
	     }
	} else {
494
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CHANNEL, "Error on sendConnect()\n");
495
496
	}

497
	/* 8. Close the socket*/
498
499
500
501
	nmxp_closeSocket(naqssock);
    }

    if(channelList) {
502
	NMXP_MEM_FREE(channelList);
503
	channelList = NULL;
504
505
506
507
508
    }

    return channelList_subset;
}

509

510
NMXP_META_CHAN_LIST *nmxp_getMetaChannelList(char * hostname, int portnum, NMXP_DATATYPE datatype, int flag_request_channelinfo, char *datas_username, char *datas_password, NMXP_CHAN_LIST **pchannelList) {
511
512
    int naqssock;
    NMXP_CHAN_PRECISLIST *precisChannelList = NULL;
513
514
515
516
    NMXP_CHAN_LIST *channelList = NULL;
    NMXP_META_CHAN_LIST *chan_list = NULL;
    NMXP_META_CHAN_LIST *iter = NULL;
    int i = 0;
517
    int32_t connection_time;
518
    int ret_sock;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
519
    int recv_errno;
520
    
521
522
    NMXP_MSG_SERVER type;
    void *buffer = NULL;
523
    int32_t length;
524
525
526
    NMXP_PRECISLISTREQUEST precisListRequestBody;
    NMXP_CHANNELINFOREQUEST channelInfoRequestBody;
    NMXP_CHANNELINFORESPONSE *channelInfo;
527
528
529
530
531
532
533

    char str_start[200], str_end[200];
    str_start[0] = 0;
    str_end[0] = 0;
    
    /* DAP Step 1: Open a socket */
    if( (naqssock = nmxp_openSocket(hostname, portnum)) == NMXP_SOCKET_ERROR) {
534
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Error opening socket!\n");
535
536
537
538
539
	return NULL;
    }

    /* DAP Step 2: Read connection time */
    if(nmxp_readConnectionTime(naqssock, &connection_time) != NMXP_SOCKET_OK) {
540
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Error reading connection time from server!\n");
541
542
543
544
545
	return NULL;
    }

    /* DAP Step 3: Send a ConnectRequest */
    if(nmxp_sendConnectRequest(naqssock, datas_username, datas_password, connection_time) != NMXP_SOCKET_OK) {
546
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Error sending connect request!\n");
547
548
549
550
551
	return NULL;
    }

    /* DAP Step 4: Wait for a Ready message */
    if(nmxp_waitReady(naqssock) != NMXP_SOCKET_OK) {
552
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Error waiting Ready message!\n");
553
554
555
	return NULL;
    }

556
557
558
559
560
561



    /* DAP Step 5: Send Data Request */
    nmxp_sendHeader(naqssock, NMXP_MSG_CHANNELLISTREQUEST, 0);
    /* DAP Step 6: Receive Data until receiving a Ready message */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
562
    ret_sock = nmxp_receiveMessage(naqssock, &type, &buffer, &length, 0, &recv_errno);
563
    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_EXTRA, "ret_sock = %d, type = %d, length = %d\n", ret_sock, type, length);
564
565
566
567
568
569
570
571

    while(ret_sock == NMXP_SOCKET_OK   &&    type != NMXP_MSG_READY) {
	channelList = buffer;

	channelList->number = ntohl(channelList->number);

	for(i = 0; i < channelList->number; i++) {
	    channelList->channel[i].key = ntohl(channelList->channel[i].key);
572
	    if(getDataTypeFromKey(channelList->channel[i].key) == datatype) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
573
		nmxp_meta_chan_add(&chan_list, channelList->channel[i].key, channelList->channel[i].name, 0, 0, NULL, NMXP_META_SORT_NAME);
574
	    }
575
576
577
	}

	/* Receive Message */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
578
	ret_sock = nmxp_receiveMessage(naqssock, &type, &buffer, &length, 0, &recv_errno);
579
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_EXTRA, "ret_sock = %d, type = %d, length = %d\n", ret_sock, type, length);
580
581
    }

582
    *pchannelList = channelList;
583

584
585
586
587
    /* DAP Step 5: Send Data Request */
    precisListRequestBody.instr_id = htonl(-1);
    precisListRequestBody.datatype = htonl(NMXP_DATA_TIMESERIES);
    precisListRequestBody.type_of_channel = htonl(-1);
588
    nmxp_sendMessage(naqssock, NMXP_MSG_PRECISLISTREQUEST, &precisListRequestBody, sizeof(NMXP_PRECISLISTREQUEST));
589

590
591

    /* DAP Step 6: Receive Data until receiving a Ready message */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
592
    ret_sock = nmxp_receiveMessage(naqssock, &type, &buffer, &length, 0, &recv_errno);
593
    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_EXTRA, "ret_sock = %d, type = %d, length = %d\n", ret_sock, type, length);
594
595
596
597
598
599
600
601
602
603
604
605
606

    while(ret_sock == NMXP_SOCKET_OK   &&    type != NMXP_MSG_READY) {
	precisChannelList = buffer;

	precisChannelList->number = ntohl(precisChannelList->number);
	for(i = 0; i < precisChannelList->number; i++) {
	    precisChannelList->channel[i].key = ntohl(precisChannelList->channel[i].key);
	    precisChannelList->channel[i].start_time = ntohl(precisChannelList->channel[i].start_time);
	    precisChannelList->channel[i].end_time = ntohl(precisChannelList->channel[i].end_time);

	    nmxp_data_to_str(str_start, precisChannelList->channel[i].start_time);
	    nmxp_data_to_str(str_end, precisChannelList->channel[i].end_time);

607
	    if(!nmxp_meta_chan_set_times(chan_list, precisChannelList->channel[i].key, precisChannelList->channel[i].start_time, precisChannelList->channel[i].end_time)) {
608
609
610
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CHANNEL, "Key %d not found for %s!\n",
			precisChannelList->channel[i].key,
			NMXP_LOG_STR(precisChannelList->channel[i].name));
611
612
613
	    }

	    /*
614
	    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_ANY, "%12d %12s %10d %10d %20s %20s\n",
615
		    precisChannelList->channel[i].key, NMXP_LOG_STR(precisChannelList->channel[i].name),
616
		    precisChannelList->channel[i].start_time, precisChannelList->channel[i].end_time,
617
		    NMXP_LOG_STR(str_start), NMXP_LOG_STR(str_end));
618
		    */
619
620
621
	}

	/* Receive Message */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
622
	ret_sock = nmxp_receiveMessage(naqssock, &type, &buffer, &length, 0, &recv_errno);
623
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_EXTRA, "ret_sock = %d, type = %d, length = %d\n", ret_sock, type, length);
624
625
626
    }


627
628
    if(flag_request_channelinfo) {
	for(iter = chan_list; iter != NULL; iter = iter->next) {
629

630
631
632
633
	    if(getChannelNumberFromKey(iter->key) == 0) {
		/* DAP Step 5: Send Data Request */
		channelInfoRequestBody.key = htonl(iter->key);
		channelInfoRequestBody.ignored = htonl(0);
634
		nmxp_sendMessage(naqssock, NMXP_MSG_CHANNELINFOREQUEST, &channelInfoRequestBody, sizeof(NMXP_CHANNELINFOREQUEST));
635

636
		/* DAP Step 6: Receive Data until receiving a Ready message */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
637
		ret_sock = nmxp_receiveMessage(naqssock, &type, &buffer, &length, 0, &recv_errno);
638
		nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_EXTRA, "ret_sock = %d, type = %d, length = %d\n", ret_sock, type, length);
639

640
641
642
		while(ret_sock == NMXP_SOCKET_OK   &&    type != NMXP_MSG_READY) {
		    channelInfo = buffer;
		    channelInfo->key = ntohl(channelInfo->key);
643

644
		    if(!nmxp_meta_chan_set_network(chan_list, channelInfo->key, channelInfo->network)) {
645
646
			nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CHANNEL, "Key %d (%d) not found for %s!\n",
				iter->key, channelInfo->key, NMXP_LOG_STR(iter->name));
647
648
		    }
		    /* Receive Message */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
649
		    ret_sock = nmxp_receiveMessage(naqssock, &type, &buffer, &length, 0, &recv_errno);
650
		    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_EXTRA, "ret_sock = %d, type = %d, length = %d\n", ret_sock, type, length);
651
652
653
		}
	    }
	}
654
655
656
657
658
659
660
661
662
663
664
665
    }



    /* DAP Step 7: Repeat steps 5 and 6 for each data request */

    /* DAP Step 8: Send a Terminate message (optional) */
    nmxp_sendTerminateSubscription(naqssock, NMXP_SHUTDOWN_NORMAL, "Bye!");

    /* DAP Step 9: Close the socket */
    nmxp_closeSocket(naqssock);

666
    return chan_list;
667
668
}

669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684

int nmxp_raw_stream_seq_no_compare(const void *a, const void *b)
{       
    int ret = 0;
    NMXP_DATA_PROCESS **ppa = (NMXP_DATA_PROCESS **) a;
    NMXP_DATA_PROCESS **ppb = (NMXP_DATA_PROCESS **) b;
    NMXP_DATA_PROCESS *pa = *ppa;
    NMXP_DATA_PROCESS *pb = *ppb;

    if(pa && pb) {
	if(pa->seq_no > pb->seq_no) {
	    ret = 1;
	} else if (pa->seq_no < pb->seq_no) {
	    ret = -1;
	}
    } else {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
685
686
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_ANY,
		"nmxp_raw_stream_seq_no_compare() pa %s NULL and pb %s NULL\n", (pa)? "!=" : "=", (pb)? "!=" : "=");
687
688
689
690
691
    }

    return ret;
}

Matteo Quintiliani's avatar
Matteo Quintiliani committed
692
void nmxp_raw_stream_init(NMXP_RAW_STREAM_DATA *raw_stream_buffer, int32_t max_tolerable_latency, int timeoutrecv) {
693
694
695
696
    int j;

    raw_stream_buffer->last_seq_no_sent = -1;
    raw_stream_buffer->last_sample_time = -1.0;
697
    raw_stream_buffer->last_latency = 0.0;
698
    /* TODO 
Matteo Quintiliani's avatar
Matteo Quintiliani committed
699
     * Suppose a packet can contain 1/4 secs of data (that is, minimum packet length is 0.25 secs of data) */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
700
701
    raw_stream_buffer->max_tolerable_latency = max_tolerable_latency;
    raw_stream_buffer->max_pdlist_items = max_tolerable_latency * 4;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
702
    raw_stream_buffer->timeoutrecv = timeoutrecv;
703
    raw_stream_buffer->n_pdlist = 0;
704
    raw_stream_buffer->pdlist = (NMXP_DATA_PROCESS **) NMXP_MEM_MALLOC(raw_stream_buffer->max_pdlist_items * sizeof(NMXP_DATA_PROCESS *));
705
706
707
708
709
710
711
712
713
    for(j=0; j<raw_stream_buffer->max_pdlist_items; j++) {
	raw_stream_buffer->pdlist[j] = NULL;
    }

}


void nmxp_raw_stream_free(NMXP_RAW_STREAM_DATA *raw_stream_buffer) {
    int j;
714
715
    if(raw_stream_buffer) {
	if(raw_stream_buffer->pdlist) {
716
	    for(j=0; j<raw_stream_buffer->n_pdlist; j++) {
717
718
		if(raw_stream_buffer->pdlist[j]) {
		    if(raw_stream_buffer->pdlist[j]->buffer) {
719
			NMXP_MEM_FREE(raw_stream_buffer->pdlist[j]->buffer);
720
721
722
			raw_stream_buffer->pdlist[j]->buffer = NULL;
		    }
		    if(raw_stream_buffer->pdlist[j]->pDataPtr) {
723
			NMXP_MEM_FREE(raw_stream_buffer->pdlist[j]->pDataPtr);
724
725
			raw_stream_buffer->pdlist[j]->pDataPtr = NULL;
		    }
726
		    NMXP_MEM_FREE(raw_stream_buffer->pdlist[j]);
727
		    raw_stream_buffer->pdlist[j] = NULL;
728
		}
729
	    }
730
	    NMXP_MEM_FREE(raw_stream_buffer->pdlist);
731
	    raw_stream_buffer->pdlist = NULL;
732
733
734
735
736
	}
    }
}


Matteo Quintiliani's avatar
Matteo Quintiliani committed
737
int nmxp_raw_stream_manage(NMXP_RAW_STREAM_DATA *p, NMXP_DATA_PROCESS *a_pd, int (*p_func_pd[NMXP_MAX_FUNC_PD]) (NMXP_DATA_PROCESS *), int n_func_pd) {
738
739
740
741
742
743
744
745
746
    int ret = 0;
    int send_again = 1;
    int seq_no_diff;
    double time_diff;
    double latency = 0.0;
    int j=0, k=0;
    int i_func_pd;
    char str_time[200];
    NMXP_DATA_PROCESS *pd = NULL;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
747
    int y, w;
748
749
    int count_null_element = 0;
    char netstachan[100];
750

751
    /* Allocate pd copy value from a_pd */
752
    if(a_pd) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
753
754
755
756
757
758
759
	/*
	if(a_pd->packet_type == 33 || a_pd->packet_type == 97) {
	    nmxp_data_log(a_pd);
	}
	*/

	/* Allocate memory for pd and copy a_pd */
760
	pd = (NMXP_DATA_PROCESS *) NMXP_MEM_MALLOC(sizeof(NMXP_DATA_PROCESS));
Matteo Quintiliani's avatar
Matteo Quintiliani committed
761
762
	memcpy(pd, a_pd, sizeof(NMXP_DATA_PROCESS));
	if(a_pd->length > 0) {
763
	    pd->buffer = NMXP_MEM_MALLOC(pd->length);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
764
765
766
767
768
	    memcpy(pd->buffer, a_pd->buffer, a_pd->length);
	} else {
	    pd->buffer = NULL;
	}
	if(a_pd->nSamp *  sizeof(int) > 0) {
769
	    pd->pDataPtr = (int *) NMXP_MEM_MALLOC(a_pd->nSamp * sizeof(int));
Matteo Quintiliani's avatar
Matteo Quintiliani committed
770
771
772
773
	    memcpy(pd->pDataPtr, a_pd->pDataPtr, a_pd->nSamp * sizeof(int));
	} else {
	    pd->pDataPtr = NULL;
	}
774
775
    } else {
	nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_RAWSTREAM,
776
		"nmxp_raw_stream_manage() passing NMXP_DATA_PROCESS pointer equal to NULL\n");
777
    }
778
    /* From here, use only pd */
779
780

    /* First time */
781
    if(p->last_seq_no_sent == -1  &&  pd) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
782
783
784
	if(p->timeoutrecv == 0) {
	    p->last_seq_no_sent = pd->seq_no - 1;
	    p->last_sample_time = pd->time;
785
	    p->last_latency = nmxp_data_latency(pd);;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
786
787
	} else {
	    p->last_seq_no_sent = 0;
788
789
	    p->last_sample_time = 0.0;
	    p->last_latency = 0.0;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
790
	}
791
	nmxp_data_to_str(str_time, pd->time);
792
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_RAWSTREAM,
793
794
795
796
797
		"%s.%s.%s [%d, %d] (%s + %.2f sec.) * First time nmxp_raw_stream_manage() * last_seq_no_sent=%d  last_sample_time=%.2f\n",
		NMXP_LOG_STR(pd->network), NMXP_LOG_STR(pd->station), NMXP_LOG_STR(pd->channel),
		pd->packet_type, pd->seq_no,
		NMXP_LOG_STR(str_time), (double) pd->nSamp / (double) pd->sampRate,
		 p->last_seq_no_sent, p->last_sample_time);
798
    }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
799

800
801
802
    if(p->n_pdlist > 0) {
	latency = nmxp_data_latency(p->pdlist[0]);
    }
803

804
805
806
    /* Add pd and sort array, in case handle the first item */
    if( (p->n_pdlist >= p->max_pdlist_items || latency >= p->max_tolerable_latency)
	    && p->timeoutrecv <= 0 ) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
807

808
	/* Supposing p->pdlist is ordered, handle the first item and over write it.  */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
809
810
811
812
813
814
	if(p->n_pdlist > 0) {
	    seq_no_diff = p->pdlist[0]->seq_no - p->last_seq_no_sent;
	    time_diff = p->pdlist[0]->time - p->last_sample_time;
	    latency = nmxp_data_latency(p->pdlist[0]);
	    nmxp_data_to_str(str_time, p->pdlist[0]->time);
	    if( seq_no_diff > 0) {
815
		nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_RAWSTREAM,
816
817
818
819
820
821
			"%s.%s.%s [%d, %d] (%s + %.2f sec.) * Force handling packet * n_pdlist=%d  seq_no_diff=%d  time_diff=%.2fs  lat. %.1fs!\n",
			NMXP_LOG_STR(p->pdlist[0]->network), NMXP_LOG_STR(p->pdlist[0]->station), NMXP_LOG_STR(p->pdlist[0]->channel),
			p->pdlist[0]->packet_type, p->pdlist[0]->seq_no,
			NMXP_LOG_STR(str_time), (double) p->pdlist[0]->nSamp / (double) p->pdlist[0]->sampRate,
			p->n_pdlist,
			seq_no_diff, time_diff, latency);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
822
823
824
825
826
		for(i_func_pd=0; i_func_pd<n_func_pd; i_func_pd++) {
		    (*p_func_pd[i_func_pd])(p->pdlist[0]);
		}
		p->last_seq_no_sent = (p->pdlist[0]->seq_no);
		p->last_sample_time = (p->pdlist[0]->time + ((double) p->pdlist[0]->nSamp / (double) p->pdlist[0]->sampRate ));
827
		p->last_latency = nmxp_data_latency(p->pdlist[0]);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
828
829
	    } else {
		/* It should not occur */
830
831
832
833
834
835
836
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_RAWSTREAM,
			"%s.%s.%s [%d, %d] (%s + %.2f sec.) * SHOULD NOT OCCUR packet discarded * n_pdlist=%d  seq_no_diff=%d  time_diff=%.2fs  lat. %.1fs!\n",
			NMXP_LOG_STR(p->pdlist[0]->network), NMXP_LOG_STR(p->pdlist[0]->station), NMXP_LOG_STR(p->pdlist[0]->channel),
			p->pdlist[0]->packet_type, p->pdlist[0]->seq_no,
			NMXP_LOG_STR(str_time), (double) p->pdlist[0]->nSamp / (double) p->pdlist[0]->sampRate,
			p->n_pdlist,
			seq_no_diff, time_diff, latency);
837
838
	    }

Matteo Quintiliani's avatar
Matteo Quintiliani committed
839
840
	    /* Free handled packet */
	    if(p->pdlist[0]->buffer) {
841
		NMXP_MEM_FREE(p->pdlist[0]->buffer);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
842
843
844
		p->pdlist[0]->buffer = NULL;
	    }
	    if(p->pdlist[0]->pDataPtr) {
845
		NMXP_MEM_FREE(p->pdlist[0]->pDataPtr);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
846
847
848
		p->pdlist[0]->pDataPtr = NULL;
	    }
	    if(p->pdlist[0]) {
849
		NMXP_MEM_FREE(p->pdlist[0]);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
850
851
852
853
854
855
856
857
858
859
		p->pdlist[0] = NULL;
	    }
	    if(pd) {
		p->pdlist[0] = pd;
	    }
	} else {
	    if(pd) {
		p->n_pdlist = 1;
		p->pdlist[0] = pd;
	    }
860
861
	}
    } else {
862
863
864
	if(pd) {
	    p->pdlist[p->n_pdlist++] = pd;
	}
865
    }
866
867
868
869
870

    /* Check if some element in pdlist is NULL and remove it */
    y=0;
    while(y < p->n_pdlist) {
	if(p->pdlist[y] == NULL) {
871
	    count_null_element++;
872
873
874
875
876
877
878
879
880
881
	    /* Shift array */
	    for(w=y+1; w < p->n_pdlist; w++) {
		p->pdlist[w-1] = p->pdlist[w];
	     }
	    p->n_pdlist--;
	} else {
	    y++;
	}
    }

882
883
884
885
886
887
888
889
890
891
892
893
894
895
    if(count_null_element > 0) {
	if(p->n_pdlist > 0) {
	    snprintf(netstachan, 100, "%s.%s.%s",
		    NMXP_LOG_STR(p->pdlist[0]->network),
		    NMXP_LOG_STR(p->pdlist[0]->station),
		    NMXP_LOG_STR(p->pdlist[0]->channel));
	} else {
	    strncpy(netstachan, "Unknown", 100);
	}
	nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_ANY,
		"nmxp_raw_stream_manage() %d NULL elements in pdlist for %s.\n",
		count_null_element, netstachan);
    }

896
    /* Sort array */
897
898
    qsort(p->pdlist, p->n_pdlist, sizeof(NMXP_DATA_PROCESS *), nmxp_raw_stream_seq_no_compare);

899
    /* TODO Check for packet duplication in pd->pdlist*/
900
901

    /* Print array, only for debugging */
902
    /*
903
904
905
    if(p->n_pdlist > 1) {
	int y = 0;
	for(y=0; y < p->n_pdlist; y++) {
906
907
908
909
910
911
912
	    nmxp_data_to_str(str_time, p->pdlist[y]->time);
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_RAWSTREAM,
		    "%s.%s.%s [%d, %d] (%s + %.2f sec.) * %02d n_pdlist=%d\n",
		    NMXP_LOG_STR(p->pdlist[y]->network), NMXP_LOG_STR(p->pdlist[y]->station), NMXP_LOG_STR(p->pdlist[y]->channel),
		    p->pdlist[y]->packet_type, p->pdlist[y]->seq_no,
		    NMXP_LOG_STR(str_time), (double) p->pdlist[y]->nSamp / (double) p->pdlist[y]->sampRate,
		    y, p->n_pdlist);
913
914
	}
    }
915
    */
916

917
    /* Condition for time-out (pd is NULL) */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
918
    if(!pd && p->n_pdlist > 0) {
919
920
921
922
923
924
925
926
927
928
929
	/* Log before changing values */
	nmxp_data_to_str(str_time, p->pdlist[0]->time);
	nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_RAWSTREAM,
		"%s.%s.%s [%d, %d] (%s + %.2f sec.) * pd is NULL and n_pdlist = %d > 0 *  last_seq_no_sent=%d, last_sample_time=%.2f\n",
		NMXP_LOG_STR(p->pdlist[0]->network), NMXP_LOG_STR(p->pdlist[0]->station), NMXP_LOG_STR(p->pdlist[0]->channel),
		p->pdlist[0]->packet_type, p->pdlist[0]->seq_no,
		NMXP_LOG_STR(str_time), (double) p->pdlist[0]->nSamp / (double) p->pdlist[0]->sampRate,
		p->n_pdlist,
		p->last_seq_no_sent, p->last_sample_time);

	/* Changing values */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
930
931
	p->last_seq_no_sent = p->pdlist[0]->seq_no - 1;
	p->last_sample_time = p->pdlist[0]->time;
932
	p->last_latency = nmxp_data_latency(p->pdlist[0]);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
933
    }
934
935
936
937
938
939
940
941
942
943
944

    /* Manage array and execute func_pd() */
    j=0;
    send_again = 1;
    while(send_again  &&  j < p->n_pdlist) {
	send_again = 0;
	seq_no_diff = p->pdlist[j]->seq_no - p->last_seq_no_sent;
	time_diff = p->pdlist[j]->time - p->last_sample_time;
	latency = nmxp_data_latency(p->pdlist[j]);
	nmxp_data_to_str(str_time, p->pdlist[j]->time);
	if(seq_no_diff <= 0) {
945
	    /* Duplicated packets: Discarded */
946
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_RAWSTREAM,
947
948
949
950
951
		    "%s.%s.%s [%d, %d] (%s + %.2f sec.) * Packet discarded * seq_no_diff=%d  time_diff=%.2fs  lat %.1fs\n",
		    NMXP_LOG_STR(p->pdlist[j]->network), NMXP_LOG_STR(p->pdlist[j]->station), NMXP_LOG_STR(p->pdlist[j]->channel),
		    p->pdlist[j]->packet_type, p->pdlist[j]->seq_no, 
		    NMXP_LOG_STR(str_time), (double) p->pdlist[j]->nSamp / (double) p->pdlist[j]->sampRate,
		    seq_no_diff, time_diff, latency);
952
953
954
	    send_again = 1;
	    j++;
	} else if(seq_no_diff == 1) {
955
	    /* Handle current packet j */
956
957
958
959
	    for(i_func_pd=0; i_func_pd<n_func_pd; i_func_pd++) {
		(*p_func_pd[i_func_pd])(p->pdlist[j]);
	    }
	    if(time_diff > TIME_TOLLERANCE || time_diff < -TIME_TOLLERANCE) {
960
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_ANY,
961
962
963
964
965
966
			"%s.%s.%s [%d, %d] (%s + %.2f sec.) * Time is not correct * last_seq_no_sent=%d  seq_no_diff=%d  time_diff=%.2fs  lat. %.1fs\n",
		    NMXP_LOG_STR(p->pdlist[j]->network), NMXP_LOG_STR(p->pdlist[j]->station), NMXP_LOG_STR(p->pdlist[j]->channel), 
		    p->pdlist[j]->packet_type, p->pdlist[j]->seq_no,
		    str_time, (double) p->pdlist[j]->nSamp /  (double) p->pdlist[j]->sampRate,
		    p->last_seq_no_sent,
		    seq_no_diff, time_diff, latency);
967
968
969
	    }
	    p->last_seq_no_sent = p->pdlist[j]->seq_no;
	    p->last_sample_time = (p->pdlist[j]->time + ((double) p->pdlist[j]->nSamp / (double) p->pdlist[j]->sampRate ));
970
	    p->last_latency = nmxp_data_latency(p->pdlist[j]);
971
972
973
	    send_again = 1;
	    j++;
	} else {
974
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_RAWSTREAM,
975
976
977
978
979
980
		    "%s.%s.%s [%d, %d] (%s + %.2f sec.) * seq_no_diff=%d > 1 * last_seq_no_sent=%d  j=%d  n_pdlist=%2d  time_diff=%.2fs  lat. %.1fs\n",
		    NMXP_LOG_STR(p->pdlist[j]->network), NMXP_LOG_STR(p->pdlist[j]->station), NMXP_LOG_STR(p->pdlist[j]->channel), 
		    p->pdlist[j]->packet_type, p->pdlist[j]->seq_no,
		    str_time, (double) p->pdlist[j]->nSamp /  (double) p->pdlist[j]->sampRate,
		    seq_no_diff, p->last_seq_no_sent, j, p->n_pdlist,
		    time_diff, latency);
981
982
983
984
985
986
987
988
	}
    }

    /* Shift and free j handled elements */
    if(j > 0) {
	for(k=0; k < p->n_pdlist; k++) {
	    if(k < j) {
		if(p->pdlist[k]->buffer) {