mirror of
https://github.com/f4exb/sdrangel.git
synced 2024-11-03 15:31:15 -05:00
435 lines
11 KiB
C++
435 lines
11 KiB
C++
#ifndef LEANSDR_GENERIC_H
|
|
#define LEANSDR_GENERIC_H
|
|
|
|
#include <sys/types.h>
|
|
|
|
#ifdef _MSC_VER
|
|
#include <stdlib.h>
|
|
#include <io.h>
|
|
#else
|
|
#include <unistd.h>
|
|
#endif
|
|
|
|
#include "leansdr/math.h"
|
|
|
|
namespace leansdr
|
|
{
|
|
|
|
//////////////////////////////////////////////////////////////////////
|
|
// Simple blocks
|
|
//////////////////////////////////////////////////////////////////////
|
|
|
|
// [file_reader] reads raw data from a file descriptor into a [pipebuf].
|
|
// If the file descriptor is seekable, data can be looped.
|
|
|
|
template<typename T>
|
|
struct file_reader: runnable
|
|
{
|
|
file_reader(scheduler *sch, int _fdin, pipebuf<T> &_out) :
|
|
runnable(sch, _out.name), loop(false), fdin(_fdin), out(_out)
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
size_t size = out.writable() * sizeof(T);
|
|
if (!size)
|
|
return;
|
|
|
|
#ifdef _MSC_VER
|
|
again: ssize_t nr = _read(fdin, out.wr(), size);
|
|
#else
|
|
again: ssize_t nr = read(fdin, out.wr(), size);
|
|
#endif
|
|
if (nr < 0)
|
|
{
|
|
fatal("leansdr::file_reader::run: read");
|
|
return;
|
|
}
|
|
if (!nr)
|
|
{
|
|
if (!loop)
|
|
return;
|
|
if (sch->debug)
|
|
fprintf(stderr, "leansdr::file_reader::run: %s looping\n", name);
|
|
#ifdef _MSC_VER
|
|
off_t res = _lseek(fdin, 0, SEEK_SET);
|
|
#else
|
|
off_t res = lseek(fdin, 0, SEEK_SET);
|
|
#endif
|
|
if (res == (off_t) -1)
|
|
{
|
|
fatal("leansdr::file_reader::run: lseek");
|
|
return;
|
|
}
|
|
goto again;
|
|
}
|
|
|
|
// Always stop at element boundary (may block)
|
|
size_t partial = nr % sizeof(T);
|
|
size_t remain = partial ? sizeof(T) - partial : 0;
|
|
while (remain)
|
|
{
|
|
if (sch->debug)
|
|
fprintf(stderr, "+");
|
|
#ifdef _MSC_VER
|
|
ssize_t nr2 = _read(fdin, (char*) out.wr() + nr, remain);
|
|
#else
|
|
ssize_t nr2 = read(fdin, (char*) out.wr() + nr, remain);
|
|
#endif
|
|
if (nr2 <= 0)
|
|
{
|
|
fatal("leansdr::file_reader::run: partial read");
|
|
return;
|
|
}
|
|
nr += nr2;
|
|
remain -= nr2;
|
|
}
|
|
|
|
out.written(nr / sizeof(T));
|
|
}
|
|
bool loop;
|
|
private:
|
|
int fdin;
|
|
pipewriter<T> out;
|
|
};
|
|
|
|
// [file_writer] writes raw data from a [pipebuf] to a file descriptor.
|
|
|
|
template<typename T>
|
|
struct file_writer: runnable
|
|
{
|
|
file_writer(scheduler *sch, pipebuf<T> &_in, int _fdout) :
|
|
runnable(sch, _in.name), in(_in), fdout(_fdout)
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
int size = in.readable() * sizeof(T);
|
|
if (!size)
|
|
return;
|
|
#ifdef _MSC_VER
|
|
int nw = _write(fdout, in.rd(), size);
|
|
#else
|
|
int nw = write(fdout, in.rd(), size);
|
|
#endif
|
|
if (!nw)
|
|
{
|
|
fatal("leansdr::file_writer::run: pipe");
|
|
return;
|
|
}
|
|
if (nw < 0)
|
|
{
|
|
fatal("leansdr::file_writer::run: write");
|
|
return;
|
|
}
|
|
if (nw % sizeof(T))
|
|
{
|
|
fatal("leansdr::file_writer::run:partial write");
|
|
return;
|
|
}
|
|
in.read(nw / sizeof(T));
|
|
}
|
|
private:
|
|
pipereader<T> in;
|
|
int fdout;
|
|
};
|
|
|
|
// [file_printer] writes data from a [pipebuf] to a file descriptor,
|
|
// with printf-style formatting and optional scaling.
|
|
|
|
template<typename T>
|
|
struct file_printer: runnable
|
|
{
|
|
file_printer(scheduler *sch, const char *_format, pipebuf<T> &_in, int _fdout, int _decimation = 1) :
|
|
runnable(sch, _in.name), scale(1), decimation(_decimation), in(_in), format(_format), fdout(_fdout), phase(0)
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
int n = in.readable();
|
|
T *pin = in.rd(), *pend = pin + n;
|
|
for (; pin < pend; ++pin)
|
|
{
|
|
if (++phase >= decimation)
|
|
{
|
|
phase -= decimation;
|
|
char buf[256];
|
|
int len = snprintf(buf, sizeof(buf), format, (*pin) * scale);
|
|
if (len < 0)
|
|
{
|
|
fatal("leansdr::file_printer::run: obsolete glibc");
|
|
return;
|
|
}
|
|
#ifdef _MSC_VER
|
|
int nw = _write(fdout, buf, len);
|
|
#else
|
|
int nw = write(fdout, buf, len);
|
|
#endif
|
|
if (nw != len)
|
|
{
|
|
fatal("leansdr::file_printer::run: partial write");
|
|
return;
|
|
}
|
|
}
|
|
}
|
|
in.read(n);
|
|
}
|
|
T scale;
|
|
int decimation;
|
|
private:
|
|
pipereader<T> in;
|
|
const char *format;
|
|
int fdout;
|
|
int phase;
|
|
};
|
|
|
|
// [file_carrayprinter] writes all data available from a [pipebuf]
|
|
// to a file descriptor on a single line.
|
|
// Special case for complex.
|
|
|
|
template<typename T>
|
|
struct file_carrayprinter: runnable
|
|
{
|
|
file_carrayprinter(scheduler *sch, const char *_head, const char *_format, const char *_sep, const char *_tail, pipebuf<complex<T> > &_in, int _fdout) :
|
|
runnable(sch, _in.name), scale(1), fixed_size(0), in(_in), head(_head), format(_format), sep(_sep), tail(_tail), fout(fdopen(_fdout, "w"))
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
int n, nmin = fixed_size ? fixed_size : 1;
|
|
while ((n = in.readable()) >= nmin)
|
|
{
|
|
if (fixed_size)
|
|
n = fixed_size;
|
|
if (fout)
|
|
{
|
|
fprintf(fout, head, n);
|
|
complex<T> *pin = in.rd();
|
|
for (int i = 0; i < n; ++i)
|
|
{
|
|
if (i)
|
|
fprintf(fout, "%s", sep);
|
|
fprintf(fout, format, pin[i].re * scale, pin[i].im * scale);
|
|
}
|
|
fprintf(fout, "%s", tail);
|
|
}
|
|
fflush(fout);
|
|
in.read(n);
|
|
}
|
|
}
|
|
T scale;
|
|
int fixed_size; // Number of elements per batch, or 0.
|
|
private:
|
|
pipereader<complex<T> > in;
|
|
const char *head, *format, *sep, *tail;
|
|
FILE *fout;
|
|
};
|
|
|
|
template<typename T, int N>
|
|
struct file_vectorprinter: runnable
|
|
{
|
|
file_vectorprinter(scheduler *sch, const char *_head, const char *_format, const char *_sep, const char *_tail, pipebuf<T[N]> &_in, int _fdout) :
|
|
runnable(sch, _in.name), scale(1), in(_in), head(_head), format(_format), sep(_sep), tail(_tail)
|
|
{
|
|
fout = fdopen(_fdout, "w");
|
|
if (!fout)
|
|
{
|
|
fatal("leansdr::file_vectorprinter::file_vectorprinter: fdopen");
|
|
}
|
|
}
|
|
void run()
|
|
{
|
|
while (in.readable() >= 1)
|
|
{
|
|
fprintf(fout, head, N);
|
|
T (*pin)[N] = in.rd();
|
|
for (int i = 0; i < N; ++i)
|
|
{
|
|
if (i)
|
|
fprintf(fout, "%s", sep);
|
|
fprintf(fout, format, (*pin)[i] * scale);
|
|
}
|
|
fprintf(fout, "%s", tail);
|
|
in.read(1);
|
|
}
|
|
fflush(fout);
|
|
}
|
|
T scale;
|
|
private:
|
|
pipereader<T[N]> in;
|
|
const char *head, *format, *sep, *tail;
|
|
FILE *fout;
|
|
};
|
|
|
|
// [itemcounter] writes the number of input items to the output [pipebuf].
|
|
// [Tout] must be a numeric type.
|
|
|
|
template<typename Tin, typename Tout>
|
|
struct itemcounter: runnable
|
|
{
|
|
itemcounter(scheduler *sch, pipebuf<Tin> &_in, pipebuf<Tout> &_out) :
|
|
runnable(sch, "itemcounter"), in(_in), out(_out)
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
if (out.writable() < 1)
|
|
return;
|
|
unsigned long count = in.readable();
|
|
if (!count)
|
|
return;
|
|
out.write(count);
|
|
in.read(count);
|
|
}
|
|
private:
|
|
pipereader<Tin> in;
|
|
pipewriter<Tout> out;
|
|
};
|
|
|
|
// [decimator] forwards 1 in N sample.
|
|
|
|
template<typename T>
|
|
struct decimator: runnable
|
|
{
|
|
unsigned int d;
|
|
|
|
decimator(scheduler *sch, int _d, pipebuf<T> &_in, pipebuf<T> &_out) :
|
|
runnable(sch, "decimator"), d(_d), in(_in), out(_out)
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
unsigned long count = min(in.readable() / d, out.writable());
|
|
T *pin = in.rd(), *pend = pin + count * d, *pout = out.wr();
|
|
for (; pin < pend; pin += d, ++pout)
|
|
*pout = *pin;
|
|
in.read(count * d);
|
|
out.written(count);
|
|
}
|
|
private:
|
|
pipereader<T> in;
|
|
pipewriter<T> out;
|
|
};
|
|
|
|
// [rate_estimator] accumulates counts of two quantities
|
|
// and periodically outputs their ratio.
|
|
|
|
template<typename T>
|
|
struct rate_estimator: runnable
|
|
{
|
|
int sample_size;
|
|
|
|
rate_estimator(scheduler *sch, pipebuf<int> &_num, pipebuf<int> &_den, pipebuf<float> &_rate) :
|
|
runnable(sch, "rate_estimator"), sample_size(10000), num(_num), den(_den), rate(_rate), acc_num(0), acc_den(0)
|
|
{
|
|
}
|
|
|
|
void run()
|
|
{
|
|
if (rate.writable() < 1)
|
|
return;
|
|
int count = min(num.readable(), den.readable());
|
|
int *pnum = num.rd(), *pden = den.rd();
|
|
for (int n = count; n--; ++pnum, ++pden)
|
|
{
|
|
acc_num += *pnum;
|
|
acc_den += *pden;
|
|
}
|
|
num.read(count);
|
|
den.read(count);
|
|
if (acc_den >= sample_size)
|
|
{
|
|
rate.write((float) acc_num / acc_den);
|
|
acc_num = acc_den = 0;
|
|
}
|
|
}
|
|
|
|
private:
|
|
pipereader<int> num, den;
|
|
pipewriter<float> rate;
|
|
T acc_num, acc_den;
|
|
};
|
|
|
|
// SERIALIZER
|
|
|
|
template<typename Tin, typename Tout>
|
|
struct serializer: runnable
|
|
{
|
|
serializer(scheduler *sch, pipebuf<Tin> &_in, pipebuf<Tout> &_out) :
|
|
nin(max((size_t) 1, sizeof(Tin) / sizeof(Tout))), nout(max((size_t) 1, sizeof(Tout) / sizeof(Tin))), in(_in), out(_out, nout)
|
|
{
|
|
(void) sch;
|
|
if (nin * sizeof(Tin) != nout * sizeof(Tout))
|
|
{
|
|
fail("serializer::serializer", "incompatible sizes");
|
|
return;
|
|
}
|
|
}
|
|
void run()
|
|
{
|
|
while (in.readable() >= nin && out.writable() >= nout)
|
|
{
|
|
memcpy(out.wr(), in.rd(), nout * sizeof(Tout));
|
|
in.read(nin);
|
|
out.written(nout);
|
|
}
|
|
}
|
|
private:
|
|
int nin, nout;
|
|
pipereader<Tin> in;
|
|
pipewriter<Tout> out;
|
|
};
|
|
// serializer
|
|
|
|
// [buffer_reader] reads from a user-supplied buffer.
|
|
|
|
template<typename T>
|
|
struct buffer_reader: runnable
|
|
{
|
|
buffer_reader(scheduler *sch, T *_data, int _count, pipebuf<T> &_out) :
|
|
runnable(sch, "buffer_reader"), data(_data), count(_count), out(_out), pos(0)
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
int n = min(out.writable(), (unsigned long) (count - pos));
|
|
memcpy(out.wr(), &data[pos], n * sizeof(T));
|
|
pos += n;
|
|
out.written(n);
|
|
}
|
|
private:
|
|
T *data;
|
|
int count;
|
|
pipewriter<T> out;
|
|
int pos;
|
|
};
|
|
// buffer_reader
|
|
|
|
// [buffer_writer] writes to a user-supplied buffer.
|
|
|
|
template<typename T>
|
|
struct buffer_writer: runnable
|
|
{
|
|
buffer_writer(scheduler *sch, pipebuf<T> &_in, T *_data, int _count) :
|
|
runnable(sch, "buffer_reader"), in(_in), data(_data), count(_count), pos(0)
|
|
{
|
|
}
|
|
void run()
|
|
{
|
|
int n = min(in.readable(), (unsigned long) (count - pos));
|
|
memcpy(&data[pos], in.rd(), n * sizeof(T));
|
|
in.read(n);
|
|
pos += n;
|
|
}
|
|
private:
|
|
pipereader<T> in;
|
|
T *data;
|
|
int count;
|
|
int pos;
|
|
};
|
|
// buffer_writer
|
|
|
|
}// namespace
|
|
|
|
#endif // LEANSDR_GENERIC_H
|