public static void SetTimeImpedances(ITripWrapper trip)
{
foreach (var time in Times[ParallelUtility.GetBatchFromThreadId()])
{
SetTimeImpedanceAndWindow(trip, time);
}
}
private static void RunTripModelSuite(ITourWrapper tour, IHalfTour halfTour, ITripWrapper trip, int batchNumber)
{
var nextTrip = GenerateIntermediateStop(halfTour, trip, batchNumber);
SetIntermediateStopDestination(trip, nextTrip, batchNumber);
SetTripModeAndTime(tour, trip, batchNumber);
if (!trip.PersonDay.IsValid)
{
return;
}
// retrieves window based on whether or not the trip's tour is home-based or work-based
var timeWindow = tour.IsHomeBasedTour ? tour.PersonDay.TimeWindow : tour.ParentTour.TimeWindow;
if (trip.IsHalfTourFromOrigin && trip.Sequence == 1)
{
// occupies minutes in window between destination and stop
timeWindow.SetBusyMinutes(tour.DestinationArrivalTime, trip.ArrivalTime);
}
else if (!trip.IsHalfTourFromOrigin && trip.Sequence == 1)
{
// occupies minutes in window between destination and stop
timeWindow.SetBusyMinutes(tour.DestinationDepartureTime, trip.ArrivalTime);
}
else
{
// occupies minutes in window from previous stop to stop
timeWindow.SetBusyMinutes(trip.GetPreviousTrip().DepartureTime, trip.ArrivalTime);
}
}
private static TWrapper CreateWrapper(ITourWrapper tourWrapper, ITripWrapper trip, int nextTripId, int intermediateStopPurpose, int destinationPurpose, IHalfTour halfTour)
{
var t = new TModel {
Id = nextTripId,
TourId = tourWrapper.Id,
HouseholdId = tourWrapper.HouseholdId,
PersonSequence = tourWrapper.PersonSequence,
Day = tourWrapper.Day,
Direction = trip.Direction,
Sequence = trip.Sequence + 1,
DestinationAddressType = trip.DestinationAddressType,
DestinationParcelId = trip.DestinationParcelId,
DestinationZoneKey = trip.DestinationZoneKey,
OriginAddressType = trip.DestinationAddressType,
OriginPurpose = intermediateStopPurpose,
DestinationPurpose = destinationPurpose,
DepartureTime = 180,
ArrivalTime = 180,
PathType = 1,
ExpansionFactor = tourWrapper.Household.ExpansionFactor
};
var type = typeof(TWrapper);
var instance = Activator.CreateInstance(type, t, tourWrapper, halfTour);
var wrapper = (TWrapper)instance;
wrapper.IsToTourOrigin = trip.IsToTourOrigin;
return(wrapper);
}
public static void SetTimeImpedances(ITripWrapper trip)
{
foreach (HTripTime time in Times[ParallelUtility.threadLocalAssignedIndex.Value])
{
SetTimeImpedanceAndWindow(trip, time);
}
}
public ITripWrapper CreateNextTrip(ITripWrapper trip, int intermediateStopPurpose, int destinationPurpose)
{
_t.PersonDay.IncrementSimulatedStops(intermediateStopPurpose);
return
(_tripWrapperCreator
.CreateWrapper(_t, trip, trip.Id + 1, intermediateStopPurpose, destinationPurpose, this));
}
public static void WriteTripForTDM(ITripWrapper trip, TDMTripListExporter tdmTripListExporter)
{
if (tdmTripListExporter == null)
{
return;
}
tdmTripListExporter.Export(trip);
}
public static void SetTimeImpedanceAndWindow(ITripWrapper trip, HTripTime time)
{
var tour = trip.Tour;
var alternativeIndex = time.Index;
var period = time.DeparturePeriod;
// set mode LOS and mode availability
if (period.End < trip.EarliestDepartureTime || period.Start > trip.LatestDepartureTime)
{
time.Available = false;
}
else
{
var pathMode = (trip.Mode >= Global.Settings.Modes.SchoolBus - 1) ? Global.Settings.Modes.Hov3 : trip.Mode;
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.Run(
trip.Household.RandomUtility,
trip.IsHalfTourFromOrigin ? trip.DestinationParcel : trip.OriginParcel,
trip.IsHalfTourFromOrigin ? trip.OriginParcel : trip.DestinationParcel,
period.Middle,
0,
tour.DestinationPurpose,
tour.CostCoefficient,
tour.TimeCoefficient,
tour.Person.IsDrivingAge,
tour.Household.VehiclesAvailable,
tour.Household.OwnsAutomatedVehicles > 0,
tour.Person.GetTransitFareDiscountFraction(),
true,
pathMode);
var pathTypeModel = pathTypeModels.First(x => x.Mode == pathMode);
time.Available = pathTypeModel.Available;
time.ModeLOS = pathTypeModel;
//set the feasible window within the small period, accounting for travel time, and recheck availability
if (time.Available)
{
time.EarliestFeasibleDepatureTime = Math.Max(period.Start,
trip.IsHalfTourFromOrigin
//JLB 20130723 replace next line
//? trip.ArrivalTimeLimit + - (int) (time.ModeLOS.PathTime + 0.5)
? trip.ArrivalTimeLimit + (int)(time.ModeLOS.PathTime + 0.5)
: trip.EarliestDepartureTime);
time.LatestFeasibleDepartureTime = Math.Min(period.End,
trip.IsHalfTourFromOrigin
? trip.LatestDepartureTime
: trip.ArrivalTimeLimit - (int)(time.ModeLOS.PathTime + 0.5));
time.Available = time.EarliestFeasibleDepatureTime < time.LatestFeasibleDepartureTime;
}
}
}
private static void SetIntermediateStopDestination(ITripWrapper trip, ITripWrapper nextTrip, int batchNumber)
{
if (nextTrip == null || trip.IsToTourOrigin || !Global.Configuration.ShouldRunIntermediateStopLocationModel)
{
if (trip.IsToTourOrigin)
{
ChoiceModelFactory.TotalTimesTripIsToTourOrigin[batchNumber]++;
}
else if (nextTrip == null)
{
ChoiceModelFactory.TotalTimesNextTripIsNull[batchNumber]++;
}
if (trip.DestinationPurpose == Global.Settings.Purposes.NoneOrHome && Global.Configuration.IsInEstimationMode && Global.Configuration.EstimationModel == "IntermediateStopLocationModel")
{
Global.PrintFile.WriteEstimationRecordExclusionMessage("ChoiceModelRunner", "SetIntermediateStopDestination", trip.Household.Id, trip.Person.Sequence, trip.Day, trip.Tour.Sequence, trip.Direction, trip.Sequence, 1);
}
return;
}
// sets the new destination for the trip
if (trip.DestinationPurpose == Global.Settings.Purposes.ChangeMode)
{
// CHANGE_MODE location is always park and ride node for tour
var parkAndRideNode = ChoiceModelFactory.ParkAndRideNodeDao.Get(trip.Tour.ParkAndRideNodeId);
if (parkAndRideNode != null)
{
trip.DestinationParcelId = parkAndRideNode.NearestParcelId;
trip.DestinationParcel = ChoiceModelFactory.Parcels[trip.DestinationParcelId];
trip.DestinationZoneKey = ChoiceModelFactory.ZoneKeys[parkAndRideNode.ZoneId];
//trip.DestinationZoneKey = parkAndRideNode.Id;
trip.DestinationAddressType = Global.Settings.AddressTypes.Other;
ChoiceModelFactory.TotalTimesChangeModeLocationSet[batchNumber]++;
}
}
else
{
ChoiceModelFactory.TotalTimesIntermediateStopLocationModelRun[batchNumber]++;
Global.ChoiceModelSession.Get <IntermediateStopLocationModel>().Run(trip, Global.Configuration.IntermediateStopLocationModelSampleSize);
}
if (Global.Configuration.IsInEstimationMode)
{
return;
}
nextTrip.OriginParcelId = trip.DestinationParcelId;
nextTrip.OriginParcel = trip.DestinationParcel;
nextTrip.OriginZoneKey = trip.DestinationZoneKey;
nextTrip.SetOriginAddressType(trip.DestinationAddressType);
}
public void Run(IHouseholdDayWrapper householdDay, ITripWrapper trip)
{
if (trip == null)
{
throw new ArgumentNullException("trip");
}
trip.PersonDay.ResetRandom(40 * (2 * trip.Tour.Sequence - 1 + trip.Direction - 1) + 50 + trip.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
}
var choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(trip.Id);
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
if (trip.DestinationParcel == null || trip.OriginParcel == null || trip.Mode <= Global.Settings.Modes.None || trip.Mode == Global.Settings.Modes.Other)
{
return;
}
RunModel(choiceProbabilityCalculator, householdDay, trip, new HTripTime(trip.DepartureTime));
choiceProbabilityCalculator.WriteObservation();
}
else
{
RunModel(choiceProbabilityCalculator, householdDay, trip);
var chosenAlternative = choiceProbabilityCalculator.SimulateChoice(trip.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", trip.PersonDay.Id);
if (!Global.Configuration.IsInEstimationMode)
{
trip.PersonDay.IsValid = false;
}
return;
}
var choice = (HTripTime)chosenAlternative.Choice;
trip.DepartureTime = choice.GetRandomFeasibleMinute(trip, choice);
}
}
public int GetDepartureTime(ITripWrapper trip)
{
if (trip == null)
{
throw new ArgumentNullException("trip");
}
var timeWindow = trip.Tour.ParentTour == null ? trip.Tour.PersonDay.TimeWindow : trip.Tour.ParentTour.TimeWindow;
var departureTime = timeWindow.GetAvailableMinute(trip.Household.RandomUtility, DeparturePeriod.Start, DeparturePeriod.End);
//if (departureTime == Constants.DEFAULT_VALUE) {
// throw new InvalidDepartureTimeException();
//}
return(departureTime);
}
public int GetRandomFeasibleMinute(ITripWrapper trip, HTripTime time)
{
if (trip == null || time == null)
{
throw new ArgumentNullException("trip time");
}
ITimeWindow timeWindow = trip.Tour.ParentTour == null ? trip.Tour.PersonDay.TimeWindow : trip.Tour.ParentTour.TimeWindow;
int departureTime = timeWindow.GetAvailableMinute(trip.Household.RandomUtility, time.EarliestFeasibleDepatureTime, time.LatestFeasibleDepartureTime);
//if (departureTime == Constants.DEFAULT_VALUE) {
// throw new InvalidDepartureTimeException();
//}
return(departureTime);
}
public int Run(ITripWrapper trip, int choice)
{
if (trip == null)
{
throw new ArgumentNullException("trip");
}
trip.PersonDay.ResetRandom(40 * (2 * trip.Tour.Sequence - 1 + trip.Direction - 1) + 20 + trip.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return(choice);
}
}
ChoiceProbabilityCalculator choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(trip.Id);
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
if (trip.OriginParcel == null)
{
return(Constants.DEFAULT_VALUE);
}
RunModel(choiceProbabilityCalculator, trip, choice);
choiceProbabilityCalculator.WriteObservation();
}
else
{
RunModel(choiceProbabilityCalculator, trip);
ChoiceProbabilityCalculator.Alternative chosenAlternative = choiceProbabilityCalculator.SimulateChoice(trip.Household.RandomUtility);
choice = (int)chosenAlternative.Choice;
}
return(choice);
}
private static ITripWrapper GenerateIntermediateStop(IHalfTour halfTour, ITripWrapper trip, int batchNumber)
{
if (!Global.Configuration.ShouldRunIntermediateStopGenerationModel)
{
return(null);
}
ITripWrapper nextTrip = null;
if (Global.Configuration.IsInEstimationMode)
{
// -- in estimation mode --
// sets the trip's destination purpose, determines whether or not a stop is generated in application mode
// uses trip instead of nextTrip, deals with subtours with tour origin at work
// need to set trip.IsToTourOrigin first
trip.IsToTourOrigin = trip.Sequence == trip.HalfTour.Trips.Count(); // last trip in half tour
var intermediateStopPurpose = trip.IsToTourOrigin ? Global.Settings.Purposes.NoneOrHome : trip.DestinationPurpose;
nextTrip = trip.GetNextTrip();
if (intermediateStopPurpose != Global.Settings.Purposes.NoneOrHome)
{
ChoiceModelFactory.TotalTimesIntermediateStopGenerated[batchNumber]++;
}
if (trip.PersonDay.GetTotalStops() > 0)
{
ChoiceModelFactory.TotalTimesIntermediateStopGenerationModelRun[batchNumber]++;
Global.ChoiceModelSession.Get <IntermediateStopGenerationModel>().Run(trip, intermediateStopPurpose);
}
}
else
{
// -- in application mode --
// sets the trip's destination purpose, determines whether or not a stop is generated
// first, if it is the first trip on a park and ride half tour, then make it a change mode stop
// TODO: this doesn't allow stops between the destination and the transit stop - can improve later
int intermediateStopPurpose;
if (trip.Sequence == 1 && trip.Tour.Mode == Global.Settings.Modes.ParkAndRide)
{
intermediateStopPurpose = Global.Settings.Purposes.ChangeMode;
ChoiceModelFactory.TotalTimesChangeModeStopGenerated[batchNumber]++;
}
else if (trip.PersonDay.GetTotalStops() == 0)
{
intermediateStopPurpose = Global.Settings.Purposes.NoneOrHome;
}
else
{
ChoiceModelFactory.TotalTimesIntermediateStopGenerationModelRun[batchNumber]++;
intermediateStopPurpose = Global.ChoiceModelSession.Get <IntermediateStopGenerationModel>().Run(trip);
}
if (intermediateStopPurpose != Global.Settings.Purposes.NoneOrHome)
{
ChoiceModelFactory.TotalTimesIntermediateStopGenerated[batchNumber]++;
var destinationPurposeForNextTrip = trip.DestinationPurpose;
// creates the next trip in the half-tour
// the next trip's destination is set to the current trip's destination
nextTrip = halfTour.CreateNextTrip(trip, intermediateStopPurpose, destinationPurposeForNextTrip);
halfTour.Trips.Add(nextTrip);
trip.DestinationAddressType = Global.Settings.AddressTypes.None;
trip.DestinationPurpose = intermediateStopPurpose;
trip.IsToTourOrigin = false;
}
else
{
trip.IsToTourOrigin = true;
}
}
return(nextTrip);
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITripWrapper trip, TripTime choice = null)
{
IPersonWrapper person = trip.Person;
IPersonDayWrapper personDay = trip.PersonDay;
ITourWrapper tour = trip.Tour;
// person inputs
int partTimeWorkerFlag = person.IsPartTimeWorker.ToFlag();
int nonworkingAdultFlag = person.IsNonworkingAdult.ToFlag();
int universityStudentFlag = person.IsUniversityStudent.ToFlag();
int retiredAdultFlag = person.IsRetiredAdult.ToFlag();
int drivingAgeStudentFlag = person.IsDrivingAgeStudent.ToFlag();
int childAge5Through15Flag = person.IsChildAge5Through15.ToFlag();
int childUnder5Flag = person.IsChildUnder5.ToFlag();
// set tour inputs
int workTourFlag = tour.IsWorkPurpose().ToFlag();
int notWorkTourFlag = (!tour.IsWorkPurpose()).ToFlag();
int notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
// set trip inputs
bool originChangeMode = trip.Sequence > 1 && trip.GetPreviousTrip().DestinationPurpose == Global.Settings.Purposes.ChangeMode;
int originWorkFlag = trip.IsWorkOriginPurpose().ToFlag();
int originSchoolFlag = trip.IsSchoolOriginPurpose().ToFlag();
int originEscortFlag = trip.IsEscortOriginPurpose().ToFlag();
int originShoppingFlag = trip.IsShoppingOriginPurpose().ToFlag();
int originPersonalBusinessFlag = trip.IsPersonalBusinessOriginPurpose().ToFlag();
int originMealFlag = trip.IsMealOriginPurpose().ToFlag();
int originSocialFlag = trip.IsSocialOriginPurpose().ToFlag();
int sovOrHovTripFlag = trip.UsesSovOrHovModes().ToFlag();
int transitTripFlag = trip.IsTransitMode().ToFlag();
int halfTourFromOriginFlag = trip.IsHalfTourFromOrigin.ToFlag();
int halfTourFromDestinationFlag = (!trip.IsHalfTourFromOrigin).ToFlag();
// set remaining inputs
ITimeWindow timeWindow = tour.IsHomeBasedTour ? personDay.TimeWindow : tour.ParentTour.TimeWindow;
ITripModeImpedance[] impedances = trip.GetTripModeImpedances();
int remainingToursCount = personDay.HomeBasedTours - personDay.GetTotalSimulatedTours();
int tripRemainingInHalfTour = (trip.DestinationParcel != null && trip.DestinationParcel != tour.OriginParcel).ToFlag(); // we don't know exact #
for (int arrivalPeriodIndex = 1; arrivalPeriodIndex < DayPeriod.SmallDayPeriods.Length; arrivalPeriodIndex++)
{
MinuteSpan arrivalPeriod = DayPeriod.SmallDayPeriods[arrivalPeriodIndex];
int previousArrivalTime = trip.GetPreviousTrip().ArrivalTime;
if (previousArrivalTime < arrivalPeriod.Start || previousArrivalTime > arrivalPeriod.End)
{
continue;
}
ITripModeImpedance arrivalImpedance = impedances[arrivalPeriod.Index]; // moved to here so not reset for every alternative
foreach (TripTime time in TripTime.Times)
{
MinuteSpan departurePeriod = time.DeparturePeriod; // moved to here so can use travel time
ITripModeImpedance departureImpedance = impedances[departurePeriod.Index];
// change availability check to include travel duration
int travelDuration = (int)Math.Round(departureImpedance.TravelTime + 0.5);
// if not the trip home, on a home-based tour, also include fastest time from the destinatinon to home
// if (trip.Tour.IsHomeBasedTour && trip.DestinationPurpose != Global.Settings.Purposes.NoneOrHome) {
// var fastestMode = Math.Min(trip.Tour.Mode, Global.Settings.Modes.Hov3);
// var pathTypeModel = PathTypeModelFactory.Singleton.Run(trip.DestinationParcel, trip.Household.ResidenceParcel, departurePeriod.Middle, 0,
// trip.Tour.DestinationPurpose, trip.Tour.CostCoefficient, trip.Tour.TimeCoefficient,
// trip.Person.IsDrivingAge, trip.Household.VehiclesAvailable, trip.Tour.Person.TransitFareDiscountFraction, false, fastestMode).First();
// travelDuration += (int) Math.Round(pathTypeModel.PathTime + 0.5);
// }
int bestArrivalTime
= trip.IsHalfTourFromOrigin
? Math.Max(departurePeriod.End - travelDuration, 1)
: Math.Min(departurePeriod.Start + travelDuration, Global.Settings.Times.MinutesInADay);
bool available =
originChangeMode
? arrivalPeriod.Index == time.DeparturePeriod.Index
: (trip.IsHalfTourFromOrigin && // if change mode, must be in same period
arrivalPeriod.Index > departurePeriod.Index &&
timeWindow.EntireSpanIsAvailable(bestArrivalTime, arrivalPeriod.Start - 1)) ||
(!trip.IsHalfTourFromOrigin &&
arrivalPeriod.Index < departurePeriod.Index &&
timeWindow.EntireSpanIsAvailable(arrivalPeriod.End, bestArrivalTime - 1)) ||
arrivalPeriod.Index == time.DeparturePeriod.Index &&
timeWindow.TotalAvailableMinutes(arrivalPeriod.Start, arrivalPeriod.End) > travelDuration;
double departurePeriodFraction = timeWindow.TotalAvailableMinutes(departurePeriod.Start, departurePeriod.End) / (departurePeriod.End - departurePeriod.Start + 1D);
int duration = Math.Abs(departurePeriod.Middle - arrivalPeriod.Middle);
available = available && departurePeriodFraction > 0;
//ensure transit path type is available in alternative
if (trip.Mode == Global.Settings.Modes.Transit && !Global.StopAreaIsEnabled)
{
double transitPathTypeInVehicleTime = ImpedanceRoster.GetValue("ivtime", trip.Mode, trip.PathType, trip.ValueOfTime, time.DeparturePeriod.Middle, trip.OriginParcel.ZoneId, trip.DestinationParcel.ZoneId).Variable;
available = available && (transitPathTypeInVehicleTime > 0);
}
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(time.Index, available, choice != null && choice.Equals(time));
// if (choice.Equals(tripTime) && !available) {
// Console.WriteLine(available);
// }
if (!alternative.Available)
{
continue;
}
alternative.Choice = time;
double departurePeriodShift = time.DeparturePeriod.Index * (48.0 / DayPeriod.SMALL_DAY_PERIOD_TOTAL_TRIP_TIMES); //adjust shift amount if period lengths change
if (trip.IsHalfTourFromOrigin)
{
// outbound "departure" (arrival) period constants
alternative.AddUtilityTerm(11, time.DeparturePeriod.Middle.IsBetween(Global.Settings.Times.ThreeAM, Global.Settings.Times.SixAM).ToFlag());
alternative.AddUtilityTerm(12, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SixAM, Global.Settings.Times.SevenAM).ToFlag());
alternative.AddUtilityTerm(13, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenAM, Global.Settings.Times.EightAM).ToFlag());
alternative.AddUtilityTerm(14, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.EightAM, Global.Settings.Times.NineAM).ToFlag());
alternative.AddUtilityTerm(15, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.NineAM, Global.Settings.Times.TenAM).ToFlag());
alternative.AddUtilityTerm(16, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.OnePM).ToFlag());
alternative.AddUtilityTerm(17, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.OnePM, Global.Settings.Times.FourPM).ToFlag());
alternative.AddUtilityTerm(18, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FourPM, Global.Settings.Times.SevenPM).ToFlag());
alternative.AddUtilityTerm(19, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenPM, Global.Settings.Times.TenPM).ToFlag());
alternative.AddUtilityTerm(20, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenPM, Global.Settings.Times.MinutesInADay).ToFlag());
}
else
{
// return departure period constants
alternative.AddUtilityTerm(21, time.DeparturePeriod.Middle.IsBetween(Global.Settings.Times.ThreeAM, Global.Settings.Times.SevenAM).ToFlag());
alternative.AddUtilityTerm(22, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenAM, Global.Settings.Times.TenAM).ToFlag());
alternative.AddUtilityTerm(23, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.OnePM).ToFlag());
alternative.AddUtilityTerm(24, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.OnePM, Global.Settings.Times.ThreePM).ToFlag());
alternative.AddUtilityTerm(124, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.ThreePM, Global.Settings.Times.FourPM).ToFlag());
alternative.AddUtilityTerm(25, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FourPM, Global.Settings.Times.FivePM).ToFlag());
alternative.AddUtilityTerm(26, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FivePM, Global.Settings.Times.SixPM).ToFlag());
alternative.AddUtilityTerm(27, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SixPM, Global.Settings.Times.SevenPM).ToFlag());
alternative.AddUtilityTerm(28, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenPM, Global.Settings.Times.NinePM).ToFlag());
alternative.AddUtilityTerm(29, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.NinePM, Global.Settings.Times.Midnight).ToFlag());
alternative.AddUtilityTerm(30, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.Midnight, Global.Settings.Times.MinutesInADay).ToFlag());
}
alternative.AddUtilityTerm(31, duration.IsRightExclusiveBetween(Global.Settings.Times.ZeroHours, Global.Settings.Times.OneHour).ToFlag()); // 0 - 1
alternative.AddUtilityTerm(32, duration.IsRightExclusiveBetween(Global.Settings.Times.OneHour, Global.Settings.Times.TwoHours).ToFlag()); // 1 - 2
alternative.AddUtilityTerm(33, duration.IsRightExclusiveBetween(Global.Settings.Times.TwoHours, Global.Settings.Times.ThreeHours).ToFlag()); // 2 - 3
alternative.AddUtilityTerm(34, duration.IsRightExclusiveBetween(Global.Settings.Times.ThreeHours, Global.Settings.Times.FiveHours).ToFlag()); // 3 - 5
alternative.AddUtilityTerm(35, duration.IsRightExclusiveBetween(Global.Settings.Times.FiveHours, Global.Settings.Times.SevenHours).ToFlag()); // 5 - 7
alternative.AddUtilityTerm(36, duration.IsRightExclusiveBetween(Global.Settings.Times.SevenHours, Global.Settings.Times.NineHours).ToFlag()); // 7 - 9
alternative.AddUtilityTerm(37, duration.IsRightExclusiveBetween(Global.Settings.Times.NineHours, Global.Settings.Times.TwelveHours).ToFlag()); // 9 - 12
alternative.AddUtilityTerm(38, duration.IsRightExclusiveBetween(Global.Settings.Times.TwelveHours, Global.Settings.Times.FourteenHours).ToFlag()); // 12 - 14
alternative.AddUtilityTerm(39, duration.IsRightExclusiveBetween(Global.Settings.Times.FourteenHours, Global.Settings.Times.EighteenHours).ToFlag()); // 14 - 18
alternative.AddUtilityTerm(40, (duration >= Global.Settings.Times.EighteenHours).ToFlag()); // 18 - 24
//these were duplicate departure shift variables before, and constrained to 0 in all the coefficient files - replaced by duration shifts
alternative.AddUtilityTerm(41, partTimeWorkerFlag * duration);
alternative.AddUtilityTerm(43, nonworkingAdultFlag * duration);
alternative.AddUtilityTerm(45, universityStudentFlag * duration);
alternative.AddUtilityTerm(47, retiredAdultFlag * duration);
alternative.AddUtilityTerm(49, drivingAgeStudentFlag * duration);
alternative.AddUtilityTerm(51, childAge5Through15Flag * duration);
alternative.AddUtilityTerm(53, childUnder5Flag * duration);
alternative.AddUtilityTerm(55, halfTourFromOriginFlag * duration);
alternative.AddUtilityTerm(131, workTourFlag * halfTourFromOriginFlag * duration);
alternative.AddUtilityTerm(133, workTourFlag * halfTourFromDestinationFlag * duration);
alternative.AddUtilityTerm(135, notWorkTourFlag * halfTourFromDestinationFlag * duration);
alternative.AddUtilityTerm(137, notHomeBasedTourFlag * duration);
alternative.AddUtilityTerm(145, originEscortFlag * duration);
alternative.AddUtilityTerm(147, originShoppingFlag * duration);
alternative.AddUtilityTerm(149, originMealFlag * duration);
alternative.AddUtilityTerm(151, originSocialFlag * duration);
alternative.AddUtilityTerm(153, originPersonalBusinessFlag * duration);
alternative.AddUtilityTerm(155, originSchoolFlag * duration);
alternative.AddUtilityTerm(42, partTimeWorkerFlag * departurePeriodShift);
alternative.AddUtilityTerm(44, nonworkingAdultFlag * departurePeriodShift);
alternative.AddUtilityTerm(46, universityStudentFlag * departurePeriodShift);
alternative.AddUtilityTerm(48, retiredAdultFlag * departurePeriodShift);
alternative.AddUtilityTerm(50, drivingAgeStudentFlag * departurePeriodShift);
alternative.AddUtilityTerm(52, childAge5Through15Flag * departurePeriodShift);
alternative.AddUtilityTerm(54, childUnder5Flag * departurePeriodShift);
alternative.AddUtilityTerm(56, halfTourFromOriginFlag * departurePeriodShift);
alternative.AddUtilityTerm(132, workTourFlag * halfTourFromOriginFlag * departurePeriodShift);
alternative.AddUtilityTerm(134, workTourFlag * halfTourFromDestinationFlag * departurePeriodShift);
alternative.AddUtilityTerm(136, notWorkTourFlag * halfTourFromDestinationFlag * departurePeriodShift);
alternative.AddUtilityTerm(138, notHomeBasedTourFlag * departurePeriodShift);
alternative.AddUtilityTerm(146, originEscortFlag * departurePeriodShift);
alternative.AddUtilityTerm(148, originShoppingFlag * departurePeriodShift);
alternative.AddUtilityTerm(150, originMealFlag * departurePeriodShift);
alternative.AddUtilityTerm(152, originSocialFlag * departurePeriodShift);
alternative.AddUtilityTerm(154, originPersonalBusinessFlag * departurePeriodShift);
alternative.AddUtilityTerm(156, originSchoolFlag * departurePeriodShift);
alternative.AddUtilityTerm(86, sovOrHovTripFlag * Math.Max(departureImpedance.GeneralizedTime, 0) * tour.TimeCoefficient);
alternative.AddUtilityTerm(88, transitTripFlag * Math.Max(departureImpedance.GeneralizedTime, 0) * tour.TimeCoefficient);
alternative.AddUtilityTerm(89, transitTripFlag * (departureImpedance.GeneralizedTime < 0).ToFlag());
alternative.AddUtilityTerm(92, halfTourFromOriginFlag * Math.Log(departurePeriodFraction));
alternative.AddUtilityTerm(92, halfTourFromDestinationFlag * Math.Log(departurePeriodFraction));
alternative.AddUtilityTerm(99, tripRemainingInHalfTour / (1D + halfTourFromOriginFlag * departureImpedance.AdjacentMinutesBefore + halfTourFromDestinationFlag * departureImpedance.AdjacentMinutesAfter));
alternative.AddUtilityTerm(97, remainingToursCount / (1D + halfTourFromOriginFlag * (arrivalImpedance.TotalMinutesAfter + departureImpedance.TotalMinutesBefore) + halfTourFromDestinationFlag * (arrivalImpedance.TotalMinutesBefore + departureImpedance.TotalMinutesAfter)));
alternative.AddUtilityTerm(98, remainingToursCount / (1D + halfTourFromOriginFlag * Math.Max(arrivalImpedance.MaxMinutesBefore, departureImpedance.MaxMinutesBefore) + halfTourFromDestinationFlag * Math.Max(arrivalImpedance.MaxMinutesBefore, departureImpedance.MaxMinutesAfter)));
}
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITripWrapper trip, int choice = Constants.DEFAULT_VALUE)
{
IHouseholdWrapper household = trip.Household;
IPersonWrapper person = trip.Person;
IPersonDayWrapper personDay = trip.PersonDay;
ITourWrapper tour = trip.Tour;
Framework.DomainModels.Models.IHalfTour halfTour = trip.HalfTour;
// household inputs
int childrenFlag = household.HasChildren.ToFlag();
int onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
// person inputs
int adultMaleFlag = person.IsAdultMale.ToFlag();
int adultFemaleFlag = person.IsAdultFemale.ToFlag();
int partTimeWorkerFlag = person.IsPartTimeWorker.ToFlag();
int retiredAdultFlag = person.IsRetiredAdult.ToFlag();
int drivingAgeStudentFlag = person.IsDrivingAgeStudent.ToFlag();
int nonworkingAdultFlag = person.IsNonworkingAdult.ToFlag();
int childAge5Through15Flag = person.IsChildAge5Through15.ToFlag();
int childUnder5Flag = person.IsChildUnder5.ToFlag();
// person-day inputs
int homeBasedTours = personDay.HomeBasedTours;
int simulatedToursFlag = personDay.SimulatedToursExist().ToFlag();
int simulatedWorkStops = personDay.SimulatedWorkStops;
int simulatedWorkStopsFlag = personDay.SimulatedWorkStopsExist().ToFlag();
int simulatedSchoolStops = personDay.SimulatedSchoolStops;
// var simulatedSchoolStopsFlag = personDay.HasSimulatedSchoolStops.ToFlag();
int simulatedEscortStops = personDay.SimulatedEscortStops;
// var simulatedEscortStopsFlag = personDay.HasSimulatedEscortStops.ToFlag();
int simulatedPersonalBusinessStops = personDay.SimulatedPersonalBusinessStops;
// var simulatedPersonalBusinessStopsFlag = personDay.HasSimulatedPersonalBusinessStops.ToFlag();
int simulatedShoppingStops = personDay.SimulatedShoppingStops;
// var simulatedShoppingStopsFlag = personDay.HasSimulatedShoppingStops.ToFlag();
int simulatedMealStops = personDay.SimulatedMealStops;
// var simulatedMealStopsFlag = personDay.HasSimulatedMealStops.ToFlag();
int simulatedSocialStops = personDay.SimulatedSocialStops;
// var simulatedSocialStopsFlag = personDay.HasSimulatedSocialStops.ToFlag();
// tour inputs
int hov2TourFlag = tour.IsHov2Mode().ToFlag();
int hov3TourFlag = tour.IsHov3Mode().ToFlag();
int transitTourFlag = tour.IsTransitMode().ToFlag();
int bikeTourFlag = tour.IsBikeMode().ToFlag();
int walkTourFlag = tour.IsWalkMode().ToFlag();
int notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
int workTourFlag = tour.IsWorkPurpose().ToFlag();
int personalBusinessOrMedicalTourFlag = tour.IsPersonalBusinessOrMedicalPurpose().ToFlag();
int socialOrRecreationTourFlag = tour.IsSocialOrRecreationPurpose().ToFlag();
int schoolTourFlag = tour.IsSchoolPurpose().ToFlag();
int escortTourFlag = tour.IsEscortPurpose().ToFlag();
int shoppingTourFlag = tour.IsShoppingPurpose().ToFlag();
int mealTourFlag = tour.IsMealPurpose().ToFlag();
// trip inputs
int oneSimulatedTripFlag = halfTour.OneSimulatedTripFlag;
int twoSimulatedTripsFlag = halfTour.TwoSimulatedTripsFlag;
int threeSimulatedTripsFlag = halfTour.ThreeSimulatedTripsFlag;
int fourSimulatedTripsFlag = halfTour.FourSimulatedTripsFlag;
int fiveSimulatedTripsFlag = halfTour.FiveSimulatedTripsFlag;
int halfTourFromOriginFlag = trip.IsHalfTourFromOrigin.ToFlag();
int halfTourFromDestinationFlag = (!trip.IsHalfTourFromOrigin).ToFlag();
int beforeMandatoryDestinationFlag = trip.IsBeforeMandatoryDestination().ToFlag();
// remaining inputs
int time = trip.IsHalfTourFromOrigin ? tour.DestinationArrivalTime : tour.DestinationDepartureTime;
int from7AMto9AMFlag = (time >= Global.Settings.Times.SevenAM && time < Global.Settings.Times.NineAM).ToFlag();
int from9AMto11AMFlag = (time >= Global.Settings.Times.NineAM && time < Global.Settings.Times.ElevenAM).ToFlag();
int from11AMto1PMFlag = (time >= Global.Settings.Times.ElevenAM && time < Global.Settings.Times.OnePM).ToFlag();
int from1PMto3PMFlag = (time >= Global.Settings.Times.OnePM && time < Global.Settings.Times.ThreePM).ToFlag();
int from3PMto5PMFlag = (time >= Global.Settings.Times.ThreePM && time < Global.Settings.Times.FivePM).ToFlag();
int from7PMto9PMFlag = (time >= Global.Settings.Times.SevenPM && time < Global.Settings.Times.NinePM).ToFlag();
int from9PMto11PMFlag = (time >= Global.Settings.Times.NinePM && time < Global.Settings.Times.ElevenPM).ToFlag();
int from11PMto7AMFlag = (time >= Global.Settings.Times.ElevenPM).ToFlag();
double foodRetailServiceMedicalQtrMileLog = household.ResidenceParcel.FoodRetailServiceMedicalLogBuffer1();
int remainingToursCount = personDay.HomeBasedTours - personDay.GetTotalSimulatedTours();
// time-of-day for trip window calculations
int startTime = trip.GetStartTime();
// time window in minutes and in hours for yet unmodeled portion of half-tour, only consider persons on this trip
int timeWindow =
tour.IsHomeBasedTour
? personDay.TimeWindow.AvailableWindow(startTime, Global.Settings.TimeDirections.Both)
: tour.ParentTour.TimeWindow.AvailableWindow(startTime, Global.Settings.TimeDirections.Both);
double duration = timeWindow / 60D;
// connectivity attributes
double c34Ratio = trip.OriginParcel.C34RatioBuffer1();
// 0 - NO MORE STOPS
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.NoneOrHome, true, choice == Global.Settings.Purposes.NoneOrHome);
alternative.Choice = Global.Settings.Purposes.NoneOrHome;
alternative.AddUtilityTerm(1, twoSimulatedTripsFlag * halfTourFromOriginFlag);
alternative.AddUtilityTerm(2, threeSimulatedTripsFlag * halfTourFromOriginFlag);
alternative.AddUtilityTerm(3, fourSimulatedTripsFlag * halfTourFromOriginFlag);
alternative.AddUtilityTerm(4, fiveSimulatedTripsFlag * halfTourFromOriginFlag);
alternative.AddUtilityTerm(5, twoSimulatedTripsFlag * halfTourFromDestinationFlag);
alternative.AddUtilityTerm(6, threeSimulatedTripsFlag * halfTourFromDestinationFlag);
alternative.AddUtilityTerm(7, fourSimulatedTripsFlag * halfTourFromDestinationFlag);
alternative.AddUtilityTerm(8, fiveSimulatedTripsFlag * halfTourFromDestinationFlag);
alternative.AddUtilityTerm(9, homeBasedTours);
alternative.AddUtilityTerm(10, simulatedToursFlag);
alternative.AddUtilityTerm(11, notHomeBasedTourFlag);
alternative.AddUtilityTerm(12, beforeMandatoryDestinationFlag);
alternative.AddUtilityTerm(16, (transitTourFlag + walkTourFlag + bikeTourFlag) * c34Ratio * foodRetailServiceMedicalQtrMileLog);
alternative.AddUtilityTerm(17, transitTourFlag);
alternative.AddUtilityTerm(40, schoolTourFlag);
// 1 - WORK STOP
if (personDay.WorkStops > 0 && tour.DestinationPurpose <= Global.Settings.Purposes.School)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Work, true, choice == Global.Settings.Purposes.Work);
alternative.Choice = Global.Settings.Purposes.Work;
alternative.AddUtilityTerm(33, workTourFlag + schoolTourFlag);
// alternative.AddUtility(47, escortTourFlag);
// alternative.AddUtility(54, personalBusinessOrMedicalTourFlag);
// alternative.AddUtility(61, shoppingTourFlag);
// alternative.AddUtility(68, mealTourFlag);
// alternative.AddUtility(75, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(82, halfTourFromOriginFlag);
alternative.AddUtilityTerm(89, simulatedWorkStops);
alternative.AddUtilityTerm(96, simulatedWorkStopsFlag);
alternative.AddUtilityTerm(103, remainingToursCount);
alternative.AddUtilityTerm(110, duration);
alternative.AddUtilityTerm(131, from9AMto11AMFlag + from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(175, adultMaleFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Work, false, choice == Global.Settings.Purposes.Work);
}
// 2 - SCHOOL STOP
if (personDay.SchoolStops > 0 && tour.DestinationPurpose <= Global.Settings.Purposes.School)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.School, true, choice == Global.Settings.Purposes.School);
alternative.Choice = Global.Settings.Purposes.School;
alternative.AddUtilityTerm(34, workTourFlag + schoolTourFlag);
// alternative.AddUtility(48, escortTourFlag);
// alternative.AddUtility(55, personalBusinessOrMedicalTourFlag);
// alternative.AddUtility(62, shoppingTourFlag);
// alternative.AddUtility(69, mealTourFlag);
// alternative.AddUtility(76, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(83, halfTourFromOriginFlag);
alternative.AddUtilityTerm(90, simulatedSchoolStops);
// alternative.AddUtility(97, simulatedSchoolStopsFlag);
alternative.AddUtilityTerm(104, remainingToursCount);
alternative.AddUtilityTerm(111, duration);
alternative.AddUtilityTerm(138, from7AMto9AMFlag + from7PMto9PMFlag + from9PMto11PMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(237, oneSimulatedTripFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.School, false, choice == Global.Settings.Purposes.School);
}
// 3 - ESCORT STOP
if (personDay.EscortStops > 0 && tour.DestinationPurpose <= Global.Settings.Purposes.Escort)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Escort, true, choice == Global.Settings.Purposes.Escort);
alternative.Choice = Global.Settings.Purposes.Escort;
alternative.AddUtilityTerm(35, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(49, escortTourFlag);
// alternative.AddUtility(56, personalBusinessOrMedicalTourFlag);
// alternative.AddUtility(63, shoppingTourFlag);
// alternative.AddUtility(70, mealTourFlag);
// alternative.AddUtility(77, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(84, halfTourFromOriginFlag);
alternative.AddUtilityTerm(91, simulatedEscortStops);
// alternative.AddUtility(98, simulatedEscortStopsFlag);
alternative.AddUtilityTerm(105, remainingToursCount);
alternative.AddUtilityTerm(112, duration);
alternative.AddUtilityTerm(146, from7AMto9AMFlag);
alternative.AddUtilityTerm(147, from9AMto11AMFlag + from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(181, childrenFlag * adultFemaleFlag);
alternative.AddUtilityTerm(183, hov2TourFlag);
alternative.AddUtilityTerm(184, hov3TourFlag);
alternative.AddUtilityTerm(238, oneSimulatedTripFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Escort, false, choice == Global.Settings.Purposes.Escort);
}
// 4 - PERSONAL BUSINESS STOP
if (personDay.PersonalBusinessStops > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.PersonalBusiness, true, choice == Global.Settings.Purposes.PersonalBusiness);
alternative.Choice = Global.Settings.Purposes.PersonalBusiness;
alternative.AddUtilityTerm(36, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(50, escortTourFlag);
alternative.AddUtilityTerm(57, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(64, shoppingTourFlag);
alternative.AddUtilityTerm(71, mealTourFlag);
alternative.AddUtilityTerm(78, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(85, halfTourFromOriginFlag);
alternative.AddUtilityTerm(92, simulatedPersonalBusinessStops);
// alternative.AddUtility(99, simulatedPersonalBusinessStopsFlag);
alternative.AddUtilityTerm(106, remainingToursCount);
alternative.AddUtilityTerm(113, duration);
alternative.AddUtilityTerm(154, from7AMto9AMFlag + from7PMto9PMFlag + from9PMto11PMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(155, from9AMto11AMFlag + from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(195, onePersonHouseholdFlag);
alternative.AddUtilityTerm(196, hov2TourFlag);
alternative.AddUtilityTerm(197, hov3TourFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.PersonalBusiness, false, choice == Global.Settings.Purposes.PersonalBusiness);
}
// 5 - SHOPPING STOP
if (personDay.ShoppingStops > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Shopping, true, choice == Global.Settings.Purposes.Shopping);
alternative.Choice = Global.Settings.Purposes.Shopping;
alternative.AddUtilityTerm(37, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(51, escortTourFlag);
alternative.AddUtilityTerm(58, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(65, shoppingTourFlag);
alternative.AddUtilityTerm(72, mealTourFlag);
alternative.AddUtilityTerm(79, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(86, halfTourFromOriginFlag);
alternative.AddUtilityTerm(93, simulatedShoppingStops);
// alternative.AddUtility(100, simulatedShoppingStopsFlag);
alternative.AddUtilityTerm(107, remainingToursCount);
alternative.AddUtilityTerm(114, duration);
alternative.AddUtilityTerm(162, from7AMto9AMFlag + from9PMto11PMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(164, from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(207, childrenFlag * adultFemaleFlag);
alternative.AddUtilityTerm(209, hov2TourFlag);
alternative.AddUtilityTerm(210, hov3TourFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Shopping, false, choice == Global.Settings.Purposes.Shopping);
}
// 6 - MEAL STOP
if (personDay.MealStops > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Meal, true, choice == Global.Settings.Purposes.Meal);
alternative.Choice = Global.Settings.Purposes.Meal;
alternative.AddUtilityTerm(38, workTourFlag);
alternative.AddUtilityTerm(46, schoolTourFlag);
alternative.AddUtilityTerm(52, escortTourFlag);
alternative.AddUtilityTerm(59, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(66, shoppingTourFlag);
alternative.AddUtilityTerm(73, mealTourFlag);
alternative.AddUtilityTerm(80, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(87, halfTourFromOriginFlag);
alternative.AddUtilityTerm(94, simulatedMealStops);
// alternative.AddUtility(101, simulatedMealStopsFlag);
alternative.AddUtilityTerm(108, remainingToursCount);
alternative.AddUtilityTerm(115, duration);
alternative.AddUtilityTerm(170, from7AMto9AMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(171, from11AMto1PMFlag + from1PMto3PMFlag);
alternative.AddUtilityTerm(172, from7PMto9PMFlag);
alternative.AddUtilityTerm(221, onePersonHouseholdFlag);
alternative.AddUtilityTerm(222, hov2TourFlag);
alternative.AddUtilityTerm(223, hov3TourFlag);
alternative.AddUtilityTerm(226, partTimeWorkerFlag + retiredAdultFlag + drivingAgeStudentFlag);
alternative.AddUtilityTerm(228, nonworkingAdultFlag + childAge5Through15Flag + childUnder5Flag);
alternative.AddUtilityTerm(239, oneSimulatedTripFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Meal, false, choice == Global.Settings.Purposes.Meal);
}
// 7 - SOCIAL (OR RECREATION) STOP
if (personDay.SocialStops > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Social, true, choice == Global.Settings.Purposes.Social);
alternative.Choice = Global.Settings.Purposes.Social;
alternative.AddUtilityTerm(39, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(53, escortTourFlag);
alternative.AddUtilityTerm(60, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(67, shoppingTourFlag);
alternative.AddUtilityTerm(74, mealTourFlag);
alternative.AddUtilityTerm(81, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(88, halfTourFromOriginFlag);
alternative.AddUtilityTerm(95, simulatedSocialStops);
// alternative.AddUtility(102, simulatedSocialStopsFlag);
alternative.AddUtilityTerm(109, remainingToursCount);
alternative.AddUtilityTerm(116, duration);
alternative.AddUtilityTerm(173, from7AMto9AMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(174, from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(235, hov2TourFlag);
alternative.AddUtilityTerm(236, hov3TourFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Social, false, choice == Global.Settings.Purposes.Social);
}
}
public int Run(ITripWrapper trip)
{
return(Run(trip, Global.Settings.Purposes.NoneOrHome));
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITripWrapper trip, IEnumerable <dynamic> pathTypeModels, IParcelWrapper originParcel, IParcelWrapper destinationParcel, int departureTime, int choice = Constants.DEFAULT_VALUE)
{
var household = trip.Household;
var householdTotals = household.HouseholdTotals;
var person = trip.Person;
var tour = trip.Tour;
var halfTour = trip.HalfTour;
// household inputs
var onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
var twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
var income0To25KFlag = household.Has0To25KIncome.ToFlag();
var income25To45KFlag = household.Has25To45KIncome.ToFlag();
var childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
var nonworkingAdults = householdTotals.NonworkingAdults;
var retiredAdults = householdTotals.RetiredAdults;
var noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(household.VehiclesAvailable);
var carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(household.VehiclesAvailable);
// person inputs
var drivingAgeStudentFlag = person.IsDrivingAgeStudent.ToFlag();
var maleFlag = person.IsMale.ToFlag();
var ageLessThan35Flag = person.AgeIsLessThan35.ToFlag();
// tour inputs
var parkAndRideTourFlag = tour.IsParkAndRideMode().ToFlag();
var transitTourFlag = tour.IsTransitMode().ToFlag();
var schoolBusTourFlag = tour.IsSchoolBusMode().ToFlag();
var hov3TourFlag = tour.IsHov3Mode().ToFlag();
var hov2TourFlag = tour.IsHov2Mode().ToFlag();
var sovTourFlag = tour.IsSovMode().ToFlag();
var bikeTourFlag = tour.IsBikeMode().ToFlag();
var walkTourFlag = tour.IsWalkMode().ToFlag();
var homeBasedWorkTourFlag = (tour.IsHomeBasedTour && tour.IsWorkPurpose()).ToFlag();
var homeBasedSchoolTourFlag = (tour.IsHomeBasedTour && tour.IsSchoolPurpose()).ToFlag();
var homeBasedEscortTourFlag = (tour.IsHomeBasedTour && tour.IsEscortPurpose()).ToFlag();
var homeBasedShoppingTourFlag = (tour.IsHomeBasedTour && tour.IsShoppingPurpose()).ToFlag();
var homeBasedMealTourFlag = (tour.IsHomeBasedTour && tour.IsMealPurpose()).ToFlag();
var homeBasedSocialTourFlag = (tour.IsHomeBasedTour && tour.IsSocialPurpose()).ToFlag();
var notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
// var homeBasedNotWorkSchoolEscortTourFlag = (tour.IsHomeBasedTour && tour.DestinationPurpose > Global.Settings.Purposes.Escort).ToFlag();
// var costco = -0.1432 * homeBasedWorkTourFlag - 0.2000 * homeBasedSchoolTourFlag - 0.3000 * homeBasedEscortTourFlag - 0.2000 * homeBasedNotWorkSchoolEscortTourFlag - 0.2000 * notHomeBasedTourFlag;
// trip inputs
var originHomeEscortFlag = (trip.IsNoneOrHomePurposeByOrigin() && trip.IsEscortPurposeByDestination()).ToFlag();
var originWorkEscortFlag = (trip.IsWorkPurposeByOrigin() && trip.IsEscortPurposeByDestination()).ToFlag();
var destinationHomeEscortFlag = (trip.IsNoneOrHomePurposeByDestination() && trip.IsEscortPurposeByOrigin()).ToFlag();
var destinationWorkEscortFlag = (trip.IsWorkPurposeByDestination() && trip.IsEscortPurposeByOrigin()).ToFlag();
// only trip on first half-tour
var onlyTripOnFirstHalfFlag = (trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && trip.IsToTourOrigin).ToFlag();
// first trip on first half-tour, not only one
var firstTripOnFirstHalfFlag = (trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && !trip.IsToTourOrigin).ToFlag();
// last trip first half-tour, not only one
var lastTripOnFirstHalfFlag = (trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips > 1 && trip.IsToTourOrigin).ToFlag();
// only trip on second half-tour
var onlyTripOnSecondHalfFlag = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && trip.IsToTourOrigin).ToFlag();
// first trip on second half-tour, not only one
var firstTripOnSecondHalfFlag = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && !trip.IsToTourOrigin).ToFlag();
// last trip second half-tour, not only one
var lastTripOnSecondHalfFlag = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips > 1 && trip.IsToTourOrigin).ToFlag();
// remaining inputs
var originMixedDensity = originParcel.MixedUse4Index1();
var originIntersectionDensity = originParcel.NetIntersectionDensity1();
var destinationParkingCost = destinationParcel.ParkingCostBuffer1(2);
var amPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.SixAM, Global.Settings.Times.TenAM).ToFlag();
var middayPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.ThreePM).ToFlag();
var pmPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.ThreePM, Global.Settings.Times.SevenPM).ToFlag();
var eveningPeriodFlag = (departureTime > Global.Settings.Times.SevenPM).ToFlag();
// availability
var tripModeAvailable = new bool[Global.Settings.Modes.TotalModes];
var isLastTripInTour = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && trip.IsToTourOrigin) || (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips > 1 && trip.IsToTourOrigin);
var frequencyPreviousTripModeIsTourMode =
trip.IsHalfTourFromOrigin
? tour.HalfTourFromOrigin.Trips.Where(x => x.Sequence < trip.Sequence).Count(x => tour.Mode == x.Mode)
: tour.HalfTourFromOrigin.Trips.Union(tour.HalfTourFromDestination.Trips.Where(x => x.Sequence < trip.Sequence)).Count(x => tour.Mode == x.Mode);
// if a park and ride tour, only car is available
if (tour.Mode == Global.Settings.Modes.ParkAndRide)
{
tripModeAvailable[Global.Settings.Modes.Sov] = tour.Household.VehiclesAvailable > 0 && tour.Person.IsDrivingAge;
tripModeAvailable[Global.Settings.Modes.Hov2] = !tripModeAvailable[Global.Settings.Modes.Sov];
}
// if the last trip of the tour and tour mode not yet used, only the tour mode is available
else if (isLastTripInTour && frequencyPreviousTripModeIsTourMode == 0)
{
tripModeAvailable[tour.Mode] = true;
}
else
{
// set availability based on tour mode
for (var mode = Global.Settings.Modes.Walk; mode <= tour.Mode; mode++)
{
tripModeAvailable[mode] = true;
}
}
// school bus is a special case - use HOV3 impedance and only available for school bus tours
var pathTypeExtra = pathTypeModels.First(x => x.Mode == Global.Settings.Modes.Hov3);
var modeExtra = Global.Settings.Modes.SchoolBus;
var availableExtra = pathTypeExtra.Available && tour.IsSchoolBusMode() && tripModeAvailable[modeExtra];
var generalizedTimeLogsumExtra = pathTypeExtra.GeneralizedTimeLogsum;
var alternative = choiceProbabilityCalculator.GetAlternative(modeExtra, availableExtra, choice == modeExtra);
alternative.Choice = modeExtra;
alternative.AddNestedAlternative(_nestedAlternativeIds[modeExtra], _nestedAlternativeIndexes[modeExtra], THETA_PARAMETER);
if (availableExtra)
{
// case Global.Settings.Modes.SchoolBus:
alternative.AddUtilityTerm(2, generalizedTimeLogsumExtra * tour.TimeCoefficient);
alternative.AddUtilityTerm(18, 1);
alternative.AddUtilityTerm(100, schoolBusTourFlag);
alternative.AddUtilityTerm(102, (schoolBusTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (schoolBusTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (schoolBusTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (schoolBusTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (schoolBusTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (schoolBusTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(112, parkAndRideTourFlag);
alternative.AddUtilityTerm(113, transitTourFlag);
}
foreach (var pathTypeModel in pathTypeModels)
{
var mode = pathTypeModel.Mode;
var available = pathTypeModel.Available && tripModeAvailable[mode];
var generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
alternative.AddNestedAlternative(_nestedAlternativeIds[mode], _nestedAlternativeIndexes[mode], THETA_PARAMETER);
if (!available)
{
continue;
}
alternative.AddUtilityTerm(2, generalizedTimeLogsum * tour.TimeCoefficient);
if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(20, 1);
alternative.AddUtilityTerm(22, carsLessThanDriversFlag);
alternative.AddUtilityTerm(100, transitTourFlag);
alternative.AddUtilityTerm(102, (transitTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (transitTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (transitTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (transitTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (transitTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (transitTourFlag * lastTripOnSecondHalfFlag));
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(34, (nonworkingAdults + retiredAdults));
alternative.AddUtilityTerm(36, onePersonHouseholdFlag);
alternative.AddUtilityTerm(37, twoPersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(100, hov3TourFlag);
alternative.AddUtilityTerm(102, (hov3TourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (hov3TourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (hov3TourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (hov3TourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (hov3TourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (hov3TourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(114, parkAndRideTourFlag);
alternative.AddUtilityTerm(115, transitTourFlag);
alternative.AddUtilityTerm(116, schoolBusTourFlag);
alternative.AddUtilityTerm(149, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(150, homeBasedSchoolTourFlag);
alternative.AddUtilityTerm(152, homeBasedEscortTourFlag);
alternative.AddUtilityTerm(153, homeBasedShoppingTourFlag);
alternative.AddUtilityTerm(154, homeBasedMealTourFlag);
alternative.AddUtilityTerm(155, homeBasedSocialTourFlag);
alternative.AddUtilityTerm(161, (destinationWorkEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(162, (originWorkEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(163, (originHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(164, (originHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(165, (originHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(166, (originHomeEscortFlag * eveningPeriodFlag));
alternative.AddUtilityTerm(167, (destinationHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(168, (destinationHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(169, (destinationHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(170, (destinationHomeEscortFlag * eveningPeriodFlag));
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(32, (childrenAge5Through15 * (1 - homeBasedEscortTourFlag)));
alternative.AddUtilityTerm(34, ((nonworkingAdults + retiredAdults) * (1 - homeBasedEscortTourFlag)));
alternative.AddUtilityTerm(38, onePersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(100, hov2TourFlag);
alternative.AddUtilityTerm(102, (hov2TourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (hov2TourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (hov2TourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (hov2TourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (hov2TourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (hov2TourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(118, transitTourFlag);
alternative.AddUtilityTerm(119, schoolBusTourFlag);
alternative.AddUtilityTerm(120, hov3TourFlag);
alternative.AddUtilityTerm(149, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(150, homeBasedSchoolTourFlag);
alternative.AddUtilityTerm(152, homeBasedEscortTourFlag);
alternative.AddUtilityTerm(153, homeBasedShoppingTourFlag);
alternative.AddUtilityTerm(154, homeBasedMealTourFlag);
alternative.AddUtilityTerm(155, homeBasedSocialTourFlag);
alternative.AddUtilityTerm(161, (destinationWorkEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(162, (originWorkEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(163, (originHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(164, (originHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(165, (originHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(166, (originHomeEscortFlag * eveningPeriodFlag));
alternative.AddUtilityTerm(167, (destinationHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(168, (destinationHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(169, (destinationHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(170, (destinationHomeEscortFlag * eveningPeriodFlag));
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient));
alternative.AddUtilityTerm(52, carsLessThanDriversFlag);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(55, income25To45KFlag);
alternative.AddUtilityTerm(59, drivingAgeStudentFlag);
alternative.AddUtilityTerm(100, sovTourFlag);
alternative.AddUtilityTerm(102, (sovTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (sovTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (sovTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (sovTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (sovTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (sovTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(121, parkAndRideTourFlag);
alternative.AddUtilityTerm(122, transitTourFlag);
alternative.AddUtilityTerm(124, hov3TourFlag);
alternative.AddUtilityTerm(125, hov2TourFlag);
}
else if (mode == Global.Settings.Modes.Bike)
{
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(62, ageLessThan35Flag);
alternative.AddUtilityTerm(65, originIntersectionDensity);
alternative.AddUtilityTerm(100, bikeTourFlag);
alternative.AddUtilityTerm(102, (bikeTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (bikeTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (bikeTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (bikeTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (bikeTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (bikeTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(127, transitTourFlag);
alternative.AddUtilityTerm(128, schoolBusTourFlag);
alternative.AddUtilityTerm(130, hov2TourFlag);
alternative.AddUtilityTerm(131, sovTourFlag);
alternative.AddUtilityTerm(147, notHomeBasedTourFlag);
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(72, ageLessThan35Flag);
alternative.AddUtilityTerm(75, originIntersectionDensity);
alternative.AddUtilityTerm(78, originMixedDensity); // origin and destination mixed use measures - geometric avg. - half mile from cell, in 1000s
alternative.AddUtilityTerm(100, walkTourFlag);
alternative.AddUtilityTerm(102, (walkTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (walkTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (walkTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (walkTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (walkTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (walkTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(141, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(142, homeBasedSchoolTourFlag);
}
}
}
public void Run(ITripWrapper trip)
{
if (trip == null)
{
throw new ArgumentNullException("trip");
}
TripTime.InitializeTripTimes();
trip.PersonDay.ResetRandom(40 * (2 * trip.Tour.Sequence - 1 + trip.Direction - 1) + 50 + trip.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
}
var choiceProbabilityCalculator = _helpers[ParallelUtility.GetBatchFromThreadId()].GetChoiceProbabilityCalculator(trip.Id);
if (_helpers[ParallelUtility.GetBatchFromThreadId()].ModelIsInEstimationMode)
{
if (trip.DestinationParcel == null || trip.OriginParcel == null || trip.Mode <= Global.Settings.Modes.None || trip.Mode == Global.Settings.Modes.Other)
{
return;
}
RunModel(choiceProbabilityCalculator, trip, new TripTime(trip.DepartureTime));
choiceProbabilityCalculator.WriteObservation();
}
else
{
RunModel(choiceProbabilityCalculator, trip);
var chosenAlternative = choiceProbabilityCalculator.SimulateChoice(trip.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", trip.PersonDay.Id);
trip.PersonDay.IsValid = false;
return;
}
var choice = (TripTime)chosenAlternative.Choice;
var departureTime = choice.GetDepartureTime(trip);
trip.DepartureTime = departureTime;
if (departureTime >= 1 && departureTime <= Global.Settings.Times.MinutesInADay)
{
trip.UpdateTripValues();
}
else
{
if (!Global.Configuration.IsInEstimationMode)
{
trip.PersonDay.IsValid = false;
}
}
}
}
ITripWrapper ITripCreator.CreateWrapper(ITourWrapper tourWrapper, ITripWrapper trip, int nextTripId, int intermediateStopPurpose, int destinationPurpose, IHalfTour halfTour)
{
return(CreateWrapper(tourWrapper, trip, nextTripId, intermediateStopPurpose, destinationPurpose, halfTour));
}
public void Run(IHouseholdDayWrapper householdDay, ITripWrapper trip)
{
if (trip == null)
{
throw new ArgumentNullException("trip");
}
trip.PersonDay.ResetRandom(40 * (2 * trip.Tour.Sequence - 1 + trip.Direction - 1) + 40 + trip.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
}
ChoiceProbabilityCalculator choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(trip.Id);
IParcelWrapper originParcel =
trip.IsHalfTourFromOrigin
? trip.DestinationParcel
: trip.OriginParcel;
// for skims - use actual travel direction, not simulation direction
IParcelWrapper destinationParcel =
trip.IsHalfTourFromOrigin
? trip.OriginParcel
: trip.DestinationParcel;
int departureTime = trip.IsHalfTourFromOrigin ? trip.LatestDepartureTime : trip.EarliestDepartureTime;
if (departureTime < 1)
{
Global.PrintFile.WriteLine("From origin / latest / earliest {0} {1} {2}", trip.IsHalfTourFromOrigin, trip.LatestDepartureTime, trip.EarliestDepartureTime);
}
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
if (destinationParcel == null || originParcel == null || trip.Mode <= Global.Settings.Modes.None || trip.Mode == Global.Settings.Modes.ParkAndRide || trip.Mode == Global.Settings.Modes.Other)
{
return;
}
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.RunAll(
trip.Household.RandomUtility,
originParcel,
destinationParcel,
departureTime,
0,
trip.Tour.DestinationPurpose,
trip.Tour.CostCoefficient,
trip.Tour.TimeCoefficient,
trip.Person.IsDrivingAge,
trip.Household.VehiclesAvailable,
trip.Person.TransitPassOwnership,
trip.Household.OwnsAutomatedVehicles > 0,
trip.Person.GetTransitFareDiscountFraction(),
false);
// there is no path type model for school bus, use HOV3
int mode = trip.Mode == Global.Settings.Modes.SchoolBus ? Global.Settings.Modes.Hov3 : trip.Mode;
IPathTypeModel pathTypeModel = pathTypeModels.First(x => x.Mode == mode);
if (!pathTypeModel.Available)
{
return;
}
RunModel(choiceProbabilityCalculator, trip, pathTypeModels, originParcel, destinationParcel, trip.Mode);
choiceProbabilityCalculator.WriteObservation();
}
else
{
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.RunAll(
trip.Household.RandomUtility,
originParcel,
destinationParcel,
departureTime,
0,
trip.Tour.DestinationPurpose,
trip.Tour.CostCoefficient,
trip.Tour.TimeCoefficient,
trip.Person.IsDrivingAge,
trip.Household.VehiclesAvailable,
trip.Person.TransitPassOwnership,
trip.Household.OwnsAutomatedVehicles > 0,
trip.Person.GetTransitFareDiscountFraction(),
false);
RunModel(choiceProbabilityCalculator, trip, pathTypeModels, originParcel, destinationParcel);
ChoiceProbabilityCalculator.Alternative chosenAlternative = choiceProbabilityCalculator.SimulateChoice(trip.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", trip.PersonDay.Id);
trip.Mode = Global.Settings.Modes.Hov3;
if (!Global.Configuration.IsInEstimationMode)
{
trip.PersonDay.IsValid = false;
}
return;
}
int choice = (int)chosenAlternative.Choice;
trip.Mode = choice;
if (choice == Global.Settings.Modes.SchoolBus)
{
trip.PathType = 0;
}
else
{
IPathTypeModel chosenPathType = pathTypeModels.First(x => x.Mode == choice);
trip.PathType = chosenPathType.PathType;
// for transit trips, overwrite origin and destination zones with stop area ids
if (Global.StopAreaIsEnabled && choice == Global.Settings.Modes.Transit &&
Global.Configuration.WriteStopAreaIDsInsteadOfZonesForTransitTrips)
{
trip.OriginZoneKey = chosenPathType.PathOriginStopAreaKey;
trip.DestinationZoneKey = chosenPathType.PathDestinationStopAreaKey;
}
}
}
}
private static void SetTripModeAndTime(ITourWrapper tour, ITripWrapper trip, int batchNumber)
{
if (Global.Configuration.ShouldRunTripModeModel)
{
// sets the trip's mode of travel to the destination
if (trip.DestinationPurpose == Global.Settings.Purposes.ChangeMode)
{
// trips to change mode destination are always by transit
ChoiceModelFactory.TotalTimesChangeModeTransitModeSet[batchNumber]++;
trip.Mode = Global.Settings.Modes.Transit;
if (Global.StopAreaIsEnabled && Global.Configuration.WriteStopAreaIDsInsteadOfZonesForTransitTrips)
{
trip.OriginZoneKey = trip.Tour.ParkAndRideOriginStopAreaKey;
trip.DestinationZoneKey = trip.Tour.ParkAndRideDestinationStopAreaKey;
}
}
else
{
ChoiceModelFactory.TotalTimesTripModeModelRun[batchNumber]++;
Global.ChoiceModelSession.Get <TripModeModel>().Run(trip);
}
if (!trip.PersonDay.IsValid)
{
return;
}
}
// sets the trip's destination arrival and departure times
if (trip.Sequence == 1)
{
if (!Global.Configuration.IsInEstimationMode)
{
trip.DepartureTime = trip.IsHalfTourFromOrigin ? tour.DestinationArrivalTime : tour.DestinationDepartureTime;
trip.UpdateTripValues();
}
}
else if (trip.OriginPurpose == Global.Settings.Purposes.ChangeMode)
{
//stay at park and ride lot assumed to be 3 minutes
if (!Global.Configuration.IsInEstimationMode)
{
int endpoint;
if (trip.IsHalfTourFromOrigin)
{
trip.DepartureTime = trip.GetPreviousTrip().ArrivalTime - 3;
endpoint = trip.DepartureTime + 1;
}
else
{
trip.DepartureTime = trip.GetPreviousTrip().ArrivalTime + 3;
endpoint = trip.DepartureTime - 1;
}
if (trip.DepartureTime >= 1 && trip.DepartureTime <= Global.Settings.Times.MinutesInADay && trip.PersonDay.TimeWindow.EntireSpanIsAvailable(endpoint, trip.DepartureTime))
{
trip.UpdateTripValues();
}
else
{
if (!Global.Configuration.IsInEstimationMode)
{
trip.PersonDay.IsValid = false;
}
}
}
}
else
{
if (Global.Configuration.ShouldRunTripTimeModel)
{
ChoiceModelFactory.TotalTimesTripTimeModelRun[batchNumber]++;
Global.ChoiceModelSession.Get <TripTimeModel>().Run(trip);
}
}
}
public void Export(ITripWrapper trip)
{
if (trip == null)
{
throw new ArgumentNullException("trip");
}
if (Global.Configuration.UseCustomTDMTripListFormat)
{
//add custom code to write data items corresponding to header here
_writer.Write(trip.Household.Id);
_writer.Write(_delimiter);
_writer.Write(trip.Person.Sequence);
_writer.Write(_delimiter);
_writer.Write(trip.Tour.Sequence);
_writer.Write(_delimiter);
_writer.Write(trip.IsHalfTourFromOrigin ? trip.Sequence : trip.Tour.HalfTour1Trips + trip.Sequence);
_writer.Write(_delimiter);
_writer.Write(trip.DepartureTime);
_writer.Write(_delimiter);
_writer.Write(trip.ArrivalTime);
_writer.Write(_delimiter);
_writer.Write(trip.ActivityEndTime);
_writer.Write(_delimiter);
_writer.Write(trip.OriginParcel.Id);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationParcel.Id);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationPurpose);
_writer.Write(_delimiter);
_writer.Write(trip.Mode);
_writer.Write(_delimiter);
_writer.Write(trip.DriverType);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationParkingNodeId);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationParkingType);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationParkingCost);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationParkingWalkTime);
_writer.Write(_delimiter);
_writer.Write(trip.Tour.CostCoefficient);
_writer.Write(_delimiter);
_writer.Write(trip.Tour.TimeCoefficient);
_writer.Write(_delimiter);
_writer.Write(trip.Person.PaidParkingAtWorkplace);
_writer.Write(_delimiter);
_writer.WriteLine(trip.ExpansionFactor);
}
else if (Global.Configuration.UseTransimsTDMTripListFormat)
{
_writer.Write(trip.Household.Id);
_writer.Write(_delimiter);
_writer.Write(trip.Person.Sequence);
_writer.Write(_delimiter);
_writer.Write(trip.Tour.Sequence);
_writer.Write(_delimiter);
_writer.Write(trip.IsHalfTourFromOrigin ? trip.Sequence : trip.Tour.HalfTour1Trips + trip.Sequence);
_writer.Write(_delimiter);
var departureTime24Hour = trip.DepartureTime.ToMinutesAfterMidnight();
if (departureTime24Hour < 180)
{
departureTime24Hour += 1440;
} // range becomes 180-1619 instead of 0-1439
_writer.Write(departureTime24Hour);
_writer.Write(_delimiter);
var arrivalTime24Hour = trip.ArrivalTime.ToMinutesAfterMidnight();
if (arrivalTime24Hour < 180)
{
arrivalTime24Hour += 1440;
} // range becomes 180-1619 instead of 0-1439
_writer.Write(arrivalTime24Hour);
_writer.Write(_delimiter);
_writer.Write(trip.ActivityEndTime - trip.ArrivalTime);
_writer.Write(_delimiter);
_writer.Write(trip.OriginParcel.LandUseCode);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationParcel.LandUseCode);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationPurpose);
_writer.Write(_delimiter);
var tsMode = trip.Mode == Global.Settings.Modes.Walk
? 1
: trip.Mode == Global.Settings.Modes.Bike
? 2
: trip.Mode >= Global.Settings.Modes.Sov && trip.Mode <= Global.Settings.Modes.Hov3 && trip.DriverType == Global.Settings.DriverTypes.Driver
? 3
: trip.Mode >= Global.Settings.Modes.Sov && trip.Mode <= Global.Settings.Modes.Hov3 && trip.DriverType != Global.Settings.DriverTypes.Driver
? 4
: trip.Mode == Global.Settings.Modes.Transit
? 5
: trip.Mode == Global.Settings.Modes.SchoolBus
? 11
: 6;
var tsPass = trip.Mode == Global.Settings.Modes.Hov2
? 1
: trip.Mode == Global.Settings.Modes.Hov3
? 2
: 0;
_writer.Write(tsMode);
_writer.Write(_delimiter);
const int zero = 0;
_writer.Write(zero);
_writer.Write(_delimiter);
_writer.Write(zero);
_writer.Write(_delimiter);
_writer.Write(trip.Household.Id * 100 + trip.Person.Sequence);
_writer.Write(_delimiter);
_writer.Write(tsPass);
_writer.Write(_delimiter);
var tripVOT = trip.ValueOfTime;
var votType = (trip.Mode <Global.Settings.Modes.Sov || trip.Mode> Global.Settings.Modes.Hov3)
? 0
: ((tripVOT < 30
? 1 + (int)tripVOT
: tripVOT < 100
? 31 + (int)((tripVOT - 30) / 5)
: 45)
+ (trip.PathType == Global.Settings.PathTypes.NoTolls ? 0 : 50));
_writer.WriteLine(votType);
}
else
{
_writer.Write(trip.Id);
_writer.Write(_delimiter);
_writer.Write(trip.OriginZoneKey);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationZoneKey);
_writer.Write(_delimiter);
_writer.Write(trip.Mode);
_writer.Write(_delimiter);
var departureTime24Hour = trip.DepartureTime.ToMinutesAfterMidnight();
_writer.Write(departureTime24Hour / 60D);
_writer.Write(_delimiter);
var arrivalTime24Hour = trip.ArrivalTime.ToMinutesAfterMidnight();
_writer.Write(arrivalTime24Hour / 60D);
_writer.Write(_delimiter);
_writer.Write(trip.ArrivalTime - trip.DepartureTime);
_writer.Write(_delimiter);
_writer.Write(trip.DestinationPurpose);
_writer.Write(_delimiter);
var valueOfTime =
trip.Household.Income < (trip.Tour.DestinationPurpose == Global.Settings.Purposes.Work ? 20000 : 40000)
? Global.Settings.ValueOfTimes.Low
: trip.Household.Income < (trip.Tour.DestinationPurpose == Global.Settings.Purposes.Work ? 45000 : 110000)
? Global.Settings.ValueOfTimes.Medium
: Global.Settings.ValueOfTimes.High;
_writer.WriteLine(valueOfTime);
}
_current++;
if (_current % 1000 == 0)
{
_writer.Flush();
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITripWrapper trip, IEnumerable <IPathTypeModel> pathTypeModels, IParcelWrapper originParcel, IParcelWrapper destinationParcel, int departureTime, int choice = Constants.DEFAULT_VALUE)
{
IHouseholdWrapper household = trip.Household;
Framework.DomainModels.Models.IHouseholdTotals householdTotals = household.HouseholdTotals;
IPersonWrapper person = trip.Person;
ITourWrapper tour = trip.Tour;
Framework.DomainModels.Models.IHalfTour halfTour = trip.HalfTour;
// household inputs
int onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
int twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
int income0To25KFlag = household.Has0To25KIncome.ToFlag();
int income25To45KFlag = household.Has25To45KIncome.ToFlag();
int incomeOver100Flag = household.Has100KPlusIncome.ToFlag();
int childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
int nonworkingAdults = householdTotals.NonworkingAdults;
int retiredAdults = householdTotals.RetiredAdults;
int noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(household.VehiclesAvailable);
int carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(household.VehiclesAvailable);
// person inputs
int drivingAgeStudentFlag = person.IsDrivingAgeStudent.ToFlag();
int maleFlag = person.IsMale.ToFlag();
int ageLessThan35Flag = person.AgeIsLessThan35.ToFlag();
// tour inputs
int parkAndRideTourFlag = tour.IsParkAndRideMode().ToFlag();
int transitTourFlag = tour.IsTransitMode().ToFlag();
int schoolBusTourFlag = tour.IsSchoolBusMode().ToFlag();
int hov3TourFlag = tour.IsHov3Mode().ToFlag();
int hov2TourFlag = tour.IsHov2Mode().ToFlag();
int sovTourFlag = tour.IsSovMode().ToFlag();
int bikeTourFlag = tour.IsBikeMode().ToFlag();
int walkTourFlag = tour.IsWalkMode().ToFlag();
int tncTourFlag = (tour.Mode == Global.Settings.Modes.PaidRideShare).ToFlag();
int homeBasedWorkTourFlag = (tour.IsHomeBasedTour && tour.IsWorkPurpose()).ToFlag();
int homeBasedSchoolTourFlag = (tour.IsHomeBasedTour && tour.IsSchoolPurpose()).ToFlag();
int homeBasedEscortTourFlag = (tour.IsHomeBasedTour && tour.IsEscortPurpose()).ToFlag();
int homeBasedShoppingTourFlag = (tour.IsHomeBasedTour && tour.IsShoppingPurpose()).ToFlag();
int homeBasedMealTourFlag = (tour.IsHomeBasedTour && tour.IsMealPurpose()).ToFlag();
int homeBasedSocialTourFlag = (tour.IsHomeBasedTour && tour.IsSocialPurpose()).ToFlag();
int notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
// var homeBasedNotWorkSchoolEscortTourFlag = (tour.IsHomeBasedTour && tour.DestinationPurpose > Global.Settings.Purposes.Escort).ToFlag();
// var costco = -0.1432 * homeBasedWorkTourFlag - 0.2000 * homeBasedSchoolTourFlag - 0.3000 * homeBasedEscortTourFlag - 0.2000 * homeBasedNotWorkSchoolEscortTourFlag - 0.2000 * notHomeBasedTourFlag;
// trip inputs
int originHomeEscortFlag = (trip.IsNoneOrHomePurposeByOrigin() && trip.IsEscortPurposeByDestination()).ToFlag();
int originWorkEscortFlag = (trip.IsWorkPurposeByOrigin() && trip.IsEscortPurposeByDestination()).ToFlag();
int destinationHomeEscortFlag = (trip.IsNoneOrHomePurposeByDestination() && trip.IsEscortPurposeByOrigin()).ToFlag();
int destinationWorkEscortFlag = (trip.IsWorkPurposeByDestination() && trip.IsEscortPurposeByOrigin()).ToFlag();
// only trip on first half-tour
int onlyTripOnFirstHalfFlag = (trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && trip.IsToTourOrigin).ToFlag();
// first trip on first half-tour, not only one
int firstTripOnFirstHalfFlag = (trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && !trip.IsToTourOrigin).ToFlag();
// last trip first half-tour, not only one
int lastTripOnFirstHalfFlag = (trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips > 1 && trip.IsToTourOrigin).ToFlag();
// only trip on second half-tour
int onlyTripOnSecondHalfFlag = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && trip.IsToTourOrigin).ToFlag();
// first trip on second half-tour, not only one
int firstTripOnSecondHalfFlag = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && !trip.IsToTourOrigin).ToFlag();
// last trip second half-tour, not only one
int lastTripOnSecondHalfFlag = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips > 1 && trip.IsToTourOrigin).ToFlag();
// remaining inputs
double originMixedDensity = originParcel.MixedUse4Index1();
double originIntersectionDensity = originParcel.NetIntersectionDensity1();
double destinationParkingCost = destinationParcel.ParkingCostBuffer1(2);
int amPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.SixAM, Global.Settings.Times.TenAM).ToFlag();
int middayPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.ThreePM).ToFlag();
int pmPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.ThreePM, Global.Settings.Times.SevenPM).ToFlag();
int eveningPeriodFlag = (departureTime > Global.Settings.Times.SevenPM).ToFlag();
// availability
bool[] tripModeAvailable = new bool[Global.Settings.Modes.TotalModes + 1];
bool isLastTripInTour = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips == 1 && trip.IsToTourOrigin) || (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips > 1 && trip.IsToTourOrigin);
int frequencyPreviousTripModeIsTourMode =
trip.IsHalfTourFromOrigin
? tour.HalfTourFromOrigin.Trips.Where(x => x.Sequence < trip.Sequence).Count(x => tour.Mode == x.Mode)
: tour.HalfTourFromOrigin.Trips.Union(tour.HalfTourFromDestination.Trips.Where(x => x.Sequence < trip.Sequence)).Count(x => tour.Mode == x.Mode);
// if a park and ride tour, only car is available
int maxAvail = tour.Mode;
if (tour.Mode == Global.Settings.Modes.ParkAndRide)
{
tripModeAvailable[Global.Settings.Modes.Sov] = tour.Household.VehiclesAvailable > 0 && tour.Person.IsDrivingAge;
tripModeAvailable[Global.Settings.Modes.Hov2] = !tripModeAvailable[Global.Settings.Modes.Sov];
}
// if the last trip of the tour and tour mode not yet used, only the tour mode is available
else if (isLastTripInTour && frequencyPreviousTripModeIsTourMode == 0 && tour.Mode <= maxAvail)
{
tripModeAvailable[tour.Mode] = true;
}
else
{
// set availability based on tour mode
for (int mode = Global.Settings.Modes.Walk; mode <= maxAvail; mode++)
{
if (mode != Global.Settings.Modes.ParkAndRide)
{
tripModeAvailable[mode] = true;
}
}
}
{
// school bus is a special case - use HOV3 impedance and only available for school bus tours
IPathTypeModel pathTypeExtra = pathTypeModels.First(x => x.Mode == Global.Settings.Modes.Hov3);
int modeExtra = Global.Settings.Modes.SchoolBus;
bool availableExtra = pathTypeExtra.Available && tour.IsSchoolBusMode() && tripModeAvailable[modeExtra];
double generalizedTimeLogsumExtra = pathTypeExtra.GeneralizedTimeLogsum;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(modeExtra, availableExtra, choice == modeExtra);
alternative.Choice = modeExtra;
alternative.AddNestedAlternative(_nestedAlternativeIds[modeExtra], _nestedAlternativeIndexes[modeExtra], THETA_PARAMETER);
if (availableExtra)
{
// case Global.Settings.Modes.SchoolBus:
alternative.AddUtilityTerm(2, generalizedTimeLogsumExtra * tour.TimeCoefficient);
alternative.AddUtilityTerm(18, 1);
alternative.AddUtilityTerm(100, schoolBusTourFlag);
alternative.AddUtilityTerm(102, (schoolBusTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (schoolBusTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (schoolBusTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (schoolBusTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (schoolBusTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (schoolBusTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(112, parkAndRideTourFlag);
alternative.AddUtilityTerm(113, transitTourFlag);
}
}
foreach (IPathTypeModel pathTypeModel in pathTypeModels)
{
int mode = pathTypeModel.Mode;
bool available = (pathTypeModel.Mode != Global.Settings.Modes.PaidRideShare || Global.Configuration.PaidRideShareModeIsAvailable) &&
pathTypeModel.Available && tripModeAvailable[mode];
double generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
alternative.AddNestedAlternative(_nestedAlternativeIds[mode], _nestedAlternativeIndexes[mode], THETA_PARAMETER);
if (!available)
{
continue;
}
double modeTimeCoefficient = (Global.Configuration.AV_IncludeAutoTypeChoice && household.OwnsAutomatedVehicles > 0 && mode >= Global.Settings.Modes.Sov && mode <= Global.Settings.Modes.Hov3) ?
tour.TimeCoefficient * (1.0 - Global.Configuration.AV_InVehicleTimeCoefficientDiscountFactor) :
(mode == Global.Settings.Modes.PaidRideShare && Global.Configuration.AV_PaidRideShareModeUsesAVs) ?
tour.TimeCoefficient * (1.0 - Global.Configuration.AV_InVehicleTimeCoefficientDiscountFactor) : tour.TimeCoefficient;
alternative.AddUtilityTerm(2, generalizedTimeLogsum * modeTimeCoefficient);
if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(20, 1);
alternative.AddUtilityTerm(22, carsLessThanDriversFlag);
alternative.AddUtilityTerm(100, transitTourFlag);
alternative.AddUtilityTerm(102, (transitTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (transitTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (transitTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (transitTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (transitTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (transitTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(136, tncTourFlag);
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(34, (nonworkingAdults + retiredAdults));
alternative.AddUtilityTerm(36, onePersonHouseholdFlag);
alternative.AddUtilityTerm(37, twoPersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(100, hov3TourFlag);
alternative.AddUtilityTerm(102, (hov3TourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (hov3TourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (hov3TourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (hov3TourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (hov3TourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (hov3TourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(114, parkAndRideTourFlag);
alternative.AddUtilityTerm(115, transitTourFlag);
alternative.AddUtilityTerm(116, schoolBusTourFlag);
alternative.AddUtilityTerm(135, tncTourFlag);
alternative.AddUtilityTerm(149, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(150, homeBasedSchoolTourFlag);
alternative.AddUtilityTerm(152, homeBasedEscortTourFlag);
alternative.AddUtilityTerm(153, homeBasedShoppingTourFlag);
alternative.AddUtilityTerm(154, homeBasedMealTourFlag);
alternative.AddUtilityTerm(155, homeBasedSocialTourFlag);
alternative.AddUtilityTerm(161, (destinationWorkEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(162, (originWorkEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(163, (originHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(164, (originHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(165, (originHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(166, (originHomeEscortFlag * eveningPeriodFlag));
alternative.AddUtilityTerm(167, (destinationHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(168, (destinationHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(169, (destinationHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(170, (destinationHomeEscortFlag * eveningPeriodFlag));
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(32, (childrenAge5Through15 * (1 - homeBasedEscortTourFlag)));
alternative.AddUtilityTerm(34, ((nonworkingAdults + retiredAdults) * (1 - homeBasedEscortTourFlag)));
alternative.AddUtilityTerm(38, onePersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(100, hov2TourFlag);
alternative.AddUtilityTerm(102, (hov2TourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (hov2TourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (hov2TourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (hov2TourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (hov2TourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (hov2TourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(118, transitTourFlag);
alternative.AddUtilityTerm(119, schoolBusTourFlag);
alternative.AddUtilityTerm(120, hov3TourFlag);
alternative.AddUtilityTerm(135, tncTourFlag);
alternative.AddUtilityTerm(149, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(150, homeBasedSchoolTourFlag);
alternative.AddUtilityTerm(152, homeBasedEscortTourFlag);
alternative.AddUtilityTerm(153, homeBasedShoppingTourFlag);
alternative.AddUtilityTerm(154, homeBasedMealTourFlag);
alternative.AddUtilityTerm(155, homeBasedSocialTourFlag);
alternative.AddUtilityTerm(161, (destinationWorkEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(162, (originWorkEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(163, (originHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(164, (originHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(165, (originHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(166, (originHomeEscortFlag * eveningPeriodFlag));
alternative.AddUtilityTerm(167, (destinationHomeEscortFlag * amPeriodFlag));
alternative.AddUtilityTerm(168, (destinationHomeEscortFlag * middayPeriodFlag));
alternative.AddUtilityTerm(169, (destinationHomeEscortFlag * pmPeriodFlag));
alternative.AddUtilityTerm(170, (destinationHomeEscortFlag * eveningPeriodFlag));
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient));
alternative.AddUtilityTerm(52, carsLessThanDriversFlag);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(55, income25To45KFlag);
alternative.AddUtilityTerm(59, drivingAgeStudentFlag);
alternative.AddUtilityTerm(100, sovTourFlag);
alternative.AddUtilityTerm(102, (sovTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (sovTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (sovTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (sovTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (sovTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (sovTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(121, parkAndRideTourFlag);
alternative.AddUtilityTerm(122, transitTourFlag);
alternative.AddUtilityTerm(124, hov3TourFlag);
alternative.AddUtilityTerm(125, hov2TourFlag);
alternative.AddUtilityTerm(134, tncTourFlag);
}
else if (mode == Global.Settings.Modes.Bike)
{
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(62, ageLessThan35Flag);
alternative.AddUtilityTerm(65, originIntersectionDensity);
alternative.AddUtilityTerm(100, bikeTourFlag);
alternative.AddUtilityTerm(102, (bikeTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (bikeTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (bikeTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (bikeTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (bikeTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (bikeTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(127, transitTourFlag);
alternative.AddUtilityTerm(128, schoolBusTourFlag);
alternative.AddUtilityTerm(130, hov2TourFlag);
alternative.AddUtilityTerm(131, sovTourFlag);
alternative.AddUtilityTerm(133, tncTourFlag);
alternative.AddUtilityTerm(147, notHomeBasedTourFlag);
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(72, ageLessThan35Flag);
alternative.AddUtilityTerm(75, originIntersectionDensity);
alternative.AddUtilityTerm(78, originMixedDensity); // origin and destination mixed use measures - geometric avg. - half mile from cell, in 1000s
alternative.AddUtilityTerm(100, walkTourFlag);
alternative.AddUtilityTerm(102, (walkTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (walkTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (walkTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (walkTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (walkTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (walkTourFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(141, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(142, homeBasedSchoolTourFlag);
}
else if (mode == Global.Settings.Modes.PaidRideShare)
{
alternative.AddUtilityTerm(100, tncTourFlag);
alternative.AddUtilityTerm(102, (tncTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (tncTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (tncTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (tncTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (tncTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (tncTourFlag * lastTripOnSecondHalfFlag));
if (Global.Configuration.PaidRideshare_UseEstimatedInsteadOfAssertedCoefficients)
{
alternative.AddUtilityTerm(80, 1.0);
alternative.AddUtilityTerm(81, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(82, carsLessThanDriversFlag);
alternative.AddUtilityTerm(83, tour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(83, tour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(84, (tour.Person.Age >= 65).ToFlag());
alternative.AddUtilityTerm(85, originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2);
alternative.AddUtilityTerm(86, destinationParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2);
alternative.AddUtilityTerm(87, income0To25KFlag);
alternative.AddUtilityTerm(88, incomeOver100Flag);
alternative.AddUtilityTerm(89, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(90, homeBasedSchoolTourFlag);
alternative.AddUtilityTerm(91, homeBasedEscortTourFlag);
alternative.AddUtilityTerm(92, homeBasedShoppingTourFlag);
}
else
{
double modeConstant = Global.Configuration.AV_PaidRideShareModeUsesAVs
? Global.Configuration.AV_PaidRideShare_ModeConstant
+ Global.Configuration.AV_PaidRideShare_DensityCoefficient * Math.Min(originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2,
(Global.Configuration.PaidRideShare_DensityMeasureCapValue > 0) ? Global.Configuration.PaidRideShare_DensityMeasureCapValue : 6000)
+ Global.Configuration.AV_PaidRideShare_AVOwnerCoefficient * (household.OwnsAutomatedVehicles > 0).ToFlag()
: Global.Configuration.PaidRideShare_ModeConstant
+ Global.Configuration.PaidRideShare_DensityCoefficient * Math.Min(originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2,
(Global.Configuration.PaidRideShare_DensityMeasureCapValue > 0) ? Global.Configuration.PaidRideShare_DensityMeasureCapValue : 6000);
alternative.AddUtilityTerm(90, modeConstant);
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_Age26to35Coefficient * tour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_Age18to25Coefficient * tour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_AgeOver65Coefficient * (tour.Person.Age >= 65).ToFlag());
}
}
}
}
public DestinationSampler(ChoiceProbabilityCalculator choiceProbabilityCalculator, int segmentZonesIndex, int sampleSize, IParcel chosenParcel, ITourWrapper tour, ITripWrapper trip) : this(choiceProbabilityCalculator, segmentZonesIndex, sampleSize, chosenParcel)
{
_tourOriginParcel = tour.OriginParcel;
_tourOriginSegmentZone = _segmentZones[_tourOriginParcel.ZoneId];
_tripOriginParcel = trip.OriginParcel;
_tripOriginSegmentZone = _segmentZones[_tripOriginParcel.ZoneId];
if (_tourOriginParcel == null || _tripOriginSegmentZone == null)
{
}
if (_choiceProbabilityCalculator.ModelIsInEstimationMode && chosenParcel == null)
{
throw new ChosenParcelNotSetInEstimationModeException();
}
}
//PCA from JLB
public DestinationSampler(ChoiceProbabilityCalculator choiceProbabilityCalculator, int segment, int sampleSize, ITourWrapper tour, ITripWrapper trip, IParcel chosenParcel)
{
_choiceProbabilityCalculator = choiceProbabilityCalculator;
_segmentZones = Global.SegmentZones[segment];
_sampleSize = sampleSize;
_tourOriginParcel = tour.OriginParcel;
_tourOriginSegmentZone = _segmentZones[_tourOriginParcel.ZoneId];
_tripOriginParcel = trip.OriginParcel;
_tripOriginSegmentZone = _segmentZones[_tripOriginParcel.ZoneId];
if (_tourOriginParcel == null || _tripOriginSegmentZone == null)
{
}
if (chosenParcel != null)
{
_chosenParcel = chosenParcel;
_chosenSegmentZone = _segmentZones[chosenParcel.ZoneId];
}
if (_choiceProbabilityCalculator.ModelIsInEstimationMode && chosenParcel == null)
{
throw new ChosenParcelNotSetInEstimationModeException();
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, IHouseholdDayWrapper householdDay, ITripWrapper trip, HTripTime choice = null)
{
var person = trip.Person;
var personDay = trip.PersonDay;
var tour = trip.Tour;
// person inputs
var partTimeWorkerFlag = person.IsPartTimeWorker.ToFlag();
var nonworkingAdultFlag = person.IsNonworkingAdult.ToFlag();
var universityStudentFlag = person.IsUniversityStudent.ToFlag();
var retiredAdultFlag = person.IsRetiredAdult.ToFlag();
var drivingAgeStudentFlag = person.IsDrivingAgeStudent.ToFlag();
var childAge5Through15Flag = person.IsChildAge5Through15.ToFlag();
var childUnder5Flag = person.IsChildUnder5.ToFlag();
// set tour inputs
var workTourFlag = tour.IsWorkPurpose().ToFlag();
var notWorkTourFlag = (!tour.IsWorkPurpose()).ToFlag();
var notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
var jointTourFlag = (tour.JointTourSequence > 0) ? 1 : 0;
var partialHalfTourFlag = (trip.IsHalfTourFromOrigin ? tour.PartialHalfTour1Sequence > 0 : tour.PartialHalfTour2Sequence > 0) ? 1 : 0;
var fullHalfTourFlag = (trip.IsHalfTourFromOrigin ? tour.FullHalfTour1Sequence > 0 : tour.FullHalfTour2Sequence > 0) ? 1 : 0;
// set trip inputs
var originChangeMode = trip.Sequence > 1 && trip.GetPreviousTrip().DestinationPurpose == Global.Settings.Purposes.ChangeMode;
var originSchoolFlag = trip.IsSchoolOriginPurpose().ToFlag();
var originEscortFlag = trip.IsEscortOriginPurpose().ToFlag();
var originShoppingFlag = trip.IsShoppingOriginPurpose().ToFlag();
var originPersonalBusinessFlag = trip.IsPersonalBusinessOriginPurpose().ToFlag();
var originMealFlag = trip.IsMealOriginPurpose().ToFlag();
var originSocialFlag = trip.IsSocialOriginPurpose().ToFlag();
var sovOrHovTripFlag = trip.UsesSovOrHovModes().ToFlag();
var transitTripFlag = trip.IsTransitMode().ToFlag();
var halfTourFromOriginFlag = trip.IsHalfTourFromOrigin.ToFlag();
var halfTourFromDestinationFlag = (!trip.IsHalfTourFromOrigin).ToFlag();
// set remaining inputs
// set remaining inputs
TimeWindow timeWindow = new TimeWindow();
if (tour.JointTourSequence > 0)
{
foreach (PersonDayWrapper pDay in householdDay.PersonDays)
{
var tInJoint = pDay.Tours.Find(t => t.JointTourSequence == tour.JointTourSequence);
if (!(tInJoint == null))
{
// set jointTour time window
timeWindow.IncorporateAnotherTimeWindow(tInJoint.PersonDay.TimeWindow);
}
}
}
else if (trip.Direction == Global.Settings.TourDirections.OriginToDestination && tour.FullHalfTour1Sequence > 0)
{
foreach (PersonDayWrapper pDay in householdDay.PersonDays)
{
var tInJoint = pDay.Tours.Find(t => t.FullHalfTour1Sequence == tour.FullHalfTour1Sequence);
if (!(tInJoint == null))
{
// set jointTour time window
timeWindow.IncorporateAnotherTimeWindow(tInJoint.PersonDay.TimeWindow);
}
}
}
else if (trip.Direction == Global.Settings.TourDirections.DestinationToOrigin && tour.FullHalfTour2Sequence > 0)
{
foreach (PersonDayWrapper pDay in householdDay.PersonDays)
{
var tInJoint = pDay.Tours.Find(t => t.FullHalfTour2Sequence == tour.FullHalfTour2Sequence);
if (!(tInJoint == null))
{
// set jointTour time window
timeWindow.IncorporateAnotherTimeWindow(tInJoint.PersonDay.TimeWindow);
}
}
}
else if (tour.ParentTour == null)
{
timeWindow.IncorporateAnotherTimeWindow(personDay.TimeWindow);
}
else
{
timeWindow.IncorporateAnotherTimeWindow(tour.ParentTour.TimeWindow);
}
//set the availability and impedances for the periods
HTripTime.SetTimeImpedances(trip);
var remainingToursCount = personDay.HomeBasedTours - personDay.GetTotalSimulatedTours();
var tripRemainingInHalfTour = (trip.DestinationParcel != null && trip.DestinationParcel != tour.OriginParcel).ToFlag(); // we don't know exact #
var previousArrivalTime = trip.IsHalfTourFromOrigin
? (trip.Sequence == 1 ? tour.DestinationDepartureTime : trip.GetPreviousTrip().ArrivalTime)
: (trip.Sequence == 1 ? tour.DestinationArrivalTime : trip.GetPreviousTrip().ArrivalTime);
var previousArrivalPeriod = new HTripTime(previousArrivalTime).DeparturePeriod;
foreach (var time in HTripTime.Times[ParallelUtility.threadLocalAssignedIndex.Value])
{
var period = time.DeparturePeriod;
var departurePeriodFraction = timeWindow.TotalAvailableMinutes(period.Start, period.End) / (period.End - period.Start + 1D);
var departureShift = period.Middle / 60.0;
var durationShift = Math.Abs(period.Middle - previousArrivalPeriod.Middle) / 60.0;
var available = time.Available && departurePeriodFraction > 0;
var alternative = choiceProbabilityCalculator.GetAlternative(time.Index, available, choice != null && choice.Equals(time));
if (!alternative.Available)
{
continue;
}
alternative.Choice = time;
var indicatedTravelTime = (int)time.ModeLOS.PathTime;
var indicatedArrivalTime = trip.IsHalfTourFromOrigin
? Math.Max(1, period.Middle - indicatedTravelTime)
: Math.Min(1440, period.Middle + indicatedTravelTime);
var totalWindowRemaining = trip.IsHalfTourFromOrigin
? timeWindow.TotalAvailableMinutesBefore(indicatedArrivalTime) + timeWindow.TotalAvailableMinutesAfter(previousArrivalTime)
: timeWindow.TotalAvailableMinutesAfter(indicatedArrivalTime) + timeWindow.TotalAvailableMinutesBefore(previousArrivalTime);
var maxWindowRemaining = trip.IsHalfTourFromOrigin
? timeWindow.MaxAvailableMinutesBefore(indicatedArrivalTime) + timeWindow.MaxAvailableMinutesAfter(previousArrivalTime)
: timeWindow.MaxAvailableMinutesAfter(indicatedArrivalTime) + timeWindow.MaxAvailableMinutesBefore(previousArrivalTime);
if (trip.IsHalfTourFromOrigin)
{
// outbound "departure" (arrival) period constants
alternative.AddUtilityTerm(11, time.DeparturePeriod.Middle.IsBetween(Global.Settings.Times.ThreeAM, Global.Settings.Times.SixAM).ToFlag());
alternative.AddUtilityTerm(12, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SixAM, Global.Settings.Times.SevenAM).ToFlag());
alternative.AddUtilityTerm(13, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenAM, Global.Settings.Times.EightAM).ToFlag());
alternative.AddUtilityTerm(14, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.EightAM, Global.Settings.Times.NineAM).ToFlag());
alternative.AddUtilityTerm(15, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.NineAM, Global.Settings.Times.TenAM).ToFlag());
alternative.AddUtilityTerm(16, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.OnePM).ToFlag());
alternative.AddUtilityTerm(17, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.OnePM, Global.Settings.Times.FourPM).ToFlag());
alternative.AddUtilityTerm(18, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FourPM, Global.Settings.Times.SevenPM).ToFlag());
alternative.AddUtilityTerm(19, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenPM, Global.Settings.Times.TenPM).ToFlag());
alternative.AddUtilityTerm(20, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenPM, Global.Settings.Times.MinutesInADay).ToFlag());
}
else
{
// return departure period constants
alternative.AddUtilityTerm(21, time.DeparturePeriod.Middle.IsBetween(Global.Settings.Times.ThreeAM, Global.Settings.Times.SevenAM).ToFlag());
alternative.AddUtilityTerm(22, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenAM, Global.Settings.Times.TenAM).ToFlag());
alternative.AddUtilityTerm(23, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.OnePM).ToFlag());
alternative.AddUtilityTerm(24, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.OnePM, Global.Settings.Times.ThreePM).ToFlag());
alternative.AddUtilityTerm(124, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.ThreePM, Global.Settings.Times.FourPM).ToFlag());
alternative.AddUtilityTerm(25, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FourPM, Global.Settings.Times.FivePM).ToFlag());
alternative.AddUtilityTerm(26, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FivePM, Global.Settings.Times.SixPM).ToFlag());
alternative.AddUtilityTerm(27, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SixPM, Global.Settings.Times.SevenPM).ToFlag());
alternative.AddUtilityTerm(28, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenPM, Global.Settings.Times.NinePM).ToFlag());
alternative.AddUtilityTerm(29, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.NinePM, Global.Settings.Times.Midnight).ToFlag());
alternative.AddUtilityTerm(30, time.DeparturePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.Midnight, Global.Settings.Times.MinutesInADay).ToFlag());
}
alternative.AddUtilityTerm(31, durationShift.IsRightExclusiveBetween(Global.Settings.Times.ZeroHours, Global.Settings.Times.OneHour).ToFlag()); // 0 - 1
alternative.AddUtilityTerm(32, durationShift.IsRightExclusiveBetween(Global.Settings.Times.OneHour, Global.Settings.Times.TwoHours).ToFlag()); // 1 - 2
alternative.AddUtilityTerm(33, durationShift.IsRightExclusiveBetween(Global.Settings.Times.TwoHours, Global.Settings.Times.ThreeHours).ToFlag()); // 2 - 3
alternative.AddUtilityTerm(34, durationShift.IsRightExclusiveBetween(Global.Settings.Times.ThreeHours, Global.Settings.Times.FiveHours).ToFlag()); // 3 - 5
alternative.AddUtilityTerm(35, durationShift.IsRightExclusiveBetween(Global.Settings.Times.FiveHours, Global.Settings.Times.SevenHours).ToFlag()); // 5 - 7
alternative.AddUtilityTerm(36, durationShift.IsRightExclusiveBetween(Global.Settings.Times.SevenHours, Global.Settings.Times.NineHours).ToFlag()); // 7 - 9
alternative.AddUtilityTerm(37, durationShift.IsRightExclusiveBetween(Global.Settings.Times.NineHours, Global.Settings.Times.TwelveHours).ToFlag()); // 9 - 12
alternative.AddUtilityTerm(38, durationShift.IsRightExclusiveBetween(Global.Settings.Times.TwelveHours, Global.Settings.Times.FourteenHours).ToFlag()); // 12 - 14
alternative.AddUtilityTerm(39, durationShift.IsRightExclusiveBetween(Global.Settings.Times.FourteenHours, Global.Settings.Times.EighteenHours).ToFlag()); // 14 - 18
alternative.AddUtilityTerm(40, (durationShift >= Global.Settings.Times.EighteenHours).ToFlag()); // 18 - 24
alternative.AddUtilityTerm(41, partTimeWorkerFlag * departureShift);
alternative.AddUtilityTerm(43, nonworkingAdultFlag * departureShift);
alternative.AddUtilityTerm(45, universityStudentFlag * departureShift);
alternative.AddUtilityTerm(47, retiredAdultFlag * departureShift);
alternative.AddUtilityTerm(49, drivingAgeStudentFlag * departureShift);
alternative.AddUtilityTerm(51, childAge5Through15Flag * departureShift);
alternative.AddUtilityTerm(53, childUnder5Flag * departureShift);
alternative.AddUtilityTerm(61, jointTourFlag * departureShift);
alternative.AddUtilityTerm(63, partialHalfTourFlag * departureShift);
alternative.AddUtilityTerm(65, fullHalfTourFlag * departureShift);
alternative.AddUtilityTerm(131, workTourFlag * halfTourFromOriginFlag * departureShift);
alternative.AddUtilityTerm(133, workTourFlag * halfTourFromDestinationFlag * departureShift);
alternative.AddUtilityTerm(135, notWorkTourFlag * halfTourFromDestinationFlag * departureShift);
alternative.AddUtilityTerm(137, notHomeBasedTourFlag * departureShift);
alternative.AddUtilityTerm(145, originEscortFlag * departureShift);
alternative.AddUtilityTerm(147, originShoppingFlag * departureShift);
alternative.AddUtilityTerm(149, originMealFlag * departureShift);
alternative.AddUtilityTerm(151, originSocialFlag * departureShift);
alternative.AddUtilityTerm(153, originPersonalBusinessFlag * departureShift);
alternative.AddUtilityTerm(155, originSchoolFlag * departureShift);
alternative.AddUtilityTerm(42, partTimeWorkerFlag * durationShift);
alternative.AddUtilityTerm(44, nonworkingAdultFlag * durationShift);
alternative.AddUtilityTerm(46, universityStudentFlag * durationShift);
alternative.AddUtilityTerm(48, retiredAdultFlag * durationShift);
alternative.AddUtilityTerm(50, drivingAgeStudentFlag * durationShift);
alternative.AddUtilityTerm(52, childAge5Through15Flag * durationShift);
alternative.AddUtilityTerm(54, childUnder5Flag * durationShift);
alternative.AddUtilityTerm(62, jointTourFlag * durationShift);
alternative.AddUtilityTerm(64, partialHalfTourFlag * durationShift);
alternative.AddUtilityTerm(66, fullHalfTourFlag * durationShift);
alternative.AddUtilityTerm(132, workTourFlag * halfTourFromOriginFlag * durationShift);
alternative.AddUtilityTerm(134, workTourFlag * halfTourFromDestinationFlag * durationShift);
alternative.AddUtilityTerm(136, notWorkTourFlag * halfTourFromDestinationFlag * durationShift);
alternative.AddUtilityTerm(138, notHomeBasedTourFlag * durationShift);
alternative.AddUtilityTerm(146, originEscortFlag * durationShift);
alternative.AddUtilityTerm(148, originShoppingFlag * durationShift);
alternative.AddUtilityTerm(150, originMealFlag * durationShift);
alternative.AddUtilityTerm(152, originSocialFlag * durationShift);
alternative.AddUtilityTerm(154, originPersonalBusinessFlag * durationShift);
alternative.AddUtilityTerm(156, originSchoolFlag * durationShift);
alternative.AddUtilityTerm(86, sovOrHovTripFlag * Math.Max(time.ModeLOS.GeneralizedTimeLogsum, 0) * tour.TimeCoefficient);
alternative.AddUtilityTerm(88, transitTripFlag * Math.Max(time.ModeLOS.GeneralizedTimeLogsum, 0) * tour.TimeCoefficient);
alternative.AddUtilityTerm(92, halfTourFromOriginFlag * Math.Log(departurePeriodFraction));
alternative.AddUtilityTerm(92, halfTourFromDestinationFlag * Math.Log(departurePeriodFraction));
alternative.AddUtilityTerm(99, tripRemainingInHalfTour / (Math.Max(1D, Math.Abs(trip.ArrivalTimeLimit - period.Middle))));
alternative.AddUtilityTerm(97, remainingToursCount / (Math.Max(1D, totalWindowRemaining)));
alternative.AddUtilityTerm(98, remainingToursCount / (Math.Max(1D, maxWindowRemaining)));
}
}
public int Run(ITripWrapper trip, IHouseholdDayWrapper householdDay)
{
return(Run(trip, householdDay, Global.Settings.Purposes.NoneOrHome));
}
public void Run(ITripWrapper trip)
{
if (trip == null)
{
throw new ArgumentNullException("trip");
}
trip.PersonDay.ResetRandom(40 * (2 * trip.Tour.Sequence - 1 + trip.Direction - 1) + 40 + trip.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
}
var choiceProbabilityCalculator = _helpers[ParallelUtility.GetBatchFromThreadId()].GetChoiceProbabilityCalculator(trip.Id);
var originParcel =
trip.IsHalfTourFromOrigin
? trip.DestinationParcel
: trip.OriginParcel;
// for skims - use actual travel direction, not simulation direction
var destinationParcel =
trip.IsHalfTourFromOrigin
? trip.OriginParcel
: trip.DestinationParcel;
// if first trip in half-tour, use tour destination time
var departureTime =
trip.IsHalfTourFromOrigin
? trip.Sequence == 1
? trip.Tour.DestinationArrivalTime
: trip.GetPreviousTrip().ArrivalTime
: trip.Sequence == 1
? trip.Tour.DestinationDepartureTime
: trip.GetPreviousTrip().ArrivalTime;
if (_helpers[ParallelUtility.GetBatchFromThreadId()].ModelIsInEstimationMode)
{
if (destinationParcel == null || originParcel == null || trip.Mode <= Global.Settings.Modes.None || trip.Mode == Global.Settings.Modes.ParkAndRide || trip.Mode == Global.Settings.Modes.Other)
{
return;
}
IEnumerable <dynamic> pathTypeModels =
PathTypeModelFactory.Model.RunAll(
trip.Household.RandomUtility,
originParcel,
destinationParcel,
departureTime,
0,
trip.Tour.DestinationPurpose,
trip.Tour.CostCoefficient,
trip.Tour.TimeCoefficient,
trip.Person.IsDrivingAge,
trip.Household.VehiclesAvailable,
trip.Person.GetTransitFareDiscountFraction(),
false);
// there is no path type model for school bus, use HOV3
var mode = trip.Mode == Global.Settings.Modes.SchoolBus ? Global.Settings.Modes.Hov3 : trip.Mode;
var pathTypeModel = pathTypeModels.First(x => x.Mode == mode);
if (!pathTypeModel.Available)
{
return;
}
RunModel(choiceProbabilityCalculator, trip, pathTypeModels, originParcel, destinationParcel, departureTime, trip.Mode);
choiceProbabilityCalculator.WriteObservation();
}
else
{
IEnumerable <dynamic> pathTypeModels =
PathTypeModelFactory.Model.RunAll(
trip.Household.RandomUtility,
originParcel,
destinationParcel,
departureTime,
0,
trip.Tour.DestinationPurpose,
trip.Tour.CostCoefficient,
trip.Tour.TimeCoefficient,
trip.Person.IsDrivingAge,
trip.Household.VehiclesAvailable,
trip.Person.GetTransitFareDiscountFraction(),
false);
RunModel(choiceProbabilityCalculator, trip, pathTypeModels, originParcel, destinationParcel, departureTime);
var chosenAlternative = choiceProbabilityCalculator.SimulateChoice(trip.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", trip.PersonDay.Id);
trip.Mode = Global.Settings.Modes.Hov3;
trip.PersonDay.IsValid = false;
return;
}
var choice = (int)chosenAlternative.Choice;
trip.Mode = choice;
if (choice == Global.Settings.Modes.SchoolBus)
{
trip.PathType = 0;
}
else
{
var chosenPathType = pathTypeModels.First(x => x.Mode == choice);
trip.PathType = chosenPathType.PathType;
// for transit trips, overwrite origin and destination zones with stop area ids
if (Global.StopAreaIsEnabled && choice == Global.Settings.Modes.Transit &&
Global.Configuration.WriteStopAreaIDsInsteadOfZonesForTransitTrips)
{
trip.OriginZoneKey = chosenPathType.PathOriginStopAreaKey;
trip.DestinationZoneKey = chosenPathType.PathDestinationStopAreaKey;
}
}
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITripWrapper trip, IHouseholdDayWrapper householdDay, int choice = Constants.DEFAULT_VALUE)
{
var household = trip.Household;
var person = trip.Person;
var personDay = trip.PersonDay;
var tour = trip.Tour;
var halfTour = trip.HalfTour;
var personDays = householdDay.PersonDays;
var isJointTour = tour.JointTourSequence > 0 ? 1 : 0;
var destinationParcel = tour.DestinationParcel;
var fullJointHalfTour = ((trip.Direction == Global.Settings.TourDirections.OriginToDestination && tour.FullHalfTour1Sequence > 0) ||
(trip.Direction == Global.Settings.TourDirections.DestinationToOrigin && tour.FullHalfTour2Sequence > 0)).ToFlag();
var partialJointHalfTour = ((trip.Direction == Global.Settings.TourDirections.OriginToDestination && tour.PartialHalfTour1Sequence > 0) ||
(trip.Direction == Global.Settings.TourDirections.DestinationToOrigin && tour.PartialHalfTour2Sequence > 0)).ToFlag();
var isJoint = (isJointTour + fullJointHalfTour + partialJointHalfTour > 0) ? 1 : 0;
var isIndividual = 1 - isJoint;
//destination parcel variables
var foodBuffer2 = 0.0;
var totEmpBuffer2 = 0.0;
var retailBuffer2 = 0.0;
if (destinationParcel != null)
{
foodBuffer2 = Math.Log(1 + destinationParcel.EmploymentFoodBuffer2);
totEmpBuffer2 = Math.Log(1 + destinationParcel.EmploymentTotalBuffer2);
retailBuffer2 = Math.Log(1 + destinationParcel.EmploymentRetailBuffer2);
}
var carOwnership = person.GetCarOwnershipSegment();
// household inputs
var onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
//var householdInc75KP = household.Has75KPlusIncome;
var votALSegment = tour.GetVotALSegment();
var transitAccessSegment = household.ResidenceParcel.TransitAccessSegment();
var totalAggregateLogsum = Global.AggregateLogsums[household.ResidenceParcel.ZoneId]
[Global.Settings.Purposes.HomeBasedComposite][carOwnership][votALSegment][transitAccessSegment];
var homeFoodBuffer2 = Math.Log(1 + household.ResidenceParcel.EmploymentFoodBuffer2);
var homeTotEmpBuffer2 = Math.Log(1 + household.ResidenceParcel.EmploymentTotalBuffer2);
var homeRetailBuffer2 = Math.Log(1 + household.ResidenceParcel.EmploymentRetailBuffer2);
// person-day inputs
var homeBasedTours = personDay.HomeBasedTours;
var simulatedToursFlag = personDay.SimulatedToursExist().ToFlag();
var simulatedWorkStops = personDay.SimulatedWorkStops;
var simulatedWorkStopsFlag = personDay.SimulatedWorkStopsExist().ToFlag();
var simulatedSchoolStops = personDay.SimulatedSchoolStops;
var simulatedEscortStops = personDay.SimulatedEscortStops;
var simulatedPersonalBusinessStops = personDay.SimulatedPersonalBusinessStops;
var simulatedShoppingStops = personDay.SimulatedShoppingStops;
var simulatedMealStops = personDay.SimulatedMealStops;
var simulatedSocialStops = personDay.SimulatedSocialStops;
var simulatedRecreationStops = personDay.SimulatedRecreationStops;
var simulatedMedicalStops = personDay.SimulatedMedicalStops;
// tour inputs
var hov2TourFlag = tour.IsHov2Mode().ToFlag();
var hov3TourFlag = tour.IsHov3Mode().ToFlag();
var transitTourFlag = tour.IsTransitMode().ToFlag();
var notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
var workTourFlag = tour.IsWorkPurpose().ToFlag();
var personalBusinessOrMedicalTourFlag = tour.IsPersonalBusinessOrMedicalPurpose().ToFlag();
var socialTourFlag = tour.IsSocialPurpose().ToFlag();
var socialOrRecreationTourFlag = tour.IsSocialOrRecreationPurpose().ToFlag();
var schoolTourFlag = tour.IsSchoolPurpose().ToFlag();
var escortTourFlag = tour.IsEscortPurpose().ToFlag();
var shoppingTourFlag = tour.IsShoppingPurpose().ToFlag();
var mealTourFlag = tour.IsMealPurpose().ToFlag();
// trip inputs
var oneSimulatedTripFlag = halfTour.OneSimulatedTripFlag;
var twoSimulatedTripsFlag = halfTour.TwoSimulatedTripsFlag;
var threeSimulatedTripsFlag = halfTour.ThreeSimulatedTripsFlag;
var fourSimulatedTripsFlag = halfTour.FourSimulatedTripsFlag;
var fivePlusSimulatedTripsFlag = halfTour.FivePlusSimulatedTripsFlag; //JLB changed 5 to 5Plus
var halfTourFromOriginFlag = trip.IsHalfTourFromOrigin.ToFlag();
var halfTourFromDestinationFlag = (!trip.IsHalfTourFromOrigin).ToFlag();
var beforeMandatoryDestinationFlag = trip.IsBeforeMandatoryDestination().ToFlag();
// remaining inputs, including joint tour variables
var time = trip.IsHalfTourFromOrigin ? tour.DestinationArrivalTime : tour.DestinationDepartureTime;
var from7AMto9AMFlag = (time >= Global.Settings.Times.SevenAM && time < Global.Settings.Times.NineAM).ToFlag();
var from9AMto11AMFlag = (time >= Global.Settings.Times.NineAM && time < Global.Settings.Times.ElevenAM).ToFlag();
var from11AMto1PMFlag = (time >= Global.Settings.Times.ElevenAM && time < Global.Settings.Times.OnePM).ToFlag();
var from1PMto3PMFlag = (time >= Global.Settings.Times.OnePM && time < Global.Settings.Times.ThreePM).ToFlag();
var from3PMto5PMFlag = (time >= Global.Settings.Times.ThreePM && time < Global.Settings.Times.FivePM).ToFlag();
var from7PMto9PMFlag = (time >= Global.Settings.Times.SevenPM && time < Global.Settings.Times.NinePM).ToFlag();
var from9PMto11PMFlag = (time >= Global.Settings.Times.NinePM && time < Global.Settings.Times.ElevenPM).ToFlag();
var from11PMto7AMFlag = (time >= Global.Settings.Times.ElevenPM).ToFlag();
var remainingToursCount = personDay.HomeBasedTours - personDay.GetTotalSimulatedTours();
// connectivity attributes
var c34Ratio = trip.OriginParcel.C34RatioBuffer1();
var adis = 0.0;
var logDist = 0.0;
var minute = DayPeriod.BigDayPeriods[DayPeriod.MIDDAY].Start;
//distance from origin to destination
if (tour.OriginParcel != null && tour.DestinationParcel != null)
{
if (trip.Direction == Global.Settings.TourDirections.OriginToDestination)
{
adis = ImpedanceRoster.GetValue("distance", Global.Settings.Modes.Sov, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, minute, tour.OriginParcel, tour.DestinationParcel).Variable;
}
else
{
adis = ImpedanceRoster.GetValue("distance", Global.Settings.Modes.Sov, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, minute, destinationParcel, tour.OriginParcel).Variable;
}
logDist = Math.Log(1 + adis);
}
var shortestTravelTime = 0D;
if (trip.Sequence > 1 && tour.OriginParcel != null && trip.OriginParcel != null)
{
if (trip.Direction == Global.Settings.TourDirections.OriginToDestination)
{
shortestTravelTime = ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov2, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, minute, tour.OriginParcel, trip.OriginParcel).Variable;
}
else
{
shortestTravelTime = ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov2, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, minute, trip.OriginParcel, tour.OriginParcel).Variable;
}
}
var windowDuration = Math.Max(1.0,
trip.IsHalfTourFromOrigin // first trip in half tour, use tour destination time
? trip.Sequence == 1
? tour.DestinationArrivalTime - tour.EarliestOriginDepartureTime - tour.IndicatedTravelTimeToDestination
: trip.GetPreviousTrip().ArrivalTime - tour.EarliestOriginDepartureTime - shortestTravelTime
: trip.Sequence == 1
? tour.LatestOriginArrivalTime - tour.DestinationDepartureTime - tour.IndicatedTravelTimeFromDestination
: tour.LatestOriginArrivalTime - trip.GetPreviousTrip().ArrivalTime - shortestTravelTime);
var duration = windowDuration / 60D;
IEnumerable <PersonDayWrapper> orderedPersonDays = householdDay.PersonDays.OrderBy(p => p.GetJointTourParticipationPriority()).ToList().Cast <PersonDayWrapper>();
int numChildrenOnJointTour = 0;
int numAdultsOnJointTour = 0;
int totHHToursJT = 0;
//int totHHStopsJT=0;
//MB use calculated variables at tour level rather than redoing time window
//TimeWindow timeWindow = new TimeWindow();
if (tour.JointTourSequence > 0)
{
foreach (PersonDayWrapper pDay in householdDay.PersonDays)
{
var tInJoint = pDay.Tours.Find(t => t.JointTourSequence == tour.JointTourSequence);
if (!(tInJoint == null))
{
// set jointTour time window
//timeWindow.IncorporateAnotherTimeWindow(tInJoint.PersonDay.TimeWindow);
totHHToursJT = personDay.HomeBasedTours + totHHToursJT;
if (pDay.Person.Age < 18)
{
numChildrenOnJointTour++;
}
if (pDay.Person.Age >= 18)
{
numAdultsOnJointTour++;
}
}
}
}
else if (trip.Direction == Global.Settings.TourDirections.OriginToDestination && tour.FullHalfTour1Sequence > 0)
{
foreach (PersonDayWrapper pDay in householdDay.PersonDays)
{
var tInJoint = pDay.Tours.Find(t => t.FullHalfTour1Sequence == tour.FullHalfTour1Sequence);
if (!(tInJoint == null))
{
// set jointTour time window
//timeWindow.IncorporateAnotherTimeWindow(tInJoint.PersonDay.TimeWindow);
totHHToursJT = personDay.HomeBasedTours + totHHToursJT;
if (pDay.Person.Age < 18)
{
numChildrenOnJointTour++;
}
if (pDay.Person.Age >= 18)
{
numAdultsOnJointTour++;
}
}
}
}
else if (trip.Direction == Global.Settings.TourDirections.DestinationToOrigin && tour.FullHalfTour2Sequence > 0)
{
foreach (PersonDayWrapper pDay in householdDay.PersonDays)
{
var tInJoint = pDay.Tours.Find(t => t.FullHalfTour2Sequence == tour.FullHalfTour2Sequence);
if (!(tInJoint == null))
{
// set jointTour time window
//timeWindow.IncorporateAnotherTimeWindow(tInJoint.PersonDay.TimeWindow);
totHHToursJT = personDay.HomeBasedTours + totHHToursJT;
if (pDay.Person.Age < 18)
{
numChildrenOnJointTour++;
}
if (pDay.Person.Age >= 18)
{
numAdultsOnJointTour++;
}
}
}
}
else if (tour.ParentTour == null)
{
//timeWindow.IncorporateAnotherTimeWindow(personDay.TimeWindow);
}
else
{
//timeWindow.IncorporateAnotherTimeWindow(tour.ParentTour.TimeWindow);
}
//timeWindow.SetBusyMinutes(Global.Settings.Times.EndOfRelevantWindow, Global.Settings.Times.MinutesInADay + 1);
// time window in minutes for yet unmodeled portion of halftour, only consider persons on this trip
//var availableWindow = timeWindow.AvailableWindow(destinationDepartureTime, Global.Settings.TimeDirections.Both);
//var duration = availableWindow/ 60D;
// 0 - NO MORE STOPS
var alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.NoneOrHome, true, choice == Global.Settings.Purposes.NoneOrHome);
alternative.Choice = Global.Settings.Purposes.NoneOrHome;
alternative.AddUtilityTerm(1, twoSimulatedTripsFlag * halfTourFromOriginFlag * isIndividual);
alternative.AddUtilityTerm(2, threeSimulatedTripsFlag * halfTourFromOriginFlag * isIndividual);
alternative.AddUtilityTerm(3, fourSimulatedTripsFlag * halfTourFromOriginFlag * isIndividual);
alternative.AddUtilityTerm(4, fivePlusSimulatedTripsFlag * halfTourFromOriginFlag * isIndividual);
alternative.AddUtilityTerm(5, twoSimulatedTripsFlag * halfTourFromDestinationFlag * isIndividual);
alternative.AddUtilityTerm(6, threeSimulatedTripsFlag * halfTourFromDestinationFlag * isIndividual);
alternative.AddUtilityTerm(7, fourSimulatedTripsFlag * halfTourFromDestinationFlag * isIndividual);
alternative.AddUtilityTerm(8, fivePlusSimulatedTripsFlag * halfTourFromDestinationFlag * isIndividual);
alternative.AddUtilityTerm(9, homeBasedTours * isIndividual);
alternative.AddUtilityTerm(10, homeBasedTours * isJointTour);
alternative.AddUtilityTerm(11, notHomeBasedTourFlag);
//alternative.AddUtilityTerm(12, beforeMandatoryDestinationFlag*isJointTour);
alternative.AddUtilityTerm(13, beforeMandatoryDestinationFlag);
alternative.AddUtilityTerm(14, numAdultsOnJointTour);
alternative.AddUtilityTerm(15, numChildrenOnJointTour);
alternative.AddUtilityTerm(16, totHHToursJT);
//alternative.AddUtilityTerm(17, totHHStopsJT);
alternative.AddUtilityTerm(22, (threeSimulatedTripsFlag + fourSimulatedTripsFlag + fivePlusSimulatedTripsFlag) * halfTourFromOriginFlag * isJoint);
alternative.AddUtilityTerm(26, threeSimulatedTripsFlag * halfTourFromDestinationFlag * isJoint);
alternative.AddUtilityTerm(27, fourSimulatedTripsFlag * halfTourFromDestinationFlag * isJoint);
alternative.AddUtilityTerm(28, fivePlusSimulatedTripsFlag * halfTourFromDestinationFlag * isJoint);
alternative.AddUtilityTerm(29, totalAggregateLogsum);
// 1 - WORK STOP
if ((personDay.WorkStops > 0 && tour.DestinationPurpose <= Global.Settings.Purposes.School) || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Work, true, choice == Global.Settings.Purposes.Work);
alternative.Choice = Global.Settings.Purposes.Work;
//alternative.AddUtilityTerm(32, isIndividualTour);
alternative.AddUtilityTerm(33, workTourFlag);
alternative.AddUtilityTerm(34, schoolTourFlag);
alternative.AddUtilityTerm(35, halfTourFromOriginFlag);
alternative.AddUtilityTerm(36, Math.Log(1 + simulatedWorkStops));
alternative.AddUtilityTerm(37, simulatedWorkStopsFlag);
alternative.AddUtilityTerm(38, (5 - simulatedWorkStops));
alternative.AddUtilityTerm(39, duration);
alternative.AddUtilityTerm(40, from9AMto11AMFlag + from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
//alternative.AddUtilityTerm(42, logDist);
//alternative.AddUtilityTerm(43, transitTourFlag);
//alternative.AddUtilityTerm(44, (person.IsPartTimeWorker).ToFlag());
alternative.AddUtilityTerm(46, totalAggregateLogsum);
//alternative.AddUtilityTerm(47,totEmpBuffer2);
alternative.AddUtilityTerm(48, hov2TourFlag + hov3TourFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Work, false, choice == Global.Settings.Purposes.Work);
}
// 2 - SCHOOL STOP
if ((personDay.SchoolStops > 0 && tour.DestinationPurpose <= Global.Settings.Purposes.School) || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.School, true, choice == Global.Settings.Purposes.School);
alternative.Choice = Global.Settings.Purposes.School;
alternative.AddUtilityTerm(51, workTourFlag);
alternative.AddUtilityTerm(52, schoolTourFlag);
//alternative.AddUtilityTerm(53, halfTourFromOriginFlag);
alternative.AddUtilityTerm(54, Math.Log(1 + simulatedSchoolStops));
alternative.AddUtilityTerm(55, (3 - simulatedSchoolStops));
//alternative.AddUtilityTerm(55, remainingToursCount);
alternative.AddUtilityTerm(56, duration);
alternative.AddUtilityTerm(57, from7AMto9AMFlag + from7PMto9PMFlag + from9PMto11PMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(58, oneSimulatedTripFlag);
//alternative.AddUtilityTerm(59, logDist);
alternative.AddUtilityTerm(61, fullJointHalfTour * numChildrenOnJointTour);
alternative.AddUtilityTerm(65, (person.Age < 12).ToFlag());
alternative.AddUtilityTerm(66, totalAggregateLogsum);
//alternative.AddUtilityTerm(66, (person.IsUniversityStudent).ToFlag());
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.School, false, choice == Global.Settings.Purposes.School);
}
// 3 - ESCORT STOP
if ((personDay.EscortStops > 0 && tour.DestinationPurpose <= Global.Settings.Purposes.Escort) || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Escort, true, choice == Global.Settings.Purposes.Escort);
alternative.Choice = Global.Settings.Purposes.Escort;
alternative.AddUtilityTerm(71, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(72, isJointTour);
alternative.AddUtilityTerm(74, escortTourFlag);
//alternative.AddUtilityTerm(75, socialOrRecreationTourFlag);
//alternative.AddUtilityTerm(76, remainingToursCount);
alternative.AddUtilityTerm(77, duration);
alternative.AddUtilityTerm(78, from7AMto9AMFlag);
//alternative.AddUtilityTerm(79, from9AMto11AMFlag + from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(81, hov2TourFlag);
alternative.AddUtilityTerm(82, hov3TourFlag);
alternative.AddUtilityTerm(83, Math.Log(1 + simulatedEscortStops * isJointTour));
alternative.AddUtilityTerm(84, Math.Log(1 + simulatedEscortStops * isIndividual));
//alternative.AddUtilityTerm(85, (3 - simulatedEscortStops));
alternative.AddUtilityTerm(85, totalAggregateLogsum);
alternative.AddUtilityTerm(86, fullJointHalfTour);
//alternative.AddUtilityTerm(88, enrollmentK8Buffer2);
alternative.AddUtilityTerm(89, numChildrenOnJointTour);
alternative.AddUtilityTerm(90, halfTourFromOriginFlag);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Escort, false, choice == Global.Settings.Purposes.Escort);
}
// 4 - PERSONAL BUSINESS STOP
if (personDay.PersonalBusinessStops > 0 || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.PersonalBusiness, true, choice == Global.Settings.Purposes.PersonalBusiness);
alternative.Choice = Global.Settings.Purposes.PersonalBusiness;
alternative.AddUtilityTerm(91, (workTourFlag + schoolTourFlag));
alternative.AddUtilityTerm(92, isJointTour);
alternative.AddUtilityTerm(93, escortTourFlag);
alternative.AddUtilityTerm(94, personalBusinessOrMedicalTourFlag * isIndividual);
alternative.AddUtilityTerm(95, shoppingTourFlag);
alternative.AddUtilityTerm(96, mealTourFlag);
alternative.AddUtilityTerm(97, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(98, halfTourFromOriginFlag);
alternative.AddUtilityTerm(99, Math.Log(1 + simulatedPersonalBusinessStops * isIndividual));
alternative.AddUtilityTerm(100, Math.Log(1 + simulatedPersonalBusinessStops * isJointTour));
alternative.AddUtilityTerm(104, (3 - simulatedPersonalBusinessStops));
alternative.AddUtilityTerm(101, duration);
alternative.AddUtilityTerm(102, (from7AMto9AMFlag + from7PMto9PMFlag + from9PMto11PMFlag + from11PMto7AMFlag));
alternative.AddUtilityTerm(103, from9AMto11AMFlag + from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(105, hov2TourFlag);
alternative.AddUtilityTerm(106, hov3TourFlag);
alternative.AddUtilityTerm(109, fullJointHalfTour);
alternative.AddUtilityTerm(110, totEmpBuffer2);
alternative.AddUtilityTerm(111, totalAggregateLogsum);
alternative.AddUtilityTerm(112, personalBusinessOrMedicalTourFlag * isJointTour);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.PersonalBusiness, false, choice == Global.Settings.Purposes.PersonalBusiness);
}
// 5 - SHOPPING STOP
if (personDay.ShoppingStops > 0 || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Shopping, true, choice == Global.Settings.Purposes.Shopping);
alternative.Choice = Global.Settings.Purposes.Shopping;
alternative.AddUtilityTerm(121, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(122, isJointTour);
alternative.AddUtilityTerm(123, escortTourFlag);
alternative.AddUtilityTerm(124, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(125, shoppingTourFlag * isIndividual);
alternative.AddUtilityTerm(126, mealTourFlag);
alternative.AddUtilityTerm(127, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(128, halfTourFromOriginFlag);
alternative.AddUtilityTerm(129, Math.Log(1 + simulatedShoppingStops * isIndividual));
alternative.AddUtilityTerm(130, Math.Log(1 + simulatedShoppingStops * isJointTour));
alternative.AddUtilityTerm(131, duration);
alternative.AddUtilityTerm(132, from7AMto9AMFlag + from9PMto11PMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(133, (from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag));
//alternative.AddUtilityTerm(134, adultFemaleOnJointTour);
alternative.AddUtilityTerm(135, hov2TourFlag);
alternative.AddUtilityTerm(136, hov3TourFlag);
alternative.AddUtilityTerm(137, Math.Log(1 + adis));
alternative.AddUtilityTerm(138, shoppingTourFlag * isJointTour);
alternative.AddUtilityTerm(140, totalAggregateLogsum);
//alternative.AddUtilityTerm(141, retailBuffer2);
//alternative.AddUtilityTerm(142, numChildrenOnJointTour);
//alternative.AddUtilityTerm(143, (household.Has100KPlusIncome).ToFlag());
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Shopping, false, choice == Global.Settings.Purposes.Shopping);
}
// 6 - MEAL STOP
if (personDay.MealStops > 0 || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Meal, true, choice == Global.Settings.Purposes.Meal);
alternative.Choice = Global.Settings.Purposes.Meal;
alternative.AddUtilityTerm(151, workTourFlag);
//alternative.AddUtilityTerm(152, isJointTour);
alternative.AddUtilityTerm(153, schoolTourFlag);
alternative.AddUtilityTerm(154, escortTourFlag);
alternative.AddUtilityTerm(155, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(156, shoppingTourFlag);
alternative.AddUtilityTerm(157, mealTourFlag);
alternative.AddUtilityTerm(158, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(159, halfTourFromOriginFlag);
alternative.AddUtilityTerm(160, Math.Log(1 + simulatedMealStops * isIndividual));
alternative.AddUtilityTerm(161, Math.Log(1 + simulatedMealStops * isJointTour));
alternative.AddUtilityTerm(162, duration);
alternative.AddUtilityTerm(164, from11AMto1PMFlag + from1PMto3PMFlag);
alternative.AddUtilityTerm(166, onePersonHouseholdFlag);
alternative.AddUtilityTerm(167, hov2TourFlag);
alternative.AddUtilityTerm(168, hov3TourFlag);
alternative.AddUtilityTerm(170, numChildrenOnJointTour);
alternative.AddUtilityTerm(171, oneSimulatedTripFlag);
alternative.AddUtilityTerm(172, Math.Log(1 + adis));
alternative.AddUtilityTerm(174, fullJointHalfTour);
alternative.AddUtilityTerm(175, foodBuffer2);
//alternative.AddUtilityTerm(176, homeFoodBuffer2);
//alternative.AddUtilityTerm(177, (household.Has100KPlusIncome).ToFlag());
alternative.AddUtilityTerm(178, numAdultsOnJointTour);
alternative.AddUtilityTerm(179, totalAggregateLogsum);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Meal, false, choice == Global.Settings.Purposes.Meal);
}
// 7 - SOCIAL STOP
if (personDay.SocialStops > 0 || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Social, true, choice == Global.Settings.Purposes.Social);
alternative.Choice = Global.Settings.Purposes.Social;
alternative.AddUtilityTerm(181, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(182, isJointTour);
alternative.AddUtilityTerm(183, escortTourFlag);
alternative.AddUtilityTerm(184, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(185, shoppingTourFlag);
alternative.AddUtilityTerm(186, mealTourFlag);
alternative.AddUtilityTerm(187, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(188, halfTourFromOriginFlag);
alternative.AddUtilityTerm(189, Math.Log(1 + simulatedSocialStops * isIndividual));
alternative.AddUtilityTerm(197, Math.Log(1 + simulatedSocialStops * isJointTour));
alternative.AddUtilityTerm(190, remainingToursCount);
alternative.AddUtilityTerm(191, duration);
alternative.AddUtilityTerm(192, from7AMto9AMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(194, hov2TourFlag);
alternative.AddUtilityTerm(195, hov3TourFlag);
alternative.AddUtilityTerm(196, logDist);
alternative.AddUtilityTerm(198, totalAggregateLogsum);
alternative.AddUtilityTerm(200, numAdultsOnJointTour);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Social, false, choice == Global.Settings.Purposes.Social);
}
// 8 - RECREATION STOP
if (personDay.RecreationStops > 0 || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Recreation, true, choice == Global.Settings.Purposes.Recreation);
alternative.Choice = Global.Settings.Purposes.Recreation;
alternative.AddUtilityTerm(211, workTourFlag + schoolTourFlag);
alternative.AddUtilityTerm(212, isJointTour);
alternative.AddUtilityTerm(213, escortTourFlag);
alternative.AddUtilityTerm(214, personalBusinessOrMedicalTourFlag);
alternative.AddUtilityTerm(215, shoppingTourFlag);
alternative.AddUtilityTerm(216, mealTourFlag);
alternative.AddUtilityTerm(217, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(218, halfTourFromOriginFlag);
alternative.AddUtilityTerm(219, Math.Log(1 + simulatedRecreationStops * isIndividual));
alternative.AddUtilityTerm(229, Math.Log(1 + simulatedRecreationStops * isJointTour));
alternative.AddUtilityTerm(220, remainingToursCount);
alternative.AddUtilityTerm(221, duration);
alternative.AddUtilityTerm(222, from7AMto9AMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(223, from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
alternative.AddUtilityTerm(225, hov3TourFlag);
alternative.AddUtilityTerm(226, numChildrenOnJointTour);
alternative.AddUtilityTerm(227, numAdultsOnJointTour);
//alternative.AddUtilityTerm(228, fullJointHalfTour);
//alternative.AddUtilityTerm(229, openSpaceBuffer2);
alternative.AddUtilityTerm(228, totalAggregateLogsum);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Recreation, false, choice == Global.Settings.Purposes.Recreation);
}
// 9 - MEDICAL STOP
if (personDay.MedicalStops > 0 || isJoint > 0)
{
alternative = choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Medical, true, choice == Global.Settings.Purposes.Medical);
alternative.Choice = Global.Settings.Purposes.Medical;
alternative.AddUtilityTerm(231, workTourFlag + schoolTourFlag + escortTourFlag);
alternative.AddUtilityTerm(232, isJointTour);
alternative.AddUtilityTerm(233, personalBusinessOrMedicalTourFlag * isIndividual);
alternative.AddUtilityTerm(234, personalBusinessOrMedicalTourFlag * isJointTour);
alternative.AddUtilityTerm(235, shoppingTourFlag);
alternative.AddUtilityTerm(236, mealTourFlag);
alternative.AddUtilityTerm(237, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(238, halfTourFromOriginFlag);
alternative.AddUtilityTerm(239, Math.Log(1 + simulatedMedicalStops * isJointTour));
alternative.AddUtilityTerm(240, Math.Log(1 + simulatedMedicalStops * isIndividual));
//alternative.AddUtilityTerm(240, fullJointHalfTour);
alternative.AddUtilityTerm(241, duration);
alternative.AddUtilityTerm(242, from7AMto9AMFlag + from11PMto7AMFlag);
alternative.AddUtilityTerm(243, from11AMto1PMFlag + from1PMto3PMFlag + from3PMto5PMFlag);
//alternative.AddUtilityTerm(248, numChildrenOnJointTour);
//alternative.AddUtilityTerm(249, adultFemaleOnJointTour);
alternative.AddUtilityTerm(244, totalAggregateLogsum);
}
else
{
choiceProbabilityCalculator.GetAlternative(Global.Settings.Purposes.Medical, false, choice == Global.Settings.Purposes.Medical);
}
}
