@glandais/vcyclist-engine 3.1.1 → 4.0.0

package/vcyclist-engine.js

@@ -39,9 +39,9 @@
   var toString = kotlin_kotlin.$_$.f7;
   var hashCode = kotlin_kotlin.$_$.m6;
   var getBooleanHashCode = kotlin_kotlin.$_$.i6;
-  var CoroutineImpl = kotlin_kotlin.$_$.x4;
   var IllegalArgumentException_init_$Create$ = kotlin_kotlin.$_$.a2;
   var Unit_instance = kotlin_kotlin.$_$.g1;
+  var CoroutineImpl = kotlin_kotlin.$_$.x4;
   var PointPerDistance_instance = kotlin_io_github_glandais_gpx.$_$.k;
   var ElevationStep_instance = kotlin_io_github_glandais_gpx.$_$.h;
   var get_COROUTINE_SUSPENDED = kotlin_kotlin.$_$.i4;
@@ -71,8 +71,9 @@
   var Pair = kotlin_kotlin.$_$.p9;
   var listOf = kotlin_kotlin.$_$.p3;
   var Default_getInstance = kotlin_kotlin.$_$.b1;
+  var isFinite = kotlin_kotlin.$_$.ca;
   var PointField_TEMPERATURE_getInstance = kotlin_io_github_glandais_gpx.$_$.c;
-  var isNaN_0 = kotlin_kotlin.$_$.da;
+  var isNaN_0 = kotlin_kotlin.$_$.ea;
   var get_entries = kotlin_io_github_glandais_gpx.$_$.u;
   var Path = kotlin_io_github_glandais_gpx.$_$.t;
   var GpxParser_instance = kotlin_io_github_glandais_gpx.$_$.d;
@@ -121,6 +122,7 @@
   initMetadataForCompanion(Companion_1);
   initMetadataForClass(Cyclist, 'Cyclist', Cyclist);
   initMetadataForClass(SimplifyPathOptions, 'SimplifyPathOptions', SimplifyPathOptions);
+  initMetadataForClass(WPrimeBalanceOptions, 'WPrimeBalanceOptions', WPrimeBalanceOptions);
   initMetadataForCompanion(Companion_2);
   initMetadataForClass(EnhanceOptions, 'EnhanceOptions', EnhanceOptions);
   initMetadataForCoroutine($enhanceCourseCOROUTINE$, CoroutineImpl);
@@ -145,17 +147,19 @@
   initMetadataForObject(MuscularPowerProvider, 'MuscularPowerProvider');
   initMetadataForObject(PowerComputer, 'PowerComputer');
   initMetadataForClass(PowerProviderConstant, 'PowerProviderConstant', VOID, CyclistPowerProviderBase);
-  initMetadataForClass(PowerProviderConstantWithTiring, 'PowerProviderConstantWithTiring', VOID, CyclistPowerProviderBase);
+  initMetadataForCompanion(Companion_6);
+  initMetadataForClass(PowerProviderDurability, 'PowerProviderDurability', VOID, CyclistPowerProviderBase);
   initMetadataForObject(PowerProviderFromData, 'PowerProviderFromData');
   initMetadataForObject(RhoProviderDefault, 'RhoProviderDefault');
   initMetadataForObject(RhoProviderEstimate, 'RhoProviderEstimate');
   initMetadataForObject(RollingResistancePowerProvider, 'RollingResistancePowerProvider');
   initMetadataForObject(VirtualizeService, 'VirtualizeService');
   initMetadataForObject(WheelBearingsPowerProvider, 'WheelBearingsPowerProvider');
-  initMetadataForCompanion(Companion_6);
+  initMetadataForCompanion(Companion_7);
   initMetadataForClass(Wind, 'Wind');
   initMetadataForObject(WindProviderNone, 'WindProviderNone');
   initMetadataForClass(WindProviderConstant, 'WindProviderConstant');
+  initMetadataForObject(WPrimeBalanceComputer, 'WPrimeBalanceComputer');
   initMetadataForLambda(enhance$slambda, CoroutineImpl, VOID, [1]);
   initMetadataForLambda(enhanceWithCourse$slambda, CoroutineImpl, VOID, [1]);
   //endregion
@@ -169,21 +173,34 @@
       new Companion();
     return Companion_instance;
   }
-  function Bike(crr, inertiaFront, inertiaRear, wheelRadiusM, efficiency) {
+  function Bike(crr, inertiaFront, inertiaRear, wheelRadiusM, efficiency, maxPedalingLeanAngleDeg) {
     Companion_getInstance_2();
     crr = crr === VOID ? 0.004 : crr;
     inertiaFront = inertiaFront === VOID ? 0.05 : inertiaFront;
     inertiaRear = inertiaRear === VOID ? 0.07 : inertiaRear;
     wheelRadiusM = wheelRadiusM === VOID ? 0.35 : wheelRadiusM;
     efficiency = efficiency === VOID ? 0.976 : efficiency;
+    maxPedalingLeanAngleDeg = maxPedalingLeanAngleDeg === VOID ? 20.0 : maxPedalingLeanAngleDeg;
     this.crr_1 = crr;
     this.inertiaFront_1 = inertiaFront;
     this.inertiaRear_1 = inertiaRear;
     this.wheelRadiusM_1 = wheelRadiusM;
     this.efficiency_1 = efficiency;
+    this.maxPedalingLeanAngleDeg_1 = maxPedalingLeanAngleDeg;
   }
+  protoOf(Bike).get_tanMaxPedalingLeanAngle_nft7tz_k$ = function () {
+    var tmp;
+    if (this.maxPedalingLeanAngleDeg_1 >= 90.0) {
+      tmp = Infinity;
+    } else {
+      // Inline function 'kotlin.math.tan' call
+      var x = this.maxPedalingLeanAngleDeg_1 * 3.141592653589793 / 180.0;
+      tmp = Math.tan(x);
+    }
+    return tmp;
+  };
   protoOf(Bike).toString = function () {
-    return 'Bike(crr=' + this.crr_1 + ', inertiaFront=' + this.inertiaFront_1 + ', inertiaRear=' + this.inertiaRear_1 + ', wheelRadiusM=' + this.wheelRadiusM_1 + ', efficiency=' + this.efficiency_1 + ')';
+    return 'Bike(crr=' + this.crr_1 + ', inertiaFront=' + this.inertiaFront_1 + ', inertiaRear=' + this.inertiaRear_1 + ', wheelRadiusM=' + this.wheelRadiusM_1 + ', efficiency=' + this.efficiency_1 + ', maxPedalingLeanAngleDeg=' + this.maxPedalingLeanAngleDeg_1 + ')';
   };
   protoOf(Bike).hashCode = function () {
     var result = getNumberHashCode(this.crr_1);
@@ -191,6 +208,7 @@
     result = imul(result, 31) + getNumberHashCode(this.inertiaRear_1) | 0;
     result = imul(result, 31) + getNumberHashCode(this.wheelRadiusM_1) | 0;
     result = imul(result, 31) + getNumberHashCode(this.efficiency_1) | 0;
+    result = imul(result, 31) + getNumberHashCode(this.maxPedalingLeanAngleDeg_1) | 0;
     return result;
   };
   protoOf(Bike).equals = function (other) {
@@ -208,6 +226,8 @@
       return false;
     if (!equals(this.efficiency_1, other.efficiency_1))
       return false;
+    if (!equals(this.maxPedalingLeanAngleDeg_1, other.maxPedalingLeanAngleDeg_1))
+      return false;
     return true;
   };
   function Course(path, cyclist, bike) {
@@ -415,6 +435,46 @@
       return false;
     return true;
   };
+  function WPrimeBalanceOptions(enabled, criticalPowerW, wPrimeJ) {
+    enabled = enabled === VOID ? true : enabled;
+    criticalPowerW = criticalPowerW === VOID ? 250.0 : criticalPowerW;
+    wPrimeJ = wPrimeJ === VOID ? 20000.0 : wPrimeJ;
+    this.enabled_1 = enabled;
+    this.criticalPowerW_1 = criticalPowerW;
+    this.wPrimeJ_1 = wPrimeJ;
+    // Inline function 'kotlin.require' call
+    if (!(this.criticalPowerW_1 > 0.0)) {
+      var message = 'criticalPowerW must be > 0, got ' + this.criticalPowerW_1;
+      throw IllegalArgumentException_init_$Create$(toString(message));
+    }
+    // Inline function 'kotlin.require' call
+    if (!(this.wPrimeJ_1 > 0.0)) {
+      var message_0 = 'wPrimeJ must be > 0, got ' + this.wPrimeJ_1;
+      throw IllegalArgumentException_init_$Create$(toString(message_0));
+    }
+  }
+  protoOf(WPrimeBalanceOptions).toString = function () {
+    return 'WPrimeBalanceOptions(enabled=' + this.enabled_1 + ', criticalPowerW=' + this.criticalPowerW_1 + ', wPrimeJ=' + this.wPrimeJ_1 + ')';
+  };
+  protoOf(WPrimeBalanceOptions).hashCode = function () {
+    var result = getBooleanHashCode(this.enabled_1);
+    result = imul(result, 31) + getNumberHashCode(this.criticalPowerW_1) | 0;
+    result = imul(result, 31) + getNumberHashCode(this.wPrimeJ_1) | 0;
+    return result;
+  };
+  protoOf(WPrimeBalanceOptions).equals = function (other) {
+    if (this === other)
+      return true;
+    if (!(other instanceof WPrimeBalanceOptions))
+      return false;
+    if (!(this.enabled_1 === other.enabled_1))
+      return false;
+    if (!equals(this.criticalPowerW_1, other.criticalPowerW_1))
+      return false;
+    if (!equals(this.wPrimeJ_1, other.wPrimeJ_1))
+      return false;
+    return true;
+  };
   function Companion_2() {
     Companion_instance_2 = this;
     this.DEFAULT_1 = new EnhanceOptions();
@@ -425,32 +485,35 @@
       new Companion_2();
     return Companion_instance_2;
   }
-  function EnhanceOptions(fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath) {
+  function EnhanceOptions(fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath, wPrimeBalance) {
     Companion_getInstance_5();
     fixElevation = fixElevation === VOID ? true : fixElevation;
     computeMaxSpeeds = computeMaxSpeeds === VOID ? true : computeMaxSpeeds;
     virtualizeTrack = virtualizeTrack === VOID ? true : virtualizeTrack;
     computeOnePointPerSecond = computeOnePointPerSecond === VOID ? true : computeOnePointPerSecond;
     simplifyPath = simplifyPath === VOID ? new SimplifyPathOptions() : simplifyPath;
+    wPrimeBalance = wPrimeBalance === VOID ? new WPrimeBalanceOptions() : wPrimeBalance;
     this.fixElevation_1 = fixElevation;
     this.computeMaxSpeeds_1 = computeMaxSpeeds;
     this.virtualizeTrack_1 = virtualizeTrack;
     this.computeOnePointPerSecond_1 = computeOnePointPerSecond;
     this.simplifyPath_1 = simplifyPath;
+    this.wPrimeBalance_1 = wPrimeBalance;
   }
-  protoOf(EnhanceOptions).copy_8buhav_k$ = function (fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath) {
-    return new EnhanceOptions(fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath);
+  protoOf(EnhanceOptions).copy_bazpxz_k$ = function (fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath, wPrimeBalance) {
+    return new EnhanceOptions(fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath, wPrimeBalance);
   };
-  protoOf(EnhanceOptions).copy$default_9zz8ww_k$ = function (fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath, $super) {
+  protoOf(EnhanceOptions).copy$default_2rfmrp_k$ = function (fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath, wPrimeBalance, $super) {
     fixElevation = fixElevation === VOID ? this.fixElevation_1 : fixElevation;
     computeMaxSpeeds = computeMaxSpeeds === VOID ? this.computeMaxSpeeds_1 : computeMaxSpeeds;
     virtualizeTrack = virtualizeTrack === VOID ? this.virtualizeTrack_1 : virtualizeTrack;
     computeOnePointPerSecond = computeOnePointPerSecond === VOID ? this.computeOnePointPerSecond_1 : computeOnePointPerSecond;
     simplifyPath = simplifyPath === VOID ? this.simplifyPath_1 : simplifyPath;
-    return $super === VOID ? this.copy_8buhav_k$(fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath) : $super.copy_8buhav_k$.call(this, fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath);
+    wPrimeBalance = wPrimeBalance === VOID ? this.wPrimeBalance_1 : wPrimeBalance;
+    return $super === VOID ? this.copy_bazpxz_k$(fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath, wPrimeBalance) : $super.copy_bazpxz_k$.call(this, fixElevation, computeMaxSpeeds, virtualizeTrack, computeOnePointPerSecond, simplifyPath, wPrimeBalance);
   };
   protoOf(EnhanceOptions).toString = function () {
-    return 'EnhanceOptions(fixElevation=' + this.fixElevation_1 + ', computeMaxSpeeds=' + this.computeMaxSpeeds_1 + ', virtualizeTrack=' + this.virtualizeTrack_1 + ', computeOnePointPerSecond=' + this.computeOnePointPerSecond_1 + ', simplifyPath=' + this.simplifyPath_1.toString() + ')';
+    return 'EnhanceOptions(fixElevation=' + this.fixElevation_1 + ', computeMaxSpeeds=' + this.computeMaxSpeeds_1 + ', virtualizeTrack=' + this.virtualizeTrack_1 + ', computeOnePointPerSecond=' + this.computeOnePointPerSecond_1 + ', simplifyPath=' + this.simplifyPath_1.toString() + ', wPrimeBalance=' + this.wPrimeBalance_1.toString() + ')';
   };
   protoOf(EnhanceOptions).hashCode = function () {
     var result = getBooleanHashCode(this.fixElevation_1);
@@ -458,6 +521,7 @@
     result = imul(result, 31) + getBooleanHashCode(this.virtualizeTrack_1) | 0;
     result = imul(result, 31) + getBooleanHashCode(this.computeOnePointPerSecond_1) | 0;
     result = imul(result, 31) + this.simplifyPath_1.hashCode() | 0;
+    result = imul(result, 31) + this.wPrimeBalance_1.hashCode() | 0;
     return result;
   };
   protoOf(EnhanceOptions).equals = function (other) {
@@ -475,6 +539,8 @@
       return false;
     if (!this.simplifyPath_1.equals(other.simplifyPath_1))
       return false;
+    if (!this.wPrimeBalance_1.equals(other.wPrimeBalance_1))
+      return false;
     return true;
   };
   function $enhanceCourseCOROUTINE$(_this__u8e3s4, course, options, elevationProvider, resultContinuation) {
@@ -532,6 +598,10 @@
               this.path0__1 = PointPerSecond_instance.computeOnePointPerSecond_7sah9a_k$(this.path0__1);
             }
 
+            if (this.options_1.wPrimeBalance_1.enabled_1) {
+              WPrimeBalanceComputer_instance.compute_ewvoaz_k$(this.path0__1, this.options_1.wPrimeBalance_1);
+            }
+
             if (this.options_1.simplifyPath_1.enabled_1) {
               this.path0__1 = PathSimplifier_instance.simplify_gxhj7b_k$(this.path0__1, this.options_1.simplifyPath_1.toleranceM_1, this.options_1.simplifyPath_1.zExaggeration_1);
             }
@@ -557,7 +627,7 @@
     return new CoursePhysics(new Course(path), RhoProviderEstimate_instance, AeroProviderConstant_instance, WindProviderNone_instance, new PowerProviderConstant(280.0));
   };
   protoOf(Enhancer).enhanceCourseDefault_jm937y_k$ = function (path, elevationProvider, options, $completion) {
-    return this.enhanceCourse_mout3j_k$(this.getDefaultCourse_gqi1xw_k$(path), options.copy$default_9zz8ww_k$(options.fixElevation_1 && !(elevationProvider == null)), elevationProvider, $completion);
+    return this.enhanceCourse_mout3j_k$(this.getDefaultCourse_gqi1xw_k$(path), options.copy$default_2rfmrp_k$(options.fixElevation_1 && !(elevationProvider == null)), elevationProvider, $completion);
   };
   protoOf(Enhancer).enhanceCourse_mout3j_k$ = function (course, options, elevationProvider, $completion) {
     var tmp = new $enhanceCourseCOROUTINE$(this, course, options, elevationProvider, $completion);
@@ -1354,9 +1424,15 @@
   function computeCorneringLimit($this, course, currentIndex, window_0) {
     var path = course.path_1;
     var radius = computeRadiusWindowed($this, path, currentIndex, window_0);
-    // Inline function 'kotlin.math.sqrt' call
-    var x = 9.80665 * radius * course.cyclist_1.get_tanMaxLeanAngle_qqiag1_k$();
-    var vMax = Math.sqrt(x);
+    var tmp;
+    if (radius >= 200.0) {
+      tmp = course.cyclist_1.get_maxSpeedMS_xz2r7y_k$();
+    } else {
+      // Inline function 'kotlin.math.sqrt' call
+      var x = 9.80665 * radius * course.cyclist_1.get_tanMaxLeanAngle_qqiag1_k$();
+      tmp = Math.sqrt(x);
+    }
+    var vMax = tmp;
     // Inline function 'kotlin.math.min' call
     var a = course.cyclist_1.get_maxSpeedMS_xz2r7y_k$();
     var result = Math.min(a, vMax);
@@ -1387,14 +1463,46 @@
     path.setRadius_p1cueq_k$(i, clamped);
     return clamped;
   }
-  function computeBrakingLimit($this, course, currentIndex, nextIndex) {
+  function computeBrakingLimit($this, course, currentIndex, nextIndex, corneringLimit) {
     var path = course.path_1;
     var vf = path.speedMax_mjo2cl_k$(nextIndex);
-    var a = course.cyclist_1.get_maxBrakeMS2_nzouxw_k$();
+    var aXMax = course.cyclist_1.get_maxBrakeMS2_nzouxw_k$();
     var d = path.distance_5s4s77_k$(nextIndex) - path.distance_5s4s77_k$(currentIndex);
+    var radius = path.radius_84q0sm_k$(currentIndex);
+    if (radius >= 200.0 || !isFinite(radius) || radius <= 0.0) {
+      // Inline function 'kotlin.math.sqrt' call
+      var x = vf * vf + 2.0 * aXMax * d;
+      // Inline function 'kotlin.math.min' call
+      var b = Math.sqrt(x);
+      return Math.min(corneringLimit, b);
+    }
+    var aYMax = 9.80665 * course.cyclist_1.get_tanMaxLeanAngle_qqiag1_k$();
+    if (computeBrakingLimit$slack(radius, aYMax, aXMax, vf, d, corneringLimit) >= 0.0)
+      return corneringLimit;
+    // Inline function 'kotlin.math.min' call
+    var low = Math.min(vf, corneringLimit);
+    var high = corneringLimit;
+    // Inline function 'kotlin.repeat' call
+    var inductionVariable = 0;
+    if (inductionVariable < 24)
+      do {
+        var index = inductionVariable;
+        inductionVariable = inductionVariable + 1 | 0;
+        var mid = (low + high) / 2.0;
+        if (computeBrakingLimit$slack(radius, aYMax, aXMax, vf, d, mid) >= 0.0)
+          low = mid;
+        else
+          high = mid;
+      }
+       while (inductionVariable < 24);
+    return low;
+  }
+  function computeBrakingLimit$slack(radius, aYMax, aXMax, vf, d, v) {
+    var lateralFraction = coerceIn(v * v / radius / aYMax, 0.0, 1.0);
     // Inline function 'kotlin.math.sqrt' call
-    var x = vf * vf + 2.0 * a * d;
-    return Math.sqrt(x);
+    var x = 1.0 - lateralFraction * lateralFraction;
+    var aX = aXMax * Math.sqrt(x);
+    return vf * vf + 2.0 * aX * d - v * v;
   }
   function MaxSpeedComputer() {
     this.MAX_RADIUS_M_1 = 200.0;
@@ -1402,6 +1510,7 @@
     this.END_SPEED_MS_1 = 2.0;
     this.BEARING_THRESHOLD_1 = 0.001;
     this.DEFAULT_WINDOW_1 = 10;
+    this.ELLIPSE_ITERATIONS_1 = 24;
   }
   protoOf(MaxSpeedComputer).computeMaxSpeeds_jc4v4o_k$ = function (course) {
     var path = course.path_1;
@@ -1415,10 +1524,7 @@
           path.setSpeedMax_7xvqir_k$(i, 2.0);
         } else {
           var cornering = computeCorneringLimit(this, course, i, 10);
-          var braking = computeBrakingLimit(this, course, i, i + 1 | 0);
-          // Inline function 'kotlin.math.min' call
-          var tmp$ret$0 = Math.min(cornering, braking);
-          path.setSpeedMax_7xvqir_k$(i, tmp$ret$0);
+          path.setSpeedMax_7xvqir_k$(i, computeBrakingLimit(this, course, i, i + 1 | 0, cornering));
         }
       }
        while (0 <= inductionVariable);
@@ -1439,12 +1545,23 @@
   function MuscularPowerProvider() {
   }
   protoOf(MuscularPowerProvider).powerAt_ejucaq_k$ = function (course, path, pointIndex) {
-    var muscular = course.cyclistPowerProvider_1.powerAt_ejucaq_k$(course, path, pointIndex);
+    var intent = course.cyclistPowerProvider_1.powerAt_ejucaq_k$(course, path, pointIndex);
+    var muscular = this.pedalsClear_b5t5mh_k$(course, path, pointIndex) ? intent : 0.0;
     path.setPCyclistProvidedMuscular_aori9q_k$(pointIndex, muscular);
     var wheel = muscular * course.get_bike_wojgoa_k$().efficiency_1;
     path.setPCyclistProvidedWheel_ofsihz_k$(pointIndex, wheel);
     return wheel;
   };
+  protoOf(MuscularPowerProvider).pedalsClear_b5t5mh_k$ = function (course, path, pointIndex) {
+    var radius = path.radius_84q0sm_k$(pointIndex);
+    if (!isFinite(radius) || radius <= 0.0)
+      return true;
+    var speed = path.speed_npav2n_k$(pointIndex);
+    if (!isFinite(speed))
+      return true;
+    var tanLean = speed * speed / (9.80665 * radius);
+    return tanLean <= course.get_bike_wojgoa_k$().get_tanMaxPedalingLeanAngle_nft7tz_k$();
+  };
   var MuscularPowerProvider_instance;
   function MuscularPowerProvider_getInstance() {
     return MuscularPowerProvider_instance;
@@ -1452,7 +1569,9 @@
   function getDtInner($this, pSum, equivalentMass, currentSpeed, dx, dt1Init, dt2Init) {
     var dt1 = dt1Init;
     var dt2 = dt2Init;
-    var tol = dx / 1.0E7;
+    // Inline function 'kotlin.math.max' call
+    var a = dx / 1.0E7;
+    var tol = Math.max(a, 1.0E-12);
     while (dt2 - dt1 >= tol) {
       var dtMid = (dt1 + dt2) / 2.0;
       var dxMid = $this.getDx_xtszre_k$(pSum, equivalentMass, currentSpeed, dtMid);
@@ -1464,6 +1583,7 @@
     return (dt1 + dt2) / 2.0;
   }
   function PowerComputer() {
+    this.MIN_DT_TOLERANCE_1 = 1.0E-12;
   }
   protoOf(PowerComputer).getNewPower_puvq9i_k$ = function (course, path, pointIndex, withCyclist) {
     var pSum = 0.0;
@@ -1477,12 +1597,14 @@
     return pSum;
   };
   protoOf(PowerComputer).getDx_xtszre_k$ = function (pSum, equivalentMass, currentSpeed, dt) {
+    // Inline function 'kotlin.math.max' call
+    var b = dt * pSum / (0.5 * equivalentMass) + currentSpeed * currentSpeed;
     // Inline function 'kotlin.math.sqrt' call
-    var x = dt * pSum / (0.5 * equivalentMass) + currentSpeed * currentSpeed;
+    var x = Math.max(0.0, b);
     var tmp0 = Math.sqrt(x);
     // Inline function 'kotlin.math.max' call
-    var b = 0.5555555555555556;
-    var newSpeed = Math.max(tmp0, b);
+    var b_0 = 0.5555555555555556;
+    var newSpeed = Math.max(tmp0, b_0);
     return (currentSpeed + newSpeed) * dt / 2.0;
   };
   protoOf(PowerComputer).getTotPower_kfbciq_k$ = function (equivalentMass, s1, s2, dt) {
@@ -1497,17 +1619,28 @@
   };
   protoOf(PowerComputer).computeCyclistPower_kgcsa3_k$ = function (course, path, equivalentMass, i) {
     if (i === 0) {
+      path.setPComputedPower_c3m6yw_k$(i, 0.0);
+      path.setPBrake_iza4qj_k$(i, 0.0);
+      return Unit_instance;
+    }
+    var dtSeconds = path.dt_17uzp4_k$(i) / 1000.0;
+    if (dtSeconds <= 0.0) {
+      path.setPComputedTotalPower_wm6di_k$(i, 0.0);
+      path.setPComputedWheelPower_qplog3_k$(i, 0.0);
+      path.setPBrake_iza4qj_k$(i, 0.0);
       path.setPComputedPower_c3m6yw_k$(i, 0.0);
       return Unit_instance;
     }
     var resistive = this.getNewPower_puvq9i_k$(course, path, i - 1 | 0, false);
     var s1 = path.speed_npav2n_k$(i - 1 | 0);
     var s2 = path.speed_npav2n_k$(i);
-    var dtSeconds = path.dt_17uzp4_k$(i) / 1000.0;
     var totPower = this.getTotPower_kfbciq_k$(equivalentMass, s1, s2, dtSeconds);
     path.setPComputedTotalPower_wm6di_k$(i, totPower);
     var powerWheel = totPower - resistive;
     path.setPComputedWheelPower_qplog3_k$(i, powerWheel);
+    // Inline function 'kotlin.math.min' call
+    var tmp$ret$0 = Math.min(0.0, powerWheel);
+    path.setPBrake_iza4qj_k$(i, tmp$ret$0);
     // Inline function 'kotlin.math.max' call
     var computed = Math.max(0.0, powerWheel) / course.get_bike_wojgoa_k$().efficiency_1;
     path.setPComputedPower_c3m6yw_k$(i, computed);
@@ -1533,24 +1666,86 @@
   protoOf(PowerProviderConstant).optimalPower_36kg19_k$ = function (course, path, pointIndex) {
     return this.power_1;
   };
-  function PowerProviderConstantWithTiring(power, useHarmonics, durationSeconds, random) {
+  function Companion_6() {
+    this.DEFAULT_DECLINE_PER_KJ_PER_KG_1 = 0.006666666666666667;
+    this.DEFAULT_MAX_DECLINE_1 = 0.2;
+  }
+  var Companion_instance_6;
+  function Companion_getInstance_9() {
+    return Companion_instance_6;
+  }
+  function PowerProviderDurability(powerW, criticalPowerW, declinePerKjPerKg, maxDecline, useHarmonics, random) {
+    criticalPowerW = criticalPowerW === VOID ? 250.0 : criticalPowerW;
+    declinePerKjPerKg = declinePerKjPerKg === VOID ? 0.006666666666666667 : declinePerKjPerKg;
+    maxDecline = maxDecline === VOID ? 0.2 : maxDecline;
     useHarmonics = useHarmonics === VOID ? false : useHarmonics;
     random = random === VOID ? Default_getInstance() : random;
     CyclistPowerProviderBase.call(this, useHarmonics, random);
-    this.power_1 = power;
-    this.durationSeconds_1 = durationSeconds;
+    this.powerW_1 = powerW;
+    this.criticalPowerW_1 = criticalPowerW;
+    this.declinePerKjPerKg_1 = declinePerKjPerKg;
+    this.maxDecline_1 = maxDecline;
     // Inline function 'kotlin.require' call
-    if (!(this.durationSeconds_1 > 0.0)) {
-      var message = 'durationSeconds must be > 0, got ' + this.durationSeconds_1;
+    if (!(this.powerW_1 > 0.0)) {
+      var message = 'powerW must be > 0, got ' + this.powerW_1;
       throw IllegalArgumentException_init_$Create$(toString(message));
     }
-  }
-  protoOf(PowerProviderConstantWithTiring).optimalPower_36kg19_k$ = function (course, path, pointIndex) {
-    var elapsedS = path.elapsed_7307rw_k$(pointIndex) / 1000.0;
-    // Inline function 'kotlin.math.max' call
-    var b = 1.0 - 0.6 * elapsedS / this.durationSeconds_1;
-    var c = Math.max(0.5, b);
-    return this.power_1 * c;
+    // Inline function 'kotlin.require' call
+    if (!(this.criticalPowerW_1 > 0.0)) {
+      var message_0 = 'criticalPowerW must be > 0, got ' + this.criticalPowerW_1;
+      throw IllegalArgumentException_init_$Create$(toString(message_0));
+    }
+    // Inline function 'kotlin.require' call
+    if (!(this.declinePerKjPerKg_1 >= 0.0)) {
+      var message_1 = 'declinePerKjPerKg must be >= 0, got ' + this.declinePerKjPerKg_1;
+      throw IllegalArgumentException_init_$Create$(toString(message_1));
+    }
+    var containsArg = this.maxDecline_1;
+    // Inline function 'kotlin.require' call
+    if (!(0.0 <= containsArg ? containsArg <= 1.0 : false)) {
+      var message_2 = 'maxDecline must be in [0, 1], got ' + this.maxDecline_1;
+      throw IllegalArgumentException_init_$Create$(toString(message_2));
+    }
+    this.lastIndex_1 = -1;
+    this.lastElapsedS_1 = 0.0;
+    this.lastPowerW_1 = 0.0;
+    this.supraCriticalWorkJ_1 = 0.0;
+    this.massKgForDose_1 = 80.0;
+  }
+  protoOf(PowerProviderDurability).get_decline_qrtukt_k$ = function () {
+    var tmp0 = this.maxDecline_1;
+    // Inline function 'kotlin.math.min' call
+    var b = this.declinePerKjPerKg_1 * this.supraCriticalWorkJ_1 / 1000.0 / this.massKgForDose_1;
+    return Math.min(tmp0, b);
+  };
+  protoOf(PowerProviderDurability).optimalPower_36kg19_k$ = function (course, path, pointIndex) {
+    this.massKgForDose_1 = course.get_cyclist_lkc0de_k$().massKg_1;
+    if (pointIndex < this.lastIndex_1) {
+      this.reset_5u6xz3_k$();
+    }
+    if (pointIndex > this.lastIndex_1) {
+      var elapsedS = path.elapsed_7307rw_k$(pointIndex) / 1000.0;
+      var dtS = elapsedS - this.lastElapsedS_1;
+      if (dtS > 0.0 && isFinite(dtS)) {
+        var tmp = this;
+        var tmp_0 = this.supraCriticalWorkJ_1;
+        // Inline function 'kotlin.math.max' call
+        var b = this.lastPowerW_1 - this.criticalPowerW_1;
+        tmp.supraCriticalWorkJ_1 = tmp_0 + Math.max(0.0, b) * dtS;
+      }
+      if (isFinite(elapsedS))
+        this.lastElapsedS_1 = elapsedS;
+      this.lastIndex_1 = pointIndex;
+    }
+    var faded = this.powerW_1 * (1.0 - this.get_decline_qrtukt_k$());
+    this.lastPowerW_1 = faded;
+    return faded;
+  };
+  protoOf(PowerProviderDurability).reset_5u6xz3_k$ = function () {
+    this.lastIndex_1 = -1;
+    this.lastElapsedS_1 = 0.0;
+    this.lastPowerW_1 = 0.0;
+    this.supraCriticalWorkJ_1 = 0.0;
   };
   function PowerProviderFromData() {
   }
@@ -1644,8 +1839,19 @@
     out.setVirtSpeedCurrent_kvg8yv_k$(0, speed);
     var i = 1;
     var iter = 0;
-    $l$loop: while (i < n) {
+    $l$loop_0: while (i < n) {
       var dx = input.distance_5s4s77_k$(i) - input.distance_5s4s77_k$(i - 1 | 0);
+      if (dx <= 0.0) {
+        copyAllFields(this, input, i, out, i);
+        out.setTime_rfmvcd_k$(i, timeMs);
+        out.setElapsed_916ug_k$(i, timeMs - startTimeMs);
+        out.setDx_uk0ky4_k$(i, 0.0);
+        out.setDt_sgarxs_k$(i, 0.0);
+        out.setSpeed_1slan1_k$(i, speed);
+        out.setVirtSpeedCurrent_kvg8yv_k$(i, speed);
+        i = i + 1 | 0;
+        continue $l$loop_0;
+      }
       var pSum = PowerComputer_instance.getNewPower_puvq9i_k$(course, out, i - 1 | 0, true);
       var dt = PowerComputer_instance.getDt_t0pcee_k$(pSum, mEq, speed, dx);
       var speedNew = 2.0 * dx / dt - speed;
@@ -1667,7 +1873,7 @@
       var _unary__edvuaz = iter;
       iter = _unary__edvuaz + 1 | 0;
       if (_unary__edvuaz > 100000)
-        break $l$loop;
+        break $l$loop_0;
     }
     var progression = out.get_indices_xelk8u_k$();
     var inductionVariable = progression.first_1;
@@ -1698,18 +1904,18 @@
   function WheelBearingsPowerProvider_getInstance() {
     return WheelBearingsPowerProvider_instance;
   }
-  function Companion_6() {
-    Companion_instance_6 = this;
+  function Companion_7() {
+    Companion_instance_7 = this;
     this.NONE_1 = new Wind(0.0, 0.0);
   }
-  var Companion_instance_6;
-  function Companion_getInstance_9() {
-    if (Companion_instance_6 == null)
-      new Companion_6();
-    return Companion_instance_6;
+  var Companion_instance_7;
+  function Companion_getInstance_10() {
+    if (Companion_instance_7 == null)
+      new Companion_7();
+    return Companion_instance_7;
   }
   function Wind(speedMS, directionRad) {
-    Companion_getInstance_9();
+    Companion_getInstance_10();
     this.speedMS_1 = speedMS;
     this.directionRad_1 = directionRad;
   }
@@ -1735,7 +1941,7 @@
   function WindProviderNone() {
   }
   protoOf(WindProviderNone).wind_umf2tf_k$ = function (course, path, pointIndex) {
-    return Companion_getInstance_9().NONE_1;
+    return Companion_getInstance_10().NONE_1;
   };
   var WindProviderNone_instance;
   function WindProviderNone_getInstance() {
@@ -1747,6 +1953,50 @@
   protoOf(WindProviderConstant).wind_umf2tf_k$ = function (course, path, pointIndex) {
     return this.wind_1;
   };
+  function WPrimeBalanceComputer() {
+  }
+  protoOf(WPrimeBalanceComputer).compute_ewvoaz_k$ = function (path, options) {
+    if (path.size_1 === 0)
+      return Unit_instance;
+    var cp = options.criticalPowerW_1;
+    var wPrime = options.wPrimeJ_1;
+    var balance = wPrime;
+    path.setWPrimeBalance_stwqkw_k$(0, balance);
+    var inductionVariable = 1;
+    var last = path.size_1;
+    if (inductionVariable < last)
+      do {
+        var i = inductionVariable;
+        inductionVariable = inductionVariable + 1 | 0;
+        var dtSeconds = (path.time_ryc845_k$(i) - path.time_ryc845_k$(i - 1 | 0)) / 1000.0;
+        var power = path.pComputedPower_g6dnn4_k$(i);
+        if (isFinite(dtSeconds) && dtSeconds > 0.0 && isFinite(power)) {
+          balance = this.step_447nu7_k$(balance, power, cp, wPrime, dtSeconds);
+        }
+        path.setWPrimeBalance_stwqkw_k$(i, balance);
+      }
+       while (inductionVariable < last);
+  };
+  protoOf(WPrimeBalanceComputer).step_447nu7_k$ = function (balance, power, cp, wPrime, dtSeconds) {
+    var tmp;
+    if (power >= cp) {
+      tmp = balance - (power - cp) * dtSeconds;
+    } else {
+      var tmp_0 = wPrime - balance;
+      // Inline function 'kotlin.math.exp' call
+      var x = -(cp - power) * dtSeconds / wPrime;
+      tmp = wPrime - tmp_0 * Math.exp(x);
+    }
+    var next = tmp;
+    // Inline function 'kotlin.math.max' call
+    // Inline function 'kotlin.math.min' call
+    var b = Math.max(0.0, next);
+    return Math.min(wPrime, b);
+  };
+  var WPrimeBalanceComputer_instance;
+  function WPrimeBalanceComputer_getInstance() {
+    return WPrimeBalanceComputer_instance;
+  }
   function parseGpxWaypoints(xml) {
     // Inline function 'kotlin.collections.map' call
     var this_0 = GpxParser_instance.parse$default_ez0f31_k$(xml).waypoints_1;
@@ -2013,7 +2263,13 @@
   function toBike(_this__u8e3s4) {
     if (_this__u8e3s4 == null)
       return new Bike();
-    return new Bike(_this__u8e3s4.crr, _this__u8e3s4.inertiaFront, _this__u8e3s4.inertiaRear, _this__u8e3s4.wheelRadiusM, _this__u8e3s4.efficiency);
+    var tmp = _this__u8e3s4.crr;
+    var tmp_0 = _this__u8e3s4.inertiaFront;
+    var tmp_1 = _this__u8e3s4.inertiaRear;
+    var tmp_2 = _this__u8e3s4.wheelRadiusM;
+    var tmp_3 = _this__u8e3s4.efficiency;
+    var tmp0_elvis_lhs = _this__u8e3s4.maxPedalingLeanAngleDeg;
+    return new Bike(tmp, tmp_0, tmp_1, tmp_2, tmp_3, tmp0_elvis_lhs == null ? 20.0 : tmp0_elvis_lhs);
   }
   function toWindProvider(_this__u8e3s4) {
     if (_this__u8e3s4 == null)
@@ -2031,13 +2287,13 @@
         var tmp2_elvis_lhs = _this__u8e3s4.useHarmonics;
         tmp = new PowerProviderConstant(tmp_0, tmp2_elvis_lhs == null ? false : tmp2_elvis_lhs);
         break;
-      case 'constant_tiring':
+      case 'durability':
         var tmp3_elvis_lhs = _this__u8e3s4.power;
         var tmp_1 = tmp3_elvis_lhs == null ? 250.0 : tmp3_elvis_lhs;
-        var tmp4_elvis_lhs = _this__u8e3s4.useHarmonics;
-        var tmp_2 = tmp4_elvis_lhs == null ? false : tmp4_elvis_lhs;
-        var tmp5_elvis_lhs = _this__u8e3s4.tiringDuration;
-        tmp = new PowerProviderConstantWithTiring(tmp_1, tmp_2, tmp5_elvis_lhs == null ? 7200.0 : tmp5_elvis_lhs);
+        var tmp4_elvis_lhs = _this__u8e3s4.criticalPower;
+        var tmp_2 = tmp4_elvis_lhs == null ? 250.0 : tmp4_elvis_lhs;
+        var tmp5_elvis_lhs = _this__u8e3s4.useHarmonics;
+        tmp = new PowerProviderDurability(tmp_1, tmp_2, VOID, VOID, tmp5_elvis_lhs == null ? false : tmp5_elvis_lhs);
         break;
       case 'from_data':
         tmp = PowerProviderFromData_instance;
@@ -2223,6 +2479,7 @@
   MaxSpeedComputer_instance = new MaxSpeedComputer();
   MuscularPowerProvider_instance = new MuscularPowerProvider();
   PowerComputer_instance = new PowerComputer();
+  Companion_instance_6 = new Companion_6();
   PowerProviderFromData_instance = new PowerProviderFromData();
   RhoProviderDefault_instance = new RhoProviderDefault();
   RhoProviderEstimate_instance = new RhoProviderEstimate();
@@ -2230,6 +2487,7 @@
   VirtualizeService_instance = new VirtualizeService();
   WheelBearingsPowerProvider_instance = new WheelBearingsPowerProvider();
   WindProviderNone_instance = new WindProviderNone();
+  WPrimeBalanceComputer_instance = new WPrimeBalanceComputer();
   //endregion
   //region block: exports
   function $jsExportAll$(_) {