/* * Copyright (c) 2010 Damien Grauser ([email protected]) * * 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 2 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, write to the Free Software Foundation, Inc., 51 * Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ #include <QMessageBox> #include "HrZones.h" #include "Colors.h" #include "TimeUtils.h" #include <QtGui> #include <QtAlgorithms> #include <qcolor.h> #include <assert.h> #include <cmath> // the infinity endpoints are indicated with extreme date ranges // but not zero dates so we can edit and compare them static const QDate date_zero(1900, 01, 01); static const QDate date_infinity(9999,12,31); // initialize default static zone parameters void HrZones::initializeZoneParameters() { static int initial_zone_default[] = { 0, 68, 83, 94, 105 }; static double initial_zone_default_trimp[] = { 0.9, 1.1, 1.2, 2.0, 5.0 }; static const QString initial_zone_default_desc[] = { tr("Active Recovery"), tr("Endurance"), tr("Tempo"), tr("Threshold"), tr("VO2Max") }; static const char *initial_zone_default_name[] = { "Z1", "Z2", "Z3", "Z4", "Z5" }; static int initial_nzones_default = sizeof(initial_zone_default) / sizeof(initial_zone_default[0]); if (run) { fileName_ = "run-hr.zones"; } else { fileName_ = "hr.zones"; } scheme.zone_default.clear(); scheme.zone_default_is_pct.clear(); scheme.zone_default_desc.clear(); scheme.zone_default_name.clear(); scheme.zone_default_trimp.clear(); scheme.nzones_default = 0; scheme.nzones_default = initial_nzones_default; for (int z = 0; z < scheme.nzones_default; z ++) { scheme.zone_default.append(initial_zone_default[z]); scheme.zone_default_is_pct.append(true); scheme.zone_default_name.append(QString(initial_zone_default_name[z])); scheme.zone_default_desc.append(QString(initial_zone_default_desc[z])); scheme.zone_default_trimp.append(initial_zone_default_trimp[z]); } } // read zone file, allowing for zones with or without end dates bool HrZones::read(QFile &file) { // // GET SET // defaults_from_user = false; scheme.zone_default.clear(); scheme.zone_default_is_pct.clear(); scheme.zone_default_name.clear(); scheme.zone_default_desc.clear(); scheme.zone_default_trimp.clear(); scheme.nzones_default = 0; ranges.clear(); // set up possible warning dialog warning = QString(); int warning_lines = 0; const int max_warning_lines = 100; // macro to append lines to the warning #define append_to_warning(s) \ if (warning_lines < max_warning_lines) \ warning.append(s); \ else if (warning_lines == max_warning_lines) \ warning.append("...\n"); \ warning_lines ++; // read using text mode takes care of end-lines if (!file.open(QIODevice::ReadOnly | QIODevice::Text)) { err = tr("can't open %1 file").arg(fileName_); return false; } QTextStream fileStream(&file); QRegExp commentrx("\\s*#.*$"); QRegExp blankrx("^[ \t]*$"); QRegExp rangerx[] = { QRegExp("^\\s*(?:from\\s+)?" // optional "from" "((\\d\\d\\d\\d)[-/](\\d{1,2})[-/](\\d{1,2})|BEGIN)" // begin date "\\s*([,:]?\\s*(LT)\\s*=\\s*(\\d+))?" // optional {LT = integer (optional %)} "\\s*([,:]?\\s*(RestHr)\\s*=\\s*(\\d+))?" // optional {RestHr = integer (optional %)} "\\s*([,:]?\\s*(MaxHr)\\s*=\\s*(\\d+))?" // optional {MaxHr = integer (optional %)} "\\s*:?\\s*$", // optional : Qt::CaseInsensitive), QRegExp("^\\s*(?:from\\s+)?" // optional "from" "((\\d\\d\\d\\d)[-/](\\d{1,2})[-/](\\d{1,2})|BEGIN)" // begin date "\\s+(?:until|to|-)\\s+" // until "((\\d\\d\\d\\d)[-/](\\d{1,2})[-/](\\d{1,2})|END)?" // end date "\\s*:?,?\\s*((LT)\\s*=\\s*(\\d+))?" // optional {LT = integer (optional %)} "\\s*:?,?\\s*((RestHr)\\s*=\\s*(\\d+))?" // optional {RestHr = integer (optional %)} "\\s*:?,?\\s*((MaxHr)\\s*=\\s*(\\d+))?" // optional {MaxHr = integer (optional %)} "\\s*:?\\s*$", // optional : Qt::CaseInsensitive) }; QRegExp zonerx("^\\s*([^ ,][^,]*),\\s*([^ ,][^,]*),\\s*" "(\\d+)\\s*(%?)\\s*(?:,\\s*(\\d+(\\.\\d+)?)\\s*)?$", Qt::CaseInsensitive);// QRegExp zonedefaultsx("^\\s*(?:zone)?\\s*defaults?\\s*:?\\s*$", Qt::CaseInsensitive); int lineno = 0; // the current range in the file // ZoneRange *range = NULL; bool in_range = false; QDate begin = date_zero, end = date_infinity; int lt = 0; int restHr = 0; int maxHr = 0; QList<HrZoneInfo> zoneInfos; // true if zone defaults are found in the file (then we need to write them) bool zones_are_defaults = false; // // READ IN hr.zones FILE // // loop through line by line while (!fileStream.atEnd()) { // starting from line 1 ++lineno; // get a line in QString line = fileStream.readLine(); int pos = commentrx.indexIn(line, 0); // strip comments if (pos != -1) line = line.left(pos); // its a blank line (we check after comments stripped) if (blankrx.indexIn(line, 0) == 0) goto next_line; // who wrote this? bleck. // check for default zone range definition (may be followed by hr zone definitions) if (zonedefaultsx.indexIn(line, 0) != -1) { zones_are_defaults = true; // defaults are allowed only at the beginning of the file if (ranges.size()) { err = tr("HR Zone defaults must be specified at head of %1 file").arg(fileName_); return false; } // only one set of defaults is allowed if (scheme.nzones_default) { err = tr("Only one set of zone defaults may be specified in %1 file").arg(fileName_); return false; } // ok move on to get defaults setup goto next_line; } // check for range specification (may be followed by zone definitions) for (int r=0; r<2; r++) { if (rangerx[r].indexIn(line, 0) != -1) { if (in_range) { // if zones are empty, then generate them HrZoneRange range(begin, end, lt, restHr, maxHr); range.zones = zoneInfos; if (range.zones.empty()) { if (range.lt > 0) setHrZonesFromLT(range); else { err = tr("line %1: read new range without reading " "any zones for previous one").arg(lineno); file.close(); return false; } } else { qSort(range.zones); } ranges.append(range); } in_range = true; zones_are_defaults = false; zoneInfos.clear(); // process the beginning date if (rangerx[r].cap(1) == "BEGIN") begin = date_zero; else { begin = QDate(rangerx[r].cap(2).toInt(), rangerx[r].cap(3).toInt(), rangerx[r].cap(4).toInt()); } // process an end date, if any, else it is null if (rangerx[r].cap(5) == "END") end = date_infinity; else if (rangerx[r].cap(6).toInt() || rangerx[r].cap(7).toInt() || rangerx[r].cap(8).toInt()) { end = QDate(rangerx[r].cap(6).toInt(), rangerx[r].cap(7).toInt(), rangerx[r].cap(8).toInt()); } else { end = QDate(); } // set up the range, capturing LT if it's specified // range = new ZoneRange(begin, end); int nLT = (r ? 11 : 7); if (rangerx[r].captureCount() >= (nLT)) lt = rangerx[r].cap(nLT).toInt(); else lt = 0; int nRestHr = (r ? 14 : 10); if (rangerx[r].captureCount() >= (nRestHr)) restHr = rangerx[r].cap(nRestHr).toInt(); else restHr = 0; int nMaxHr = (r ? 17 : 13); if (rangerx[r].captureCount() >= (nRestHr)) maxHr = rangerx[r].cap(nMaxHr).toInt(); else maxHr = 0; // bleck goto next_line; } } // check for zone definition if (zonerx.indexIn(line, 0) != -1) { if (! (in_range || zones_are_defaults)) { err = tr("line %1: read zone without preceding date range").arg(lineno); file.close(); return false; } int lo = zonerx.cap(3).toInt(); double trimp = zonerx.cap(5).toDouble(); // allow for zone specified as % of LT bool lo_is_pct = false; if (zonerx.cap(4) == "%") { if (zones_are_defaults) lo_is_pct = true; else if (lt > 0) lo = int(lo * lt / 100); else { err = tr("attempt to set zone based on % of LT without setting LT in line number %1.\n"). arg(lineno); file.close(); return false; } } int hi = -1; // signal an undefined number double tr = zonerx.cap(5).toDouble(); if (zones_are_defaults) { scheme.nzones_default ++; scheme.zone_default_is_pct.append(lo_is_pct); scheme.zone_default.append(lo); scheme.zone_default_name.append(zonerx.cap(1)); scheme.zone_default_desc.append(zonerx.cap(2)); scheme.zone_default_trimp.append(trimp); defaults_from_user = true; } else { HrZoneInfo zone(zonerx.cap(1), zonerx.cap(2), lo, hi, tr); zoneInfos.append(zone); } } next_line: ; } // did we drop out mid way through ? if (in_range) { HrZoneRange range(begin, end, lt, restHr, maxHr); range.zones = zoneInfos; if (range.zones.empty()) { if (range.lt > 0) setHrZonesFromLT(range); else { err = tr("file ended without reading any zones for last range"); file.close(); return false; } } else { qSort(range.zones); } ranges.append(range); } // reading done file.close(); // sort the ranges qSort(ranges); // // POST-PROCESS / FIX-UP ZONES // // set the default zones if not in file if (!scheme.nzones_default) { // do we have a zone which is explicitly set? for (int i=0; i<ranges.count(); i++) { if (ranges[i].hrZonesSetFromLT == false) { // set the defaults using this one! scheme.nzones_default = ranges[i].zones.count(); for (int j=0; j<scheme.nzones_default; j++) { scheme.zone_default.append(((double)ranges[i].zones[j].lo / (double)ranges[i].lt) * 100.00); scheme.zone_default_is_pct.append(true); scheme.zone_default_name.append(ranges[i].zones[j].name); scheme.zone_default_desc.append(ranges[i].zones[j].desc); scheme.zone_default_trimp.append(ranges[i].zones[j].trimp); } } } // still not set then reset to defaults as usual if (!scheme.nzones_default) initializeZoneParameters(); } // resolve undefined endpoints in ranges and zones for (int nr = 0; nr < ranges.size(); nr ++) { // clean up gaps or overlaps in zone ranges if (ranges[nr].end.isNull()) { ranges[nr].end = (nr < ranges.size() - 1) ? ranges[nr + 1].begin : date_infinity; } else if ((nr < ranges.size() - 1) && (ranges[nr + 1].begin != ranges[nr].end)) { append_to_warning(tr("Setting end date of range %1 to start date of range %2.\n").arg(nr + 1).arg(nr + 2)); ranges[nr].end = ranges[nr + 1].begin; } else if ((nr == ranges.size() - 1) && (ranges[nr].end < QDate::currentDate())) { append_to_warning(tr("Extending final range %1 to infinite to include present date.\n").arg(nr + 1)); ranges[nr].end = date_infinity; } if (ranges[nr].lt <= 0) { err = tr("LT must be greater than zero in zone range %1 of %2").arg(nr + 1).arg(fileName_); return false; } if (ranges[nr].zones.size()) { // check that the first zone starts with zero // ranges[nr].zones[0].lo = 0; // there is no reason we should enforce this, so removing it. // resolve zone end powers for (int nz = 0; nz < ranges[nr].zones.size(); nz ++) { if (ranges[nr].zones[nz].hi == -1) ranges[nr].zones[nz].hi = (nz < ranges[nr].zones.size() - 1) ? ranges[nr].zones[nz + 1].lo : INT_MAX; else if ((nz < ranges[nr].zones.size() - 1) && (ranges[nr].zones[nz].hi != ranges[nr].zones[nz + 1].lo)) { if (abs(ranges[nr].zones[nz].hi - ranges[nr].zones[nz + 1].lo) > 4) { append_to_warning(tr("Range %1: matching top of zone %2 (%3) to bottom of zone %4 (%5).\n") .arg(nr+1).arg(ranges[nr].zones[nz].name).arg(ranges[nr].zones[nz].hi) .arg(ranges[nr].zones[nz + 1].name) .arg(ranges[nr].zones[nz + 1].lo)); } ranges[nr].zones[nz].hi = ranges[nr].zones[nz + 1].lo; } else if ((nz == ranges[nr].zones.size() - 1) && (ranges[nr].zones[nz].hi < INT_MAX)) { append_to_warning(tr("Range %1: setting top of zone %2 from %3 to MAX.\n") .arg(nr + 1).arg(ranges[nr].zones[nz].name).arg(ranges[nr].zones[nz].hi)); ranges[nr].zones[nz].hi = INT_MAX; } } } } // mark zones as modified so pages which depend on zones can be updated modificationTime = QDateTime::currentDateTime(); return true; } // note empty dates are treated as automatic matches for begin or // end of range int HrZones::whichRange(const QDate &date) const { for (int rnum = 0; rnum < ranges.size(); ++rnum) { const HrZoneRange &range = ranges[rnum]; if (((date >= range.begin) || (range.begin.isNull())) && ((date < range.end) || (range.end.isNull()))) return rnum; } return -1; } int HrZones::numZones(int rnum) const { if (rnum < 0 || rnum >= ranges.size()) return 0; return ranges[rnum].zones.size(); } int HrZones::whichZone(int rnum, double value) const { if (rnum < 0 || rnum > ranges.size()) return 0; const HrZoneRange &range = ranges[rnum]; for (int j = 0; j < range.zones.size(); ++j) { const HrZoneInfo &info = range.zones[j]; // note: the "end" of range is actually in the next zone if ((value >= info.lo) && (value < info.hi)) return j; } // if we got here either it is negative, nan, inf or way high return -1; } void HrZones::zoneInfo(int rnum, int znum, QString &name, QString &description, int &low, int &high, double &trimp) const { assert(rnum < ranges.size()); const HrZoneRange &range = ranges[rnum]; assert(znum < range.zones.size()); const HrZoneInfo &zone = range.zones[znum]; name = zone.name; description = zone.desc; low = zone.lo; high = zone.hi; trimp= zone.trimp; } int HrZones::getLT(int rnum) const { if (rnum < 0 || rnum > ranges.size()) return 0; return ranges[rnum].lt; } void HrZones::setLT(int rnum, int lt) { ranges[rnum].lt = lt; modificationTime = QDateTime::currentDateTime(); } // generate a list of zones from LT int HrZones::lowsFromLT(QList <int> *lows, int lt) const { lows->clear(); for (int z = 0; z < scheme.nzones_default; z++) lows->append(scheme.zone_default_is_pct[z] ? scheme.zone_default[z] * lt / 100 : scheme.zone_default[z]); return scheme.nzones_default; } int HrZones::getRestHr(int rnum) const { if (rnum < 0 || rnum > ranges.size()) return 0; return ranges[rnum].restHr; } void HrZones::setRestHr(int rnum, int restHr) { ranges[rnum].restHr = restHr; modificationTime = QDateTime::currentDateTime(); } int HrZones::getMaxHr(int rnum) const { if (rnum < 0 || rnum > ranges.size()) return 0; return ranges[rnum].maxHr; } void HrZones::setMaxHr(int rnum, int maxHr) { ranges[rnum].maxHr = maxHr; modificationTime = QDateTime::currentDateTime(); } // access the zone name QString HrZones::getDefaultZoneName(int z) const { return scheme.zone_default_name[z]; } // access the zone description QString HrZones::getDefaultZoneDesc(int z) const { return scheme.zone_default_desc[z]; } // set the zones from the LT value (the cp variable) void HrZones::setHrZonesFromLT(HrZoneRange &range) { range.zones.clear(); if (scheme.nzones_default == 0) initializeZoneParameters(); for (int i = 0; i < scheme.nzones_default; i++) { int lo = scheme.zone_default_is_pct[i] ? scheme.zone_default[i] * range.lt / 100 : scheme.zone_default[i]; int hi = lo; double trimp = scheme.zone_default_trimp[i]; HrZoneInfo zone(scheme.zone_default_name[i], scheme.zone_default_desc[i], lo, hi, trimp); range.zones.append(zone); } // sort the zones (some may be pct, others absolute, so zones need to be sorted, // rather than the defaults qSort(range.zones); // set zone end dates for (int i = 0; i < range.zones.size(); i ++) range.zones[i].hi = (i < scheme.nzones_default - 1) ? range.zones[i + 1].lo : INT_MAX; // mark that the zones were set from LT, so if zones are subsequently // written, only LT is saved range.hrZonesSetFromLT = true; } void HrZones::setHrZonesFromLT(int rnum) { assert((rnum >= 0) && (rnum < ranges.size())); setHrZonesFromLT(ranges[rnum]); } // return the list of starting values of zones for a given range QList <int> HrZones::getZoneLows(int rnum) const { if (rnum >= ranges.size()) return QList <int>(); const HrZoneRange &range = ranges[rnum]; QList <int> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].lo); return return_values; } // return the list of ending values of zones for a given range QList <int> HrZones::getZoneHighs(int rnum) const { if (rnum >= ranges.size()) return QList <int>(); const HrZoneRange &range = ranges[rnum]; QList <int> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].hi); return return_values; } // return the list of zone names QList <QString> HrZones::getZoneNames(int rnum) const { if (rnum >= ranges.size()) return QList <QString>(); const HrZoneRange &range = ranges[rnum]; QList <QString> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].name); return return_values; } // return the list of zone trimp coef QList <double> HrZones::getZoneTrimps(int rnum) const { if (rnum >= ranges.size()) return QList <double>(); const HrZoneRange &range = ranges[rnum]; QList <double> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].trimp); return return_values; } QString HrZones::summarize(int rnum, QVector<double> &time_in_zone, QColor color) const { assert(rnum < ranges.size()); const HrZoneRange &range = ranges[rnum]; if (time_in_zone.size() < range.zones.size()) return ""; QString summary; if(range.lt > 0){ summary += "<table align=\"center\" width=\"70%\" border=\"0\">"; summary += "<tr><td align=\"center\">"; summary += tr("Threshold (bpm): %1").arg(range.lt); summary += "</td></tr></table>"; } summary += "<table align=\"center\" width=\"70%\" "; summary += "border=\"0\">"; summary += "<tr>"; summary += tr("<td align=\"center\">Zone</td>"); summary += tr("<td align=\"center\">Description</td>"); summary += tr("<td align=\"center\">Low (bpm)</td>"); summary += tr("<td align=\"center\">High (bpm)</td>"); summary += tr("<td align=\"center\">Time</td>"); summary += tr("<td align=\"center\">%</td>"); summary += "</tr>"; double duration = 0; foreach(double v, time_in_zone) { duration += v; } for (int zone = 0; zone < time_in_zone.size(); ++zone) { if (time_in_zone[zone] > 0.0) { QString name, desc; int lo, hi; double trimp; zoneInfo(rnum, zone, name, desc, lo, hi, trimp); if (zone % 2 == 0) summary += "<tr bgcolor='" + color.name() + "'>"; else summary += "<tr>"; summary += QString("<td align=\"center\">%1</td>").arg(name); summary += QString("<td align=\"center\">%1</td>").arg(desc); summary += QString("<td align=\"center\">%1</td>").arg(lo); if (hi == INT_MAX) summary += "<td align=\"center\">MAX</td>"; else summary += QString("<td align=\"center\">%1</td>").arg(hi); summary += QString("<td align=\"center\">%1</td>") .arg(time_to_string((unsigned) round(time_in_zone[zone]))); summary += QString("<td align=\"center\">%1</td>") .arg((double)time_in_zone[zone]/duration * 100, 0, 'f', 0); summary += "</tr>"; } } summary += "</table>"; return summary; } #define USE_SHORT_POWER_ZONES_FORMAT true /* whether a less redundent format should be used */ void HrZones::write(QDir home) { QString strzones; // always write the defaults (config pane can adjust) strzones += QString("DEFAULTS:\n"); for (int z = 0 ; z < scheme.nzones_default; z ++) strzones += QString("%1,%2,%3%4,%5\n"). arg(scheme.zone_default_name[z]). arg(scheme.zone_default_desc[z]). arg(scheme.zone_default[z]). arg(scheme.zone_default_is_pct[z]?"%":""). arg(scheme.zone_default_trimp[z]); strzones += QString("\n"); for (int i = 0; i < ranges.size(); i++) { int lt = getLT(i); int restHr = getRestHr(i); int maxHr = getMaxHr(i); // print header for range // note this explicitly sets the first and last ranges such that all time is spanned // note: BEGIN is not needed anymore // since it becomes Jan 01 1900 strzones += QString("%1: LT=%2, RestHr=%3, MaxHr=%4").arg(getStartDate(i).toString("yyyy/MM/dd")).arg(lt).arg(restHr).arg(maxHr); strzones += QString("\n"); // step through and print the zones if they've been explicitly set if (! ranges[i].hrZonesSetFromLT) { for (int j = 0; j < ranges[i].zones.size(); j ++) { const HrZoneInfo &zi = ranges[i].zones[j]; strzones += QString("%1,%2,%3,%4\n").arg(zi.name).arg(zi.desc).arg(zi.lo).arg(zi.trimp); } strzones += QString("\n"); } } QFile file(home.canonicalPath() + "/" + fileName_); if (file.open(QFile::WriteOnly)) { QTextStream stream(&file); stream << strzones; file.close(); } else { QMessageBox msgBox; msgBox.setIcon(QMessageBox::Critical); msgBox.setText(tr("Problem Saving Heartrate Zones")); msgBox.setInformativeText(tr("File: %1 cannot be opened for 'Writing'. Please check file properties.").arg(home.canonicalPath() + "/" + fileName_)); msgBox.exec(); return; } }
/* * Copyright (c) 2010 Damien Grauser ([email protected]) * * 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 2 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, write to the Free Software Foundation, Inc., 51 * Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ #include <QMessageBox> #include "HrZones.h" #include "Colors.h" #include "TimeUtils.h" #include <QtGui> #include <QtAlgorithms> #include <qcolor.h> #include <assert.h> #include <math.h> #include <boost/crc.hpp> // the infinity endpoints are indicated with extreme date ranges // but not zero dates so we can edit and compare them static const QDate date_zero(1900, 01, 01); static const QDate date_infinity(9999,12,31); // initialize default static zone parameters void HrZones::initializeZoneParameters() { static int initial_zone_default[] = { 0, 68, 83, 94, 105 }; static double initial_zone_default_trimp[] = { 0.9, 1.1, 1.2, 2.0, 5.0 }; static const QString initial_zone_default_desc[] = { tr("Active Recovery"), tr("Endurance"), tr("Tempo"), tr("Threshold"), tr("VO2Max") }; static const char *initial_zone_default_name[] = { "Z1", "Z2", "Z3", "Z4", "Z5" }; static int initial_nzones_default = sizeof(initial_zone_default) / sizeof(initial_zone_default[0]); scheme.zone_default.clear(); scheme.zone_default_is_pct.clear(); scheme.zone_default_desc.clear(); scheme.zone_default_name.clear(); scheme.zone_default_trimp.clear(); scheme.nzones_default = 0; scheme.nzones_default = initial_nzones_default; for (int z = 0; z < scheme.nzones_default; z ++) { scheme.zone_default.append(initial_zone_default[z]); scheme.zone_default_is_pct.append(true); scheme.zone_default_name.append(QString(initial_zone_default_name[z])); scheme.zone_default_desc.append(QString(initial_zone_default_desc[z])); scheme.zone_default_trimp.append(initial_zone_default_trimp[z]); } } // read zone file, allowing for zones with or without end dates bool HrZones::read(QFile &file) { defaults_from_user = false; scheme.zone_default.clear(); scheme.zone_default_is_pct.clear(); scheme.zone_default_name.clear(); scheme.zone_default_desc.clear(); scheme.zone_default_trimp.clear(); scheme.nzones_default = 0; ranges.clear(); // set up possible warning dialog warning = QString(); int warning_lines = 0; const int max_warning_lines = 100; // macro to append lines to the warning #define append_to_warning(s) \ if (warning_lines < max_warning_lines) \ warning.append(s); \ else if (warning_lines == max_warning_lines) \ warning.append("...\n"); \ warning_lines ++; // read using text mode takes care of end-lines if (!file.open(QIODevice::ReadOnly | QIODevice::Text)) { err = "can't open file"; return false; } QTextStream fileStream(&file); QRegExp commentrx("\\s*#.*$"); QRegExp blankrx("^[ \t]*$"); QRegExp rangerx[] = { QRegExp("^\\s*(?:from\\s+)?" // optional "from" "((\\d\\d\\d\\d)[-/](\\d{1,2})[-/](\\d{1,2})|BEGIN)" // begin date "\\s*([,:]?\\s*(LT)\\s*=\\s*(\\d+))?" // optional {LT = integer (optional %)} "\\s*([,:]?\\s*(RestHr)\\s*=\\s*(\\d+))?" // optional {RestHr = integer (optional %)} "\\s*([,:]?\\s*(MaxHr)\\s*=\\s*(\\d+))?" // optional {MaxHr = integer (optional %)} "\\s*:?\\s*$", // optional : Qt::CaseInsensitive), QRegExp("^\\s*(?:from\\s+)?" // optional "from" "((\\d\\d\\d\\d)[-/](\\d{1,2})[-/](\\d{1,2})|BEGIN)" // begin date "\\s+(?:until|to|-)\\s+" // until "((\\d\\d\\d\\d)[-/](\\d{1,2})[-/](\\d{1,2})|END)?" // end date "\\s*:?,?\\s*((LT)\\s*=\\s*(\\d+))?" // optional {LT = integer (optional %)} "\\s*:?,?\\s*((RestHr)\\s*=\\s*(\\d+))?" // optional {RestHr = integer (optional %)} "\\s*:?,?\\s*((MaxHr)\\s*=\\s*(\\d+))?" // optional {MaxHr = integer (optional %)} "\\s*:?\\s*$", // optional : Qt::CaseInsensitive) }; QRegExp zonerx("^\\s*([^ ,][^,]*),\\s*([^ ,][^,]*),\\s*" "(\\d+)\\s*(%?)\\s*(?:,\\s*(\\d+(\\.\\d+)?)\\s*)?$", Qt::CaseInsensitive);// QRegExp zonedefaultsx("^\\s*(?:zone)?\\s*defaults?\\s*:?\\s*$", Qt::CaseInsensitive); int lineno = 0; // the current range in the file // ZoneRange *range = NULL; bool in_range = false; QDate begin = date_zero, end = date_infinity; int lt = 0; int restHr = 0; int maxHr = 0; QList<HrZoneInfo> zoneInfos; // true if zone defaults are found in the file (then we need to write them) bool zones_are_defaults = false; while (! fileStream.atEnd() ) { ++lineno; QString line = fileStream.readLine(); int pos = commentrx.indexIn(line, 0); if (pos != -1) line = line.left(pos); if (blankrx.indexIn(line, 0) == 0) goto next_line; // check for default zone range definition (may be followed by hr zone definitions) if (zonedefaultsx.indexIn(line, 0) != -1) { zones_are_defaults = true; // defaults are allowed only at the beginning of the file if (ranges.size()) { err = "HR Zone defaults must be specified at head of hr.zones file"; return false; } // only one set of defaults is allowed if (scheme.nzones_default) { err = "Only one set of zone defaults may be specified in hr.zones file"; return false; } goto next_line; } // check for range specification (may be followed by zone definitions) for (int r=0; r<2; r++) { if (rangerx[r].indexIn(line, 0) != -1) { if (in_range) { // if zones are empty, then generate them HrZoneRange range(begin, end, lt, restHr, maxHr); range.zones = zoneInfos; if (range.zones.empty()) { if (range.lt > 0) setHrZonesFromLT(range); else { err = tr("line %1: read new range without reading " "any zones for previous one").arg(lineno); file.close(); return false; } } else { qSort(range.zones); } ranges.append(range); } in_range = true; zones_are_defaults = false; zoneInfos.clear(); // process the beginning date if (rangerx[r].cap(1) == "BEGIN") begin = date_zero; else { begin = QDate(rangerx[r].cap(2).toInt(), rangerx[r].cap(3).toInt(), rangerx[r].cap(4).toInt()); } // process an end date, if any, else it is null if (rangerx[r].cap(5) == "END") end = date_infinity; else if (rangerx[r].cap(6).toInt() || rangerx[r].cap(7).toInt() || rangerx[r].cap(8).toInt()) { end = QDate(rangerx[r].cap(6).toInt(), rangerx[r].cap(7).toInt(), rangerx[r].cap(8).toInt()); } else { end = QDate(); } // set up the range, capturing LT if it's specified // range = new ZoneRange(begin, end); int nLT = (r ? 11 : 7); if (rangerx[r].numCaptures() >= (nLT)) lt = rangerx[r].cap(nLT).toInt(); else lt = 0; int nRestHr = (r ? 14 : 10); if (rangerx[r].numCaptures() >= (nRestHr)) restHr = rangerx[r].cap(nRestHr).toInt(); else restHr = 0; int nMaxHr = (r ? 17 : 13); if (rangerx[r].numCaptures() >= (nRestHr)) maxHr = rangerx[r].cap(nMaxHr).toInt(); else maxHr = 0; // bleck goto next_line; } } // check for zone definition if (zonerx.indexIn(line, 0) != -1) { if (! (in_range || zones_are_defaults)) { err = tr("line %1: read zone without " "preceeding date range").arg(lineno); file.close(); return false; } int lo = zonerx.cap(3).toInt(); double trimp = zonerx.cap(5).toDouble(); // allow for zone specified as % of LT bool lo_is_pct = false; if (zonerx.cap(4) == "%") { if (zones_are_defaults) lo_is_pct = true; else if (lt > 0) lo = int(lo * lt / 100); else { err = tr("attempt to set zone based on % of " "LT without setting LT in line number %1.\n"). arg(lineno); file.close(); return false; } } int hi = -1; // signal an undefined number double tr = zonerx.cap(5).toDouble(); if (zones_are_defaults) { scheme.nzones_default ++; scheme.zone_default_is_pct.append(lo_is_pct); scheme.zone_default.append(lo); scheme.zone_default_name.append(zonerx.cap(1)); scheme.zone_default_desc.append(zonerx.cap(2)); scheme.zone_default_trimp.append(trimp); defaults_from_user = true; } else { HrZoneInfo zone(zonerx.cap(1), zonerx.cap(2), lo, hi, tr); zoneInfos.append(zone); } } next_line: {} } if (in_range) { HrZoneRange range(begin, end, lt, restHr, maxHr); range.zones = zoneInfos; if (range.zones.empty()) { if (range.lt > 0) setHrZonesFromLT(range); else { err = tr("file ended without reading any zones for last range"); file.close(); return false; } } else { qSort(range.zones); } ranges.append(range); } file.close(); // sort the ranges qSort(ranges); // set the default zones if not in file if (!scheme.nzones_default) { // do we have a zone which is explicitly set? for (int i=0; i<ranges.count(); i++) { if (ranges[i].hrZonesSetFromLT == false) { // set the defaults using this one! scheme.nzones_default = ranges[i].zones.count(); for (int j=0; j<scheme.nzones_default; j++) { scheme.zone_default.append(((double)ranges[i].zones[j].lo / (double)ranges[i].lt) * 100.00); scheme.zone_default_is_pct.append(true); scheme.zone_default_name.append(ranges[i].zones[j].name); scheme.zone_default_desc.append(ranges[i].zones[j].desc); scheme.zone_default_trimp.append(ranges[i].zones[j].trimp); } } } // still not set then reset to defaults as usual if (!scheme.nzones_default) initializeZoneParameters(); } // resolve undefined endpoints in ranges and zones for (int nr = 0; nr < ranges.size(); nr ++) { // clean up gaps or overlaps in zone ranges if (ranges[nr].end.isNull()) ranges[nr].end = (nr < ranges.size() - 1) ? ranges[nr + 1].begin : date_infinity; else if ((nr < ranges.size() - 1) && (ranges[nr + 1].begin != ranges[nr].end)) { append_to_warning(tr("Setting end date of range %1 " "to start date of range %2.\n"). arg(nr + 1). arg(nr + 2) ); ranges[nr].end = ranges[nr + 1].begin; } else if ((nr == ranges.size() - 1) && (ranges[nr].end < QDate::currentDate())) { append_to_warning(tr("Extending final range %1 to infinite " "to include present date.\n").arg(nr + 1)); ranges[nr].end = date_infinity; } if (ranges[nr].lt <= 0) { err = tr("LT must be greater than zero in zone " "range %1 of hr.zones").arg(nr + 1); return false; } if (ranges[nr].zones.size()) { // check that the first zone starts with zero ranges[nr].zones[0].lo = 0; // resolve zone end powers for (int nz = 0; nz < ranges[nr].zones.size(); nz ++) { if (ranges[nr].zones[nz].hi == -1) ranges[nr].zones[nz].hi = (nz < ranges[nr].zones.size() - 1) ? ranges[nr].zones[nz + 1].lo : INT_MAX; else if ((nz < ranges[nr].zones.size() - 1) && (ranges[nr].zones[nz].hi != ranges[nr].zones[nz + 1].lo)) { if (abs(ranges[nr].zones[nz].hi - ranges[nr].zones[nz + 1].lo) > 4) { append_to_warning(tr("Range %1: matching top of zone %2 " "(%3) to bottom of zone %4 (%5).\n"). arg(nr + 1). arg(ranges[nr].zones[nz].name). arg(ranges[nr].zones[nz].hi). arg(ranges[nr].zones[nz + 1].name). arg(ranges[nr].zones[nz + 1].lo) ); } ranges[nr].zones[nz].hi = ranges[nr].zones[nz + 1].lo; } else if ((nz == ranges[nr].zones.size() - 1) && (ranges[nr].zones[nz].hi < INT_MAX)) { append_to_warning(tr("Range %1: setting top of zone %2 from %3 to MAX.\n"). arg(nr + 1). arg(ranges[nr].zones[nz].name). arg(ranges[nr].zones[nz].hi) ); ranges[nr].zones[nz].hi = INT_MAX; } } } } // mark zones as modified so pages which depend on zones can be updated modificationTime = QDateTime::currentDateTime(); return true; } // note empty dates are treated as automatic matches for begin or // end of range int HrZones::whichRange(const QDate &date) const { for (int rnum = 0; rnum < ranges.size(); ++rnum) { const HrZoneRange &range = ranges[rnum]; if (((date >= range.begin) || (range.begin.isNull())) && ((date < range.end) || (range.end.isNull()))) return rnum; } return -1; } int HrZones::numZones(int rnum) const { if (rnum < 0 || rnum >= ranges.size()) return 0; return ranges[rnum].zones.size(); } int HrZones::whichZone(int rnum, double value) const { if (rnum < 0 || rnum > ranges.size()) return 0; const HrZoneRange &range = ranges[rnum]; for (int j = 0; j < range.zones.size(); ++j) { const HrZoneInfo &info = range.zones[j]; // note: the "end" of range is actually in the next zone if ((value >= info.lo) && (value < info.hi)) return j; } // if we got here either it is negative, nan, inf or way high if (value < 0 || isnan(value)) return 0; else return range.zones.size()-1; } void HrZones::zoneInfo(int rnum, int znum, QString &name, QString &description, int &low, int &high, double &trimp) const { assert(rnum < ranges.size()); const HrZoneRange &range = ranges[rnum]; assert(znum < range.zones.size()); const HrZoneInfo &zone = range.zones[znum]; name = zone.name; description = zone.desc; low = zone.lo; high = zone.hi; trimp= zone.trimp; } int HrZones::getLT(int rnum) const { if (rnum < 0 || rnum > ranges.size()) return 0; return ranges[rnum].lt; } void HrZones::setLT(int rnum, int lt) { ranges[rnum].lt = lt; modificationTime = QDateTime::currentDateTime(); } // generate a list of zones from LT int HrZones::lowsFromLT(QList <int> *lows, int lt) const { lows->clear(); for (int z = 0; z < scheme.nzones_default; z++) lows->append(scheme.zone_default_is_pct[z] ? scheme.zone_default[z] * lt / 100 : scheme.zone_default[z]); return scheme.nzones_default; } int HrZones::getRestHr(int rnum) const { if (rnum < 0 || rnum > ranges.size()) return 0; return ranges[rnum].restHr; } void HrZones::setRestHr(int rnum, int restHr) { ranges[rnum].restHr = restHr; modificationTime = QDateTime::currentDateTime(); } int HrZones::getMaxHr(int rnum) const { if (rnum < 0 || rnum > ranges.size()) return 0; return ranges[rnum].maxHr; } void HrZones::setMaxHr(int rnum, int maxHr) { ranges[rnum].maxHr = maxHr; modificationTime = QDateTime::currentDateTime(); } // access the zone name QString HrZones::getDefaultZoneName(int z) const { return scheme.zone_default_name[z]; } // access the zone description QString HrZones::getDefaultZoneDesc(int z) const { return scheme.zone_default_desc[z]; } // set the zones from the LT value (the cp variable) void HrZones::setHrZonesFromLT(HrZoneRange &range) { range.zones.clear(); if (scheme.nzones_default == 0) initializeZoneParameters(); for (int i = 0; i < scheme.nzones_default; i++) { int lo = scheme.zone_default_is_pct[i] ? scheme.zone_default[i] * range.lt / 100 : scheme.zone_default[i]; int hi = lo; double trimp = scheme.zone_default_trimp[i]; HrZoneInfo zone(scheme.zone_default_name[i], scheme.zone_default_desc[i], lo, hi, trimp); range.zones.append(zone); } // sort the zones (some may be pct, others absolute, so zones need to be sorted, // rather than the defaults qSort(range.zones); // set zone end dates for (int i = 0; i < range.zones.size(); i ++) range.zones[i].hi = (i < scheme.nzones_default - 1) ? range.zones[i + 1].lo : INT_MAX; // mark that the zones were set from LT, so if zones are subsequently // written, only LT is saved range.hrZonesSetFromLT = true; } void HrZones::setHrZonesFromLT(int rnum) { assert((rnum >= 0) && (rnum < ranges.size())); setHrZonesFromLT(ranges[rnum]); } // return the list of starting values of zones for a given range QList <int> HrZones::getZoneLows(int rnum) const { if (rnum >= ranges.size()) return QList <int>(); const HrZoneRange &range = ranges[rnum]; QList <int> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].lo); return return_values; } // return the list of ending values of zones for a given range QList <int> HrZones::getZoneHighs(int rnum) const { if (rnum >= ranges.size()) return QList <int>(); const HrZoneRange &range = ranges[rnum]; QList <int> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].hi); return return_values; } // return the list of zone names QList <QString> HrZones::getZoneNames(int rnum) const { if (rnum >= ranges.size()) return QList <QString>(); const HrZoneRange &range = ranges[rnum]; QList <QString> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].name); return return_values; } // return the list of zone trimp coef QList <double> HrZones::getZoneTrimps(int rnum) const { if (rnum >= ranges.size()) return QList <double>(); const HrZoneRange &range = ranges[rnum]; QList <double> return_values; for (int i = 0; i < range.zones.size(); i ++) return_values.append(ranges[rnum].zones[i].trimp); return return_values; } QString HrZones::summarize(int rnum, QVector<double> &time_in_zone) const { assert(rnum < ranges.size()); const HrZoneRange &range = ranges[rnum]; if (time_in_zone.size() < range.zones.size()) return ""; QString summary; if(range.lt > 0){ summary += "<table align=\"center\" width=\"70%\" border=\"0\">"; summary += "<tr><td align=\"center\">"; summary += tr("Threshold (bpm): %1").arg(range.lt); summary += "</td></tr></table>"; } summary += "<table align=\"center\" width=\"70%\" "; summary += "border=\"0\">"; summary += "<tr>"; summary += tr("<td align=\"center\">Zone</td>"); summary += tr("<td align=\"center\">Description</td>"); summary += tr("<td align=\"center\">Low (bpm)</td>"); summary += tr("<td align=\"center\">High (bpm)</td>"); summary += tr("<td align=\"center\">Time</td>"); summary += tr("<td align=\"center\">%</td>"); summary += "</tr>"; QColor color = QApplication::palette().alternateBase().color(); color = QColor::fromHsv(color.hue(), color.saturation() * 2, color.value()); double duration = 0; foreach(double v, time_in_zone) { duration += v; } for (int zone = 0; zone < time_in_zone.size(); ++zone) { if (time_in_zone[zone] > 0.0) { QString name, desc; int lo, hi; double trimp; zoneInfo(rnum, zone, name, desc, lo, hi, trimp); if (zone % 2 == 0) summary += "<tr bgcolor='" + color.name() + "'>"; else summary += "<tr>"; summary += QString("<td align=\"center\">%1</td>").arg(name); summary += QString("<td align=\"center\">%1</td>").arg(desc); summary += QString("<td align=\"center\">%1</td>").arg(lo); if (hi == INT_MAX) summary += "<td align=\"center\">MAX</td>"; else summary += QString("<td align=\"center\">%1</td>").arg(hi); summary += QString("<td align=\"center\">%1</td>") .arg(time_to_string((unsigned) round(time_in_zone[zone]))); summary += QString("<td align=\"center\">%1</td>") .arg((double)time_in_zone[zone]/duration * 100, 0, 'f', 0); summary += "</tr>"; } } summary += "</table>"; return summary; } #define USE_SHORT_POWER_ZONES_FORMAT true /* whether a less redundent format should be used */ void HrZones::write(QDir home) { QString strzones; // always write the defaults (config pane can adjust) strzones += QString("DEFAULTS:\n"); for (int z = 0 ; z < scheme.nzones_default; z ++) strzones += QString("%1,%2,%3%4,%5\n"). arg(scheme.zone_default_name[z]). arg(scheme.zone_default_desc[z]). arg(scheme.zone_default[z]). arg(scheme.zone_default_is_pct[z]?"%":""). arg(scheme.zone_default_trimp[z]); strzones += QString("\n"); for (int i = 0; i < ranges.size(); i++) { int lt = getLT(i); int restHr = getRestHr(i); int maxHr = getMaxHr(i); // print header for range // note this explicitly sets the first and last ranges such that all time is spanned // note: BEGIN is not needed anymore // since it becomes Jan 01 1900 strzones += QString("%1: LT=%2, RestHr=%3, MaxHr=%4").arg(getStartDate(i).toString("yyyy/MM/dd")).arg(lt).arg(restHr).arg(maxHr); strzones += QString("\n"); // step through and print the zones if they've been explicitly set if (! ranges[i].hrZonesSetFromLT) { for (int j = 0; j < ranges[i].zones.size(); j ++) { const HrZoneInfo &zi = ranges[i].zones[j]; strzones += QString("%1,%2,%3,%4\n").arg(zi.name).arg(zi.desc).arg(zi.lo).arg(zi.trimp); } strzones += QString("\n"); } } QFile file(home.absolutePath() + "/hr.zones"); if (file.open(QFile::WriteOnly)) { QTextStream stream(&file); stream << strzones; file.close(); } }
