Demand/Dynamic User Assignment

The tool  /tools/assign/duaIterate.py  can be used to compute the (approximate) dynamic user equilibrium.

python duaIterate.py -n  -t  -l 

= General behavior = This script tries to calculate a user equilibrium, that is, it tries to find a route for each vehicle (each trip from the trip-file above) such that each vehicle cannot reduce its travel cost (usually the travel time) by using a different route. It does so iteratively (hence the name) by first routing the vehicles in a network with default edge costs (usually the travel time if driving with maximum speed) then simulating the calculated routes resulting in "real" edge costs which are fed back into the router for calculating new routes and so forth. The number of iterations is fixed by a parameter see above. In order to ensure convergence there are different methods employed to calculate the route choice probability from the route cost (so the vehicle does not always choose the "cheapest" route). In general, new routes will be added by the router to the route set of each vehicle in each iteration (at least if none of the present routes is the "cheapest") and may be chosen according to the route choice mechanisms described below.

The two methods which are implemented are called Gawron and Logit (reference needed!!!) in the following. The input for each of the methods is a weight or cost function $$w$$ on the edges of the net, coming from the simulation or default costs (in the first step or for edges which have not been traveled yet), and a set of routes $$R$$ where each route $$r$$ has an old cost $$c_r$$ and an old probability $$p_r$$ (from the last iteration) and needs a new cost $$c_r'$$ and a new probability $$p_r'$$.

= Logit = The Logit mechanism applies a fixed formula to each route to calculate the new probability. It ignores old costs and old probabilities and takes the route cost directly as the sum of the edge costs from the last simulation.

$$c_r' = \sum_{e\in r}w(e)$$

The probabilities are calculated from an exponential function with parameter $$\theta$$ scaled by the sum over all route values:

$$p_r' = \frac{\exp(\theta c_r')}{\sum_{s\in R}\exp(\theta c_s')}$$

= Gawron = = Parameters =