nmelone/GoldenCheetah
1
0
Fork 0
mirror of https://github.com/GoldenCheetah/GoldenCheetah.git synced 2026-02-16 01:19:57 +00:00
Code Issues Packages Projects Releases Wiki Activity

Anaerobic TISS

Browse Source 
Added Anaerobic TISS from Dr Skiba to complement
the Aerobic TISS added recently.

Added to;
* Metrics
* Ride Plot
* PMC metrics
This commit is contained in:
Mark Liversedge
2014-03-22 10:47:44 +00:00
parent bc47e0a2d1
commit 191219bfa4
12 changed files with 322 additions and 29 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
src/RideFile.cpp
View File

@@ -668,6 +668,7 @@ RideFile::isDataPresent(SeriesType series)
case watts : return dataPresent.watts; break;
case aPower : return dataPresent.apower; break;
case aTISS : return dataPresent.atiss; break;
case anTISS : return dataPresent.antiss; break;
case alt : return dataPresent.alt; break;
case lon : return dataPresent.lon; break;
case lat : return dataPresent.lat; break;
@@ -733,6 +734,7 @@ RideFilePoint::value(RideFile::SeriesType series) const
case RideFile::xPower : return xp; break;
case RideFile::aPower : return apower; break;
case RideFile::aTISS : return atiss; break;
case RideFile::anTISS : return antiss; break;
default:
case RideFile::none : break;
@@ -1041,9 +1043,15 @@ RideFile::recalculateDerivedSeries()
static const double b = -7.5095428451195f;
static const double c = -0.86118031563782f;
static const double t = 2;
// anTISS
static const double an = 0.238923886004611f;
static const double bn = -12.2066385296127f;
static const double cn = -1.73549567522521f;
int CP = 0;
int WPRIME = 0;
double aTISS = 0.0f;
double anTISS = 0.0f;
// set WPrime and CP
if (context->athlete->zones()) {
@@ -1188,12 +1196,14 @@ RideFile::recalculateDerivedSeries()
APtotal += p->apower;
APcount++;
// Aerobic TISS
// Anaerobic and Aerobic TISS
if (CP && dataPresent.watts) {
// a * exp (b * exp (c * fraction of cp) )
aTISS += recIntSecs_ * (a * exp(b * exp(c * (double(p->watts) / double(CP)))));
anTISS += recIntSecs_ * (an * exp(bn * exp(cn * (double(p->watts) / double(CP)))));
p->atiss = aTISS;
p->antiss = anTISS;
}
// last point