chris/OSCAR-code
1
0
Fork 0
mirror of https://gitlab.com/pholy/OSCAR-code.git synced 2025-04-04 18:20:42 +00:00
Code Issues Packages Projects Releases Wiki Activity
d5d77647d3
OSCAR-code/oscar/SleepLib/day.cpp

1682 lines
45 KiB
C++
Raw Blame History

/* SleepLib Day Class Implementation
*
* Copyright (c) 2019-2024 The OSCAR Team
* Copyright (c) 2011-2018 Mark Watkins <mark@jedimark.net>
*
* This file is subject to the terms and conditions of the GNU General Public
* License. See the file COPYING in the main directory of the source code
* for more details. */
#define TEST_MACROS_ENABLEDoff
#include <test_macros.h>
#include <QMultiMap>
#include <algorithm>
#include <cmath>
#include <limits>
#include <QDebug>
#include "day.h"
#include "profiles.h"
Day::Day()
{
d_firstsession = true;
d_summaries_open = false;
d_events_open = false;
d_invalidate = true;
}
Day::~Day()
{
for (auto & sess : sessions) {
delete sess;
}
}
void Day::updateCPAPCache()
{
d_count.clear();
d_sum.clear();
OpenSummary();
QList<ChannelID> channels = getSortedMachineChannels(MT_CPAP, schema::FLAG | schema::MINOR_FLAG | schema::SPAN);
for (const auto code : channels) {
d_count[code] = count(code);
d_sum[code] = count(code);
d_machhours[MT_CPAP] = hours(MT_CPAP);
}
}
Session * Day::firstSession(MachineType type)
{
for (auto & sess : sessions) {
if (!sess->enabled()) continue;
if (sess->type() == type) {
return sess;
}
}
return nullptr;
}
bool Day::addMachine(Machine *mach)
{
invalidate();
if (!machines.contains(mach->type())) {
machines[mach->type()] = mach;
return true;
}
return false;
}
Machine *Day::machine(MachineType type)
{
auto it = machines.find(type);
if (it != machines.end())
return it.value();
return nullptr;
}
QList<Session *> Day::getSessions(MachineType type, bool ignore_enabled)
{
QList<Session *> newlist;
for (auto & sess : sessions) {
if (!ignore_enabled && !sess->enabled())
continue;
if (sess->type() == type)
newlist.append(sess);
}
return newlist;
}
Session *Day::find(SessionID sessid)
{
for (auto & sess : sessions) {
if (sess->session() == sessid) {
return sess;
}
}
return nullptr;
}
Session *Day::find(SessionID sessid, MachineType mt)
{
for (auto & sess : sessions) {
if ((sess->session() == sessid) && (sess->s_machtype == mt)) {
return sess;
}
}
return nullptr;
}
void Day::addSession(Session *s)
{
if (s == nullptr) {
qDebug() << "addSession called with null session pointer";
return;
}
invalidate();
auto mi = machines.find(s->type());
if (mi != machines.end()) {
if (mi.value() != s->machine()) {
qDebug() << "OSCAR can't add session" << s->session()
<< "["+QDateTime::fromTime_t(s->session()).toString("MMM dd, yyyy hh:mm:ss")+"]"
<< "from machine" << mi.value()->serial() << "to machine" << s->machine()->serial()
<< "to this day record, as it already contains a different machine of the same MachineType" << s->type();
return;
}
} else {
machines[s->type()] = s->machine();
}
if (s->first() == 0) {
qWarning() << "Day::addSession discarding session" << s->session()
<< "from machine" << s->machine()->serial() << "with first=0";
return;
}
for (auto & sess : sessions) {
if (sess->session() == s->session() && sess->type() == s->type()) {
// This usually indicates a problem in purging or cleanup somewhere,
// unless there's a problem with a parser.
qCritical() << "Day object" << this->date().toString() << "adding duplicate session" << s->session()
<< "["+QDateTime::fromTime_t(s->session()).toString("MMM dd, yyyy hh:mm:ss")+"]";
// Don't skip this one, since it might have replaced the original elsewhere already.
//return;
}
}
sessions.push_back(s);
}
EventDataType Day::calcMiddle(ChannelID code)
{
int c = p_profile->general->prefCalcMiddle();
if (c == 0) {
return percentile(code, 0.5); // Median
} else if (c == 1 ) {
return wavg(code); // Weighted Average
} else {
return avg(code); // Average
}
}
EventDataType Day::calcMax(ChannelID code)
{
return p_profile->general->prefCalcMax() ? percentile(code, 0.995f) : Max(code);
}
EventDataType Day::calcPercentile(ChannelID code)
{
double p = p_profile->general->prefCalcPercentile() / 100.0;
return percentile(code, p);
}
QString Day::calcMiddleLabel(ChannelID code)
{
int c = p_profile->general->prefCalcMiddle();
if (c == 0) {
return QString("%1 %2").arg(STR_TR_Median).arg(schema::channel[code].label());
} else if (c == 1) {
return QString("%1 %2").arg(STR_TR_WAvg).arg(schema::channel[code].label());
} else {
return QString("%1 %2").arg(STR_TR_Avg).arg(schema::channel[code].label());
}
}
QString Day::calcMaxLabel(ChannelID code)
{
return QString("%1 %2").arg(p_profile->general->prefCalcMax() ? QString("99.5%") : STR_TR_Max).arg(schema::channel[code].label());
}
QString Day::calcPercentileLabel(ChannelID code)
{
return QString("%1% %2").arg(p_profile->general->prefCalcPercentile(),0, 'f',0).arg(schema::channel[code].label());
}
EventDataType Day::countInsideSpan(ChannelID span, ChannelID code)
{
int count = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
count += sess->countInsideSpan(span, code);
}
}
return count;
}
EventDataType Day::lookupValue(ChannelID code, qint64 time, bool square)
{
for (auto & sess : sessions) {
if (sess->enabled()) {
// drift is handled by session first/last/SearchValue
if ((time > sess->first()) && (time < sess->last())) {
if (sess->channelExists(code)) {
return sess->SearchValue(code,time,square);
}
}
}
}
return 0;
}
EventDataType Day::timeAboveThreshold(ChannelID code, EventDataType threshold)
{
EventDataType val = 0;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_availableChannels.contains(code)) {
val += sess->timeAboveThreshold(code,threshold);
}
}
return val;
}
EventDataType Day::timeBelowThreshold(ChannelID code, EventDataType threshold)
{
EventDataType val = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
val += sess->timeBelowThreshold(code,threshold);
}
}
return val;
}
EventDataType Day::settings_sum(ChannelID code)
{
EventDataType val = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
auto set = sess->settings.find(code);
if (set != sess->settings.end()) {
val += set.value().toDouble();
}
}
}
return val;
}
EventDataType Day::settings_max(ChannelID code)
{
EventDataType min = -std::numeric_limits<EventDataType>::max();
EventDataType max = min;
EventDataType value;
for (auto & sess : sessions) {
if (sess->enabled() && sess->s_machtype != MT_JOURNAL) {
value = sess->settings.value(code, min).toFloat();
if (value > max) {
max = value;
}
}
}
return max;
}
EventDataType Day::settings_min(ChannelID code)
{
EventDataType max = std::numeric_limits<EventDataType>::max();
EventDataType min = max;
EventDataType value;
for (auto & sess : sessions) {
if (sess->enabled() && sess->s_machtype != MT_JOURNAL) {
value = sess->settings.value(code, max).toFloat();
if (value < min) {
min = value;
}
}
}
return min;
}
EventDataType Day::settings_avg(ChannelID code)
{
EventDataType val = 0;
int cnt = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
auto set = sess->settings.find(code);
if (set != sess->settings.end()) {
val += set.value().toDouble();
cnt++;
}
}
}
val = (cnt > 0) ? (val / EventDataType(cnt)) : val;
return val;
}
EventDataType Day::settings_wavg(ChannelID code)
{
double s0 = 0, s1 = 0, s2 = 0, tmp;
for (auto & sess : sessions) {
if (sess->enabled()) {
auto set = sess->settings.find(code);
if (set != sess->settings.end()) {
if (code == CPAP_Mode && sess->type() != MT_CPAP) {
// There used to be a bug in gLineChart::SetDay that inserted a CPAP_Mode
// setting in any session that didn't already have one. That shouldn't
// happen any more, but leave this diagnostic message here in case it does.
qWarning() << sess->session() << "non-CPAP session with CPAP mode setting";
continue;
}
s0 = sess->hours();
tmp = set.value().toDouble();
s1 += tmp * s0;
s2 += s0;
}
}
}
if (s2 == 0) { return 0; }
tmp = (s1 / s2);
return tmp;
}
EventDataType Day::percentile(ChannelID code, EventDataType percentile)
{
// Cache this calculation?
// QHash<ChannelID, QHash<EventDataType, EventDataType> >::iterator pi;
// pi=perc_cache.find(code);
// if (pi!=perc_cache.end()) {
// QHash<EventDataType, EventDataType> & hsh=pi.value();
// QHash<EventDataType, EventDataType>::iterator hi=hsh.find(
// if (hi!=pi.value().end()) {
// return hi.value();
// }
// }
QHash<EventStoreType, qint64> wmap; // weight map
QHash<EventStoreType, qint64>::iterator wmapit;
qint64 SN = 0;
EventDataType lastgain = 0, gain = 0;
// First Calculate count of all events
bool timeweight;
for (auto & sess : sessions) {
if (!sess->enabled()) { continue; }
auto ei = sess->m_valuesummary.find(code);
if (ei == sess->m_valuesummary.end()) { continue; }
auto tei = sess->m_timesummary.find(code);
timeweight = (tei != sess->m_timesummary.end());
gain = sess->m_gain[code];
// Here's assuming gains don't change accross a days sessions
// Can't assume this in any multi day calculations..
if (lastgain > 0) {
if (gain != lastgain) {
qDebug() << "Gains differ across sessions: " << gain << lastgain;
}
}
lastgain = gain;
qint64 weight;
//qint64 tval;
if (timeweight) {
wmap.reserve(wmap.size() + tei.value().size());
for (auto it = tei.value().begin(), teival_end=tei.value().end(); it != teival_end; ++it) {
weight = it.value();
SN += weight;
wmap[it.key()] += weight;
}
} else {
wmap.reserve(wmap.size() + ei.value().size());
for (auto it = ei.value().begin(), eival_end=ei.value().end(); it != eival_end; ++it) {
weight = it.value();
SN += weight;
wmap[it.key()] += weight;
}
}
}
QVector<ValueCount> valcnt;
valcnt.resize(wmap.size());
// Build sorted list of value/counts
auto wmap_end = wmap.end();
int ii=0;
for (auto it = wmap.begin(); it != wmap_end; ++it) {
valcnt[ii++]=ValueCount(EventDataType(it.key()) * gain, it.value(), 0);
}
// sort by weight, then value
//qSort(valcnt);
std::sort(valcnt.begin(), valcnt.end());
//double SN=100.0/double(N); // 100% / overall sum
double p = 100.0 * percentile;
double nth = double(SN) * percentile; // index of the position in the unweighted set would be
double nthi = floor(nth);
qint64 sum1 = 0, sum2 = 0;
qint64 w1, w2 = 0;
double v1 = 0, v2;
int N = valcnt.size();
int k = 0;
for (k = 0; k < N; k++) {
v1 = valcnt.at(k).value;
w1 = valcnt.at(k).count;
sum1 += w1;
if (sum1 > nthi) {
return v1;
}
if (sum1 == nthi) {
break; // boundary condition
}
}
if (k >= N) {
return v1;
}
if (valcnt.size() == 1) {
return valcnt[0].value;
}
v2 = valcnt[k + 1].value;
w2 = valcnt[k + 1].count;
sum2 = sum1 + w2;
// value lies between v1 and v2
double px = 100.0 / double(SN); // Percentile represented by one full value
// calculate percentile ranks
double p1 = px * (double(sum1) - (double(w1) / 2.0));
double p2 = px * (double(sum2) - (double(w2) / 2.0));
// calculate linear interpolation
double v = v1 + ((p - p1) / (p2 - p1)) * (v2 - v1);
return v;
// p1.....p.............p2
// 37 55 70
}
EventDataType Day::p90(ChannelID code)
{
return percentile(code, 0.90F);
}
EventDataType Day::rangeCount(ChannelID code, qint64 st, qint64 et)
{
int cnt = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
cnt += sess->rangeCount(code, st, et);
}
}
return cnt;
}
EventDataType Day::rangeSum(ChannelID code, qint64 st, qint64 et)
{
double val = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
val += sess->rangeSum(code, st, et);
}
}
return val;
}
EventDataType Day::rangeAvg(ChannelID code, qint64 st, qint64 et)
{
double val = 0;
int cnt = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
val += sess->rangeSum(code, st, et);
cnt += sess->rangeCount(code, st,et);
}
}
if (cnt == 0) { return 0; }
val /= double(cnt);
return val;
}
EventDataType Day::rangeWavg(ChannelID code, qint64 st, qint64 et)
{
double sum = 0;
double cnt = 0;
qint64 lasttime, time;
double data, duration;
for (auto & sess : sessions) {
auto EVEC = sess->eventlist.find(code);
if (EVEC == sess->eventlist.end()) continue;
for (auto & el : EVEC.value()) {
if (el->count() < 1) continue;
lasttime = el->time(0);
if (lasttime < st)
lasttime = st;
for (unsigned i=1; i<el->count(); i++) {
data = el->data(i);
time = el->time(i);
if (time < st) {
lasttime = st;
continue;
}
if (time > et) {
time = et;
}
duration = double(time - lasttime) / 1000.0;
sum += data * duration;
cnt += duration;
if (time >= et) break;
lasttime = time;
}
}
}
if (cnt < 0.000001)
return 0;
return sum / cnt;
}
// Boring non weighted percentile
EventDataType Day::rangePercentile(ChannelID code, float p, qint64 st, qint64 et)
{
int count = rangeCount(code, st,et);
QVector<EventDataType> list;
list.resize(count);
int idx = 0;
qint64 time;
for (auto & sess : sessions) {
auto EVEC = sess->eventlist.find(code);
if (EVEC == sess->eventlist.end()) continue;
for (auto & el : EVEC.value()) {
for (unsigned i=0; i<el->count(); i++) {
time = el->time(i);
if ((time < st) || (time > et)) continue;
list[idx++] = el->data(i);
}
}
}
// TODO: use nth_element instead..
//qSort(list);
std::sort(list.begin(), list.end());
float b = float(idx) * p;
int a = floor(b);
int c = ceil(b);
if ((a == c) || (c >= idx)) {
return list[a];
}
EventDataType v1 = list[a];
EventDataType v2 = list[c];
EventDataType diff = v2 - v1; // the whole == C-A
double ba = b - float(a); // A....B...........C == B-A
double val = v1 + diff * ba;
return val;
}
EventDataType Day::avg(ChannelID code)
{
double val = 0;
// Cache this?
int cnt = 0;
for (auto & sess : sessions) {
if (sess->enabled()) {
val += sess->sum(code);
cnt += sess->count(code);
}
}
if (cnt == 0) { return 0; }
val /= double(cnt);
return val;
}
EventDataType Day::sum(ChannelID code)
{
// Cache this?
EventDataType val = 0;
if (code == AllAhiChannels) {
for (int i = 0; i < ahiChannels.size(); i++)
val += sum(ahiChannels.at(i));
return val;
}
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_sum.contains(code)) {
val += sess->sum(code);
}
}
return val;
}
EventDataType Day::wavg(ChannelID code)
{
double s0 = 0, s1 = 0, s2 = 0;
qint64 d;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_wavg.contains(code)) {
d = sess->length();
s0 = double(d) / 3600000.0;
if (s0 > 0) {
s1 += sess->wavg(code) * s0;
s2 += s0;
}
}
}
if (s2 == 0) {
return 0;
}
return (s1 / s2);
}
// Total session time in milliseconds
qint64 Day::total_time()
{
qint64 d_totaltime = 0;
QMultiMap<qint64, bool> range;
//range.reserve(size()*2);
// Remember sessions may overlap..
qint64 first, last;
for (auto & sess : sessions) {
int slicesize = sess->m_slices.size();
if (sess->enabled() && (sess->type() != MT_JOURNAL)) {
first = sess->first();
last = sess->last();
if (slicesize == 0) {
// This algorithm relies on non zero length, and correctly ordered sessions
if (last > first) {
range.insert(first, 0);
range.insert(last, 1);
d_totaltime += sess->length();
if (sess->length() == 0) {
qWarning() << sess->s_session << "0 length session";
}
}
} else {
for (auto & slice : sess->m_slices) {
if (slice.status == MaskOn) {
range.insert(slice.start, 0);
range.insert(slice.end, 1);
d_totaltime += slice.end - slice.start;
if (slice.end - slice.start == 0) {
qWarning() << sess->s_session << "0 length slice";
}
}
}
}
}
}
bool b;
int nest = 0;
qint64 ti = 0;
qint64 total = 0;
// This is my implementation of a typical "brace counting" algorithm mentioned here:
// http://stackoverflow.com/questions/7468948/problem-calculating-overlapping-date-ranges
auto rend = range.end();
for (auto rit = range.begin(); rit != rend; ++rit) {
b = rit.value();
if (!b) {
if (!ti) {
ti = rit.key();
}
nest++;
} else {
if (--nest <= 0) {
total += rit.key() - ti;
ti = 0;
}
}
}
if (total != d_totaltime) {
// They can overlap.. tough.
// qDebug() << "Sessions Times overlaps!" << total << d_totaltime;
}
return total; //d_totaltime;
}
// Total session time in milliseconds, only considering device type
qint64 Day::total_time(MachineType type)
{
qint64 d_totaltime = 0;
QMultiMap<qint64, bool> range;
//range.reserve(size()*2);
// Remember sessions may overlap..
qint64 first, last;
for (auto & sess : sessions) {
int slicesize = sess->m_slices.size();
if (sess->enabled() && (sess->type() == type)) {
first = sess->first();
last = sess->last();
// This algorithm relies on non zero length, and correctly ordered sessions
if (slicesize == 0) {
if (last > first) {
range.insert(first, 0);
range.insert(last, 1);
d_totaltime += sess->length();
if (sess->length() == 0) {
qWarning() << sess->s_session << "0 length session";
}
}
} else {
for (const auto & slice : sess->m_slices) {
if (slice.status == MaskOn) {
range.insert(slice.start, 0);
range.insert(slice.end, 1);
d_totaltime += slice.end - slice.start;
if (slice.end - slice.start == 0) {
qWarning() << sess->s_session << "0 length slice";
}
}
}
}
}
}
bool b;
int nest = 0;
qint64 ti = 0;
qint64 total = 0;
// This is my implementation of a typical "brace counting" algorithm mentioned here:
// http://stackoverflow.com/questions/7468948/problem-calculating-overlapping-date-ranges
auto rend = range.end();
for (auto rit = range.begin(); rit != rend; ++rit) {
b = rit.value();
if (!b) {
if (!ti) {
ti = rit.key();
}
nest++;
} else {
if (--nest <= 0) {
total += rit.key() - ti;
ti = 0;
}
}
}
if (total != d_totaltime) {
// They can overlap.. tough.
// qDebug() << "Sessions Times overlaps!" << total << d_totaltime;
}
return total; //d_totaltime;
}
ChannelID Day::getPressureChannelID() {
// TODO: This is an awful hack that depends on the enum ordering of the generic CPAP_Mode channel.
// See the comment in getCPAPModeStr().
// Determined the preferred pressure channel (CPAP_IPAP or CPAP_Pressure)
CPAPMode cpapmode = (CPAPMode)(int)settings_max(CPAP_Mode);
// TODO: PRS1 ventilators in CPAP mode report IPAP rather than pressure...but their pressure setting channel is CPAP_Pressure,
// so this currently gets fixed in the welcome screen manually.
// And why would ASV or AVAPS have Pressure channels?
QList<ChannelID> preferredIDs = { CPAP_Pressure, CPAP_PressureSet, CPAP_IPAP, CPAP_IPAPSet };
if (cpapmode == MODE_ASV || cpapmode == MODE_ASV_VARIABLE_EPAP || cpapmode == MODE_AVAPS ||
cpapmode == MODE_BILEVEL_FIXED || cpapmode == MODE_BILEVEL_AUTO_FIXED_PS || cpapmode == MODE_BILEVEL_AUTO_VARIABLE_PS || cpapmode == MODE_TRILEVEL_AUTO_VARIABLE_PDIFF) {
preferredIDs = { CPAP_IPAP, CPAP_IPAPSet, CPAP_Pressure, CPAP_PressureSet };
}
for (auto & preferredID : preferredIDs) {
// If preferred channel has data, return it
if (channelHasData(preferredID)) {
// if ( cpapmode == MODE_AVAPS )
// qDebug() << QString("Found pressure channel 0x%1").arg(preferredID, 4, 16, QChar('0'));
return preferredID;
}
}
qDebug() << "No pressure channel for " << getCPAPModeStr();
return NoChannel;
}
bool Day::hasEnabledSessions()
{
for (auto & sess : sessions) {
if (sess->enabled()) {
return true;
}
}
return false;
}
bool Day::hasEnabledSessions(MachineType type)
{
for (auto & sess : sessions) {
if ((sess->type() == type) && sess->enabled()) {
return true;
}
}
return false;
}
/*EventDataType Day::percentile(ChannelID code,double percent)
{
double val=0;
int cnt=0;
for (auto & sess : sessions) {
if (sess->eventlist.find(code)!=sess->eventlist.end()) {
val+=sess->percentile(code,percent);
cnt++;
}
}
if (cnt==0) return 0;
return EventDataType(val/cnt);
}*/
qint64 Day::first(ChannelID code)
{
qint64 date = 0;
qint64 tmp;
for (auto & sess : sessions) {
if (sess->enabled()) {
tmp = sess->first(code);
if (!tmp) { continue; }
if (!date) {
date = tmp;
} else {
if (tmp < date) { date = tmp; }
}
}
}
return date;
}
qint64 Day::last(ChannelID code)
{
qint64 date = 0;
qint64 tmp;
for (auto & sess : sessions) {
if (sess->enabled()) {
tmp = sess->last(code);
if (!tmp) { continue; }
if (!date) {
date = tmp;
} else {
if (tmp > date) { date = tmp; }
}
}
}
return date;
}
EventDataType Day::Min(ChannelID code)
{
EventDataType min = 0;
EventDataType tmp;
bool first = true;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_min.contains(code)) {
tmp = sess->Min(code);
if (first) {
min = tmp;
first = false;
} else {
if (tmp < min) { min = tmp; }
}
}
}
return min;
}
EventDataType Day::physMin(ChannelID code)
{
EventDataType min = 0;
EventDataType tmp;
bool first = true;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_min.contains(code)) {
tmp = sess->physMin(code);
if (first) {
min = tmp;
first = false;
} else {
if (tmp < min) { min = tmp; }
}
}
}
return min;
}
bool Day::hasData(ChannelID code, SummaryType type)
{
bool has = false;
for (auto & sess : sessions) {
if (sess->type() == MT_JOURNAL) continue;
if (sess->enabled()) {
switch (type) {
// case ST_90P:
// has=sess->m_90p.contains(code);
// break;
case ST_PERC:
has = sess->m_valuesummary.contains(code);
break;
case ST_MIN:
has = sess->m_min.contains(code);
break;
case ST_MAX:
has = sess->m_max.contains(code);
break;
case ST_CNT:
has = sess->m_cnt.contains(code);
break;
case ST_AVG:
has = sess->m_avg.contains(code);
break;
case ST_WAVG:
has = sess->m_wavg.contains(code);
break;
case ST_CPH:
has = sess->m_cph.contains(code);
break;
case ST_SPH:
has = sess->m_sph.contains(code);
break;
case ST_FIRST:
has = sess->m_firstchan.contains(code);
break;
case ST_LAST:
has = sess->m_lastchan.contains(code);
break;
case ST_SUM:
has = sess->m_sum.contains(code);
break;
default:
break;
}
if (has) { break; }
}
}
return has;
}
EventDataType Day::Max(ChannelID code)
{
EventDataType max = 0;
EventDataType tmp;
bool first = true;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_max.contains(code)) {
tmp = sess->Max(code);
if (first) {
max = tmp;
first = false;
} else {
if (tmp > max) { max = tmp; }
}
}
}
return max;
}
EventDataType Day::physMax(ChannelID code)
{
EventDataType max = 0;
EventDataType tmp;
bool first = true;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_max.contains(code)) {
tmp = sess->physMax(code);
if (first) {
max = tmp;
first = false;
} else {
if (tmp > max) { max = tmp; }
}
}
}
return max;
}
EventDataType Day::cph(ChannelID code)
{
double sum = 0;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_cnt.contains(code)) {
sum += sess->count(code);
}
}
sum /= hours();
return sum;
}
EventDataType Day::sph(ChannelID code)
{
EventDataType sum = 0;
EventDataType h = 0;
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_sum.contains(code)) {
sum += sess->sum(code) / 3600.0;
}
}
h = hours();
sum = (100.0 / h) * sum;
return sum;
}
EventDataType Day::count(ChannelID code)
{
EventDataType total = 0;
if (code == AllAhiChannels) {
for (int i = 0; i < ahiChannels.size(); i++)
total += count(ahiChannels.at(i));
return total;
}
for (auto & sess : sessions) {
if (sess->enabled() && sess->m_cnt.contains(code)) {
total += sess->count(code);
}
}
return total;
}
bool Day::noSettings(Machine * mach)
{
for (auto & sess : sessions) {
if ((mach == nullptr) && sess->noSettings()) {
// If this day generally has just summary data.
return true;
} else if ((mach == sess->machine()) && sess->noSettings()) {
// Focus only on device match
return true;
}
}
return false;
}
bool Day::summaryOnly(Machine * mach)
{
for (auto & sess : sessions) {
if ((mach == nullptr) && sess->summaryOnly()) {
// If this day generally has just summary data.
return true;
} else if ((mach == sess->machine()) && sess->summaryOnly()) {
// Focus only on device match
return true;
}
}
return false;
}
bool Day::settingExists(ChannelID id)
{
for (auto & sess : sessions) {
if (sess->enabled()) {
auto set = sess->settings.find(id);
if (set != sess->settings.end()) {
return true;
}
}
}
return false;
}
bool Day::eventsLoaded()
{
for (auto & sess : sessions) {
if (sess->eventsLoaded()) {
return true;
}
}
return false;
}
bool Day::channelExists(ChannelID id)
{
for (auto & sess : sessions) {
if (sess->enabled() && sess->eventlist.contains(id)) {
return true;
}
}
return false;
}
bool Day::hasEvents() {
for (auto & sess : sessions) {
if (sess->eventlist.size() > 0) return true;
}
return false;
}
bool Day::channelHasData(ChannelID id)
{
for (auto & sess : sessions) {
if (sess->enabled()) {
if (sess->m_cnt.contains(id)) {
return true;
}
if (sess->eventlist.contains(id)) {
return true;
}
if (sess->m_valuesummary.contains(id)) {
return true;
}
}
}
return false;
}
void Day::OpenEvents()
{
for (auto & sess : sessions) {
if (sess->type() != MT_JOURNAL)
sess->OpenEvents();
}
d_events_open = true;
}
void Day::OpenSummary()
{
if (d_summaries_open) return;
for (auto & sess : sessions) {
sess->LoadSummary();
}
d_summaries_open = true;
}
void Day::CloseEvents()
{
for (auto & sess : sessions) {
sess->TrashEvents();
}
d_events_open = false;
}
QList<ChannelID> Day::getSortedMachineChannels(MachineType type, quint32 chantype)
{
QList<ChannelID> available;
auto mi_end = machines.end();
for (auto mi = machines.begin(); mi != mi_end; mi++) {
if (mi.key() != type) continue;
available.append(mi.value()->availableChannels(chantype));
}
QMultiMap<int, ChannelID> order;
for (const auto code : available) {
order.insert(schema::channel[code].order(), code);
}
QList<ChannelID> channels;
for (auto it = order.begin(); it != order.end(); ++it) {
ChannelID code = it.value();
channels.append(code);
}
return channels;
}
QList<ChannelID> Day::getSortedMachineChannels(quint32 chantype)
{
QList<ChannelID> available;
auto mi_end = machines.end();
for (auto mi = machines.begin(); mi != mi_end; mi++) {
if (mi.key() == MT_JOURNAL) continue;
available.append(mi.value()->availableChannels(chantype));
}
QMultiMap<int, ChannelID> order;
for (auto code : available) {
order.insert(schema::channel[code].order(), code);
}
QList<ChannelID> channels;
for (auto it = order.begin(); it != order.end(); ++it) {
ChannelID code = it.value();
channels.append(code);
}
return channels;
}
qint64 Day::first(MachineType type)
{
qint64 date = 0;
qint64 tmp;
for (auto & sess : sessions) {
if ((sess->type() == type) && sess->enabled()) {
tmp = sess->first();
if (!tmp) { continue; }
if (!date) {
date = tmp;
} else {
if (tmp < date) { date = tmp; }
}
}
}
return date;
}
qint64 Day::first()
{
qint64 date = 0;
qint64 tmp;
for (auto & sess : sessions) {
if (sess->type() == MT_JOURNAL) continue;
if (sess->enabled()) {
tmp = sess->first();
if (!tmp) { continue; }
if (!date) {
date = tmp;
} else {
if (tmp < date) { date = tmp; }
}
}
}
return date;
}
//! \brief Returns the last session time of this day
qint64 Day::last()
{
qint64 date = 0;
qint64 tmp;
for (auto & sess : sessions) {
if (sess->type() == MT_JOURNAL) continue;
if (sess->enabled()) {
tmp = sess->last();
if (!tmp) { continue; }
if (!date) {
date = tmp;
} else {
if (tmp > date) { date = tmp; }
}
}
}
return date;
}
qint64 Day::last(MachineType type)
{
qint64 date = 0;
qint64 tmp;
for (auto & sess : sessions) {
if ((sess->type() == type) && sess->enabled()) {
tmp = sess->last();
if (!tmp) { continue; }
if (!date) {
date = tmp;
} else {
if (tmp > date) { date = tmp; }
}
}
}
return date;
}
bool Day::removeSession(Session *sess)
{
sess->machine()->sessionlist.remove(sess->session());
MachineType mt = sess->type();
bool b = sessions.removeAll(sess) > 0;
if (!searchMachine(mt)) {
machines.remove(mt);
}
return b;
}
bool Day::searchMachine(MachineType mt) {
for (auto & sess : sessions) {
if (sess->type() == mt)
return true;
}
return false;
}
bool Day::hasMachine(Machine * mach)
{
for (auto & sess : sessions) {
if (sess->machine() == mach)
return true;
}
return false;
}
void Day::removeMachine(Machine * mach)
{
// Yell about and fix any dangling references rather than possibly crashing later.
//
// This has no functional use and can be removed when the data structures are cleaned up
// with better encapsulation and fewer unnecessary references between each other.
QList<Session*> list = sessions; // make a copy so the iterator doesn't get broken by removals
for (auto & sess : list) {
if (sess->machine() == mach) {
// This indicates a problem with the device class not tracking all of its sessions, for
// example if there's a duplicate session ID.
qCritical() << "Day object" << this->date().toString()
<< "session" << sess->session() << "refers to machine" << mach->serial();
removeSession(sess);
}
}
for (auto & m : machines.keys()) {
if (machines[m] == mach) {
// This indicates a problem internal to the Day class, since removeSession should remove
// machines from this list if there are no longer any sessions pointing to it.
qCritical() << "Day object" << this->date().toString()
<< "refers to machine" << mach->serial();
machines.remove(m);
}
}
}
int Day::getCPAPMode()
{
// NOTE: This needs to return the generic mode, unlike getCPAPModeStr().
// This function is used only to determine whether to use advanced graphs,
// which refer to the generic mode.
/*
Machine * mach = machine(MT_CPAP);
if (!mach) return 0;
CPAPLoader * loader = qobject_cast<CPAPLoader *>(mach->loader());
ChannelID modechan = loader->CPAPModeChannel();
*/
ChannelID modechan = CPAP_Mode;
// schema::Channel & chan = schema::channel[modechan];
// TODO: This is an awful hack that depends on the enum ordering of the device-specific CPAP mode.
// See the comment in getCPAPModeStr().
int mode = (CPAPMode)(int)qRound(settings_wavg(modechan));
return mode;
}
QString Day::getCPAPModeStr()
{
Machine * mach = machine(MT_CPAP);
if (!mach) return STR_MessageBox_Error;
CPAPLoader * loader = qobject_cast<CPAPLoader *>(mach->loader());
ChannelID modechan = loader->CPAPModeChannel();
schema::Channel & chan = schema::channel[modechan];
// TODO: This is an awful hack that depends on the enum ordering of the device-specific CPAP mode.
// Instead, we should calculate how long each mode was in operation and
// determine the one that was running the longest, along with the settings
// while that mode was in operation.
int mode = (CPAPMode)(int)qRound(settings_wavg(modechan));
return chan.option(mode);
// if (mode == MODE_CPAP) {
// return QObject::tr("Fixed");
// } else if (mode == MODE_APAP) {
// return QObject::tr("Auto");
// } else if (mode == MODE_BILEVEL_FIXED ) {
// return QObject::tr("Fixed Bi-Level");
// } else if (mode == MODE_BILEVEL_AUTO_FIXED_PS) {
// return QObject::tr("Auto Bi-Level (Fixed PS)");
// } else if (mode == MODE_BILEVEL_AUTO_VARIABLE_PS) {
// return QObject::tr("Auto Bi-Level (Variable PS)");
// } else if (mode == MODE_TRILEVEL_AUTO_VARIABLE_PDIFF) {
// return QObject::tr("Auto TriLevel (Variable PDIFF)");
// } else if (mode == MODE_ASV) {
// return QObject::tr("ASV Fixed EPAP");
// } else if (mode == MODE_ASV_VARIABLE_EPAP) {
// return QObject::tr("ASV Variable EPAP");
// }
// return STR_TR_Unknown;
}
QString Day::getPressureRelief()
{
Machine * mach = machine(MT_CPAP);
if (!mach) return STR_MessageBox_Error;
CPAPLoader * loader = qobject_cast<CPAPLoader *>(mach->loader());
if (!loader) return STR_MessageBox_Error;
QString pr_str;
ChannelID pr_level_chan = loader->PresReliefLevel();
ChannelID pr_mode_chan = loader->PresReliefMode();
// Separate calculation for SleepStyle devices
if (mach->info.loadername == "SleepStyle") {
pr_str = loader->PresReliefLabel();
int pr_level = -1;
if (pr_level_chan != NoChannel && settingExists(pr_level_chan)) {
pr_level = qRound(settings_wavg(pr_level_chan));
}
if (pr_level == -1)
return STR_TR_None;
if ((pr_mode_chan != NoChannel) && settingExists(pr_mode_chan)) {
schema::Channel & chan = schema::channel[pr_level_chan];
QString level = chan.option(pr_level);
if (level.isEmpty()) {
level = QString().number(pr_level) + " " + chan.units();
if (settings_min(pr_level_chan) != settings_max(pr_level_chan))
level = QObject::tr("varies");
}
pr_str += QString(" %1").arg(level);
}
return pr_str;
}
if ((pr_mode_chan != NoChannel) && settingExists(pr_mode_chan)) {
// TODO: This is an awful hack that depends on the enum ordering of the pressure relief mode.
// See the comment in getCPAPModeStr().
int pr_mode = qRound(settings_wavg(pr_mode_chan));
pr_str = QString("%1%2").arg(loader->PresReliefLabel()).arg(schema::channel[pr_mode_chan].option(pr_mode));
int pr_level = -1;
if (pr_level_chan != NoChannel && settingExists(pr_level_chan)) {
pr_level = qRound(settings_wavg(pr_level_chan));
}
if (pr_level >= 0) {
// TODO: Ideally the formatting of LOOKUP datatypes should be done in only one place.
schema::Channel & chan = schema::channel[pr_level_chan];
QString level = chan.option(pr_level);
if (level.isEmpty()) {
level = QString().number(pr_level) + " " + chan.units();;
}
pr_str += QString(" %1").arg(level);
}
} else pr_str = STR_TR_None;
return pr_str;
}
QString Day::validPressure(float pressure) {
if (fabsf(pressure) == std::numeric_limits<EventDataType>::max())
return QObject::tr("n/a");
return QString("%1").arg(pressure, 0, 'f', 1);
}
QString Day::getPressureSettings()
{
if (machine(MT_CPAP) == nullptr) {
qCritical("getPressureSettings called with no CPAP machine record");
return QString();
}
// TODO: This is an awful hack that depends on the enum ordering of the generic CPAP_Mode channel.
// See the comment in getCPAPModeStr().
CPAPMode mode = (CPAPMode)(int)settings_max(CPAP_Mode);
QString units = schema::channel[CPAP_Pressure].units();
if (mode == MODE_CPAP) {
return QObject::tr("Fixed %1 (%2)").arg(validPressure(settings_min(CPAP_Pressure))).
arg(units);
} else if (mode == MODE_APAP) {
return QObject::tr("Min %1 Max %2 (%3)").arg(validPressure(settings_min(CPAP_PressureMin))).
arg(validPressure(settings_max(CPAP_PressureMax))).
arg(units);
} else if (mode == MODE_BILEVEL_FIXED ) {
return QObject::tr("EPAP %1 IPAP %2 (%3)").arg(validPressure(settings_min(CPAP_EPAP))).
arg(validPressure(settings_max(CPAP_IPAP))).
arg(units);
} else if (mode == MODE_BILEVEL_AUTO_FIXED_PS) {
return QObject::tr("PS %1 over %2-%3 (%4)").arg(validPressure(settings_max(CPAP_PS))).
arg(validPressure(settings_min(CPAP_EPAPLo))).
arg(validPressure(settings_max(CPAP_IPAPHi))).
arg(units);
} else if (mode == MODE_BILEVEL_AUTO_VARIABLE_PS) {
return QObject::tr("Min EPAP %1 Max IPAP %2 PS %3-%4 (%5)").arg(validPressure(settings_min(CPAP_EPAPLo))).
arg(validPressure(settings_max(CPAP_IPAPHi))).
arg(validPressure(settings_min(CPAP_PSMin))).
arg(validPressure(settings_max(CPAP_PSMax))).arg(units);
} else if (mode == MODE_TRILEVEL_AUTO_VARIABLE_PDIFF) {
return QObject::tr("Min EEPAP %1 Max EEPAP %2 PDIFF %3-%4 (%5)").arg(validPressure(settings_min(CPAP_EEPAPLo))).
arg(validPressure(settings_max(CPAP_EEPAPHi))).
arg(validPressure(settings_min(CPAP_PSMin))).
arg(validPressure(settings_max(CPAP_PSMax))).arg(units);
} else if (mode == MODE_ASV) {
return QObject::tr("EPAP %1 PS %2-%3 (%4)").arg(validPressure(settings_min(CPAP_EPAP))).
arg(validPressure(settings_min(CPAP_PSMin))).
arg(validPressure(settings_max(CPAP_PSMax))).arg(units);
} else if (mode == MODE_ASV_VARIABLE_EPAP) {
return QObject::tr("Min EPAP %1 Max IPAP %2 PS %3-%4 (%5)").
arg(validPressure(settings_min(CPAP_EPAPLo))).
arg(validPressure(settings_max(CPAP_IPAPHi))).
arg(validPressure(settings_max(CPAP_PSMin))).
arg(validPressure(settings_min(CPAP_PSMax))).
arg(units);
} else if (mode == MODE_AVAPS) {
// qDebug() << "AVAPS: EPAP" << settings_min(CPAP_EPAP) << "IPAP min" << settings_max(CPAP_IPAPLo) <<
// "IPAP max" << settings_max(CPAP_IPAPHi);
QString retStr;
if (settings_min(CPAP_EPAPLo) == settings_max(CPAP_EPAPHi))
retStr = QObject::tr("EPAP %1 IPAP %2-%3 (%4)").
arg(validPressure(settings_min(CPAP_EPAP))).
arg(validPressure(settings_max(CPAP_IPAPLo))).
arg(validPressure(settings_max(CPAP_IPAPHi))).
arg(units);
else
retStr = QObject::tr("EPAP %1-%2 IPAP %3-%4 (%5)").
arg(validPressure(settings_min(CPAP_EPAPLo))).
arg(validPressure(settings_min(CPAP_EPAPHi))).
arg(validPressure(settings_max(CPAP_IPAPLo))).
arg(validPressure(settings_max(CPAP_IPAPHi))).
arg(units);
// qDebug() << "AVAPS mode:" << retStr;
return retStr;
}
return STR_TR_Unknown;
}
EventDataType Day::calc(ChannelID code, ChannelCalcType type)
{
EventDataType value;
switch(type) {
case Calc_Min:
value = Min(code);
break;
case Calc_Middle:
value = calcMiddle(code);
break;
case Calc_Perc:
value = calcPercentile(code);
break;
case Calc_Max:
value = calcMax(code);
break;
case Calc_UpperThresh:
value = schema::channel[code].upperThreshold();
break;
case Calc_LowerThresh:
value = schema::channel[code].lowerThreshold();
break;
case Calc_Zero:
default:
value = 0;
break;
};
return value;
}
Reference in New Issue View Git Blame Copy Permalink