/*
 * Copyright (c) 2010 Damien Grauser (Damien.Grauser@pev-geneve.ch)
 *
 * 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 <QRegExp>
#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 (sport_.isEmpty() || sport_ == "Bike") {
        fileName_ = "hr.zones";
    } else {
        fileName_ = sport_.toLower() + "-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);
#if QT_VERSION < 0x060000
    fileStream.setCodec("UTF-8");
#endif

    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*((AeT)\\s*=\\s*(\\d+))?"                // optional {AeT = 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*((AeT)\\s*=\\s*(\\d+))?"                // optional {AeT = 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 aet = 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, aet, 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 {
                        std::sort(range.zones.begin(), range.zones.end());
                    }
                    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;

                int nAeT = (r ? 20 : 16);
                if (rangerx[r].captureCount() >= (nAeT)) aet = rangerx[r].cap(nAeT).toInt();
                else aet = 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, aet, 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 {

            std::sort(range.zones.begin(), range.zones.end());
        }
        ranges.append(range);
    }

    // reading done
    file.close();

    // sort the ranges
    std::sort(ranges.begin(), ranges.end());

    //
    // 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::getAeT(int rnum) const
{
    if (rnum < 0 || rnum > ranges.size()) return 0;
    if (ranges[rnum].aet > 0) return ranges[rnum].aet;
    return round(0.9 * ranges[rnum].lt); // Default to 90% LTHR when not set
}

void HrZones::setAeT(int rnum, int aet)
{
    ranges[rnum].aet = aet;
    modificationTime = QDateTime::currentDateTime();
}

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
    std::sort(range.zones.begin(), range.zones.end());

    // 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 aet = getAeT(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, AeT=%5").arg(getStartDate(i).toString("yyyy/MM/dd")).arg(lt).arg(restHr).arg(maxHr).arg(aet);
        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);
#if QT_VERSION < 0x060000
        stream.setCodec("UTF-8");
#endif
        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;
    }
}

void HrZones::addHrZoneRange(QDate _start, QDate _end, int _lt, int _aet, int _restHr, int _maxHr)
{
    ranges.append(HrZoneRange(_start, _end, _lt, _aet, _restHr, _maxHr));
}

// insert a new zone range using the current scheme
// return the range number
int HrZones::addHrZoneRange(QDate _start, int _lt, int _aet, int _restHr, int _maxHr)
{
    int rnum;

    // where to add this range?
    for(rnum=0; rnum < ranges.count(); rnum++) if (ranges[rnum].begin > _start) break;

    // at the end ?
    if (rnum == ranges.count()) ranges.append(HrZoneRange(_start, date_infinity, _lt, _aet, _restHr, _maxHr));
    else ranges.insert(rnum, HrZoneRange(_start, ranges[rnum].begin, _lt, _aet, _restHr, _maxHr));

    // modify previous end date
    if (rnum) ranges[rnum-1].end = _start;

    // set zones from LT
    if (_lt > 0) {
        setLT(rnum, _lt);
        setHrZonesFromLT(rnum);
    }

    return rnum;
}

void HrZones::addHrZoneRange()
{
    ranges.append(HrZoneRange(date_zero, date_infinity));
}

void HrZones::setEndDate(int rnum, QDate endDate)
{
    ranges[rnum].end = endDate;
    modificationTime = QDateTime::currentDateTime();
}
void HrZones::setStartDate(int rnum, QDate startDate)
{
    ranges[rnum].begin = startDate;
    modificationTime = QDateTime::currentDateTime();
}

QDate HrZones::getStartDate(int rnum) const
{
    assert(rnum >= 0);
    return ranges[rnum].begin;
}

QString HrZones::getStartDateString(int rnum) const
{
    assert(rnum >= 0);
    QDate d = ranges[rnum].begin;
    return (d.isNull() ? "BEGIN" : d.toString());
}

QDate HrZones::getEndDate(int rnum) const
{
    assert(rnum >= 0);
    return ranges[rnum].end;
}

QString HrZones::getEndDateString(int rnum) const
{
    assert(rnum >= 0);
    QDate d = ranges[rnum].end;
    return (d.isNull() ? "END" : d.toString());
}

int HrZones::getRangeSize() const
{
    return ranges.size();
}

// generate a zone color with a specific number of zones
QColor hrZoneColor(int z, int) {
    switch(z) {

    case 0  : return GColor(CHZONE1); break;
    case 1  : return GColor(CHZONE2); break;
    case 2  : return GColor(CHZONE3); break;
    case 3  : return GColor(CHZONE4); break;
    case 4  : return GColor(CHZONE5); break;
    case 5  : return GColor(CHZONE6); break;
    case 6  : return GColor(CHZONE7); break;
    case 7  : return GColor(CHZONE8); break;
    case 8  : return GColor(CHZONE9); break;
    case 9  : return GColor(CHZONE10); break;
    default: return QColor(128,128,128); break;
    }
}

// delete a range, extend an adjacent (prior if available, otherwise next)
// range to cover the same time period, then return the number of the new range
// covering the date range of the deleted range or -1 if none left
int HrZones::deleteRange(int rnum) {

    // check bounds - silently fail, don't assert
    assert (rnum < ranges.count() && rnum >= 0);

    // extend the previous range to the end of this range
    // but only if we have a previous range
    if (rnum > 0) setEndDate(rnum-1, getEndDate(rnum));

    // delete this range then
    ranges.removeAt(rnum);

    return rnum-1;
}

// insert a new range starting at the given date extending to the end of the zone currently
// containing that date.  If the start date of that zone is prior to the specified start
// date, then that zone range is shorted.
int HrZones::insertRangeAtDate(QDate date, int lt) {
    assert(date.isValid());
    int rnum;

    if (ranges.empty()) {
        addHrZoneRange();
        rnum = 0;
    }
    else {
        rnum = whichRange(date);
        assert(rnum >= 0);
        QDate date1 = getStartDate(rnum);

        // if the old range has dates before the specified, then truncate
        // the old range and shift up the existing ranges
        if (date > date1) {
            QDate endDate = getEndDate(rnum);
            setEndDate(rnum, date);
            ranges.insert(++ rnum, HrZoneRange(date, endDate));
        }
    }

    if (lt > 0) {
        setLT(rnum, lt);
        setHrZonesFromLT(rnum);
    }

    return rnum;
}

quint16
HrZones::getFingerprint() const
{
    quint64 x = 0;
    for (int i=0; i<ranges.size(); i++) {

        // from
        x += ranges[i].begin.toJulianDay();

        // to
        x += ranges[i].end.toJulianDay();

        // LT
        x += ranges[i].lt;

        // AeT
        x += ranges[i].aet;

        // Rest HR
        x += ranges[i].restHr;

        // Max HR
        x += ranges[i].maxHr;

        // each zone definition (manual edit/default changed)
        for (int j=0; j<ranges[i].zones.count(); j++) {
            x += ranges[i].zones[j].lo;
            x += int(double(100.0f * ranges[i].zones[j].trimp));
        }
    }
    QByteArray ba = QByteArray::number(x);
#if QT_VERSION < 0x060000
    return qChecksum(ba, ba.length());
#else
    return qChecksum(ba);
#endif
}

quint16
HrZones::getFingerprint(QDate forDate) const
{
    quint64 x = 0;

    int i = whichRange(forDate);
    if (i >= 0) {

        // zone parameters...
        x += ranges[i].lt;
        x += ranges[i].aet;
        x += ranges[i].restHr;
        x += ranges[i].maxHr;

        // each zone definition (manual edit/default changed)
        for (int j=0; j<ranges[i].zones.count(); j++) {
            x += ranges[i].zones[j].lo;
            x += int(double(100.0f * ranges[i].zones[j].trimp));
        }
    }
    QByteArray ba = QByteArray::number(x);
#if QT_VERSION < 0x060000
    return qChecksum(ba, ba.length());
#else
    return qChecksum(ba);
#endif
}
