Fabio/brouter
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions 1

Merge pull request #127 from Phyks/travelTime

Browse source 
Few fixes to travel time computation
This commit is contained in:
abrensch 2018-12-12 22:22:33 +01:00 committed by GitHub
commit cab8ca3dc8
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
2 changed files with 11 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
brouter-core/src/main/java/btools/router/RoutingContext.java
View file

@ -175,8 +175,8 @@ public final class RoutingContext
turnInstructionRoundabouts = expctxGlobal.getVariableValue( "turnInstructionRoundabouts", 1.f ) != 0.f;
// Speed computation model (for bikes)
// Mass of the biker + bike + luggages, in kg
bikeMass = expctxGlobal.getVariableValue( "bikeMass", 90.f );
// Total mass (biker + bike + luggages or hiker), in kg
totalMass = expctxGlobal.getVariableValue( "totalMass", 90.f );
// Max speed (before braking), in km/h in profile and m/s in code
maxSpeed = expctxGlobal.getVariableValue( "maxSpeed", 45.f ) / 3.6;
// Equivalent surface for wind, S * C_x, F = -1/2 * S * C_x * v^2 = - S_C_x * v^2
@ -224,7 +224,7 @@ public final class RoutingContext
public boolean turnInstructionRoundabouts;
// Speed computation model (for bikes)
public double bikeMass;
public double totalMass;
public double maxSpeed;
public double S_C_x;
public double defaultC_r;

11
brouter-core/src/main/java/btools/router/StdPath.java
View file

@ -243,16 +243,21 @@ final class StdPath extends OsmPath
double incline = calcIncline( dist );
double speed; // Travel speed
double f_roll = rc.bikeMass * GRAVITY * ( rc.defaultC_r + incline );
if (rc.footMode )
double energy;
if (rc.footMode || rc.expctxWay.getCostfactor() > 5)
{
// Use Tobler's hiking function for walking sections
speed = 6 * exp(-3.5 * Math.abs( incline + 0.05)) / 3.6;
energy = rc.totalMass * GRAVITY * incline * dist;
}
else if (rc.bikeMode)
{
double f_roll = rc.totalMass * GRAVITY * ( rc.defaultC_r + incline );
speed = solveCubic( rc.S_C_x, f_roll, rc.bikerPower );
speed = Math.min(speed, rc.maxSpeed);
// Don't compute energy assuming constant biker power, as speed is capped
// to a max value.
energy = dist*(rc.S_C_x*speed*speed + f_roll);
}
else
{
@ -260,7 +265,7 @@ final class StdPath extends OsmPath
}
float dt = (float) ( dist / speed );
totalTime += dt;
totalEnergy += dist*(rc.S_C_x*speed*speed + f_roll);
totalEnergy += energy;
}
private static double solveCubic( double a, double c, double d )