private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, IHouseholdDayWrapper householdDay, int sampleSize, IParcelWrapper choice = null)
{
timesStartedRunModel++;
//Console.WriteLine("Started {0} Finished {1}", timesStarted,timesStartedRunModel);
var household = tour.Household;
var person = tour.Person;
var personDay = tour.PersonDay;
var householdTotals = household.HouseholdTotals;
var originParcel = tour.OriginParcel;
int childUnder16Flag = (person.Age < 16).ToFlag();
int noCarsInHouseholdFlag = (household.VehiclesAvailable == 0).ToFlag();
int fewerCarsThanDriversFlag = household.GetFlagForCarsLessThanDrivers(household.VehiclesAvailable);
int lowIncomeFlag = (household.Income >= 0 && household.Income <= 40000).ToFlag();
int highIncomeFlag = (household.Income >= 80000).ToFlag();
int missingIncomeFlag = (household.Income < 0).ToFlag();
int originTransitBand1 = (originParcel.GetDistanceToTransit() >= 0 && originParcel.GetDistanceToTransit() <= 0.25).ToFlag();
int originTransitBand3 = (originParcel.GetDistanceToTransit() > 0.5).ToFlag();
TimeWindow timeWindow = new TimeWindow();
if (tour.JointTourSequence > 0)
{
foreach (IPersonDayWrapper pDay in householdDay.PersonDays)
{
ITourWrapper tInJoint = (ITourWrapper)pDay.Tours.Find(t => t.JointTourSequence == tour.JointTourSequence);
if (!(tInJoint == null))
{
// set jointTour time window
timeWindow.IncorporateAnotherTimeWindow(tInJoint.PersonDay.TimeWindow);
}
}
}
else if (tour.ParentTour == null)
{
timeWindow.IncorporateAnotherTimeWindow(personDay.TimeWindow);
}
timeWindow.SetBusyMinutes(Global.Settings.Times.EndOfRelevantWindow, Global.Settings.Times.MinutesInADay + 1);
var maxAvailableMinutes =
(tour.JointTourSequence > 0 || tour.ParentTour == null)
? timeWindow.MaxAvailableMinutesAfter(Global.Settings.Times.FiveAM)
: tour.ParentTour.DestinationDepartureTime - tour.ParentTour.DestinationArrivalTime;
var householdHasChildren = household.HasChildren;
var householdHasNoChildren = householdHasChildren ? false : true;
var fastestAvailableTimeOfDay =
tour.IsHomeBasedTour || tour.ParentTour == null
? 1
: tour.ParentTour.DestinationArrivalTime + (tour.ParentTour.DestinationDepartureTime - tour.ParentTour.DestinationArrivalTime) / 2;
var tourCategory = tour.GetTourCategory();
var secondaryFlag = ChoiceModelUtility.GetSecondaryFlag(tourCategory);
var workOrSchoolPatternFlag = personDay.GetIsWorkOrSchoolPattern().ToFlag();
var otherPatternFlag = personDay.GetIsOtherPattern().ToFlag();
int jointTourFlag = (tour.JointTourSequence > 0).ToFlag();
ChoiceModelUtility.DrawRandomTourTimePeriods(tour, tourCategory);
var segment = Global.ContainerDaySim.GetInstance <SamplingWeightsSettingsFactory>().SamplingWeightsSettings.GetTourDestinationSegment(tour.DestinationPurpose, tour.IsHomeBasedTour ? Global.Settings.TourPriorities.HomeBasedTour : Global.Settings.TourPriorities.WorkBasedTour, Global.Settings.Modes.Sov, person.PersonType);
var destinationSampler = new DestinationSampler(choiceProbabilityCalculator, segment, sampleSize, choice, tour.OriginParcel);
var tourDestinationUtilities = new TourDestinationUtilities(tour, sampleSize, secondaryFlag, personDay.GetIsWorkOrSchoolPattern().ToFlag(), personDay.GetIsOtherPattern().ToFlag(), fastestAvailableTimeOfDay, maxAvailableMinutes);
// get destination sample and perform code that used to be in SetUtilities below
var sampleItems = destinationSampler.SampleAndReturnTourDestinations(tourDestinationUtilities);
int observedMode = 0;
if (tour.Mode == Global.Settings.Modes.Bike || tour.Mode == Global.Settings.Modes.Walk)
{
observedMode = 0;
}
else if (tour.Mode == Global.Settings.Modes.Sov)
{
observedMode = 1;
}
else if (tour.Mode == Global.Settings.Modes.Hov2 || tour.Mode == Global.Settings.Modes.Hov3)
{
observedMode = 2;
}
else if (tour.Mode == Global.Settings.Modes.Transit)
{
observedMode = 3;
}
int purpose = tour.DestinationPurpose;
if (tour.ParentTour != null)
{
purpose = 0;
}
int prefix = 0;
int index = 0;
int destindex = 0;
int indChosen = -1;
foreach (var sampleItem in sampleItems)
{
bool parcelAvailable = sampleItem.Key.Available;
bool isChosen = sampleItem.Key.IsChosen;
double adjustmentFactor = sampleItem.Key.AdjustmentFactor;
var destinationParcel = ChoiceModelFactory.Parcels[sampleItem.Key.ParcelId];
destindex++;
if (isChosen)
{
indChosen = destindex;
}
if (!parcelAvailable)
{
continue;
}
int destinationTransitBand1 = (destinationParcel.GetDistanceToTransit() >= 0 && destinationParcel.GetDistanceToTransit() <= 0.25).ToFlag();
int destinationTransitBand3 = (destinationParcel.GetDistanceToTransit() > 0.5).ToFlag();
// var destinationSubzone = _zoneSubzones[destinationParcel.ZoneId][destinationTransitBand3];
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.Run(
tour.Household.RandomUtility,
originParcel.ZoneId,
destinationParcel.ZoneId,
tour.DestinationArrivalTime,
tour.DestinationDepartureTime,
tour.DestinationPurpose,
tour.CostCoefficient,
tour.TimeCoefficient,
tour.Person.IsDrivingAge,
tour.Household.VehiclesAvailable,
(tour.Household.OwnsAutomatedVehicles > 0),
tour.Person.GetTransitFareDiscountFraction(),
false,
Global.Settings.Modes.Walk, Global.Settings.Modes.Sov, Global.Settings.Modes.Hov2, Global.Settings.Modes.Transit);
for (var modeIndex = 0; modeIndex <= 3; modeIndex++)
{
var pathTypeModel = modeIndex == 0 ? pathTypeModels.First(x => x.Mode == Global.Settings.Modes.Walk)
: modeIndex == 1 ? pathTypeModels.First(x => x.Mode == Global.Settings.Modes.Sov)
: modeIndex == 2 ? pathTypeModels.First(x => x.Mode == Global.Settings.Modes.Hov2)
: pathTypeModels.First(x => x.Mode == Global.Settings.Modes.Transit);
var modeAvailable = pathTypeModel.Available;
var chosen = modeIndex == observedMode && destinationParcel == tour.DestinationParcel;
//if (chosen) Global.PrintFile.WriteLine("Sequence {0}: Chosen parcel {1} Mode {2} Available {3}",timesStartedRunModel,destinationParcel.Id,observedMode,modeAvailable);
var alternative = choiceProbabilityCalculator.GetAlternative(index++, modeAvailable, chosen);
alternative.Choice = destinationParcel;
if (!modeAvailable)
{
continue;
}
alternative.AddUtilityTerm(1, tour.Id);
alternative.AddUtilityTerm(9, purpose);
if (modeIndex == 0)
{
// WALK
alternative.AddUtilityTerm(prefix + 3, pathTypeModel.GeneralizedTimeLogsum);
alternative.AddUtilityTerm(prefix + 4, destinationParcel.NetIntersectionDensity2());
alternative.AddUtilityTerm(prefix + 5, destinationParcel.MixedUse4Index2());
}
else if (modeIndex == 1)
{
// SOV
alternative.AddUtilityTerm(prefix + 3, pathTypeModel.GeneralizedTimeLogsum);
alternative.AddUtilityTerm(prefix + 11, 1.0);
alternative.AddUtilityTerm(prefix + 12, destinationParcel.ParkingOffStreetPaidHourlyPriceBuffer2);
alternative.AddUtilityTerm(prefix + 14, fewerCarsThanDriversFlag);
alternative.AddUtilityTerm(prefix + 15, lowIncomeFlag);
alternative.AddUtilityTerm(prefix + 16, highIncomeFlag);
alternative.AddUtilityTerm(prefix + 17, missingIncomeFlag);
}
else if (modeIndex == 2)
{
// HOV
alternative.AddUtilityTerm(prefix + 3, pathTypeModel.GeneralizedTimeLogsum);
alternative.AddUtilityTerm(prefix + 21, 1.0);
alternative.AddUtilityTerm(prefix + 22, childUnder16Flag);
alternative.AddUtilityTerm(prefix + 23, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(prefix + 24, fewerCarsThanDriversFlag);
alternative.AddUtilityTerm(prefix + 25, lowIncomeFlag);
alternative.AddUtilityTerm(prefix + 26, highIncomeFlag);
alternative.AddUtilityTerm(prefix + 27, missingIncomeFlag);
}
else
{
// TRANSIT
alternative.AddUtilityTerm(prefix + 3, pathTypeModel.GeneralizedTimeLogsum);
alternative.AddUtilityTerm(prefix + 31, 1.0);
alternative.AddUtilityTerm(prefix + 32, childUnder16Flag);
alternative.AddUtilityTerm(prefix + 33, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(prefix + 34, fewerCarsThanDriversFlag);
alternative.AddUtilityTerm(prefix + 35, lowIncomeFlag);
alternative.AddUtilityTerm(prefix + 36, highIncomeFlag);
alternative.AddUtilityTerm(prefix + 37, missingIncomeFlag);
alternative.AddUtilityTerm(prefix + 38, originTransitBand1);
alternative.AddUtilityTerm(prefix + 39, originTransitBand3);
//alternative.AddUtilityTerm(prefix + 40, destinationTransitBand3);
}
alternative.AddUtilityTerm(2, adjustmentFactor);
// Size terms
alternative.AddUtilityTerm(prefix + 50, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(prefix + 51, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(prefix + 52, destinationParcel.EmploymentGovernment);
//alternative.AddUtilityTerm(prefix + 53, destinationParcel.EmploymentIndustrial);
alternative.AddUtilityTerm(prefix + 54, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(prefix + 55, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(prefix + 56, destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(prefix + 57, destinationParcel.EmploymentFood);
//alternative.AddUtilityTerm(prefix + 58, destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(prefix + 59, destinationParcel.Households);
//alternative.AddUtilityTerm(prefix + 60, destinationParcel.StudentsK8 + destinationParcel.StudentsHighSchool);
//alternative.AddUtilityTerm(prefix + 61, destinationParcel.StudentsUniversity);
}
}
//Global.PrintFile.WriteLine("Sequence {0}: Chosen parcel {1} is sample item {2} of {3}",timesStartedRunModel,tour.DestinationParcelId,indChosen,sampleItems.Count);
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper subtour, IEnumerable <IPathTypeModel> pathTypeModels, IParcelWrapper destinationParcel, int parentTourMode, int choice = Constants.DEFAULT_VALUE)
{
var household = subtour.Household;
var person = subtour.Person;
var personDay = subtour.PersonDay;
// household inputs
var income0To25KFlag = household.Has0To25KIncome.ToFlag();
var income25To50KFlag = household.Has25To50KIncome.ToFlag();
// person inputs
var maleFlag = person.IsMale.ToFlag();
// tour inputs
var sovTourFlag = (parentTourMode == Global.Settings.Modes.Sov).ToFlag();
var hov2TourFlag = (parentTourMode == Global.Settings.Modes.Hov2).ToFlag();
var bikeTourFlag = (parentTourMode == Global.Settings.Modes.Bike).ToFlag();
var walkTourFlag = (parentTourMode == Global.Settings.Modes.Walk).ToFlag();
// remaining inputs
var originParcel = subtour.OriginParcel;
var parkingDuration = ChoiceModelUtility.GetParkingDuration(person.IsFulltimeWorker);
var destinationParkingCost = destinationParcel.ParkingCostBuffer1(parkingDuration);
double escortPercentage;
double nonEscortPercentage;
ChoiceModelUtility.SetEscortPercentages(personDay, out escortPercentage, out nonEscortPercentage, true);
// paidRideShare is a special case
if (Global.Configuration.PaidRideShareModeIsAvailable)
{
var pathTypeExtra = pathTypeModels.First(x => x.Mode == Global.Settings.Modes.PaidRideShare);
var modeExtra = Global.Settings.Modes.PaidRideShare;
var availableExtra = pathTypeExtra.Available;
var generalizedTimeLogsumExtra = pathTypeExtra.GeneralizedTimeLogsum;
var distanceExtra = pathTypeExtra.PathDistance;
var alternative = choiceProbabilityCalculator.GetAlternative(modeExtra - 2, availableExtra, choice == modeExtra);
alternative.Choice = modeExtra;
alternative.AddNestedAlternative(_nestedAlternativeIds[modeExtra], _nestedAlternativeIndexes[modeExtra], THETA_PARAMETER);
if (availableExtra)
{
var extraCostPerMile = Global.Configuration.AV_PaidRideShareModeUsesAVs ?
Global.Configuration.AV_PaidRideShare_ExtraCostPerDistanceUnit : Global.Configuration.PaidRideShare_ExtraCostPerDistanceUnit;
var fixedCostPerRide = Global.Configuration.AV_PaidRideShareModeUsesAVs ?
Global.Configuration.AV_PaidRideShare_FixedCostPerRide : Global.Configuration.PaidRideShare_FixedCostPerRide;
var autoTimeCoefficient = Global.Configuration.AV_PaidRideShareModeUsesAVs ?
subtour.TimeCoefficient * (1.0 - Global.Configuration.AV_InVehicleTimeCoefficientDiscountFactor) : subtour.TimeCoefficient;
alternative.AddUtilityTerm(2, generalizedTimeLogsumExtra * autoTimeCoefficient);
alternative.AddUtilityTerm(2, distanceExtra * extraCostPerMile * subtour.CostCoefficient);
alternative.AddUtilityTerm(2, fixedCostPerRide * subtour.CostCoefficient);
var modeConstant = Global.Configuration.AV_PaidRideShareModeUsesAVs
? Global.Configuration.AV_PaidRideShare_ModeConstant
+ Global.Configuration.AV_PaidRideShare_DensityCoefficient * Math.Min(originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2, 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, 6000);
alternative.AddUtilityTerm(90, modeConstant);
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_Age26to35Coefficient * subtour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_Age18to25Coefficient * subtour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_AgeOver65Coefficient * (subtour.Person.Age >= 65).ToFlag());
}
}
foreach (var pathTypeModel in pathTypeModels)
{
var mode = pathTypeModel.Mode;
var available = pathTypeModel.Mode != Global.Settings.Modes.ParkAndRide && pathTypeModel.Available;
var generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
var alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
alternative.AddNestedAlternative(_nestedAlternativeIds[pathTypeModel.Mode], _nestedAlternativeIndexes[pathTypeModel.Mode], THETA_PARAMETER);
if (!available)
{
continue;
}
var modeTimeCoefficient = (Global.Configuration.AV_IncludeAutoTypeChoice && household.OwnsAutomatedVehicles > 0 && mode >= Global.Settings.Modes.Sov && mode <= Global.Settings.Modes.Hov3) ?
subtour.TimeCoefficient * (1.0 - Global.Configuration.AV_InVehicleTimeCoefficientDiscountFactor) : subtour.TimeCoefficient;
alternative.AddUtilityTerm(2, generalizedTimeLogsum * modeTimeCoefficient);
if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(20, 1);
// alternative.AddUtility(129, destinationParcel.MixedUse2Index1());
// alternative.AddUtility(128, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(127, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(126, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(125, originParcel.TotalEmploymentDensity1());
// alternative.AddUtility(124, originParcel.MixedUse2Index1());
// alternative.AddUtility(123, Math.Log(destinationParcel.StopsTransitBuffer1+1));
// alternative.AddUtility(122, Math.Log(originParcel.StopsTransitBuffer1+1));
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * subtour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(88, sovTourFlag);
alternative.AddUtilityTerm(89, hov2TourFlag);
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * subtour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(88, sovTourFlag);
alternative.AddUtilityTerm(89, hov2TourFlag);
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * subtour.CostCoefficient));
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(55, income25To50KFlag);
alternative.AddUtilityTerm(58, sovTourFlag);
alternative.AddUtilityTerm(59, hov2TourFlag);
}
else if (mode == Global.Settings.Modes.Bike)
{
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(69, bikeTourFlag);
// alternative.AddUtility(169, destinationParcel.MixedUse4Index1());
// alternative.AddUtility(168, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(167, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(166, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(165, originParcel.TotalEmploymentDensity1());
// alternative.AddUtility(164, originParcel.MixedUse4Index1());
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(70, 1);
alternative.AddUtilityTerm(79, walkTourFlag);
// alternative.AddUtility(179, destinationParcel.MixedUse4Index1());
// alternative.AddUtility(178, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(177, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(176, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(175, originParcel.TotalEmploymentDensity1());
// alternative.AddUtility(174, originParcel.MixedUse4Index1());
}
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, IHouseholdDayWrapper householdDay, ITourWrapper tour,
IParcelWrapper destinationParcel, int householdCars,
int constrainedMode, int constrainedArrivalTime, int constrainedDepartureTime, HTourModeTime choice = null)
{
IHouseholdWrapper household = tour.Household;
IPersonWrapper person = tour.Person;
IPersonDayWrapper personDay = tour.PersonDay;
IHouseholdTotals householdTotals = household.HouseholdTotals;
// household inputs
int childrenUnder5 = householdTotals.ChildrenUnder5;
int childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
int nonworkingAdults = householdTotals.NonworkingAdults;
int retiredAdults = householdTotals.RetiredAdults;
int onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
int twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
//var householdCars = household.VehiclesAvailable; MABADD now an input parameter
int noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(householdCars);
int carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(householdCars);
int carsLessThanWorkersFlag = household.GetFlagForCarsLessThanWorkers(householdCars);
int income0To25KFlag = household.Has0To25KIncome.ToFlag();
int income100KPlusFlag = household.Has100KPlusIncome.ToFlag();
// 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 fulltimeWorkerFlag = person.IsFulltimeWorker.ToFlag();
int childAge5Through15Flag = person.IsChildAge5Through15.ToFlag();
int childUnder5Flag = person.IsChildUnder5.ToFlag();
int maleFlag = person.IsMale.ToFlag();
int ageUnder30Flag = person.AgeIsLessThan30.ToFlag();
int ageBetween51And98Flag = person.AgeIsBetween51And98.ToFlag();
int adultFlag = person.IsAdult.ToFlag();
// person-day inputs
int homeBasedToursOnlyFlag = (personDay == null) ? 1 : personDay.OnlyHomeBasedToursExist().ToFlag();
int firstSimulatedHomeBasedTourFlag = (personDay == null) ? 1 : personDay.IsFirstSimulatedHomeBasedTour().ToFlag();
int laterSimulatedHomeBasedTourFlag = (personDay == null) ? 0 : personDay.IsLaterSimulatedHomeBasedTour().ToFlag();
int totalStops = (personDay == null) ? 0 : personDay.GetTotalStops();
int totalSimulatedStops = (personDay == null) ? 0 : personDay.GetTotalSimulatedStops();
int escortStops = (personDay == null) ? 0 : personDay.EscortStops;
int homeBasedTours = (personDay == null) ? 1 : personDay.HomeBasedTours;
int simulatedHomeBasedTours = (personDay == null) ? 0 : personDay.SimulatedHomeBasedTours;
// tour inputs
int escortTourFlag = tour.IsEscortPurpose().ToFlag();
int shoppingTourFlag = tour.IsShoppingPurpose().ToFlag();
int mealTourFlag = tour.IsMealPurpose().ToFlag();
int socialTourFlag = tour.IsSocialPurpose().ToFlag();
int personalBusinessTourFlag = tour.IsPersonalBusinessPurpose().ToFlag();
int recreationTourFlag = tour.IsRecreationPurpose().ToFlag();
int medicalTourFlag = tour.IsMedicalPurpose().ToFlag();
IParcelWrapper originParcel = tour.OriginParcel;
//var destinationParcel = tour.DestinationParcel; MABADD now an input parameter
int jointTourFlag = (tour.JointTourSequence > 0) ? 1 : 0;
int partialHalfTour1Flag = (tour.PartialHalfTour1Sequence > 0) ? 1 : 0;
int partialHalfTour2Flag = (tour.PartialHalfTour2Sequence > 0) ? 1 : 0;
bool partialHalfTour = (tour.PartialHalfTour1Sequence > 0 || tour.PartialHalfTour2Sequence > 0);
int fullHalfTour1Flag = (tour.FullHalfTour1Sequence > 0) ? 1 : 0;
int fullHalfTour2Flag = (tour.FullHalfTour2Sequence > 0) ? 1 : 0;
int parentTourMode = tour.ParentTour == null ? 0 : tour.ParentTour.Mode;
// remaining inputs
// Higher priority tour of 2+ tours for the same purpose
int highPrioritySameFlag = (personDay == null) ? 1 : (tour.GetTotalToursByPurpose() > tour.GetTotalSimulatedToursByPurpose() && tour.GetTotalSimulatedToursByPurpose() == 1).ToFlag();
// Lower priority tour(s) of 2+ tours for the same purpose
int lowPrioritySameFlag = (personDay == null) ? 0 : (tour.GetTotalSimulatedToursByPurpose() > 1).ToFlag();
// Higher priority tour of 2+ tours for different purposes
int highPriorityDifferentFlag = (personDay == null) ? 0 : (personDay.IsFirstSimulatedHomeBasedTour() && personDay.HomeBasedToursExist()).ToFlag() * (1 - highPrioritySameFlag);
// Lower priority tour of 2+ tours for different purposes
int lowPriorityDifferentFlag = (personDay == null) ? 0 : (personDay.IsLaterSimulatedHomeBasedTour() && personDay.HomeBasedToursExist()).ToFlag() * (1 - lowPrioritySameFlag);
ITimeWindow timeWindow = (householdDay == null) ? new TimeWindow() : tour.GetRelevantTimeWindow(householdDay);
int totalMinutesAvailableInDay = timeWindow.TotalAvailableMinutes(1, 1440);
int bigPeriodCount = DayPeriod.H_BIG_DAY_PERIOD_TOTAL_TOUR_TIMES;
int nPeriodCombs = bigPeriodCount * (bigPeriodCount + 1) / 2;
bool useTimeComponents = Global.Configuration.IsInEstimationMode || constrainedArrivalTime == 0 || constrainedDepartureTime == 0;
int componentIndex = 0;
int periodComb = -1;
//set components
if (useTimeComponents)
{
for (int arrivalPeriodIndex = 0; arrivalPeriodIndex < bigPeriodCount; arrivalPeriodIndex++)
{
MinuteSpan arrivalPeriod = DayPeriod.HBigDayPeriods[arrivalPeriodIndex];
int arrivalPeriodAvailableMinutes = timeWindow.TotalAvailableMinutes(arrivalPeriod.Start, arrivalPeriod.End);
for (int departurePeriodIndex = arrivalPeriodIndex; departurePeriodIndex < bigPeriodCount; departurePeriodIndex++)
{
MinuteSpan departurePeriod = DayPeriod.HBigDayPeriods[departurePeriodIndex];
int departurePeriodAvailableMinutes = timeWindow.TotalAvailableMinutes(departurePeriod.Start, departurePeriod.End);
if (arrivalPeriod == departurePeriod)
{
componentIndex = arrivalPeriodIndex;
choiceProbabilityCalculator.CreateUtilityComponent(componentIndex);
ChoiceProbabilityCalculator.Component arrivalComponent = choiceProbabilityCalculator.GetUtilityComponent(componentIndex);
if (arrivalPeriodAvailableMinutes > 0)
{
double hoursArrival = arrivalPeriod.Middle / 60.0;
int firstCoef = 300;
arrivalComponent.AddUtilityTerm(300, Math.Log(arrivalPeriodAvailableMinutes));
//arrival shift variables
arrivalComponent.AddUtilityTerm(firstCoef + 2, partTimeWorkerFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 3, nonworkingAdultFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 4, universityStudentFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 5, retiredAdultFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 6, drivingAgeStudentFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 7, childAge5Through15Flag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 8, childUnder5Flag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 9, escortTourFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 10, shoppingTourFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 11, mealTourFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 12, socialTourFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 13, personalBusinessTourFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 14, recreationTourFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 15, medicalTourFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 16, income0To25KFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 17, income100KPlusFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 18, highPrioritySameFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 19, lowPrioritySameFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 20, highPriorityDifferentFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 21, lowPriorityDifferentFlag * hoursArrival);
// arrivalComponent.AddUtilityTerm(firstCoef + 22, jointTourFlag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 23, partialHalfTour1Flag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 24, fullHalfTour1Flag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 25, partialHalfTour2Flag * hoursArrival);
arrivalComponent.AddUtilityTerm(firstCoef + 26, fullHalfTour2Flag * hoursArrival);
}
componentIndex = bigPeriodCount + departurePeriodIndex;
choiceProbabilityCalculator.CreateUtilityComponent(componentIndex);
ChoiceProbabilityCalculator.Component departureComponent = choiceProbabilityCalculator.GetUtilityComponent(componentIndex);
if (departurePeriodAvailableMinutes > 0)
{
departureComponent.AddUtilityTerm(300, Math.Log(departurePeriodAvailableMinutes));
}
}
// set period combination component
periodComb++;
componentIndex = 2 * bigPeriodCount + periodComb;
choiceProbabilityCalculator.CreateUtilityComponent(componentIndex);
ChoiceProbabilityCalculator.Component combinationComponent = choiceProbabilityCalculator.GetUtilityComponent(componentIndex);
if (arrivalPeriodAvailableMinutes > 0 && departurePeriodAvailableMinutes > 0)
{
double hoursDuration = (departurePeriod.Middle - arrivalPeriod.Middle) / 60.0;
int firstCoef = 700;
//combination constants
combinationComponent.AddUtilityTerm(firstCoef + periodComb, 1.0);
// duration shift variables
// combinationComponent.AddUtilityTerm(firstCoef + 31, escortTourFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 32, shoppingTourFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 33, mealTourFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 34, socialTourFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 35, personalBusinessTourFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 36, recreationTourFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 37, medicalTourFlag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 38, highPrioritySameFlag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 39, lowPrioritySameFlag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 40, highPriorityDifferentFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 41, lowPriorityDifferentFlag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 42, partTimeWorkerFlag * hoursDuration);
// combinationComponent.AddUtilityTerm(firstCoef + 43, jointTourFlag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 44, partialHalfTour1Flag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 45, fullHalfTour1Flag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 46, partialHalfTour2Flag * hoursDuration);
combinationComponent.AddUtilityTerm(firstCoef + 47, fullHalfTour2Flag * hoursDuration);
// peak-to-peak variables
if (arrivalPeriod.Index == DayPeriod.AM_PEAK && departurePeriod.Index == DayPeriod.PM_PEAK)
{
combinationComponent.AddUtilityTerm(firstCoef + 48, fulltimeWorkerFlag);
combinationComponent.AddUtilityTerm(firstCoef + 49, income0To25KFlag);
combinationComponent.AddUtilityTerm(firstCoef + 50, income100KPlusFlag);
}
}
}
}
}
for (int mode = Global.Settings.Modes.Walk; mode <= Global.Settings.Modes.WalkRideBike; mode++)
{
componentIndex = 2 * bigPeriodCount + nPeriodCombs + mode - 1;
choiceProbabilityCalculator.CreateUtilityComponent(componentIndex);
ChoiceProbabilityCalculator.Component modeComponent = choiceProbabilityCalculator.GetUtilityComponent(componentIndex);
if (mode == Global.Settings.Modes.Walk)
{
modeComponent.AddUtilityTerm(70, 1.0);
modeComponent.AddUtilityTerm(71, maleFlag);
// modeComponent.AddUtilityTerm(72, ageUnder30Flag);
// modeComponent.AddUtilityTerm(73, ageBetween51And98Flag);
// modeComponent.AddUtilityTerm(74, income0To25KFlag);
modeComponent.AddUtilityTerm(75, income100KPlusFlag);
// modeComponent.AddUtilityTerm(76, noCarsInHouseholdFlag);
modeComponent.AddUtilityTerm(77, carsLessThanDriversFlag);
modeComponent.AddUtilityTerm(179, destinationParcel.MixedUse4Index1());
modeComponent.AddUtilityTerm(179, destinationParcel.TotalEmploymentDensity1() / 5000.0);
modeComponent.AddUtilityTerm(179, destinationParcel.NetIntersectionDensity1() / 50.0);
modeComponent.AddUtilityTerm(179, originParcel.NetIntersectionDensity1() / 50.0);
modeComponent.AddUtilityTerm(179, originParcel.HouseholdDensity1() / 1000.0);
modeComponent.AddUtilityTerm(179, originParcel.MixedUse4Index1());
}
else if (mode == Global.Settings.Modes.Bike)
{
modeComponent.AddUtilityTerm(60, 1);
modeComponent.AddUtilityTerm(61, maleFlag);
modeComponent.AddUtilityTerm(62, ageUnder30Flag);
modeComponent.AddUtilityTerm(63, ageBetween51And98Flag);
modeComponent.AddUtilityTerm(64, income0To25KFlag);
// modeComponent.AddUtilityTerm(65, income100KPlusFlag);
// modeComponent.AddUtilityTerm(66, noCarsInHouseholdFlag);
modeComponent.AddUtilityTerm(67, carsLessThanDriversFlag);
modeComponent.AddUtilityTerm(169, destinationParcel.MixedUse4Index2());
modeComponent.AddUtilityTerm(169, destinationParcel.TotalEmploymentDensity2() / 20000.0);
modeComponent.AddUtilityTerm(169, destinationParcel.NetIntersectionDensity2() / 200.0);
modeComponent.AddUtilityTerm(164, originParcel.NetIntersectionDensity2() / 200.0);
modeComponent.AddUtilityTerm(164, originParcel.HouseholdDensity2() / 4000.0);
modeComponent.AddUtilityTerm(164, originParcel.MixedUse4Index2());
}
else if (mode == Global.Settings.Modes.Sov)
{
// modeComponent.AddUtilityTerm(50, 1);
modeComponent.AddUtilityTerm(54, income0To25KFlag);
// modeComponent.AddUtilityTerm(55, income100KPlusFlag);
modeComponent.AddUtilityTerm(57, carsLessThanWorkersFlag);
}
else if (mode == Global.Settings.Modes.HovDriver)
{
modeComponent.AddUtilityTerm(31, childrenUnder5);
modeComponent.AddUtilityTerm(32, childrenAge5Through15);
modeComponent.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
modeComponent.AddUtilityTerm(36, noCarsInHouseholdFlag);
modeComponent.AddUtilityTerm(37, carsLessThanDriversFlag);
modeComponent.AddUtilityTerm(40, 1);
modeComponent.AddUtilityTerm(41, onePersonHouseholdFlag);
}
else if (mode == Global.Settings.Modes.HovPassenger)
{
modeComponent.AddUtilityTerm(30, 1);
modeComponent.AddUtilityTerm(31, childrenUnder5);
modeComponent.AddUtilityTerm(32, childrenAge5Through15);
modeComponent.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
modeComponent.AddUtilityTerm(38, onePersonHouseholdFlag);
modeComponent.AddUtilityTerm(39, twoPersonHouseholdFlag);
modeComponent.AddUtilityTerm(36, noCarsInHouseholdFlag);
modeComponent.AddUtilityTerm(37, carsLessThanDriversFlag);
}
else if (mode == Global.Settings.Modes.Transit)
{
modeComponent.AddUtilityTerm(20, 1);
modeComponent.AddUtilityTerm(21, maleFlag);
modeComponent.AddUtilityTerm(22, ageUnder30Flag);
// modeComponent.AddUtilityTerm(23, ageBetween51And98Flag);
modeComponent.AddUtilityTerm(24, income0To25KFlag);
modeComponent.AddUtilityTerm(25, income100KPlusFlag);
modeComponent.AddUtilityTerm(26, noCarsInHouseholdFlag);
modeComponent.AddUtilityTerm(27, carsLessThanDriversFlag);
modeComponent.AddUtilityTerm(129, destinationParcel.MixedUse2Index1());
modeComponent.AddUtilityTerm(129, destinationParcel.TotalEmploymentDensity1() / 5000.0);
modeComponent.AddUtilityTerm(129, destinationParcel.NetIntersectionDensity1() / 50.0);
modeComponent.AddUtilityTerm(124, originParcel.NetIntersectionDensity1() / 50.0);
modeComponent.AddUtilityTerm(124, originParcel.HouseholdDensity1() / 1000.0);
modeComponent.AddUtilityTerm(124, originParcel.MixedUse2Index1());
modeComponent.AddUtilityTerm(123, Math.Log(destinationParcel.StopsTransitBuffer1 + 1));
modeComponent.AddUtilityTerm(122, Math.Log(originParcel.StopsTransitBuffer1 + 1));
}
else if (mode == Global.Settings.Modes.CarParkRideWalk)
{
modeComponent.AddUtilityTerm(10, 1);
modeComponent.AddUtilityTerm(16, noCarsInHouseholdFlag);
modeComponent.AddUtilityTerm(17, carsLessThanWorkersFlag);
modeComponent.AddUtilityTerm(129, destinationParcel.MixedUse2Index1());
modeComponent.AddUtilityTerm(129, destinationParcel.TotalEmploymentDensity1() / 5000.0);
modeComponent.AddUtilityTerm(129, destinationParcel.NetIntersectionDensity1() / 50.0);
modeComponent.AddUtilityTerm(123, Math.Log(destinationParcel.StopsTransitBuffer1 + 1));
}
// JLB 201508 Need to provide mode-specific constants and add any other needed mode-specific terms
else if (mode == Global.Settings.Modes.CarKissRideWalk)
{
modeComponent.AddUtilityTerm(10, 1);
}
else if (mode == Global.Settings.Modes.BikeParkRideWalk)
{
modeComponent.AddUtilityTerm(10, 1);
}
else if (mode == Global.Settings.Modes.BikeParkRideBike)
{
modeComponent.AddUtilityTerm(10, 1);
}
else if (mode == Global.Settings.Modes.BikeOnTransit)
{
modeComponent.AddUtilityTerm(10, 1);
}
else if (mode == Global.Settings.Modes.CarParkRideBike)
{
modeComponent.AddUtilityTerm(10, 1);
}
else if (mode == Global.Settings.Modes.WalkRideBike)
{
modeComponent.AddUtilityTerm(10, 1);
}
// JLB 201508 may need to enhance next block to accommodate new transit modes
if (mode == Global.Settings.Modes.Walk || mode == Global.Settings.Modes.Bike || mode == Global.Settings.Modes.Hov2 ||
mode == Global.Settings.Modes.Hov3 || mode == Global.Settings.Modes.Transit)
{
int firstCoef = 200 + 10 * mode;
// modeComponent.AddUtilityTerm(firstCoef + 0, escortTourFlag);
// modeComponent.AddUtilityTerm(firstCoef + 1, shoppingTourFlag);
// modeComponent.AddUtilityTerm(firstCoef + 2, mealTourFlag);
// modeComponent.AddUtilityTerm(firstCoef + 3, socialTourFlag);
// modeComponent.AddUtilityTerm(firstCoef + 4, personalBusinessTourFlag);
// modeComponent.AddUtilityTerm(firstCoef + 5, recreationTourFlag);
// modeComponent.AddUtilityTerm(firstCoef + 6, medicalTourFlag);
// modeComponent.AddUtilityTerm(firstCoef + 7, jointTourFlag);
modeComponent.AddUtilityTerm(firstCoef + 8, Math.Min(partialHalfTour1Flag + partialHalfTour2Flag, 1.0));
modeComponent.AddUtilityTerm(firstCoef + 9, Math.Min(fullHalfTour1Flag + fullHalfTour2Flag, 1.0));
// modeComponent.AddUtilityTerm(290+mode, mode == parentTourMode ? 1:0 );
}
// modeComponent.AddUtilityTerm(298, mode>=Global.Settings.Modes.Sov && mode<=Global.Settings.Modes.Hov3 && parentTourMode==Global.Settings.Modes.Sov ? 1:0 );
// modeComponent.AddUtilityTerm(299, mode>=Global.Settings.Modes.Sov && mode<=Global.Settings.Modes.Hov3 && parentTourMode>=Global.Settings.Modes.Hov2 && parentTourMode<=Global.Settings.Modes.Hov3 ? 1:0 );
}
//loop on all alternatives, using modeTimes objects
{
foreach (HTourModeTime modeTimes in HTourModeTime.ModeTimes[ParallelUtility.threadLocalAssignedIndex.Value])
{
MinuteSpan arrivalPeriod = modeTimes.ArrivalPeriod;
int arrivalPeriodAvailableMinutes = timeWindow.TotalAvailableMinutes(arrivalPeriod.Start, arrivalPeriod.End);
MinuteSpan departurePeriod = modeTimes.DeparturePeriod;
int departurePeriodAvailableMinutes = timeWindow.TotalAvailableMinutes(departurePeriod.Start, departurePeriod.End);
periodComb = modeTimes.PeriodCombinationIndex;
int mode = modeTimes.Mode;
int altIndex = modeTimes.Index;
//set availabillity based on time window variables and any constrained choices
bool available = (modeTimes.LongestFeasibleWindow != null) && (mode > 0)
//&& (mode <= Global.Settings.Modes.Transit) //JLB 201508 commented out to make transit combo modes available
&& (person.Age >= 18 || (modeTimes.Mode != Global.Settings.Modes.Sov && modeTimes.Mode != Global.Settings.Modes.HovDriver))
//TODO JLB 201508 may need to adjust the followign line to allow transit combo or other modes modes on partial half tour mode choice: line currently prohibits them as well as SOV and HOVPassenger
&& (constrainedMode > 0 || mode == Global.Settings.Modes.Walk || mode == Global.Settings.Modes.Bike || mode == Global.Settings.Modes.HovDriver || mode == Global.Settings.Modes.Transit || !partialHalfTour) &&
(constrainedMode <= 0 || constrainedMode == mode) &&
(constrainedArrivalTime <= 0 || (constrainedArrivalTime >= arrivalPeriod.Start && constrainedArrivalTime <= arrivalPeriod.End)) &&
(constrainedDepartureTime <= 0 || (constrainedDepartureTime >= departurePeriod.Start && constrainedDepartureTime <= departurePeriod.End));
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(altIndex, available,
choice != null && choice.Index == altIndex);
alternative.Choice = modeTimes; // JLB added 20130420
alternative.AddNestedAlternative(HTourModeTime.TotalTourModeTimes + periodComb + 1, periodComb, THETA_PARAMETER);
if (Global.Configuration.IsInEstimationMode && choice != null && altIndex == choice.Index)
{
Global.PrintFile.WriteLine("Aper Dper Mode {0} {1} {2} Travel Times {3} {4} Window {5} {6}",
arrivalPeriod.Index, departurePeriod.Index, mode,
modeTimes.ModeAvailableToDestination ? modeTimes.TravelTimeToDestination : -1,
modeTimes.ModeAvailableFromDestination ? modeTimes.TravelTimeFromDestination : -1,
modeTimes.LongestFeasibleWindow != null ? modeTimes.LongestFeasibleWindow.Start : -1,
modeTimes.LongestFeasibleWindow != null ? modeTimes.LongestFeasibleWindow.End : -1);
}
/* if (altIndex == 0)
* {
* alternative.AddUtilityTerm(991, tour.Household.Id);
* alternative.AddUtilityTerm(992, tour.Person.Id);
* alternative.AddUtilityTerm(993, tour.PersonDay.Day);
* alternative.AddUtilityTerm(994, tour.Sequence);
* alternative.AddUtilityTerm(995, constrainedMode);
* alternative.AddUtilityTerm(996, constrainedArrivalTime);
* alternative.AddUtilityTerm(997, constrainedDepartureTime);
* alternative.AddUtilityTerm(998, tour.DestinationPurpose);
* alternative.AddUtilityTerm(999, (tour.ParentTour == null) ? 0 : 1);
* }
*/
//if in application mode and combination is not available, can skip the rest
if (!Global.Configuration.IsInEstimationMode && !alternative.Available)
{
continue;
}
if (useTimeComponents)
{
// arrival period utility component
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(arrivalPeriod.Index));
// departure period utility component
alternative.AddUtilityComponent(
choiceProbabilityCalculator.GetUtilityComponent(bigPeriodCount + departurePeriod.Index));
// period combination utility component
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(2 * bigPeriodCount + periodComb));
}
// mode utility component
alternative.AddUtilityComponent(
choiceProbabilityCalculator.GetUtilityComponent(2 * bigPeriodCount + nPeriodCombs + mode - 1));
//even in estimation mode, do not need the rest of the code if not available
if (!alternative.Available)
{
continue;
}
// set parking cost for period combination
double parkingDuration = (departurePeriod == arrivalPeriod
? (arrivalPeriod.End - arrivalPeriod.Start) / 2.0
: departurePeriod.Middle - arrivalPeriod.Middle) / 60.0;
// parking at work is free if no paid parking at work and tour goes to usual workplace
//TODO JLB 20150827 Check to see if parkign cost is included in HTourModeTime impedance. If so, drop it here. otherwise fix this to use TOD-based parking cost.
double destinationParkingCost = (!Global.Configuration.IsInEstimationMode &&
Global.Configuration.ShouldRunPayToParkAtWorkplaceModel &&
tour.Person.UsualWorkParcel != null &&
destinationParcel == tour.Person.UsualWorkParcel &&
person.PaidParkingAtWorkplace == 0)
? 0.0
: destinationParcel.ParkingCostBuffer1(parkingDuration);
double parkingCostFraction = (mode == Global.Settings.Modes.Sov)
? 1.0
: (mode == Global.Settings.Modes.Hov2)
? 1.0 / Global.Configuration.Coefficients_HOV2CostDivisor_Work
: (mode == Global.Settings.Modes.Hov3)
? 1.0 / Global.Configuration.Coefficients_HOV3CostDivisor_Work
: 0.0;
double minimumTimeNeeded = modeTimes.TravelTimeToDestination + modeTimes.TravelTimeFromDestination +
Global.Settings.Times.MinimumActivityDuration;
alternative.AddUtilityTerm(1, modeTimes.GeneralizedTimeToDestination + modeTimes.GeneralizedTimeFromDestination);
alternative.AddUtilityTerm(2, destinationParkingCost * parkingCostFraction);
//alternative.AddUtilityTerm(3,
// Math.Log(Math.Min(1140, modeTimes.LongestFeasibleWindow.End - modeTimes.LongestFeasibleWindow.Start -
// minimumTimeNeeded + 1.0)));
// JLB 20140204 replaced coeff 3 with a different time window formulation: time pressure
// instead of having positive utility for increasing time window, have negative utility for decreasing time window
alternative.AddUtilityTerm(3,
Math.Log(Math.Max(Constants.EPSILON, 1 -
Math.Pow(minimumTimeNeeded / (Math.Min(1140, modeTimes.LongestFeasibleWindow.End - modeTimes.LongestFeasibleWindow.Start)), 0.4)
)));
//alternative.AddUtilityTerm(4, Math.Log((totalMinutesAvailableInDay + 1.0)/(minimumTimeNeeded + 1.0)));
//alternative.AddUtilityTerm(4,
// Math.Log(Math.Max(Constants.EPSILON, 1 - minimumTimeNeeded/(Math.Min(1140, totalMinutesAvailableInDay)))));
alternative.AddUtilityTerm(5,
(maleFlag == 0 && mode == Global.Settings.Modes.Walk &&
arrivalPeriod.Index >= DayPeriod.EVENING)
? 1
: 0);
alternative.AddUtilityTerm(5,
(maleFlag == 0 && mode == Global.Settings.Modes.Walk &&
departurePeriod.Index >= DayPeriod.EVENING)
? 1
: 0);
}
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, IEnumerable <IPathTypeModel> pathTypeModels, IParcelWrapper destinationParcel, int householdCars, int choice = Constants.DEFAULT_VALUE)
{
IHouseholdWrapper household = tour.Household;
Framework.DomainModels.Models.IHouseholdTotals householdTotals = household.HouseholdTotals;
IPersonDayWrapper personDay = tour.PersonDay;
IPersonWrapper person = tour.Person;
// household inputs
int childrenUnder5 = householdTotals.ChildrenUnder5;
int childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
// var nonworkingAdults = householdTotals.NonworkingAdults;
// var retiredAdults = householdTotals.RetiredAdults;
int onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
int twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
int noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(householdCars);
int carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(householdCars);
int carsLessThanWorkersFlag = household.GetFlagForCarsLessThanWorkers(householdCars);
int income0To25KFlag = household.Has0To25KIncome.ToFlag();
int income75kPlusFlag = household.Has75KPlusIncome.ToFlag();
int incomeOver100Flag = household.Has100KPlusIncome.ToFlag();
// person inputs
int maleFlag = person.IsMale.ToFlag();
int ageBetween51And98Flag = person.AgeIsBetween51And98.ToFlag();
int univStudentFlag = person.IsUniversityStudent.ToFlag();
IParcelWrapper originParcel = tour.OriginParcel;
int parkingDuration = ChoiceModelUtility.GetParkingDuration(person.IsFulltimeWorker);
// parking at work is free if no paid parking at work and tour goes to usual workplace
double destinationParkingCost = (Global.Configuration.ShouldRunPayToParkAtWorkplaceModel && tour.Person.UsualWorkParcel != null &&
destinationParcel == tour.Person.UsualWorkParcel && person.PaidParkingAtWorkplace == 0) ? 0.0 : destinationParcel.ParkingCostBuffer1(parkingDuration);
ChoiceModelUtility.SetEscortPercentages(personDay, out double escortPercentage, out double nonEscortPercentage);
foreach (IPathTypeModel pathTypeModel in pathTypeModels)
{
int mode = pathTypeModel.Mode;
double generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
bool available = (pathTypeModel.Mode != Global.Settings.Modes.PaidRideShare || Global.Configuration.PaidRideShareModeIsAvailable) &&
pathTypeModel.Available;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
alternative.AddNestedAlternative(_nestedAlternativeIds[pathTypeModel.Mode], _nestedAlternativeIndexes[pathTypeModel.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.ParkAndRide)
{
alternative.AddUtilityTerm(10, 1);
alternative.AddUtilityTerm(11, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(13, carsLessThanWorkersFlag);
// alternative.AddUtility(129, destinationParcel.MixedUse2Index1());
alternative.AddUtilityTerm(130, Math.Log(destinationParcel.TotalEmploymentDensity1() + 1) * 2553.0 / Math.Log(2553.0));
alternative.AddUtilityTerm(128, destinationParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(127, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(123, Math.Log(destinationParcel.StopsTransitBuffer1 + 1));
//alternative.AddUtilityTerm(120, Math.Log(destinationParcel.TotalEmploymentDensity1() + 1) * 2553.0 / Math.Log(2553.0));
alternative.AddUtilityTerm(118, destinationParcel.TotalEmploymentDensity1());
//alternative.AddUtilityTerm(117, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(113, Math.Log(destinationParcel.StopsTransitBuffer1 + 1));
alternative.AddUtilityTerm(112, Math.Log(originParcel.StopsTransitBuffer1 + 1));
alternative.AddUtilityTerm(115, Math.Log(pathTypeModel.PathParkAndRideNodeCapacity + 1));
}
else if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(20, 1);
alternative.AddUtilityTerm(21, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(22, carsLessThanDriversFlag); //for calibration
// alternative.AddUtility(129, destinationParcel.MixedUse2Index1());
alternative.AddUtilityTerm(130, Math.Log(destinationParcel.TotalEmploymentDensity1() + 1) * 2553.0 / Math.Log(2553.0));
alternative.AddUtilityTerm(128, destinationParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(127, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(126, originParcel.NetIntersectionDensity1());
//alternative.AddUtilityTerm(125, originParcel.HouseholdDensity1());
//alternative.AddUtilityTerm(124, originParcel.MixedUse2Index1());
alternative.AddUtilityTerm(123, Math.Log(destinationParcel.StopsTransitBuffer1 + 1));
alternative.AddUtilityTerm(122, Math.Log(originParcel.StopsTransitBuffer1 + 1));
alternative.AddUtilityTerm(180, univStudentFlag);
alternative.AddUtilityTerm(100, income75kPlusFlag);
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / Global.Configuration.Coefficients_HOV3CostDivisor_Work));
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(31, childrenUnder5);
alternative.AddUtilityTerm(32, childrenAge5Through15);
// alternative.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
alternative.AddUtilityTerm(35, pathTypeModel.PathDistance.AlmostEquals(0) ? 0 : Math.Log(pathTypeModel.PathDistance));
alternative.AddUtilityTerm(36, noCarsInHouseholdFlag); // for calibration of hov3 vs hov2
alternative.AddUtilityTerm(37, carsLessThanDriversFlag); // for calibration of hov3 vs hov2
alternative.AddUtilityTerm(38, onePersonHouseholdFlag);
alternative.AddUtilityTerm(39, twoPersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(42, carsLessThanDriversFlag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / Global.Configuration.Coefficients_HOV2CostDivisor_Work));
alternative.AddUtilityTerm(31, childrenUnder5);
alternative.AddUtilityTerm(32, childrenAge5Through15);
// alternative.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
alternative.AddUtilityTerm(35, pathTypeModel.PathDistance.AlmostEquals(0) ? 0 : Math.Log(pathTypeModel.PathDistance));
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(42, carsLessThanDriversFlag);
alternative.AddUtilityTerm(48, onePersonHouseholdFlag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(1, (destinationParkingCost) * tour.CostCoefficient);
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(53, carsLessThanWorkersFlag);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(131, escortPercentage);
alternative.AddUtilityTerm(132, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Bike)
{
double class1Dist
= Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions
? ImpedanceRoster.GetValue("class1distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable
: 0;
double class2Dist =
Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions
? ImpedanceRoster.GetValue("class2distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable
: 0;
// double worstDist = Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions ?
// ImpedanceRoster.GetValue("worstdistance", mode, Global.Settings.PathTypes.FullNetwork,
// Global.Settings.VotGroups.Medium, tour.DestinationArrivalTime,originParcel, destinationParcel).Variable : 0;
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(63, ageBetween51And98Flag);
alternative.AddUtilityTerm(67, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(68, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(169, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(168, destinationParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(167, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(166, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(165, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(164, originParcel.MixedUse4Index1());
alternative.AddUtilityTerm(162, (class1Dist > 0).ToFlag());
alternative.AddUtilityTerm(162, (class2Dist > 0).ToFlag());
// alternative.AddUtility(163, (worstDist > 0).ToFlag());
alternative.AddUtilityTerm(170, 1.0 * destinationParcel.MixedUse4Index1()
+ 0.00002 * destinationParcel.TotalEmploymentDensity1()
+ 0.001 * destinationParcel.NetIntersectionDensity1()
+ 0.001 * originParcel.NetIntersectionDensity1()
+ 0.0002 * originParcel.HouseholdDensity1()
+ 1.0 * originParcel.MixedUse4Index1());
//alternative.AddUtilityTerm(261, originParcel.PCA_DensityTerm_Buffer1());
//alternative.AddUtilityTerm(262, originParcel.PCA_WalkabilityTerm_Buffer1());
//alternative.AddUtilityTerm(263, originParcel.PCA_MixedUseTerm_Buffer1());
//alternative.AddUtilityTerm(264, originParcel.PCA_TransitAccessTerm_Buffer1());
//alternative.AddUtilityTerm(261, destinationParcel.PCA_DensityTerm_Buffer1());
//alternative.AddUtilityTerm(262, destinationParcel.PCA_WalkabilityTerm_Buffer1());
//alternative.AddUtilityTerm(263, destinationParcel.PCA_MixedUseTerm_Buffer1());
//alternative.AddUtilityTerm(264, destinationParcel.PCA_TransitAccessTerm_Buffer1());
//alternative.AddUtilityTerm(265, destinationParcel.PCA_DensityTerm_Buffer1());
//alternative.AddUtilityTerm(266, destinationParcel.PCA_WalkabilityTerm_Buffer1());
//alternative.AddUtilityTerm(267, destinationParcel.PCA_MixedUseTerm_Buffer1());
//alternative.AddUtilityTerm(268, destinationParcel.PCA_TransitAccessTerm_Buffer1());
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(70, 1); //for calibration
alternative.AddUtilityTerm(71, maleFlag);
alternative.AddUtilityTerm(73, ageBetween51And98Flag);
alternative.AddUtilityTerm(75, pathTypeModel.PathDistance * pathTypeModel.PathDistance);
alternative.AddUtilityTerm(77, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(78, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(179, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(178, destinationParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(177, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(176, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(175, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(179, originParcel.MixedUse4Index1());
alternative.AddUtilityTerm(181, 1.0 * destinationParcel.MixedUse4Index1()
+ 0.00001 * destinationParcel.TotalEmploymentDensity1()
+ 0.001 * destinationParcel.NetIntersectionDensity1()
+ 0.001 * originParcel.NetIntersectionDensity1()
+ 0.0001 * originParcel.HouseholdDensity1()
+ 1.0 * originParcel.MixedUse4Index1());
//alternative.AddUtilityTerm(271, originParcel.PCA_DensityTerm_Buffer1());
//alternative.AddUtilityTerm(272, originParcel.PCA_WalkabilityTerm_Buffer1());
//alternative.AddUtilityTerm(273, originParcel.PCA_MixedUseTerm_Buffer1());
//alternative.AddUtilityTerm(274, originParcel.PCA_TransitAccessTerm_Buffer1());
//alternative.AddUtilityTerm(271, destinationParcel.PCA_DensityTerm_Buffer1());
//alternative.AddUtilityTerm(272, destinationParcel.PCA_WalkabilityTerm_Buffer1());
//alternative.AddUtilityTerm(273, destinationParcel.PCA_MixedUseTerm_Buffer1());
//alternative.AddUtilityTerm(274, destinationParcel.PCA_TransitAccessTerm_Buffer1());
//alternative.AddUtilityTerm(275, destinationParcel.PCA_DensityTerm_Buffer1());
//alternative.AddUtilityTerm(276, destinationParcel.PCA_WalkabilityTerm_Buffer1());
//alternative.AddUtilityTerm(277, destinationParcel.PCA_MixedUseTerm_Buffer1());
//alternative.AddUtilityTerm(278, destinationParcel.PCA_TransitAccessTerm_Buffer1());
}
else if (mode == Global.Settings.Modes.PaidRideShare)
{
if (Global.Configuration.PaidRideshare_UseEstimatedInsteadOfAssertedCoefficients)
{
alternative.AddUtilityTerm(80, 1.0);
alternative.AddUtilityTerm(81, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(82, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(83, tour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(83, tour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(84, (tour.Person.Age >= 65).ToFlag()); //alternative.AddUtilityTerm(81, tour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(85, originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2);
alternative.AddUtilityTerm(86, destinationParcel.HouseholdsBuffer2 + destinationParcel.StudentsUniversityBuffer2 + destinationParcel.EmploymentTotalBuffer2);
alternative.AddUtilityTerm(87, income0To25KFlag);
alternative.AddUtilityTerm(88, incomeOver100Flag);
}
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());
}
}
RegionSpecificCustomizations(alternative, tour, pathTypeModel.PathType, mode, destinationParcel);
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper subtour, IEnumerable <IPathTypeModel> pathTypeModels, IParcelWrapper destinationParcel, int parentTourMode, int choice = Constants.DEFAULT_VALUE)
{
IHouseholdWrapper household = subtour.Household;
IPersonWrapper person = subtour.Person;
IPersonDayWrapper personDay = subtour.PersonDay;
// household inputs
int income0To25KFlag = household.Has0To25KIncome.ToFlag();
int income25To50KFlag = household.Has25To50KIncome.ToFlag();
int incomeOver100Flag = household.Has100KPlusIncome.ToFlag();
// person inputs
int maleFlag = person.IsMale.ToFlag();
// tour inputs
int sovTourFlag = (parentTourMode == Global.Settings.Modes.Sov).ToFlag();
int hov2TourFlag = (parentTourMode == Global.Settings.Modes.Hov2).ToFlag();
int bikeTourFlag = (parentTourMode == Global.Settings.Modes.Bike).ToFlag();
int walkTourFlag = (parentTourMode == Global.Settings.Modes.Walk).ToFlag();
int tncTourFlag = (parentTourMode == Global.Settings.Modes.PaidRideShare).ToFlag();
// remaining inputs
IParcelWrapper originParcel = subtour.OriginParcel;
int parkingDuration = ChoiceModelUtility.GetParkingDuration(person.IsFulltimeWorker);
double destinationParkingCost = destinationParcel.ParkingCostBuffer1(parkingDuration);
ChoiceModelUtility.SetEscortPercentages(personDay, out double escortPercentage, out double nonEscortPercentage, true);
foreach (IPathTypeModel pathTypeModel in pathTypeModels)
{
int mode = pathTypeModel.Mode;
bool available = (pathTypeModel.Mode != Global.Settings.Modes.PaidRideShare || Global.Configuration.PaidRideShareModeIsAvailable) &&
pathTypeModel.Available;
double generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
alternative.AddNestedAlternative(_nestedAlternativeIds[pathTypeModel.Mode], _nestedAlternativeIndexes[pathTypeModel.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) ?
subtour.TimeCoefficient * (1.0 - Global.Configuration.AV_InVehicleTimeCoefficientDiscountFactor) :
(mode == Global.Settings.Modes.PaidRideShare && Global.Configuration.AV_PaidRideShareModeUsesAVs) ?
subtour.TimeCoefficient * (1.0 - Global.Configuration.AV_InVehicleTimeCoefficientDiscountFactor) : subtour.TimeCoefficient;
alternative.AddUtilityTerm(2, generalizedTimeLogsum * modeTimeCoefficient);
if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(20, 1);
// alternative.AddUtility(129, destinationParcel.MixedUse2Index1());
// alternative.AddUtility(128, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(127, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(126, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(125, originParcel.TotalEmploymentDensity1());
// alternative.AddUtility(124, originParcel.MixedUse2Index1());
// alternative.AddUtility(123, Math.Log(destinationParcel.StopsTransitBuffer1+1));
// alternative.AddUtility(122, Math.Log(originParcel.StopsTransitBuffer1+1));
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * subtour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(88, sovTourFlag);
alternative.AddUtilityTerm(89, hov2TourFlag);
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * subtour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(88, sovTourFlag);
alternative.AddUtilityTerm(89, hov2TourFlag);
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * subtour.CostCoefficient));
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(55, income25To50KFlag);
alternative.AddUtilityTerm(58, sovTourFlag);
alternative.AddUtilityTerm(59, hov2TourFlag);
}
else if (mode == Global.Settings.Modes.Bike)
{
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(69, bikeTourFlag);
alternative.AddUtilityTerm(169, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(168, destinationParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(167, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(166, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(165, originParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(164, originParcel.MixedUse4Index1());
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(70, 1);
alternative.AddUtilityTerm(79, walkTourFlag);
alternative.AddUtilityTerm(179, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(178, destinationParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(177, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(176, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(175, originParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(174, originParcel.MixedUse4Index1());
}
else if (mode == Global.Settings.Modes.PaidRideShare)
{
if (Global.Configuration.PaidRideshare_UseEstimatedInsteadOfAssertedCoefficients)
{
alternative.AddUtilityTerm(80, 1.0);
alternative.AddUtilityTerm(81, sovTourFlag);
alternative.AddUtilityTerm(82, walkTourFlag + bikeTourFlag);
alternative.AddUtilityTerm(83, subtour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(83, subtour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(84, tncTourFlag);
alternative.AddUtilityTerm(85, originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2);
alternative.AddUtilityTerm(86, income0To25KFlag);
alternative.AddUtilityTerm(87, incomeOver100Flag);
}
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 * subtour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_Age18to25Coefficient * subtour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(90, Global.Configuration.PaidRideShare_AgeOver65Coefficient * (subtour.Person.Age >= 65).ToFlag());
}
}
}
}
ITripWrapper ITripCreator.CreateWrapper(ITourWrapper tourWrapper, int nextTripId, int direction, int sequence, bool isToTourOrigin, IHalfTour halfTour)
{
return(CreateWrapper(tourWrapper, nextTripId, direction, sequence, isToTourOrigin, halfTour));
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, TripWrapper trip,
IEnumerable <IPathTypeModel> pathTypeModels, IParcelWrapper originParcel,
IParcelWrapper destinationParcel,
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 HHwithChildrenFlag = household.HasChildren.ToFlag();
int HHwithSmallChildrenFlag = household.HasChildrenUnder5.ToFlag();
int childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
int HHwithLowIncomeFlag = (household.Income >= 300000 && household.Income < 600000).ToFlag();
int HHwithMidleIncomeFlag = (household.Income >= 600000 && household.Income < 900000).ToFlag();
int HHwithHighIncomeFlag = (household.Income >= 900000).ToFlag();
int nonworkingAdults = householdTotals.NonworkingAdults;
int retiredAdults = householdTotals.RetiredAdults;
int twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
int noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(household.VehiclesAvailable);
int carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(household.VehiclesAvailable);
// person inputs
//var drivingAgeStudentFlag = person.IsDrivingAgeStudent.ToFlag();
int maleFlag = person.IsMale.ToFlag();
int femaleFlag = person.IsFemale.ToFlag();
int PTpass = person.TransitPassOwnership;
//var ageLessThan35Flag = person.AgeIsLessThan35.ToFlag();
// tour inputs
int bikeTourFlag = tour.IsBikeMode().ToFlag();
int walkTourFlag = tour.IsWalkMode().ToFlag();
int CarDrivAloneFlag = tour.IsSovMode().ToFlag();
int CarDrivNotAloneFlag = tour.IsHov2Mode().ToFlag();
int CarPassengerFlag = tour.IsHov3Mode().ToFlag();
int transitTourFlag = tour.IsTransitMode().ToFlag();
int paidRideShareTourFlag = ((IActumTourWrapper)tour).IsPaidRideShareMode().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 homeBasedSocialTourFlag = (tour.IsHomeBasedTour && tour.IsSocialPurpose()).ToFlag();
int notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
int homeBasedNotWorkSchoolEscortTourFlag =
(tour.IsHomeBasedTour && tour.DestinationPurpose > Global.Settings.Purposes.Escort).ToFlag();
int jointTourFlag = (tour.JointTourSequence > 0) ? 1 : 0;
int partialHalfTourFlag = (trip.IsHalfTourFromOrigin
? tour.PartialHalfTour1Sequence > 0
: tour.PartialHalfTour2Sequence > 0)
? 1
: 0;
int fullHalfTourFlag = (trip.IsHalfTourFromOrigin ? tour.FullHalfTour1Sequence > 0 : tour.FullHalfTour2Sequence > 0)
? 1
: 0;
// 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
int departureTime = trip.IsHalfTourFromOrigin ? trip.LatestDepartureTime : trip.EarliestDepartureTime;
double originMixedDensity = originParcel.MixedUse4Index1();
double originIntersectionDensity = originParcel.NetIntersectionDensity1();
double destinationParkingCost = destinationParcel.ParkingCostBuffer1(2);
int amPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.SixAM, Global.Settings.Times.NineAM).ToFlag();
//GV changed to 6-9 am
int middayPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.NineAM, Global.Settings.Times.ThreePM).ToFlag();
int pmPeriodFlag = departureTime.IsLeftExclusiveBetween(Global.Settings.Times.ThreePM, Global.Settings.Times.SixPM).ToFlag();
//GV changed to 3-6 pm
int eveningPeriodFlag = (departureTime > Global.Settings.Times.SixPM).ToFlag(); //GV changed to 6 pm
// availability
bool[] tripModeAvailable = new bool[Global.Settings.Modes.TotalModes];
bool isLastTripInTour = (!trip.IsHalfTourFromOrigin && halfTour.SimulatedTrips >= 1 && trip.IsToTourOrigin);
int frequencyPreviousTripModeIsTourMode = 0;
if (trip.IsHalfTourFromOrigin)
{
frequencyPreviousTripModeIsTourMode +=
tour.HalfTourFromOrigin.Trips.Where(x => x.Sequence < trip.Sequence).Count(x => tour.Mode == x.Mode);
}
else
{
if (tour.HalfTourFromOrigin != null)
{
frequencyPreviousTripModeIsTourMode +=
tour.HalfTourFromOrigin.Trips.Where(x => x.Sequence > 0).Count(x => tour.Mode == x.Mode);
}
frequencyPreviousTripModeIsTourMode +=
tour.HalfTourFromDestination.Trips.Where(x => x.Sequence < trip.Sequence).Count(x => tour.Mode == x.Mode);
}
// GV commented out park and ride for COMPAS1; JOHN restored it for COMPAS2
//// 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 && person.Age >= 18;
tripModeAvailable[Global.Settings.Modes.HovDriver] = tour.Household.VehiclesAvailable > 0 && person.Age >= 18;
tripModeAvailable[Global.Settings.Modes.HovPassenger] = !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 (int mode = Global.Settings.Modes.Walk; mode <= tour.Mode; mode++)
{
tripModeAvailable[mode] = true;
}
}
if (person.Age < 18)
{
tripModeAvailable[Global.Settings.Modes.Sov] = false;
tripModeAvailable[Global.Settings.Modes.HovDriver] = false;
}
// GV commented out School Bus
// 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);
//const int modeExtra = Global.Settings.Modes.SchoolBus;
//var availableExtra = pathTypeExtra.Available && tour.IsSchoolBusMode && tripModeAvailable[modeExtra]
// && (trip.IsHalfTourFromOrigin
// ? trip.LatestDepartureTime - pathTypeExtra.PathTime >= trip.ArrivalTimeLimit
// : trip.EarliestDepartureTime + pathTypeExtra.PathTime <= trip.ArrivalTimeLimit);
//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);
//}
ChoiceProbabilityCalculator.Alternative alternative;
foreach (IPathTypeModel pathTypeModel in pathTypeModels)
{
int mode = pathTypeModel.Mode;
bool available = pathTypeModel.Available && tripModeAvailable[mode] &&
(trip.IsHalfTourFromOrigin
? trip.LatestDepartureTime - pathTypeModel.PathTime >= trip.ArrivalTimeLimit
: trip.EarliestDepartureTime + pathTypeModel.PathTime <= trip.ArrivalTimeLimit);
double generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
alternative = choiceProbabilityCalculator.GetAlternative(mode - 1, available, choice == mode);
alternative.Choice = mode;
// GV for tree Logit uncomment this
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);
//GV: income is not sign. 22. june 2016
//alternative.AddUtilityTerm(54, HHwithLowIncomeFlag);
//alternative.AddUtilityTerm(55, HHwithMidleIncomeFlag);
//alternative.AddUtilityTerm(56, HHwithHighIncomeFlag);
alternative.AddUtilityTerm(57, PTpass);
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(181, jointTourFlag);
alternative.AddUtilityTerm(182, fullHalfTourFlag + partialHalfTourFlag);
}
else if (mode == Global.Settings.Modes.HovPassenger)
{
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(33, femaleFlag);
alternative.AddUtilityTerm(34, (nonworkingAdults + retiredAdults));
alternative.AddUtilityTerm(36, onePersonHouseholdFlag);
alternative.AddUtilityTerm(37, twoPersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(100, CarPassengerFlag);
alternative.AddUtilityTerm(102, (CarPassengerFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (CarPassengerFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (CarPassengerFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (CarPassengerFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (CarPassengerFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (CarPassengerFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(115, transitTourFlag);
alternative.AddUtilityTerm(149, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(150, homeBasedSchoolTourFlag);
alternative.AddUtilityTerm(152, homeBasedEscortTourFlag);
alternative.AddUtilityTerm(153, homeBasedShoppingTourFlag);
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));
alternative.AddUtilityTerm(183, jointTourFlag);
alternative.AddUtilityTerm(184, fullHalfTourFlag + partialHalfTourFlag);
}
else if (mode == Global.Settings.Modes.HovDriver)
{
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(42, (childrenAge5Through15 * (1 - homeBasedEscortTourFlag)));
alternative.AddUtilityTerm(43, (femaleFlag * (1 - homeBasedEscortTourFlag)));
alternative.AddUtilityTerm(44, ((nonworkingAdults + retiredAdults) * (1 - homeBasedEscortTourFlag)));
alternative.AddUtilityTerm(38, onePersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(100, CarDrivNotAloneFlag);
alternative.AddUtilityTerm(102, (CarDrivNotAloneFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (CarDrivNotAloneFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (CarDrivNotAloneFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (CarDrivNotAloneFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (CarDrivNotAloneFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (CarDrivNotAloneFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(118, transitTourFlag);
alternative.AddUtilityTerm(120, CarPassengerFlag);
alternative.AddUtilityTerm(149, homeBasedWorkTourFlag);
alternative.AddUtilityTerm(150, homeBasedSchoolTourFlag);
alternative.AddUtilityTerm(152, homeBasedEscortTourFlag);
alternative.AddUtilityTerm(153, homeBasedShoppingTourFlag);
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(100, CarDrivAloneFlag);
alternative.AddUtilityTerm(102, (CarDrivAloneFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (CarDrivAloneFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (CarDrivAloneFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (CarDrivAloneFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (CarDrivAloneFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (CarDrivAloneFlag * lastTripOnSecondHalfFlag));
alternative.AddUtilityTerm(122, transitTourFlag);
alternative.AddUtilityTerm(124, CarPassengerFlag);
alternative.AddUtilityTerm(125, CarDrivNotAloneFlag);
alternative.AddUtilityTerm(126, maleFlag);
alternative.AddUtilityTerm(185, jointTourFlag);
alternative.AddUtilityTerm(186, fullHalfTourFlag + partialHalfTourFlag);
}
else if (mode == Global.Settings.Modes.Bike)
{
alternative.AddUtilityTerm(60, 1);
//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(130, CarDrivNotAloneFlag);
alternative.AddUtilityTerm(131, CarDrivAloneFlag);
alternative.AddUtilityTerm(132, maleFlag);
alternative.AddUtilityTerm(147, notHomeBasedTourFlag);
alternative.AddUtilityTerm(187, jointTourFlag);
alternative.AddUtilityTerm(188, fullHalfTourFlag + partialHalfTourFlag);
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(70, 1);
//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);
alternative.AddUtilityTerm(187, jointTourFlag);
alternative.AddUtilityTerm(188, fullHalfTourFlag + partialHalfTourFlag);
}
else if (mode == Global.Settings.Modes.PaidRideShare)
{
double modeConstant = Global.Configuration.AV_PaidRideShareModeUsesAVs
? Global.Configuration.AV_PaidRideShare_ModeConstant
+ Global.Configuration.AV_PaidRideShare_DensityCoefficient * Math.Min(originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2, 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, 6000);
alternative.AddUtilityTerm(80, modeConstant);
alternative.AddUtilityTerm(80, Global.Configuration.PaidRideShare_Age26to35Coefficient * tour.Person.AgeIsBetween26And35.ToFlag());
alternative.AddUtilityTerm(80, Global.Configuration.PaidRideShare_Age18to25Coefficient * tour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(80, Global.Configuration.PaidRideShare_AgeOver65Coefficient * (tour.Person.Age >= 65).ToFlag());
alternative.AddUtilityTerm(100, paidRideShareTourFlag);
alternative.AddUtilityTerm(102, (paidRideShareTourFlag * onlyTripOnFirstHalfFlag));
alternative.AddUtilityTerm(103, (paidRideShareTourFlag * onlyTripOnSecondHalfFlag));
alternative.AddUtilityTerm(104, (paidRideShareTourFlag * firstTripOnFirstHalfFlag));
alternative.AddUtilityTerm(105, (paidRideShareTourFlag * firstTripOnSecondHalfFlag));
alternative.AddUtilityTerm(106, (paidRideShareTourFlag * lastTripOnFirstHalfFlag));
alternative.AddUtilityTerm(107, (paidRideShareTourFlag * lastTripOnSecondHalfFlag));
}
}
}
ITourWrapper ITourCreator.CreateWrapper(ITour subtour, ITourWrapper tourWrapper, int purpose, bool suppressRandomVOT)
{
return(CreateWrapper(subtour, tourWrapper, purpose, suppressRandomVOT));
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, IEnumerable <IPathTypeModel> pathTypeModels, IParcelWrapper destinationParcel, int choice = Constants.DEFAULT_VALUE)
{
IHouseholdWrapper household = tour.Household;
Framework.DomainModels.Models.IHouseholdTotals householdTotals = household.HouseholdTotals;
IPersonWrapper person = tour.Person;
// household inputs
int childrenUnder5 = householdTotals.ChildrenUnder5;
int childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
int drivingAgeStudents = householdTotals.DrivingAgeStudents;
int noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(household.VehiclesAvailable);
int carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(household.VehiclesAvailable);
// person inputs
int ageBetween51And98Flag = person.AgeIsBetween51And98.ToFlag();
// other inputs
// in estimation mode, add an unavailable alternative 0 so Alogit will estimate it
if (Global.Configuration.IsInEstimationMode)
{
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(0, false, choice == 0);
alternative.Choice = 0;
alternative.AddUtilityTerm(30, 0);
}
foreach (IPathTypeModel pathTypeModel in pathTypeModels)
{
int mode = pathTypeModel.Mode;
bool available = (mode <= Global.Settings.Modes.Hov3 && mode >= Global.Settings.Modes.Walk) && pathTypeModel.Available;
double generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
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) : tour.TimeCoefficient;
if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(2, generalizedTimeLogsum * modeTimeCoefficient);
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(31, childrenUnder5);
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(33, drivingAgeStudents);
alternative.AddUtilityTerm(36, noCarsInHouseholdFlag); // for calibration of hov3 vs hov2
alternative.AddUtilityTerm(37, carsLessThanDriversFlag); // for calibration of hov3 vs hov2
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(2, generalizedTimeLogsum * modeTimeCoefficient);
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(42, carsLessThanDriversFlag);
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(52, carsLessThanDriversFlag); // for calibration
}
else if (mode == Global.Settings.Modes.Bike)
{
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(67, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(68, carsLessThanDriversFlag); //for calibration
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(2, generalizedTimeLogsum * modeTimeCoefficient);
alternative.AddUtilityTerm(70, 1);
alternative.AddUtilityTerm(73, ageBetween51And98Flag);
alternative.AddUtilityTerm(76, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(77, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(78, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(81, childrenUnder5);
alternative.AddUtilityTerm(82, childrenAge5Through15);
alternative.AddUtilityTerm(83, drivingAgeStudents);
}
}
}
public void Run(ITourWrapper tour, IHouseholdDayWrapper householdDay, int sampleSize, IParcelWrapper constrainedParcel = null)
{
if (tour == null)
{
throw new ArgumentNullException("tour");
}
tour.PersonDay.ResetRandom(20 + tour.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
if (!TourDestinationUtilities.ShouldRunInEstimationModeForModel(tour))
{
return;
}
// JLB 20140421 add the following to keep from estimatign twice for the same tour
if (tour.DestinationModeAndTimeHaveBeenSimulated)
{
return;
}
}
else if (tour.DestinationPurpose == Global.Settings.Purposes.School)
{
// the following lines were redundant. Usual destination properties are set in GetMandatoryTourSimulatedData();
// sets the destination for the school tour
//tour.DestinationParcelId = tour.Person.UsualSchoolParcelId;
//tour.DestinationParcel = tour.Person.UsualSchoolParcel;
//tour.DestinationZoneKey = tour.Person.UsualSchoolZoneKey;
//tour.DestinationAddressType = Global.Settings.AddressTypes.UsualSchool;
return;
}
else if (tour.DestinationPurpose == Global.Settings.Purposes.Work &&
tour.DestinationAddressType == Global.Settings.AddressTypes.UsualWorkplace)
{
return;
}
else if (constrainedParcel != null)
{
tour.DestinationParcel = constrainedParcel;
tour.DestinationParcelId = constrainedParcel.Id;
tour.DestinationZoneKey = ChoiceModelFactory.ZoneKeys[constrainedParcel.ZoneId]; // mb fixed
tour.DestinationAddressType = Global.Settings.AddressTypes.Other;
return;
}
var choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(tour.Id);
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
if (constrainedParcel == null)
{
RunModel(choiceProbabilityCalculator, tour, householdDay, sampleSize, tour.DestinationParcel);
choiceProbabilityCalculator.WriteObservation();
}
}
else
{
RunModel(choiceProbabilityCalculator, tour, householdDay, sampleSize);
var chosenAlternative = choiceProbabilityCalculator.SimulateChoice(tour.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", tour.PersonDay.Id);
if (!Global.Configuration.IsInEstimationMode)
{
tour.PersonDay.IsValid = false;
tour.PersonDay.HouseholdDay.IsValid = false;
}
return;
}
var choice = (ParcelWrapper)chosenAlternative.Choice;
tour.DestinationParcelId = choice.Id;
tour.DestinationParcel = choice;
tour.DestinationZoneKey = ChoiceModelFactory.ZoneKeys[choice.ZoneId];
tour.DestinationAddressType = Global.Settings.AddressTypes.Other;
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, IEnumerable <IPathTypeModel> pathTypeModels, IParcelWrapper destinationParcel, int choice = Constants.DEFAULT_VALUE)
{
IHouseholdWrapper household = tour.Household;
Framework.DomainModels.Models.IHouseholdTotals householdTotals = household.HouseholdTotals;
IPersonWrapper person = tour.Person;
IPersonDayWrapper personDay = tour.PersonDay;
// household inputs
int income0To25KFlag = household.Has0To25KIncome.ToFlag();
int incomeOver100Flag = household.Has100KPlusIncome.ToFlag();
int childrenUnder5 = householdTotals.ChildrenUnder5;
int childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
int nonworkingAdults = householdTotals.NonworkingAdults;
int retiredAdults = householdTotals.RetiredAdults;
int onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
int twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
int noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(household.VehiclesAvailable);
int carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(household.VehiclesAvailable);
int carsLessThanWorkersFlag = household.GetFlagForCarsLessThanWorkers(household.VehiclesAvailable);
// person inputs
int maleFlag = person.IsMale.ToFlag();
int ageBetween51And98Flag = person.AgeIsBetween51And98.ToFlag();
int univStudentFlag = person.IsUniversityStudent.ToFlag();
// tour inputs
int shoppingTourFlag = (tour.DestinationPurpose == Global.Settings.Purposes.Shopping).ToFlag();
int mealTourFlag = (tour.DestinationPurpose == Global.Settings.Purposes.Meal).ToFlag();
int socialOrRecreationTourFlag = (tour.DestinationPurpose == Global.Settings.Purposes.Social).ToFlag();
// remaining inputs
IParcelWrapper originParcel = tour.OriginParcel;
int parkingDuration = ChoiceModelUtility.GetParkingDuration(person.IsFulltimeWorker);
double destinationParkingCost = destinationParcel.ParkingCostBuffer1(parkingDuration);
ChoiceModelUtility.SetEscortPercentages(personDay, out double escortPercentage, out double nonEscortPercentage, true);
//foreach (var pathTypeModel in pathTypeModels) {
foreach (IPathTypeModel pathTypeModel in pathTypeModels)
{
IPathTypeModel ipathTypeModel = pathTypeModel;
int mode = pathTypeModel.Mode;
bool available = (pathTypeModel.Mode != Global.Settings.Modes.ParkAndRide || Global.Configuration.IncludeParkAndRideInOtherHomeBasedTourModeModel) &&
(pathTypeModel.Mode != Global.Settings.Modes.PaidRideShare || Global.Configuration.PaidRideShareModeIsAvailable) &&
pathTypeModel.Available;
double generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
alternative.AddNestedAlternative(_nestedAlternativeIds[pathTypeModel.Mode], _nestedAlternativeIndexes[pathTypeModel.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.ParkAndRide)
{
alternative.AddUtilityTerm(5, 1);
alternative.AddUtilityTerm(6, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(7, carsLessThanWorkersFlag);
alternative.AddUtilityTerm(128, destinationParcel.TotalEmploymentDensity1());
}
else if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(20, 1);
alternative.AddUtilityTerm(21, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(22, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(120, shoppingTourFlag);
alternative.AddUtilityTerm(121, mealTourFlag);
// alternative.AddUtility(129, destinationParcel.MixedUse2Index1());
alternative.AddUtilityTerm(128, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(127, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(126, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(125, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(124, originParcel.MixedUse2Index1());
// alternative.AddUtility(123, Math.Log(destinationParcel.StopsTransitBuffer1+1));
// alternative.AddUtility(122, Math.Log(originParcel.StopsTransitBuffer1+1));
alternative.AddUtilityTerm(180, univStudentFlag);
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(31, childrenUnder5);
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
alternative.AddUtilityTerm(35, pathTypeModel.PathDistance.AlmostEquals(0) ? 0 : Math.Log(pathTypeModel.PathDistance));
alternative.AddUtilityTerm(36, noCarsInHouseholdFlag); // for calibration of hov3 vs hov2
alternative.AddUtilityTerm(37, carsLessThanWorkersFlag); // for calibration of hov3 vs hov2
alternative.AddUtilityTerm(38, onePersonHouseholdFlag);
alternative.AddUtilityTerm(39, twoPersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(43, carsLessThanWorkersFlag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
alternative.AddUtilityTerm(136, shoppingTourFlag);
alternative.AddUtilityTerm(137, mealTourFlag);
alternative.AddUtilityTerm(138, socialOrRecreationTourFlag);
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(31, childrenUnder5);
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
alternative.AddUtilityTerm(35, pathTypeModel.PathDistance.AlmostEquals(0) ? 0 : Math.Log(pathTypeModel.PathDistance));
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(43, carsLessThanWorkersFlag);
alternative.AddUtilityTerm(48, onePersonHouseholdFlag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
alternative.AddUtilityTerm(136, shoppingTourFlag);
alternative.AddUtilityTerm(137, mealTourFlag);
alternative.AddUtilityTerm(138, socialOrRecreationTourFlag);
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient));
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(52, carsLessThanDriversFlag);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(131, escortPercentage);
alternative.AddUtilityTerm(132, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Bike)
{
double class1Dist = Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions ?
ImpedanceRoster.GetValue("class1distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable : 0;
double class2Dist = Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions ?
ImpedanceRoster.GetValue("class2distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable : 0;
double worstDist = Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions ?
ImpedanceRoster.GetValue("worstdistance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable : 0;
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(63, ageBetween51And98Flag);
alternative.AddUtilityTerm(67, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(68, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(160, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(169, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(168, destinationParcel.TotalEmploymentDensity1());
alternative.AddUtilityTerm(167, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(166, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(165, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(164, originParcel.MixedUse4Index1());
alternative.AddUtilityTerm(161, (class1Dist > 0).ToFlag());
alternative.AddUtilityTerm(162, (class2Dist > 0).ToFlag());
alternative.AddUtilityTerm(163, (worstDist > 0).ToFlag());
alternative.AddUtilityTerm(170, 1.0 * destinationParcel.MixedUse4Index1()
+ 0.00002 * destinationParcel.TotalEmploymentDensity1()
+ 0.001 * destinationParcel.NetIntersectionDensity1()
+ 0.001 * originParcel.NetIntersectionDensity1()
+ 0.0002 * originParcel.HouseholdDensity1()
+ 1.0 * originParcel.MixedUse4Index1());
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(70, 1);
alternative.AddUtilityTerm(73, ageBetween51And98Flag);
alternative.AddUtilityTerm(77, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(78, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(171, mealTourFlag);
alternative.AddUtilityTerm(172, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(179, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(178, destinationParcel.HouseholdDensity1());
alternative.AddUtilityTerm(177, destinationParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(176, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(175, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(174, originParcel.MixedUse4Index1());
alternative.AddUtilityTerm(181, 1.0 * destinationParcel.MixedUse4Index1()
+ 0.00001 * destinationParcel.TotalEmploymentDensity1()
+ 0.001 * destinationParcel.NetIntersectionDensity1()
+ 0.001 * originParcel.NetIntersectionDensity1()
+ 0.0001 * originParcel.HouseholdDensity1()
+ 1.0 * originParcel.MixedUse4Index1());
}
else if (mode == Global.Settings.Modes.PaidRideShare)
{
if (Global.Configuration.PaidRideshare_UseEstimatedInsteadOfAssertedCoefficients)
{
alternative.AddUtilityTerm(80, 1.0);
//alternative.AddUtilityTerm(81, tour.Person.AgeIsBetween26And35.ToFlag());
//alternative.AddUtilityTerm(82, tour.Person.AgeIsBetween18And25.ToFlag());
alternative.AddUtilityTerm(81, noCarsInHouseholdFlag); //for calibration
alternative.AddUtilityTerm(82, carsLessThanDriversFlag); //for calibration
alternative.AddUtilityTerm(83, (tour.Person.Age >= 65).ToFlag());
alternative.AddUtilityTerm(84, originParcel.HouseholdsBuffer2 + originParcel.StudentsUniversityBuffer2 + originParcel.EmploymentTotalBuffer2);
alternative.AddUtilityTerm(85, destinationParcel.HouseholdsBuffer2 + destinationParcel.StudentsUniversityBuffer2 + destinationParcel.EmploymentTotalBuffer2);
alternative.AddUtilityTerm(86, income0To25KFlag);
alternative.AddUtilityTerm(87, incomeOver100Flag);
alternative.AddUtilityTerm(88, mealTourFlag);
alternative.AddUtilityTerm(89, shoppingTourFlag);
}
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());
}
}
RegionSpecificCustomizations(alternative, tour, pathTypeModel.PathType, mode, destinationParcel);
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, IHouseholdDayWrapper householdDay, int sampleSize, IParcelWrapper choice = null)
{
// var household = tour.Household;
var person = tour.Person;
var personDay = tour.PersonDay;
// var totalAvailableMinutes =
// tour.ParentTour == null
// ? personDay.TimeWindow.TotalAvailableMinutes(1, Global.Settings.Times.MinutesInADay)
// : tour.ParentTour.TimeWindow.TotalAvailableMinutes(1, Global.Settings.Times.MinutesInADay);
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 (tour.ParentTour == null)
{
timeWindow.IncorporateAnotherTimeWindow(personDay.TimeWindow);
}
timeWindow.SetBusyMinutes(Global.Settings.Times.EndOfRelevantWindow, Global.Settings.Times.MinutesInADay + 1);
var maxAvailableMinutes =
(tour.JointTourSequence > 0 || tour.ParentTour == null)
? timeWindow.MaxAvailableMinutesAfter(Global.Settings.Times.FiveAM)
: tour.ParentTour.DestinationDepartureTime - tour.ParentTour.DestinationArrivalTime;
// var hoursAvailableInverse =
// tour.IsHomeBasedTour
// ? (personDay.HomeBasedTours - personDay.SimulatedHomeBasedTours + 1) / (Math.Max(totalAvailableMinutes - 360, 30) / 60D)
// : 1 / (Math.Max(totalAvailableMinutes, 1) / 60D);
var fastestAvailableTimeOfDay =
tour.IsHomeBasedTour || tour.ParentTour == null
? 1
: tour.ParentTour.DestinationArrivalTime + (tour.ParentTour.DestinationDepartureTime - tour.ParentTour.DestinationArrivalTime) / 2;
var tourCategory = tour.GetTourCategory();
// var primaryFlag = ChoiceModelUtility.GetPrimaryFlag(tourCategory);
var secondaryFlag = ChoiceModelUtility.GetSecondaryFlag(tourCategory);
ChoiceModelUtility.DrawRandomTourTimePeriods(tour, tourCategory);
var segment = Global.ContainerDaySim.GetInstance <SamplingWeightsSettingsFactory>().SamplingWeightsSettings.GetTourDestinationSegment(tour.DestinationPurpose, tour.IsHomeBasedTour ? Global.Settings.TourPriorities.HomeBasedTour : Global.Settings.TourPriorities.WorkBasedTour, Global.Settings.Modes.Sov, person.PersonType);
var destinationSampler = new DestinationSampler(choiceProbabilityCalculator, segment, sampleSize, choice, tour.OriginParcel);
var tourDestinationUtilities = new TourDestinationUtilities(tour, sampleSize, secondaryFlag, personDay.GetIsWorkOrSchoolPattern().ToFlag(), personDay.GetIsOtherPattern().ToFlag(), fastestAvailableTimeOfDay, maxAvailableMinutes);
destinationSampler.SampleTourDestinations(tourDestinationUtilities);
}
protected override void RegionSpecificOtherTourDistrictCoefficients(ChoiceProbabilityCalculator.Alternative alternative, ITourWrapper _tour, IParcelWrapper destinationParcel)
{
//add any region-specific new terms in region-specific class, using coefficient numbers 114-120, or other unused variable #
//Global.PrintFile.WriteLine("Fresno_OtherTourDestinationModel.RegionSpecificOtherTourDistrictCoefficients called");
int origdist = _tour.OriginParcel.District;
//var destZone = destinationParcel.ZoneId;
int destParcel = destinationParcel.Id;
//var destRiverPark = (destZone == 514 || destZone == 515 || destZone == 516 || destZone == 1486 || destZone == 1487) ? 1 : 0;
int destRiverPark = (destParcel == 19599 || destParcel == 19600 || destParcel == 19601 || destParcel == 19598 || destParcel == 12156 || destParcel == 12158 || destParcel == 19603) ? 1 : 0;
int origDistDestRiverPark = (origdist == 4 || origdist == 7 || origdist == 10) && (destParcel == 19599 || destParcel == 19600 || destParcel == 19601 || destParcel == 19598 || destParcel == 12156 || destParcel == 12158 || destParcel == 19603) ? 1 : 0;
int origEastDestRiverPark = (origdist == 5) && (destParcel == 19599 || destParcel == 19600 || destParcel == 19601 || destParcel == 19598 || destParcel == 12156 || destParcel == 12158 || destParcel == 19603) ? 1 : 0;
if (_tour.DestinationPurpose == Global.Settings.Purposes.Shopping || _tour.DestinationPurpose == Global.Settings.Purposes.Meal)
{
//Global.PrintFile.WriteLine("Destination is River Park Mall");
alternative.AddUtilityTerm(115, destRiverPark);
alternative.AddUtilityTerm(116, origDistDestRiverPark);
alternative.AddUtilityTerm(117, origEastDestRiverPark);
}
//alternative.AddUtilityTerm(115, origEastDestCBD);
//alternative.AddUtilityTerm(116, origKitDestTRP);
//alternative.AddUtilityTerm(117, origTacDestKit);
//alternative.AddUtilityTerm(118, origKitDestNotKit);
//alternative.AddUtilityTerm(119, origSTacWorkCBD);
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, IHouseholdDayWrapper householdDay, ITripWrapper trip, HTripTime 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();
int jointTourFlag = (tour.JointTourSequence > 0) ? 1 : 0;
int partialHalfTourFlag = (trip.IsHalfTourFromOrigin ? tour.PartialHalfTour1Sequence > 0 : tour.PartialHalfTour2Sequence > 0) ? 1 : 0;
int fullHalfTourFlag = (trip.IsHalfTourFromOrigin ? tour.FullHalfTour1Sequence > 0 : tour.FullHalfTour2Sequence > 0) ? 1 : 0;
// set trip inputs
bool originChangeMode = trip.Sequence > 1 && trip.GetPreviousTrip().DestinationPurpose == Global.Settings.Purposes.ChangeMode;
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
// set remaining inputs
TimeWindow timeWindow = new TimeWindow();
if (tour.JointTourSequence > 0)
{
foreach (PersonDayWrapper pDay in householdDay.PersonDays)
{
ITourWrapper 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)
{
ITourWrapper 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)
{
ITourWrapper 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);
int remainingToursCount = personDay.HomeBasedTours - personDay.GetTotalSimulatedTours();
int tripRemainingInHalfTour = (trip.DestinationParcel != null && trip.DestinationParcel != tour.OriginParcel).ToFlag(); // we don't know exact #
int previousArrivalTime = trip.IsHalfTourFromOrigin
? (trip.Sequence == 1 ? tour.DestinationDepartureTime : trip.GetPreviousTrip().ArrivalTime)
: (trip.Sequence == 1 ? tour.DestinationArrivalTime : trip.GetPreviousTrip().ArrivalTime);
MinuteSpan previousArrivalPeriod = new HTripTime(previousArrivalTime).DeparturePeriod;
foreach (HTripTime time in HTripTime.Times[ParallelUtility.threadLocalAssignedIndex.Value])
{
MinuteSpan period = time.DeparturePeriod;
double departurePeriodFraction = timeWindow.TotalAvailableMinutes(period.Start, period.End) / (period.End - period.Start + 1D);
double departureShift = period.Middle / 60.0;
double durationShift = Math.Abs(period.Middle - previousArrivalPeriod.Middle) / 60.0;
bool available = time.Available && departurePeriodFraction > 0;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(time.Index, available, choice != null && choice.Equals(time));
if (!alternative.Available)
{
continue;
}
alternative.Choice = time;
int indicatedTravelTime = (int)time.ModeLOS.PathTime;
int indicatedArrivalTime = trip.IsHalfTourFromOrigin
? Math.Max(1, period.Middle - indicatedTravelTime)
: Math.Min(1440, period.Middle + indicatedTravelTime);
int totalWindowRemaining = trip.IsHalfTourFromOrigin
? timeWindow.TotalAvailableMinutesBefore(indicatedArrivalTime) + timeWindow.TotalAvailableMinutesAfter(previousArrivalTime)
: timeWindow.TotalAvailableMinutesAfter(indicatedArrivalTime) + timeWindow.TotalAvailableMinutesBefore(previousArrivalTime);
int 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 void Run(ITourWrapper tour)
{
if (tour == null)
{
throw new ArgumentNullException("tour");
}
tour.PersonDay.ResetRandom(40 + tour.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
}
var choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(tour.Id);
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
if (tour.DestinationParcel == null || tour.Mode <= Global.Settings.Modes.None || (tour.Mode > Global.Settings.Modes.Transit && tour.Mode != Global.Settings.Modes.SchoolBus))
{
return;
}
var pathTypeModels =
PathTypeModel.RunAll(
tour.Household.RandomUtility,
tour.OriginParcel,
tour.DestinationParcel,
tour.DestinationArrivalTime,
tour.DestinationDepartureTime,
tour.DestinationPurpose,
tour.CostCoefficient,
tour.TimeCoefficient,
tour.Person.IsDrivingAge,
tour.Household.VehiclesAvailable,
tour.Person.GetTransitFareDiscountFraction(),
false);
var mode = (tour.Mode == Global.Settings.Modes.SchoolBus) ? Global.Settings.Modes.Hov3 : tour.Mode; // use HOV3 for school bus impedance
var pathTypeModel = pathTypeModels.First(x => x.Mode == mode);
if (!pathTypeModel.Available)
{
return;
}
RunModel(choiceProbabilityCalculator, tour, pathTypeModels, tour.DestinationParcel, tour.Household.VehiclesAvailable, tour.Mode);
choiceProbabilityCalculator.WriteObservation();
}
else
{
var pathTypeModels =
PathTypeModel.RunAll(
tour.Household.RandomUtility,
tour.OriginParcel,
tour.DestinationParcel,
tour.DestinationArrivalTime,
tour.DestinationDepartureTime,
tour.DestinationPurpose,
tour.CostCoefficient,
tour.TimeCoefficient,
tour.Person.IsDrivingAge,
tour.Household.VehiclesAvailable,
tour.Person.GetTransitFareDiscountFraction(),
false);
RunModel(choiceProbabilityCalculator, tour, pathTypeModels, tour.DestinationParcel, tour.Household.VehiclesAvailable);
var chosenAlternative = choiceProbabilityCalculator.SimulateChoice(tour.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", tour.PersonDay.Id);
tour.Mode = Global.Settings.Modes.Hov3;
if (!Global.Configuration.IsInEstimationMode)
{
tour.PersonDay.IsValid = false;
}
return;
}
var choice = (int)chosenAlternative.Choice;
tour.Mode = choice;
if (choice == Global.Settings.Modes.SchoolBus)
{
tour.PathType = 0;
}
else
{
var chosenPathType = pathTypeModels.First(x => x.Mode == choice);
tour.PathType = chosenPathType.PathType;
tour.ParkAndRideNodeId = choice == Global.Settings.Modes.ParkAndRide ? chosenPathType.PathParkAndRideNodeId : 0;
}
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, TourTime choice = null)
{
IHouseholdWrapper household = tour.Household;
IPersonWrapper person = tour.Person;
IPersonDayWrapper personDay = tour.PersonDay;
// household inputs
int income0To25KFlag = household.Has0To25KIncome.ToFlag();
int income100KPlusFlag = household.Has100KPlusIncome.ToFlag();
// 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 fulltimeWorkerFlag = person.IsFulltimeWorker.ToFlag();
int childAge5Through15Flag = person.IsChildAge5Through15.ToFlag();
int childUnder5Flag = person.IsChildUnder5.ToFlag();
// person-day inputs
int homeBasedToursOnlyFlag = personDay.OnlyHomeBasedToursExist().ToFlag();
int firstSimulatedHomeBasedTourFlag = personDay.IsFirstSimulatedHomeBasedTour().ToFlag();
int laterSimulatedHomeBasedTourFlag = personDay.IsLaterSimulatedHomeBasedTour().ToFlag();
int totalStopPurposes = personDay.GetTotalStopPurposes(); // JLB 20150401 change from Stops to StopPurposes to make definition consistent in apply and estimate
int totalSimulatedStops = personDay.GetTotalSimulatedStops();
int escortStops = personDay.EscortStops;
int homeBasedTours = personDay.HomeBasedTours;
int simulatedHomeBasedTours = personDay.SimulatedHomeBasedTours;
// tour inputs
int escortTourFlag = tour.IsEscortPurpose().ToFlag();
int shoppingTourFlag = tour.IsShoppingPurpose().ToFlag();
int mealTourFlag = tour.IsMealPurpose().ToFlag();
int socialTourFlag = tour.IsSocialPurpose().ToFlag();
int personalBusinessTourFlag = tour.IsPersonalBusinessPurpose().ToFlag();
int sovOrHovTourFlag = tour.IsAnAutoMode().ToFlag();
int transitTourFlag = tour.UsesTransitModes().ToFlag();
// remaining inputs
// Higher priority tour of 2+ tours for the same purpose
int highPrioritySameFlag = (tour.GetTotalToursByPurpose() > tour.GetTotalSimulatedToursByPurpose() && tour.GetTotalSimulatedToursByPurpose() == 1).ToFlag();
// Lower priority tour(s) of 2+ tours for the same purpose
int lowPrioritySameFlag = (tour.GetTotalSimulatedToursByPurpose() > 1).ToFlag();
// Higher priority tour of 2+ tours for different purposes
int highPriorityDifferentFlag = (personDay.IsFirstSimulatedHomeBasedTour() && personDay.HomeBasedToursExist()).ToFlag() * (1 - highPrioritySameFlag);
// Lower priority tour of 2+ tours for different purposes
int lowPriorityDifferentFlag = (personDay.IsLaterSimulatedHomeBasedTour() && personDay.HomeBasedToursExist()).ToFlag() * (1 - lowPrioritySameFlag);
ITimeWindow timeWindow = tour.ParentTour == null ? personDay.TimeWindow : tour.ParentTour.TimeWindow;
ITourModeImpedance[] impedances = tour.GetTourModeImpedances();
int period1Middle = 15;
int smallPeriodCount = DayPeriod.SmallDayPeriods.Count();
for (int periodIndex = 0; periodIndex < smallPeriodCount; periodIndex++)
{
// set arrival period component
int arrivalPeriodIndex = periodIndex;
double arrivalPeriodShift = arrivalPeriodIndex * (48.0 / DayPeriod.SMALL_DAY_PERIOD_TOTAL_TRIP_TIMES); //adjust shift amount if period lengths change
MinuteSpan arrivalPeriod = DayPeriod.SmallDayPeriods[arrivalPeriodIndex];
int earlyArriveFlag = arrivalPeriod.Middle.IsBetween(Global.Settings.Times.ThreeAM, Global.Settings.Times.SevenAM).ToFlag();
ITourModeImpedance arrivalImpedance = impedances[arrivalPeriodIndex];
int arrivalPeriodAvailableMinutes = timeWindow.TotalAvailableMinutes(arrivalPeriod.Start, arrivalPeriod.End);
choiceProbabilityCalculator.CreateUtilityComponent(3 * periodIndex + 0);
ChoiceProbabilityCalculator.Component arrivalComponent = choiceProbabilityCalculator.GetUtilityComponent(3 * periodIndex + 0);
if (arrivalPeriodAvailableMinutes > 0 || Global.Configuration.IsInEstimationMode)
{
arrivalComponent.AddUtilityTerm(11, arrivalPeriod.Middle.IsBetween(Global.Settings.Times.ThreeAM, Global.Settings.Times.SixAM).ToFlag());
arrivalComponent.AddUtilityTerm(12, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SixAM, Global.Settings.Times.SevenAM).ToFlag());
arrivalComponent.AddUtilityTerm(13, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenAM, Global.Settings.Times.EightAM).ToFlag());
arrivalComponent.AddUtilityTerm(14, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.EightAM, Global.Settings.Times.NineAM).ToFlag());
arrivalComponent.AddUtilityTerm(15, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.NineAM, Global.Settings.Times.TenAM).ToFlag());
arrivalComponent.AddUtilityTerm(16, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.OnePM).ToFlag());
arrivalComponent.AddUtilityTerm(17, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.OnePM, Global.Settings.Times.FourPM).ToFlag());
arrivalComponent.AddUtilityTerm(18, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FourPM, Global.Settings.Times.SevenPM).ToFlag());
arrivalComponent.AddUtilityTerm(19, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenPM, Global.Settings.Times.TenPM).ToFlag());
arrivalComponent.AddUtilityTerm(20, arrivalPeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenPM, Global.Settings.Times.MinutesInADay).ToFlag());
arrivalComponent.AddUtilityTerm(41, partTimeWorkerFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(43, nonworkingAdultFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(45, universityStudentFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(47, retiredAdultFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(49, drivingAgeStudentFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(139, fulltimeWorkerFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(141, childAge5Through15Flag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(143, childUnder5Flag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(145, escortTourFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(147, shoppingTourFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(149, mealTourFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(151, socialTourFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(153, personalBusinessTourFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(51, income0To25KFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(53, income100KPlusFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(55, highPrioritySameFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(57, lowPrioritySameFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(155, highPriorityDifferentFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(157, lowPriorityDifferentFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(59, homeBasedToursOnlyFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(61, totalStopPurposes * homeBasedToursOnlyFlag * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(63, totalStopPurposes * (1 - homeBasedToursOnlyFlag) * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(65, (totalStopPurposes - totalSimulatedStops) * (1 - homeBasedToursOnlyFlag) * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(67, escortStops * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(69, tour.TotalSubtours * arrivalPeriodShift);
arrivalComponent.AddUtilityTerm(72, income100KPlusFlag * earlyArriveFlag);
arrivalComponent.AddUtilityTerm(175, escortTourFlag * earlyArriveFlag);
arrivalComponent.AddUtilityTerm(176, shoppingTourFlag * earlyArriveFlag);
arrivalComponent.AddUtilityTerm(177, mealTourFlag * earlyArriveFlag);
arrivalComponent.AddUtilityTerm(85, sovOrHovTourFlag * arrivalImpedance.GeneralizedTimeFromOrigin * tour.TimeCoefficient);
arrivalComponent.AddUtilityTerm(87, transitTourFlag * Math.Max(0, arrivalImpedance.GeneralizedTimeFromOrigin) * tour.TimeCoefficient);
arrivalComponent.AddUtilityTerm(89, transitTourFlag * (arrivalImpedance.GeneralizedTimeFromOrigin < 0).ToFlag());
arrivalComponent.AddUtilityTerm(91, Math.Log(Math.Max(1, arrivalPeriodAvailableMinutes)));
arrivalComponent.AddUtilityTerm(93, arrivalImpedance.AdjacentMinutesBefore * firstSimulatedHomeBasedTourFlag);
arrivalComponent.AddUtilityTerm(95, arrivalImpedance.AdjacentMinutesBefore * laterSimulatedHomeBasedTourFlag);
arrivalComponent.AddUtilityTerm(99, (totalStopPurposes - totalSimulatedStops) / (1D + arrivalImpedance.AdjacentMinutesBefore));
}
// set departure period component
int departurePeriodIndex = periodIndex;
MinuteSpan departurePeriod = DayPeriod.SmallDayPeriods[departurePeriodIndex];
int lateDepartFlag = departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.NinePM, Global.Settings.Times.MinutesInADay).ToFlag();
ITourModeImpedance departureImpedance = impedances[departurePeriodIndex];
int departurePeriodAvailableMinutes = timeWindow.TotalAvailableMinutes(departurePeriod.Start, departurePeriod.End);
choiceProbabilityCalculator.CreateUtilityComponent(3 * periodIndex + 1);
ChoiceProbabilityCalculator.Component departureComponent = choiceProbabilityCalculator.GetUtilityComponent(3 * periodIndex + 1);
if (departurePeriodAvailableMinutes > 0 || Global.Configuration.IsInEstimationMode)
{
departureComponent.AddUtilityTerm(21, departurePeriod.Middle.IsBetween(Global.Settings.Times.ThreeAM, Global.Settings.Times.SevenAM).ToFlag());
departureComponent.AddUtilityTerm(22, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenAM, Global.Settings.Times.TenAM).ToFlag());
departureComponent.AddUtilityTerm(23, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.TenAM, Global.Settings.Times.OnePM).ToFlag());
departureComponent.AddUtilityTerm(24, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.OnePM, Global.Settings.Times.ThreePM).ToFlag());
departureComponent.AddUtilityTerm(124, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.ThreePM, Global.Settings.Times.FourPM).ToFlag());
departureComponent.AddUtilityTerm(25, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FourPM, Global.Settings.Times.FivePM).ToFlag());
departureComponent.AddUtilityTerm(26, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.FivePM, Global.Settings.Times.SixPM).ToFlag());
departureComponent.AddUtilityTerm(27, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SixPM, Global.Settings.Times.SevenPM).ToFlag());
departureComponent.AddUtilityTerm(28, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.SevenPM, Global.Settings.Times.NinePM).ToFlag());
departureComponent.AddUtilityTerm(29, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.NinePM, Global.Settings.Times.Midnight).ToFlag());
departureComponent.AddUtilityTerm(30, departurePeriod.Middle.IsLeftExclusiveBetween(Global.Settings.Times.Midnight, Global.Settings.Times.MinutesInADay).ToFlag());
departureComponent.AddUtilityTerm(73, income100KPlusFlag * lateDepartFlag);
departureComponent.AddUtilityTerm(178, escortTourFlag * lateDepartFlag);
departureComponent.AddUtilityTerm(179, shoppingTourFlag * lateDepartFlag);
departureComponent.AddUtilityTerm(180, mealTourFlag * lateDepartFlag);
departureComponent.AddUtilityTerm(86, sovOrHovTourFlag * departureImpedance.GeneralizedTimeFromDestination * tour.TimeCoefficient);
departureComponent.AddUtilityTerm(88, transitTourFlag * Math.Max(0, departureImpedance.GeneralizedTimeFromDestination) * tour.TimeCoefficient);
departureComponent.AddUtilityTerm(89, transitTourFlag * (departureImpedance.GeneralizedTimeFromDestination < 0).ToFlag());
departureComponent.AddUtilityTerm(92, Math.Log(Math.Max(1, departurePeriodAvailableMinutes)));
departureComponent.AddUtilityTerm(94, departureImpedance.AdjacentMinutesAfter * firstSimulatedHomeBasedTourFlag);
departureComponent.AddUtilityTerm(96, departureImpedance.AdjacentMinutesAfter * laterSimulatedHomeBasedTourFlag);
departureComponent.AddUtilityTerm(100, (totalStopPurposes - totalSimulatedStops) / (1D + departureImpedance.AdjacentMinutesAfter));
}
// set duration component (relative to period 1)
double periodSpan = periodIndex * (48.0 / DayPeriod.SMALL_DAY_PERIOD_TOTAL_TRIP_TIMES); //adjust shift amount if period lengths change;
if (arrivalPeriodIndex == 0)
{
period1Middle = arrivalPeriod.Middle;
}
int duration = departurePeriod.Middle - period1Middle;
int durationUnder1HourFlag = ChoiceModelUtility.GetDurationUnder1HourFlag(duration);
int duration1To2HoursFlag = ChoiceModelUtility.GetDuration1To2HoursFlag(duration);
int durationUnder4HoursFlag = ChoiceModelUtility.GetDurationUnder4HoursFlag(duration);
int durationUnder8HoursFlag = ChoiceModelUtility.GetDurationUnder8HoursFlag(duration);
int durationUnder9HoursFlag = ChoiceModelUtility.GetDurationUnder9HoursFlag(duration);
choiceProbabilityCalculator.CreateUtilityComponent(3 * periodIndex + 2);
ChoiceProbabilityCalculator.Component durationComponent = choiceProbabilityCalculator.GetUtilityComponent(3 * periodIndex + 2);
if (tour.IsWorkPurpose() || tour.IsSchoolPurpose())
{
durationComponent.AddUtilityTerm(31, duration.IsRightExclusiveBetween(Global.Settings.Times.ZeroHours, Global.Settings.Times.ThreeHours).ToFlag());
durationComponent.AddUtilityTerm(32, duration.IsRightExclusiveBetween(Global.Settings.Times.ThreeHours, Global.Settings.Times.FiveHours).ToFlag());
durationComponent.AddUtilityTerm(33, duration.IsRightExclusiveBetween(Global.Settings.Times.FiveHours, Global.Settings.Times.SevenHours).ToFlag());
durationComponent.AddUtilityTerm(34, duration.IsRightExclusiveBetween(Global.Settings.Times.SevenHours, Global.Settings.Times.NineHours).ToFlag());
durationComponent.AddUtilityTerm(35, duration.IsRightExclusiveBetween(Global.Settings.Times.NineHours, Global.Settings.Times.TenHours).ToFlag());
durationComponent.AddUtilityTerm(36, duration.IsRightExclusiveBetween(Global.Settings.Times.TenHours, Global.Settings.Times.ElevenHours).ToFlag());
durationComponent.AddUtilityTerm(37, duration.IsRightExclusiveBetween(Global.Settings.Times.ElevenHours, Global.Settings.Times.TwelveHours).ToFlag());
durationComponent.AddUtilityTerm(38, duration.IsRightExclusiveBetween(Global.Settings.Times.TwelveHours, Global.Settings.Times.FourteenHours).ToFlag());
durationComponent.AddUtilityTerm(39, duration.IsRightExclusiveBetween(Global.Settings.Times.FourteenHours, Global.Settings.Times.EighteenHours).ToFlag());
durationComponent.AddUtilityTerm(40, (duration >= Global.Settings.Times.EighteenHours).ToFlag());
}
else
{
durationComponent.AddUtilityTerm(31, duration.IsRightExclusiveBetween(Global.Settings.Times.ZeroHours, Global.Settings.Times.OneHour).ToFlag());
durationComponent.AddUtilityTerm(32, duration.IsRightExclusiveBetween(Global.Settings.Times.OneHour, Global.Settings.Times.TwoHours).ToFlag());
durationComponent.AddUtilityTerm(33, duration.IsRightExclusiveBetween(Global.Settings.Times.TwoHours, Global.Settings.Times.ThreeHours).ToFlag());
durationComponent.AddUtilityTerm(34, duration.IsRightExclusiveBetween(Global.Settings.Times.ThreeHours, Global.Settings.Times.FiveHours).ToFlag());
durationComponent.AddUtilityTerm(35, duration.IsRightExclusiveBetween(Global.Settings.Times.FiveHours, Global.Settings.Times.SevenHours).ToFlag());
durationComponent.AddUtilityTerm(36, duration.IsRightExclusiveBetween(Global.Settings.Times.SevenHours, Global.Settings.Times.NineHours).ToFlag());
durationComponent.AddUtilityTerm(37, duration.IsRightExclusiveBetween(Global.Settings.Times.NineHours, Global.Settings.Times.TwelveHours).ToFlag());
durationComponent.AddUtilityTerm(38, duration.IsRightExclusiveBetween(Global.Settings.Times.TwelveHours, Global.Settings.Times.FourteenHours).ToFlag());
durationComponent.AddUtilityTerm(39, duration.IsRightExclusiveBetween(Global.Settings.Times.FourteenHours, Global.Settings.Times.EighteenHours).ToFlag());
durationComponent.AddUtilityTerm(40, (duration >= Global.Settings.Times.EighteenHours).ToFlag());
}
durationComponent.AddUtilityTerm(42, partTimeWorkerFlag * periodSpan);
durationComponent.AddUtilityTerm(44, nonworkingAdultFlag * periodSpan);
durationComponent.AddUtilityTerm(46, universityStudentFlag * periodSpan);
durationComponent.AddUtilityTerm(48, retiredAdultFlag * periodSpan);
durationComponent.AddUtilityTerm(50, drivingAgeStudentFlag * periodSpan);
durationComponent.AddUtilityTerm(140, fulltimeWorkerFlag * periodSpan);
durationComponent.AddUtilityTerm(142, childAge5Through15Flag * periodSpan);
durationComponent.AddUtilityTerm(144, childUnder5Flag * periodSpan);
durationComponent.AddUtilityTerm(146, escortTourFlag * periodSpan);
durationComponent.AddUtilityTerm(148, shoppingTourFlag * periodSpan);
durationComponent.AddUtilityTerm(150, mealTourFlag * periodSpan);
durationComponent.AddUtilityTerm(152, socialTourFlag * periodSpan);
durationComponent.AddUtilityTerm(154, personalBusinessTourFlag * periodSpan);
durationComponent.AddUtilityTerm(52, income0To25KFlag * periodSpan);
durationComponent.AddUtilityTerm(54, income100KPlusFlag * periodSpan);
durationComponent.AddUtilityTerm(56, highPrioritySameFlag * periodSpan);
durationComponent.AddUtilityTerm(58, lowPrioritySameFlag * periodSpan);
durationComponent.AddUtilityTerm(156, highPriorityDifferentFlag * periodSpan);
durationComponent.AddUtilityTerm(158, lowPriorityDifferentFlag * periodSpan);
durationComponent.AddUtilityTerm(60, homeBasedToursOnlyFlag * periodSpan);
durationComponent.AddUtilityTerm(62, totalStopPurposes * homeBasedToursOnlyFlag * periodSpan);
durationComponent.AddUtilityTerm(64, totalStopPurposes * (1 - homeBasedToursOnlyFlag) * periodSpan);
durationComponent.AddUtilityTerm(66, (totalStopPurposes - totalSimulatedStops) * (1 - homeBasedToursOnlyFlag) * periodSpan);
durationComponent.AddUtilityTerm(68, escortStops * periodSpan);
durationComponent.AddUtilityTerm(70, tour.TotalSubtours * periodSpan);
durationComponent.AddUtilityTerm(71, fulltimeWorkerFlag * durationUnder9HoursFlag);
durationComponent.AddUtilityTerm(169, escortTourFlag * durationUnder1HourFlag);
durationComponent.AddUtilityTerm(170, shoppingTourFlag * durationUnder1HourFlag);
durationComponent.AddUtilityTerm(171, mealTourFlag * durationUnder1HourFlag);
durationComponent.AddUtilityTerm(172, escortTourFlag * duration1To2HoursFlag);
durationComponent.AddUtilityTerm(173, shoppingTourFlag * duration1To2HoursFlag);
durationComponent.AddUtilityTerm(174, mealTourFlag * duration1To2HoursFlag);
durationComponent.AddUtilityTerm(81, highPrioritySameFlag * durationUnder8HoursFlag);
durationComponent.AddUtilityTerm(82, lowPrioritySameFlag * durationUnder8HoursFlag);
durationComponent.AddUtilityTerm(83, highPriorityDifferentFlag * durationUnder8HoursFlag);
durationComponent.AddUtilityTerm(84, lowPriorityDifferentFlag * durationUnder4HoursFlag);
}
foreach (TourTime time in TourTime.Times)
{
bool available =
(time.ArrivalPeriod.Index < time.DeparturePeriod.Index &&
timeWindow.EntireSpanIsAvailable(time.ArrivalPeriod.End, time.DeparturePeriod.Start)) ||
(time.ArrivalPeriod.Index == time.DeparturePeriod.Index &&
timeWindow.TotalAvailableMinutes(time.ArrivalPeriod.Start, time.ArrivalPeriod.End) > 1);
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(time.Index, available, choice != null && choice.Equals(time));
if (!alternative.Available && !Global.Configuration.IsInEstimationMode)
{
continue;
}
int arrivalPeriodIndex = time.ArrivalPeriod.Index;
ITourModeImpedance arrivalImpedance = impedances[arrivalPeriodIndex];
int departurePeriodIndex = time.DeparturePeriod.Index;
ITourModeImpedance departureImpedance = impedances[departurePeriodIndex];
alternative.Choice = time;
// arrival period utility component
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(3 * time.ArrivalPeriod.Index + 0));
// departure period utility component
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(3 * time.DeparturePeriod.Index + 1));
// duration utility components
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(3 * (time.DeparturePeriod.Index - time.ArrivalPeriod.Index) + 2));
alternative.AddUtilityTerm(97, Math.Min(0.3, (homeBasedTours - simulatedHomeBasedTours) / (1.0 + arrivalImpedance.TotalMinutesBefore + departureImpedance.TotalMinutesAfter)));
alternative.AddUtilityTerm(98, Math.Min(0.3, (homeBasedTours - simulatedHomeBasedTours) / (1.0 + Math.Max(arrivalImpedance.MaxMinutesBefore, departureImpedance.MaxMinutesAfter))));
}
}
ITripWrapper ITripCreator.CreateWrapper(ITrip trip, ITourWrapper tourWrapper, IHalfTour halfTour)
{
return(CreateWrapper(trip, tourWrapper, halfTour));
}
protected virtual void RegionSpecificOtherTourDistrictCoefficients(ChoiceProbabilityCalculator.Alternative alternative, ITourWrapper tour, IParcelWrapper destinationParcel, IPersonWrapper person)
{
//see PSRC_OtherTourDestinationModel for example
//Global.PrintFile.WriteLine("Generic OtherTourDestinationModel.RegionSpecificOtherTourDistrictCoefficients being called so must not be overridden by CustomizationDll");
}
ITripWrapper ITripCreator.CreateWrapper(ITourWrapper tourWrapper, ITripWrapper trip, int nextTripId, int intermediateStopPurpose, int destinationPurpose, IHalfTour halfTour)
{
return(CreateWrapper(tourWrapper, trip, nextTripId, intermediateStopPurpose, destinationPurpose, halfTour));
}
public static bool ShouldRunInEstimationModeForModel(ITourWrapper tour)
{
// determine validity and need, then characteristics
// detect and skip invalid trip records (error = true) and those that trips that don't require stop location choice (need = false)
int excludeReason = 0;
// if (_maxZone == -1) {
// // TODO: Verify / Optimize
// _maxZone = ChoiceModelRunner.ZoneKeys.Max(z => z.Key);
// }
//
// if (_maxParcel == -1) {
// // TODO: Optimize
// _maxParcel = ChoiceModelRunner.Parcels.Values.Max(parcel => parcel.Id);
// }
if (Global.Configuration.IsInEstimationMode)
{
// if (tour.OriginParcelId > _maxParcel) {
// excludeReason = 3;
// }
if (tour.OriginParcelId <= 0)
{
excludeReason = 4;
}
// else if (tour.DestinationAddressType > _maxParcel) {
// excludeReason = 5;
// }
else if (tour.DestinationParcelId <= 0)
{
excludeReason = 6;
}
// else if (tour.OriginParcelId > _maxParcel) {
// excludeReason = 7;
// }
// else if (tour.OriginParcelId <= 0) {
// excludeReason = 8;
// }
else if (tour.OriginParcelId == tour.DestinationParcelId && Global.Configuration.DestinationScale == 0)
{
excludeReason = 9;
}
else if (tour.OriginParcel.ZoneId == -1)
{
// TODO: Verify this condition... it used to check that the zone was == null.
// I'm not sure what the appropriate condition should be though.
excludeReason = 10;
}
if (excludeReason > 0)
{
Global.PrintFile.WriteEstimationRecordExclusionMessage(CHOICE_MODEL_NAME, "ShouldRunInEstimationModeForModel", tour.Household.Id, tour.Person.Sequence, 0, tour.Sequence, 0, 0, excludeReason);
}
}
bool shouldRun = (excludeReason == 0);
return(shouldRun);
}
protected virtual void RegionSpecificCustomizations(ChoiceProbabilityCalculator.Alternative alternative, ITourWrapper tour, int pathType, int mode, IParcelWrapper destinationParcel)
{
//see PSRC customization dll for example
//Global.PrintFile.WriteLine("Generic Default WorkTourModeModel.RegionSpecificCustomizations being called so must not be overridden by CustomizationDll");
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, int sampleSize, IParcelWrapper choice = null)
{
// var household = tour.Household;
IPersonWrapper person = tour.Person;
IPersonDayWrapper personDay = tour.PersonDay;
// var totalAvailableMinutes =
// tour.ParentTour == null
// ? personDay.TimeWindow.TotalAvailableMinutes(1, Global.Settings.Times.MinutesInADay)
// : tour.ParentTour.TimeWindow.TotalAvailableMinutes(1, Global.Settings.Times.MinutesInADay);
int maxAvailableMinutes =
tour.ParentTour == null
? personDay.TimeWindow.MaxAvailableMinutesAfter(121)
: tour.ParentTour.DestinationDepartureTime - tour.ParentTour.DestinationArrivalTime;
// var hoursAvailableInverse =
// tour.IsHomeBasedTour
// ? (personDay.HomeBasedTours - personDay.SimulatedHomeBasedTours + 1) / (Math.Max(totalAvailableMinutes - 360, 30) / 60D)
// : 1 / (Math.Max(totalAvailableMinutes, 1) / 60D);
int fastestAvailableTimeOfDay =
tour.IsHomeBasedTour || tour.ParentTour == null
? 1
: tour.ParentTour.DestinationArrivalTime + (tour.ParentTour.DestinationDepartureTime - tour.ParentTour.DestinationArrivalTime) / 2;
int tourCategory = tour.GetTourCategory();
// var primaryFlag = ChoiceModelUtility.GetPrimaryFlag(tourCategory);
int secondaryFlag = ChoiceModelUtility.GetSecondaryFlag(tourCategory);
ChoiceModelUtility.DrawRandomTourTimePeriods(tour, tourCategory);
int segment = Global.ContainerDaySim.GetInstance <SamplingWeightsSettingsFactory>().SamplingWeightsSettings.GetTourDestinationSegment(tour.DestinationPurpose, tour.IsHomeBasedTour ? Global.Settings.TourPriorities.HomeBasedTour : Global.Settings.TourPriorities.WorkBasedTour, Global.Settings.Modes.Sov, person.PersonType);
DestinationSampler destinationSampler = new DestinationSampler(choiceProbabilityCalculator, segment, sampleSize, choice, tour.OriginParcel);
TourDestinationUtilities tourDestinationUtilities = new TourDestinationUtilities(this, tour, sampleSize, secondaryFlag, personDay.GetIsWorkOrSchoolPattern().ToFlag(), personDay.GetIsOtherPattern().ToFlag(), fastestAvailableTimeOfDay, maxAvailableMinutes);
destinationSampler.SampleTourDestinations(tourDestinationUtilities);
}
public void Run(ITourWrapper tour)
{
if (tour == null)
{
throw new ArgumentNullException("tour");
}
tour.PersonDay.ResetRandom(40 + tour.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
}
ChoiceProbabilityCalculator choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(tour.Id);
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
if (tour.DestinationParcel == null ||
tour.Mode <= Global.Settings.Modes.None ||
tour.Mode == Global.Settings.Modes.SchoolBus ||
(tour.Mode == Global.Settings.Modes.PaidRideShare && !Global.Configuration.PaidRideShareModeIsAvailable) ||
tour.Mode > Global.Settings.Modes.PaidRideShare)
{
return;
}
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.RunAllPlusParkAndRide(
tour.Household.RandomUtility,
tour.OriginParcel,
tour.DestinationParcel,
tour.DestinationArrivalTime,
tour.DestinationDepartureTime,
tour.DestinationPurpose,
tour.CostCoefficient,
tour.TimeCoefficient,
tour.Person.IsDrivingAge,
tour.Household.VehiclesAvailable,
tour.Person.TransitPassOwnership,
tour.Household.OwnsAutomatedVehicles > 0,
tour.Person.GetTransitFareDiscountFraction(),
false);
IPathTypeModel pathTypeModel = pathTypeModels.First(x => x.Mode == tour.Mode);
if (!pathTypeModel.Available)
{
return;
}
RunModel(choiceProbabilityCalculator, tour, pathTypeModels, tour.DestinationParcel, tour.Household.VehiclesAvailable, tour.Mode);
choiceProbabilityCalculator.WriteObservation();
}
else
{
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.RunAllPlusParkAndRide(
tour.Household.RandomUtility,
tour.OriginParcel,
tour.DestinationParcel,
tour.DestinationArrivalTime,
tour.DestinationDepartureTime,
tour.DestinationPurpose,
tour.CostCoefficient,
tour.TimeCoefficient,
tour.Person.IsDrivingAge,
tour.Household.VehiclesAvailable,
tour.Person.TransitPassOwnership,
tour.Household.OwnsAutomatedVehicles > 0,
tour.Person.GetTransitFareDiscountFraction(),
false);
RunModel(choiceProbabilityCalculator, tour, pathTypeModels, tour.DestinationParcel, tour.Household.VehiclesAvailable);
ChoiceProbabilityCalculator.Alternative chosenAlternative = choiceProbabilityCalculator.SimulateChoice(tour.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", tour.PersonDay.Id);
tour.Mode = Global.Settings.Modes.Hov3;
tour.PersonDay.IsValid = false;
return;
}
int choice = (int)chosenAlternative.Choice;
tour.Mode = choice;
if (choice == Global.Settings.Modes.SchoolBus || choice == Global.Settings.Modes.PaidRideShare)
{
tour.PathType = 0;
}
else
{
IPathTypeModel chosenPathType = pathTypeModels.First(x => x.Mode == choice);
tour.PathType = chosenPathType.PathType;
if (choice == Global.Settings.Modes.ParkAndRide)
{
tour.ParkAndRideNodeId = chosenPathType.PathParkAndRideNodeId;
tour.ParkAndRidePathType = chosenPathType.PathType;
tour.ParkAndRideTransitTime = chosenPathType.PathParkAndRideTransitTime;
tour.ParkAndRideTransitDistance = chosenPathType.PathParkAndRideTransitDistance;
tour.ParkAndRideTransitCost = chosenPathType.PathParkAndRideTransitCost;
tour.ParkAndRideWalkAccessEgressTime = chosenPathType.PathParkAndRideWalkAccessEgressTime;
tour.ParkAndRideTransitGeneralizedTime = chosenPathType.PathParkAndRideTransitGeneralizedTime;
if (Global.StopAreaIsEnabled)
{
tour.ParkAndRideOriginStopAreaKey = chosenPathType.PathOriginStopAreaKey;
tour.ParkAndRideDestinationStopAreaKey = chosenPathType.PathDestinationStopAreaKey;
}
}
}
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITripWrapper trip, IHouseholdDayWrapper householdDay, 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;
List <IPersonDayWrapper> personDays = householdDay.PersonDays;
int isJointTour = tour.JointTourSequence > 0 ? 1 : 0;
IParcelWrapper destinationParcel = tour.DestinationParcel;
int fullJointHalfTour = ((trip.Direction == Global.Settings.TourDirections.OriginToDestination && tour.FullHalfTour1Sequence > 0) ||
(trip.Direction == Global.Settings.TourDirections.DestinationToOrigin && tour.FullHalfTour2Sequence > 0)).ToFlag();
int partialJointHalfTour = ((trip.Direction == Global.Settings.TourDirections.OriginToDestination && tour.PartialHalfTour1Sequence > 0) ||
(trip.Direction == Global.Settings.TourDirections.DestinationToOrigin && tour.PartialHalfTour2Sequence > 0)).ToFlag();
int isJoint = (isJointTour + fullJointHalfTour + partialJointHalfTour > 0) ? 1 : 0;
int isIndividual = 1 - isJoint;
//destination parcel variables
double foodBuffer2 = 0.0;
double totEmpBuffer2 = 0.0;
double retailBuffer2 = 0.0;
if (destinationParcel != null)
{
foodBuffer2 = Math.Log(1 + destinationParcel.EmploymentFoodBuffer2);
totEmpBuffer2 = Math.Log(1 + destinationParcel.EmploymentTotalBuffer2);
retailBuffer2 = Math.Log(1 + destinationParcel.EmploymentRetailBuffer2);
}
int carOwnership = person.GetCarOwnershipSegment();
// household inputs
int onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
//var householdInc75KP = household.Has75KPlusIncome;
int votALSegment = tour.GetVotALSegment();
int transitAccessSegment = household.ResidenceParcel.TransitAccessSegment();
double totalAggregateLogsum = Global.AggregateLogsums[household.ResidenceParcel.ZoneId]
[Global.Settings.Purposes.HomeBasedComposite][carOwnership][votALSegment][transitAccessSegment];
double homeFoodBuffer2 = Math.Log(1 + household.ResidenceParcel.EmploymentFoodBuffer2);
double homeTotEmpBuffer2 = Math.Log(1 + household.ResidenceParcel.EmploymentTotalBuffer2);
double homeRetailBuffer2 = Math.Log(1 + household.ResidenceParcel.EmploymentRetailBuffer2);
// 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;
int simulatedEscortStops = personDay.SimulatedEscortStops;
int simulatedPersonalBusinessStops = personDay.SimulatedPersonalBusinessStops;
int simulatedShoppingStops = personDay.SimulatedShoppingStops;
int simulatedMealStops = personDay.SimulatedMealStops;
int simulatedSocialStops = personDay.SimulatedSocialStops;
int simulatedRecreationStops = personDay.SimulatedRecreationStops;
int simulatedMedicalStops = personDay.SimulatedMedicalStops;
// tour inputs
int hov2TourFlag = tour.IsHov2Mode().ToFlag();
int hov3TourFlag = tour.IsHov3Mode().ToFlag();
int transitTourFlag = tour.IsTransitMode().ToFlag();
int notHomeBasedTourFlag = (!tour.IsHomeBasedTour).ToFlag();
int workTourFlag = tour.IsWorkPurpose().ToFlag();
int personalBusinessOrMedicalTourFlag = tour.IsPersonalBusinessOrMedicalPurpose().ToFlag();
int socialTourFlag = tour.IsSocialPurpose().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 fivePlusSimulatedTripsFlag = halfTour.FivePlusSimulatedTripsFlag; //JLB changed 5 to 5Plus
int halfTourFromOriginFlag = trip.IsHalfTourFromOrigin.ToFlag();
int halfTourFromDestinationFlag = (!trip.IsHalfTourFromOrigin).ToFlag();
int beforeMandatoryDestinationFlag = trip.IsBeforeMandatoryDestination().ToFlag();
// remaining inputs, including joint tour variables
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();
int remainingToursCount = personDay.HomeBasedTours - personDay.GetTotalSimulatedTours();
// connectivity attributes
double c34Ratio = trip.OriginParcel.C34RatioBuffer1();
double adis = 0.0;
double logDist = 0.0;
int 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);
}
double 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;
}
}
double 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);
double 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)
{
ITourWrapper 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)
{
ITourWrapper 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)
{
ITourWrapper 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
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 * 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);
}
}
protected override void RegionSpecificCustomizations(ChoiceProbabilityCalculator.Alternative alternative, ITourWrapper tour, int pathType, int mode, IParcelWrapper destinationParcel)
{
//Global.PrintFile.WriteLine("Default PSRC_WorkTourModeModel.RegionSpecificCustomizations called");
if (mode == Global.Settings.Modes.Transit && pathType != Global.Settings.PathTypes.LightRail && pathType != Global.Settings.PathTypes.CommuterRail && pathType != Global.Settings.PathTypes.Ferry)
{
alternative.AddUtilityTerm(200 + tour.OriginParcel.District, 1); //district specific transit calibration constant
alternative.AddUtilityTerm(300 + destinationParcel.District, 1); //district specific transit calibration constant
}
if (mode == Global.Settings.Modes.ParkAndRide)
{
alternative.AddUtilityTerm(250, pathType == 3 ? 1 : 0);
alternative.AddUtilityTerm(251, pathType == 4 ? 1 : 0);
alternative.AddUtilityTerm(252, pathType == 5 ? 1 : 0);
alternative.AddUtilityTerm(253, pathType == 6 ? 1 : 0);
alternative.AddUtilityTerm(254, pathType == 7 ? 1 : 0);
}
else if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(255, pathType == 3 ? 1 : 0);
alternative.AddUtilityTerm(256, pathType == 4 ? 1 : 0);
alternative.AddUtilityTerm(257, pathType == 5 ? 1 : 0);
alternative.AddUtilityTerm(258, pathType == 6 ? 1 : 0);
alternative.AddUtilityTerm(259, pathType == 7 ? 1 : 0);
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, IEnumerable <PathTypeModel> pathTypeModels, IParcelWrapper destinationParcel, int choice = Constants.DEFAULT_VALUE)
{
var household = tour.Household;
var householdTotals = household.HouseholdTotals;
var person = tour.Person;
var personDay = tour.PersonDay;
// household inputs
var income0To25KFlag = household.Has0To25KIncome.ToFlag();
var childrenUnder5 = householdTotals.ChildrenUnder5;
var childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
var nonworkingAdults = householdTotals.NonworkingAdults;
var retiredAdults = householdTotals.RetiredAdults;
var onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
var twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
var noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(household.VehiclesAvailable);
var carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(household.VehiclesAvailable);
var carsLessThanWorkersFlag = household.GetFlagForCarsLessThanWorkers(household.VehiclesAvailable);
// person inputs
var maleFlag = person.IsMale.ToFlag();
var ageBetween51And98Flag = person.AgeIsBetween51And98.ToFlag();
// tour inputs
var shoppingTourFlag = (tour.DestinationPurpose == Global.Settings.Purposes.Shopping).ToFlag();
var mealTourFlag = (tour.DestinationPurpose == Global.Settings.Purposes.Meal).ToFlag();
var socialOrRecreationTourFlag = (tour.DestinationPurpose == Global.Settings.Purposes.Social).ToFlag();
// remaining inputs
var originParcel = tour.OriginParcel;
var parkingDuration = ChoiceModelUtility.GetParkingDuration(person.IsFulltimeWorker);
var destinationParkingCost = destinationParcel.ParkingCostBuffer1(parkingDuration);
double escortPercentage;
double nonEscortPercentage;
ChoiceModelUtility.SetEscortPercentages(personDay, out escortPercentage, out nonEscortPercentage, true);
foreach (var pathTypeModel in pathTypeModels)
{
var mode = pathTypeModel.Mode;
var available = pathTypeModel.Mode != Global.Settings.Modes.ParkAndRide && pathTypeModel.Available;
var generalizedTimeLogsum = pathTypeModel.GeneralizedTimeLogsum;
var alternative = choiceProbabilityCalculator.GetAlternative(mode, available, choice == mode);
alternative.Choice = mode;
alternative.AddNestedAlternative(_nestedAlternativeIds[pathTypeModel.Mode], _nestedAlternativeIndexes[pathTypeModel.Mode], THETA_PARAMETER);
if (!available)
{
continue;
}
alternative.AddUtilityTerm(2, generalizedTimeLogsum * tour.TimeCoefficient);
if (mode == Global.Settings.Modes.Transit)
{
alternative.AddUtilityTerm(20, 1);
alternative.AddUtilityTerm(21, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(120, shoppingTourFlag);
// alternative.AddUtility(121, mealTourFlag);
// alternative.AddUtility(129, destinationParcel.MixedUse2Index1());
alternative.AddUtilityTerm(128, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(127, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(126, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(125, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(124, originParcel.MixedUse2Index1());
// alternative.AddUtility(123, Math.Log(destinationParcel.StopsTransitBuffer1+1));
// alternative.AddUtility(122, Math.Log(originParcel.StopsTransitBuffer1+1));
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(31, childrenUnder5);
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
alternative.AddUtilityTerm(35, pathTypeModel.PathDistance.AlmostEquals(0) ? 0 : Math.Log(pathTypeModel.PathDistance));
alternative.AddUtilityTerm(38, onePersonHouseholdFlag);
alternative.AddUtilityTerm(39, twoPersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(43, carsLessThanWorkersFlag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
alternative.AddUtilityTerm(136, shoppingTourFlag);
alternative.AddUtilityTerm(137, mealTourFlag);
alternative.AddUtilityTerm(138, socialOrRecreationTourFlag);
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(31, childrenUnder5);
alternative.AddUtilityTerm(32, childrenAge5Through15);
alternative.AddUtilityTerm(34, nonworkingAdults + retiredAdults);
alternative.AddUtilityTerm(35, pathTypeModel.PathDistance.AlmostEquals(0) ? 0 : Math.Log(pathTypeModel.PathDistance));
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(43, carsLessThanWorkersFlag);
alternative.AddUtilityTerm(48, onePersonHouseholdFlag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
alternative.AddUtilityTerm(136, shoppingTourFlag);
alternative.AddUtilityTerm(137, mealTourFlag);
alternative.AddUtilityTerm(138, socialOrRecreationTourFlag);
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient));
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(52, carsLessThanDriversFlag);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(131, escortPercentage);
alternative.AddUtilityTerm(132, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Bike)
{
var class1Dist = Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions ?
ImpedanceRoster.GetValue("class1distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable : 0;
var class2Dist = Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions ?
ImpedanceRoster.GetValue("class2distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable : 0;
var worstDist = Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions ?
ImpedanceRoster.GetValue("worstdistance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable : 0;
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(63, ageBetween51And98Flag);
alternative.AddUtilityTerm(160, socialOrRecreationTourFlag);
alternative.AddUtilityTerm(169, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(168, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(167, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(166, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(165, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(164, originParcel.MixedUse4Index1());
alternative.AddUtilityTerm(161, (class1Dist > 0).ToFlag());
alternative.AddUtilityTerm(162, (class2Dist > 0).ToFlag());
alternative.AddUtilityTerm(163, (worstDist > 0).ToFlag());
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(73, ageBetween51And98Flag);
alternative.AddUtilityTerm(171, mealTourFlag);
alternative.AddUtilityTerm(172, socialOrRecreationTourFlag);
// alternative.AddUtility(179, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(178, destinationParcel.HouseholdDensity1());
// alternative.AddUtility(177, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(176, originParcel.NetIntersectionDensity1());
alternative.AddUtilityTerm(175, originParcel.HouseholdDensity1());
// alternative.AddUtility(174, originParcel.MixedUse4Index1());
}
}
}
public void Run(ITourWrapper subtour)
{
if (subtour == null)
{
throw new ArgumentNullException("subtour");
}
subtour.PersonDay.ResetRandom(40 + subtour.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
// JLB 20150331 added following to exclude records with missing coordinates from estimation
if (subtour.OriginParcelId == 0 || subtour.DestinationParcelId == 0)
{
return;
}
}
var choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(subtour.Id);
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
if (subtour.DestinationParcel == null || subtour.Mode <= Global.Settings.Modes.None || subtour.Mode > Global.Settings.Modes.Transit)
{
return;
}
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.RunAll(
subtour.Household.RandomUtility,
subtour.OriginParcel,
subtour.DestinationParcel,
subtour.DestinationArrivalTime,
subtour.DestinationDepartureTime,
subtour.DestinationPurpose,
subtour.CostCoefficient,
subtour.TimeCoefficient,
subtour.Person.IsDrivingAge,
subtour.Household.VehiclesAvailable,
subtour.Household.OwnsAutomatedVehicles > 0,
subtour.Person.GetTransitFareDiscountFraction(),
false);
var pathTypeModel = pathTypeModels.First(x => x.Mode == subtour.Mode);
if (!pathTypeModel.Available)
{
return;
}
RunModel(choiceProbabilityCalculator, subtour, pathTypeModels, subtour.DestinationParcel, subtour.ParentTour.Mode, subtour.Mode);
choiceProbabilityCalculator.WriteObservation();
}
else
{
IEnumerable <IPathTypeModel> pathTypeModels =
PathTypeModelFactory.Singleton.RunAll(
subtour.Household.RandomUtility,
subtour.OriginParcel,
subtour.DestinationParcel,
subtour.DestinationArrivalTime,
subtour.DestinationDepartureTime,
subtour.DestinationPurpose,
subtour.CostCoefficient,
subtour.TimeCoefficient,
subtour.Person.IsDrivingAge,
subtour.Household.VehiclesAvailable,
subtour.Household.OwnsAutomatedVehicles > 0,
subtour.Person.GetTransitFareDiscountFraction(),
false);
RunModel(choiceProbabilityCalculator, subtour, pathTypeModels, subtour.DestinationParcel, subtour.ParentTour.Mode);
var chosenAlternative = choiceProbabilityCalculator.SimulateChoice(subtour.Household.RandomUtility);
if (chosenAlternative == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", subtour.PersonDay.Id);
subtour.Mode = Global.Settings.Modes.Hov3;
subtour.PersonDay.IsValid = false;
return;
}
var choice = (int)chosenAlternative.Choice;
subtour.Mode = choice;
if (choice == Global.Settings.Modes.SchoolBus || choice == Global.Settings.Modes.PaidRideShare)
{
subtour.PathType = 0;
}
else
{
var chosenPathType = pathTypeModels.First(x => x.Mode == choice);
subtour.PathType = chosenPathType.PathType;
subtour.ParkAndRideNodeId = choice == Global.Settings.Modes.ParkAndRide ? chosenPathType.PathParkAndRideNodeId : 0;
}
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, ITourWrapper tour, List <PathTypeModel> pathTypeModels, IParcelWrapper destinationParcel, int householdCars, int choice = Constants.DEFAULT_VALUE)
{
var household = tour.Household;
var person = tour.Person;
var personDay = tour.PersonDay;
// household inputs
var income0To25KFlag = household.Has0To25KIncome.ToFlag();
var income25To50KFlag = household.Has25To50KIncome.ToFlag();
var income75KPlusFlag = household.Has75KPlusIncome.ToFlag();
var onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
var twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
var noCarsInHouseholdFlag = household.GetFlagForNoCarsInHousehold(householdCars);
var carsLessThanDriversFlag = household.GetFlagForCarsLessThanDrivers(householdCars);
// person inputs
var childUnder5Flag = person.IsChildUnder5.ToFlag();
var adultFlag = person.IsAdult.ToFlag();
var drivingAgeStudentFlag = person.IsDrivingAgeStudent.ToFlag();
var maleFlag = person.IsMale.ToFlag();
// remaining inputs
var originParcel = tour.OriginParcel;
var destinationParkingCost = destinationParcel.ParkingCostBuffer1(6);
double escortPercentage;
double nonEscortPercentage;
ChoiceModelUtility.SetEscortPercentages(personDay, out escortPercentage, out nonEscortPercentage);
// school bus is a special case - use HOV3 impedance
var pathTypeExtra = pathTypeModels.First(x => x.Mode == Global.Settings.Modes.Hov3);
var modeExtra = Global.Settings.Modes.SchoolBus;
var availableExtra = pathTypeExtra.Available;
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(10, 1);
alternative.AddUtilityTerm(17, childUnder5Flag);
alternative.AddUtilityTerm(18, adultFlag);
}
foreach (var pathTypeModel in pathTypeModels)
{
var mode = pathTypeModel.Mode;
var available = pathTypeModel.Available;
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(21, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(22, carsLessThanDriversFlag);
alternative.AddUtilityTerm(27, childUnder5Flag);
alternative.AddUtilityTerm(28, adultFlag);
alternative.AddUtilityTerm(29, drivingAgeStudentFlag);
alternative.AddUtilityTerm(129, destinationParcel.MixedUse2Index1());
// alternative.AddUtility(128, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(127, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(126, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(125, originParcel.HouseholdDensity1());
// alternative.AddUtility(124, originParcel.MixedUse2Index1());
// alternative.AddUtility(123, Math.Log(destinationParcel.StopsTransitBuffer1+1));
alternative.AddUtilityTerm(122, Math.Log(originParcel.StopsTransitBuffer1 + 1));
}
else if (mode == Global.Settings.Modes.Hov3)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT3));
alternative.AddUtilityTerm(30, 1);
alternative.AddUtilityTerm(37, twoPersonHouseholdFlag);
alternative.AddUtilityTerm(37, onePersonHouseholdFlag);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(44, income0To25KFlag);
alternative.AddUtilityTerm(45, income25To50KFlag);
alternative.AddUtilityTerm(47, childUnder5Flag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Hov2)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient / ChoiceModelUtility.CPFACT2));
alternative.AddUtilityTerm(38, onePersonHouseholdFlag);
alternative.AddUtilityTerm(40, 1);
alternative.AddUtilityTerm(41, noCarsInHouseholdFlag);
alternative.AddUtilityTerm(44, income0To25KFlag);
alternative.AddUtilityTerm(45, income25To50KFlag);
alternative.AddUtilityTerm(47, childUnder5Flag);
alternative.AddUtilityTerm(133, escortPercentage);
alternative.AddUtilityTerm(134, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Sov)
{
alternative.AddUtilityTerm(1, (destinationParkingCost * tour.CostCoefficient));
alternative.AddUtilityTerm(50, 1);
alternative.AddUtilityTerm(52, carsLessThanDriversFlag);
alternative.AddUtilityTerm(54, income0To25KFlag);
alternative.AddUtilityTerm(56, income75KPlusFlag);
alternative.AddUtilityTerm(59, drivingAgeStudentFlag);
alternative.AddUtilityTerm(131, escortPercentage);
alternative.AddUtilityTerm(132, nonEscortPercentage);
}
else if (mode == Global.Settings.Modes.Bike)
{
var class1Dist =
Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions
? ImpedanceRoster.GetValue("class1distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable
: 0;
var class2Dist =
Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions
? ImpedanceRoster.GetValue("class2distance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable
: 0;
var worstDist =
Global.Configuration.PathImpedance_BikeUseTypeSpecificDistanceFractions
? ImpedanceRoster.GetValue("worstdistance", mode, Global.Settings.PathTypes.FullNetwork,
Global.Settings.ValueOfTimes.DefaultVot, tour.DestinationArrivalTime, originParcel, destinationParcel).Variable
: 0;
alternative.AddUtilityTerm(60, 1);
alternative.AddUtilityTerm(61, maleFlag);
alternative.AddUtilityTerm(69, adultFlag);
alternative.AddUtilityTerm(169, destinationParcel.MixedUse4Index1());
alternative.AddUtilityTerm(168, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(167, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(166, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(165, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(164, originParcel.MixedUse4Index1());
alternative.AddUtilityTerm(161, (class1Dist > 0).ToFlag());
alternative.AddUtilityTerm(162, (class2Dist > 0).ToFlag());
alternative.AddUtilityTerm(163, (worstDist > 0).ToFlag());
}
else if (mode == Global.Settings.Modes.Walk)
{
alternative.AddUtilityTerm(79, adultFlag);
alternative.AddUtilityTerm(179, destinationParcel.MixedUse4Index1());
// alternative.AddUtility(178, destinationParcel.TotalEmploymentDensity1());
// alternative.AddUtility(177, destinationParcel.NetIntersectionDensity1());
// alternative.AddUtility(176, originParcel.NetIntersectionDensity1());
// alternative.AddUtility(175, originParcel.HouseholdDensity1());
alternative.AddUtilityTerm(179, originParcel.MixedUse4Index1());
}
}
}
public void Run(IHouseholdDayWrapper householdDay, ITourWrapper tour,
int constrainedMode, int constrainedArrivalTime, int constrainedDepartureTime)
{
if (tour == null)
{
throw new ArgumentNullException("tour");
}
//HTourModeTime.InitializeTourModeTimes();
tour.PersonDay.ResetRandom(50 + tour.Sequence - 1);
if (Global.Configuration.IsInEstimationMode)
{
if (Global.Configuration.EstimationModel != CHOICE_MODEL_NAME)
{
return;
}
// JLB 20140421 add the following to keep from estimatign twice for the same tour
if (tour.DestinationModeAndTimeHaveBeenSimulated)
{
return;
}
if (tour.DestinationParcel == null || tour.OriginParcel == null || tour.Mode < Global.Settings.Modes.Walk || tour.Mode > Global.Settings.Modes.WalkRideBike)
{
return;
}
}
// set remaining inputs
HTourModeTime.SetModeTimeImpedances(householdDay, tour, constrainedMode, constrainedArrivalTime, constrainedDepartureTime, -1, -1.0);
ChoiceProbabilityCalculator choiceProbabilityCalculator = _helpers[ParallelUtility.threadLocalAssignedIndex.Value].GetChoiceProbabilityCalculator(tour.Id);
if (_helpers[ParallelUtility.threadLocalAssignedIndex.Value].ModelIsInEstimationMode)
{
HTourModeTime observedChoice = new HTourModeTime(tour.Mode, tour.DestinationArrivalTime, tour.DestinationDepartureTime);
RunModel(choiceProbabilityCalculator, householdDay, tour, tour.DestinationParcel, tour.Household.VehiclesAvailable,
constrainedMode, constrainedArrivalTime, constrainedDepartureTime,
observedChoice);
choiceProbabilityCalculator.WriteObservation();
}
else if (Global.Configuration.TestEstimationModelInApplicationMode)
{
Global.Configuration.IsInEstimationMode = false;
RunModel(choiceProbabilityCalculator, householdDay, tour, tour.DestinationParcel, tour.Household.VehiclesAvailable,
constrainedMode, constrainedArrivalTime, constrainedDepartureTime);
HTourModeTime observedChoice = new HTourModeTime(tour.Mode, tour.DestinationArrivalTime, tour.DestinationDepartureTime);
ChoiceProbabilityCalculator.Alternative simulatedChoice = choiceProbabilityCalculator.SimulateChoice(tour.Household.RandomUtility, tour.Id, observedChoice.Index);
Global.Configuration.IsInEstimationMode = true;
}
else
{
HTourModeTime choice;
if (constrainedMode > 0 && constrainedArrivalTime > 0 && constrainedDepartureTime > 0)
{
choice = new HTourModeTime(constrainedMode, constrainedArrivalTime, constrainedDepartureTime);
}
else
{
RunModel(choiceProbabilityCalculator, householdDay, tour, tour.DestinationParcel, tour.Household.VehiclesAvailable,
constrainedMode, constrainedArrivalTime, constrainedDepartureTime);
ChoiceProbabilityCalculator.Alternative simulatedChoice = choiceProbabilityCalculator.SimulateChoice(tour.Household.RandomUtility);
if (simulatedChoice == null)
{
Global.PrintFile.WriteNoAlternativesAvailableWarning(CHOICE_MODEL_NAME, "Run", tour.PersonDay.Id);
if (!Global.Configuration.IsInEstimationMode)
{
tour.PersonDay.IsValid = false;
tour.PersonDay.HouseholdDay.IsValid = false;
}
return;
}
choice = (HTourModeTime)simulatedChoice.Choice;
}
tour.Mode = choice.Mode;
if (Global.StopAreaIsEnabled && tour.Mode == Global.Settings.Modes.ParkAndRide)
{
tour.ParkAndRideOriginStopAreaKey = choice.ParkAndRideOriginStopAreaKey;
tour.ParkAndRideDestinationStopAreaKey = choice.ParkAndRideDestinationStopAreaKey;
}
MinuteSpan arrivalPeriod = choice.ArrivalPeriod;
MinuteSpan departurePeriod = choice.DeparturePeriod;
//use constrained times to set temporary arrival and departure times with minimum duration of stay for time window calculations
if (constrainedArrivalTime > 0 || constrainedDepartureTime > 0)
{
if (constrainedArrivalTime > 0)
{
tour.DestinationArrivalTime = constrainedArrivalTime;
}
else
{
tour.DestinationArrivalTime = Math.Min(arrivalPeriod.End, constrainedDepartureTime - Global.Settings.Times.MinimumActivityDuration);
}
if (constrainedDepartureTime > 0)
{
tour.DestinationDepartureTime = constrainedDepartureTime;
}
else
{
tour.DestinationDepartureTime = Math.Max(departurePeriod.Start, constrainedArrivalTime + Global.Settings.Times.MinimumActivityDuration);
}
}
//or if times aren't constrained use periods to set temporary arrival and departure times with minimum duration of stay for time window calculations
else if (arrivalPeriod == departurePeriod)
{
int departureTime = Math.Max(choice.GetRandomDepartureTime(householdDay, tour), departurePeriod.Start + Global.Settings.Times.MinimumActivityDuration);
tour.DestinationArrivalTime = departureTime - Global.Settings.Times.MinimumActivityDuration;
tour.DestinationDepartureTime = departureTime;
}
else if (arrivalPeriod.End == departurePeriod.Start - 1)
{
tour.DestinationArrivalTime = arrivalPeriod.End;
tour.DestinationDepartureTime = arrivalPeriod.End + Global.Settings.Times.MinimumActivityDuration;
}
else
{
tour.DestinationArrivalTime = arrivalPeriod.End;
tour.DestinationDepartureTime = departurePeriod.Start;
}
}
}
protected override void RegionSpecificCustomizations(ChoiceProbabilityCalculator.Alternative alternative, ITourWrapper tour, int pathType, int mode, IParcelWrapper destinationParcel)
{
//Global.PrintFile.WriteLine("Default PSRC_SchoolTourModeModel.RegionSpecificCustomizations called");
if (mode == Global.Settings.Modes.Transit && pathType != Global.Settings.PathTypes.LightRail && pathType != Global.Settings.PathTypes.CommuterRail && pathType != Global.Settings.PathTypes.Ferry)
{
alternative.AddUtilityTerm(200 + tour.OriginParcel.District, 1); //district specific transit calibration constant
alternative.AddUtilityTerm(300 + destinationParcel.District, 1); //district specific transit calibration constant
}
}
