voltlogger_analyzer/interpolator.c

/*
    voltlogger_analyzer
    
    Copyright (C) 2015 Dmitry Yu Okunev <dyokunev@ut.mephi.ru> 0x8E30679C
    
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
    
    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
    
    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */

#define _GNU_SOURCE

#include "configuration.h"
#include "macros.h"

#include <stdio.h>	/* getline()	*/
#include <stdlib.h>	/* atof()	*/
#include <unistd.h>	/* getopt()	*/
#include <math.h>
#include <assert.h>

#include "error.h"
#include "analyzer.h"

enum functype {
	FT_SIN,
	FT_REPEAT,
};
typedef enum functype functype_t;

static char *const functypes[] = {
	[FT_REPEAT] = "repeat",
	NULL
};



int main(int argc, char *argv[]) {
	functype_t functype = FT_REPEAT;
	char realtime = 0;
	float frequency = 50;
	int concurrency = 1, addperiods = 0;
	int channelsNum = 1;
	char *checkpointpath = NULL;
	FILE *input  = stdin;
	FILE *output = stdout;

	// Initializing output subsystem
	{
		int output_method     = OM_STDERR;
		int output_quiet      = 0;
		int output_verbosity  = 9;
		int output_debuglevel = 9;

		error_init(&output_method, &output_quiet, &output_verbosity, &output_debuglevel);
	}

	// Parsing arguments
	char c;
	while ((c = getopt (argc, argv, "ri:N:f:F:c:C:e:s:a:")) != -1) {
		char *arg;
		arg = optarg;

		switch (c)
		{
			case 'f': {
				char *value;
				functype = getsubopt(&arg, functypes, &value);
				break;
			}
			case 'F':
				frequency = atof(optarg);
				break;
			case 'r':
				realtime = 1;
				break;
			case 'N':
				channelsNum = atoi(optarg);
				break;
			case 'c':
				concurrency = atoi(optarg);
				break;
			case 'C':
				checkpointpath = arg;
				break;
			case 'i':
				assert ((input = fopen(arg, "r")) != NULL);
				break;
			case 'a':
				addperiods	= atoi(optarg);
				break;
			default:
				abort ();
		}
	}

	switch (functype) {
		case FT_REPEAT:
			vl_interpolate_repeat(input, output, channelsNum, 0, addperiods, checkpointpath, concurrency, frequency, realtime);
			break;
		default:
			fprintf(stderr, "Unknown approximation function\n");
			abort ();
	}

	return 0;
}