compute_speeds_2(densities, flows, frflows, celldata, outflow, ts)

function speeds = compute_speeds_2(densities, flows, frflows, celldata, outflow, ts)
% COMPUTE_SPEEDS_2 - compute speeds from densities and flows.
%
% Call:   speeds = compute_speeds_2(densities, flows, frflows, celldata, outflow, ts)
%
% Parameters:
%             densities - vector of densities;
%             flows     - vector of flows, must have the same size as 'densities';
%             frflows   - vector of off-ramp flows, must have the same size as 'densities';
%             celldata  - array of cell data structures,
%                         whose length must be the same as size of 'densities';
%             outflow   - how much can be allowed through the right boundary;
%             ts        - sampling period.
%
% Returns:   speeds - vector of speeds.
%
% Last modified:   10/06/2006.

%
% Alex Kurzhanskiy   <[e-mail]>
%

K = size(densities, 1);
L = size(flows, 1);
M = size(frflows, 1);
N = size(celldata, 2);  % row vector

if K ~= L
  error('COMPUTE_SPEEDS_2: number of densities does not match number of flows.');
end

if K ~= M
  error('COMPUTE_SPEEDS_2: number of densities does not match number of off-ramp flows.');
end

if K ~= N
  error('COMPUTE_SPEEDS_2: number of densities does not match number of cells.');
end

speeds = [];

for i = 1:N
  if densities(i, 1) <= 1e-7
    v = celldata(i).FDfmax / celldata(i).FDrhocrit;
  else
    if i == N
      v = get_rbflow_2(densities(N, 1), flows(N, 1), frflows(N, 1), celldata(N), outflow, ts) / densities(i, 1);
    else
      v = (flows(i+1, 1) + frflows(i, 1)) / densities(i, 1);
    end
  end

  v = min([v (celldata(i).FDfmax / celldata(i).FDrhocrit)]);

  speeds = [speeds; v];
end

return;