srmio/split.c at master · rclasen/srmio · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
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
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
/*
 * Copyright (c) 2008 Rainer Clasen
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms described in the file LICENSE included in this
 * distribution.
 *
 */

#include "common.h"

/*
 * split data
 * - when recint changes
 * - on negative time steps
 * - on gaps
 * - when data overlaps significantly
 *
 */
srmio_data_t *srmio_data_split( srmio_data_t src, srmio_time_t gap,
	srmio_time_t overlap, srmio_error_t *err )
{
	srmio_data_t *list = NULL;
	unsigned alloc = 0;
	unsigned used = 0;
	unsigned c, delta;

	assert( src );

	srmio_error_set( err, "no data" );

	for( c = 0; c < src->cused; ++c ){
		int split = 0;

		if( c == 0 ){
			++split;

		} else {
			srmio_time_t end = src->chunks[c-1]->time
				+ src->chunks[c-1]->dur;

			if( src->chunks[c]->dur != src->chunks[c-1]->dur ){
				DPRINTF( "found recint change" );
				++split;

			} else if( src->chunks[c]->time > end ){
				if( src->chunks[c]->time - end > gap ){
					DPRINTF( "found gap" );
					++split;
				}

			} else if( end > src->chunks[c]->time ){
				if( end - src->chunks[c]->time > overlap ){
					DPRINTF( "found overlap" );
					++split;
				}

			}
		}

		if( split ){
			DPRINTF( "starting new data %u", used );
			if( used +1 >= alloc ){
				srmio_data_t *new;

				if( NULL == (new = realloc( list,
					sizeof(srmio_data_t) * ( alloc +1 +10)  ) )){
					srmio_error_errno( err, "realloc");
					goto clean;
				}

				alloc += 10;
				list = new;
			}

			if( NULL == (list[used] = srmio_data_header( src, err)))
				goto clean;

			list[++used] = NULL;
		}

		srmio_data_add_chunk( list[used-1], src->chunks[c], err );
	}

	delta = 0;
	for( c = 0; c < used; ++c ){
		unsigned last = delta + list[c]->cused -1;
		unsigned m;

		DPRINTF( "file %d: %d - %d", c, delta, last );

		for( m = 0; m < src->mused; ++m ){
			srmio_marker_t mk = src->marker[m];
			srmio_marker_t nm;

			if( mk->last < delta || mk->first > last )
				continue;

			DPRINTF( "found marker %u: %d - %d", m,
				mk->first, mk->last );

			if( NULL == (nm = srmio_marker_clone( mk, err )))
				goto clean;

			if( nm->first < delta )
				nm->first = 0;
			else
				nm->first -= delta;

			// nm->last < delta was handled, before
			if( nm->last > last )
				nm->last = list[c]->cused -1;
			else
				nm->last -= delta;

			srmio_data_add_markerp( list[c], nm, err );
		}

		delta += list[c]->cused;
	}

	return list;

clean:
	for( c = 0; c < used;  ++c )
		srmio_data_free( list[c] );
	free( list );

	return NULL;
}