private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, PersonDayWrapper personDay, HouseholdDayWrapper householdDay, int choice = Constants.DEFAULT_VALUE)
{
Framework.DomainModels.Wrappers.IHouseholdWrapper household = personDay.Household;
Framework.DomainModels.Wrappers.IPersonWrapper person = personDay.Person;
IEnumerable <PersonDayWrapper> personTypeOrderedPersonDays = householdDay.PersonDays.OrderBy(p => p.Person.PersonType).ToList().Cast <PersonDayWrapper>();
int mandatoryCount = 0;
int nonMandatoryCount = 0;
int homeCount = 0;
int i = 0;
foreach (PersonDayWrapper pDay in personTypeOrderedPersonDays)
{
i++;
if (i <= 5)
{
if (pDay.PatternType == Global.Settings.PatternTypes.Mandatory)
{
mandatoryCount++;
}
else if (pDay.PatternType == Global.Settings.PatternTypes.Optional)
{
nonMandatoryCount++;
}
else
{
homeCount++;
}
}
}
bool mandatoryAvailableFlag = true;
if (personDay.Person.IsNonworkingAdult || personDay.Person.IsRetiredAdult ||
(!personDay.Person.IsWorker && !personDay.Person.IsStudent) ||
(!Global.Configuration.IsInEstimationMode && !personDay.Person.IsWorker && personDay.Person.UsualSchoolParcel == null)
)
{
mandatoryAvailableFlag = false;
}
int carOwnership =
household.VehiclesAvailable == 0
? Global.Settings.CarOwnerships.NoCars
: household.VehiclesAvailable < household.HouseholdTotals.DrivingAgeMembers
? Global.Settings.CarOwnerships.LtOneCarPerAdult
: Global.Settings.CarOwnerships.OneOrMoreCarsPerAdult;
int noCarsFlag = FlagUtility.GetNoCarsFlag(carOwnership);
int carCompetitionFlag = FlagUtility.GetCarCompetitionFlag(carOwnership);
int votALSegment = Global.Settings.VotALSegments.Medium; // TODO: calculate a VOT segment that depends on household income
int transitAccessSegment = household.ResidenceParcel.TransitAccessSegment();
double personalBusinessAggregateLogsum = Global.AggregateLogsums[household.ResidenceParcel.ZoneId]
[Global.Settings.Purposes.PersonalBusiness][carOwnership][votALSegment][transitAccessSegment];
double shoppingAggregateLogsum = Global.AggregateLogsums[household.ResidenceParcel.ZoneId]
[Global.Settings.Purposes.Shopping][carOwnership][votALSegment][transitAccessSegment];
double mealAggregateLogsum = Global.AggregateLogsums[household.ResidenceParcel.ZoneId]
[Global.Settings.Purposes.Meal][carOwnership][votALSegment][transitAccessSegment];
double socialAggregateLogsum = Global.AggregateLogsums[household.ResidenceParcel.ZoneId]
[Global.Settings.Purposes.Social][carOwnership][votALSegment][transitAccessSegment];
//var compositeLogsum = Global.AggregateLogsums[household.ResidenceZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votALSegment][transitAccessSegment];
double compositeLogsum = Global.AggregateLogsums[household.ResidenceZoneId][Global.Settings.Purposes.HomeBasedComposite][Global.Settings.CarOwnerships.NoCars][votALSegment][transitAccessSegment];
// Pattern Type Mandatory on tour (at least one work or school tour)
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(0, mandatoryAvailableFlag, choice == 1);
alternative.Choice = 1;
alternative.AddUtilityTerm(1, 1);
alternative.AddUtilityTerm(2, person.IsChildUnder5.ToFlag());
alternative.AddUtilityTerm(3, person.IsChildAge5Through15.ToFlag());
alternative.AddUtilityTerm(4, person.IsFulltimeWorker.ToFlag());
alternative.AddUtilityTerm(5, person.IsMale.ToFlag());
//alternative.AddUtilityTerm(4, person.IsPartTimeWorker.ToFlag());
alternative.AddUtilityTerm(7, householdDay.Household.HasChildrenUnder5.ToFlag());
alternative.AddUtilityTerm(8, householdDay.Household.HasChildrenAge5Through15.ToFlag());
alternative.AddUtilityTerm(10, (householdDay.AdultsInSharedHomeStay == 2 && householdDay.Household.HouseholdTotals.FullAndPartTimeWorkers >= 2).ToFlag());
//alternative.AddUtilityTerm(14, (householdDay.AdultsInSharedHomeStay == 2).ToFlag());
//alternative.AddUtilityTerm(15, (householdDay.Household.HouseholdTotals.FullAndPartTimeWorkers >= 2).ToFlag());
alternative.AddUtilityTerm(11, (householdDay.AdultsInSharedHomeStay == 1 && householdDay.Household.HasChildren).ToFlag());
//alternative.AddUtilityTerm(12, (householdDay.Household.Size == 2).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(13, (householdDay.Household.Size == 3).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(14, (householdDay.Household.Size >= 4).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(12, (householdDay.Household.VehiclesAvailable == 1).ToFlag());
//alternative.AddUtilityTerm(13, (householdDay.Household.VehiclesAvailable >= 2).ToFlag());
//GV: not sign. 10. juni 2016
//alternative.AddUtilityTerm(15, (householdDay.Household.VehiclesAvailable == 1 && household.Has2Drivers).ToFlag());
alternative.AddUtilityTerm(16, (householdDay.Household.VehiclesAvailable >= 2 && household.Has2Drivers).ToFlag());
//GV: introduced again - 10. june 2016
//GV: logsum for mandatory - wrong sign
//alternative.AddUtilityTerm(17, compositeLogsum);
//alternative.AddUtilityTerm(17, (householdDay.Household.Income >= 300000 && householdDay.Household.Income < 600000).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(18, (householdDay.Household.Income >= 600000 && householdDay.Household.Income < 900000).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(19, (householdDay.Household.Income >= 900000).ToFlag()); //GV; 16. april 2013, not significant
alternative.AddUtilityTerm(20, householdDay.PrimaryPriorityTimeFlag);
//alternative.AddUtilityTerm(19, (mandatoryCount == 0)? 1 : 0); //GV - goes to infinity
// PatternType NonMandatory on tour (tours, but none for work or school)
alternative = choiceProbabilityCalculator.GetAlternative(1, true, choice == 2);
alternative.Choice = 2;
alternative.AddUtilityTerm(22, person.IsRetiredAdult.ToFlag());
alternative.AddUtilityTerm(23, person.IsNonworkingAdult.ToFlag());
//GV: not sign. 10. june 2016
//alternative.AddUtilityTerm(24, householdDay.Household.HasChildrenUnder5.ToFlag());
//alternative.AddUtilityTerm(25, householdDay.Household.HasChildrenAge5Through15.ToFlag());
//alternative.AddUtilityTerm(31, (householdDay.AdultsInSharedHomeStay == 2).ToFlag());
//alternative.AddUtilityTerm(33, (householdDay.Household.HouseholdTotals.FullAndPartTimeWorkers >= 2).ToFlag());
alternative.AddUtilityTerm(26, (householdDay.AdultsInSharedHomeStay == 2 && householdDay.Household.HouseholdTotals.FullAndPartTimeWorkers >= 2).ToFlag());
alternative.AddUtilityTerm(27, (householdDay.AdultsInSharedHomeStay == 1 && householdDay.Household.HasChildren).ToFlag());
//alternative.AddUtilityTerm(28, (householdDay.Household.Size == 2).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(29, (householdDay.Household.Size == 3).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(30, (householdDay.Household.Size >= 4).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(27, (householdDay.Household.VehiclesAvailable == 1).ToFlag());
//alternative.AddUtilityTerm(28, (householdDay.Household.VehiclesAvailable >= 2).ToFlag());
alternative.AddUtilityTerm(31, (householdDay.Household.VehiclesAvailable == 1 && household.Has2Drivers).ToFlag());
alternative.AddUtilityTerm(32, (householdDay.Household.VehiclesAvailable >= 2 && household.Has2Drivers).ToFlag());
//alternative.AddUtilityTerm(33, compositeLogsum); //GV: logsum for non-mandatory
//alternative.AddUtilityTerm(33, (householdDay.Household.Income >= 300000 && householdDay.Household.Income < 600000).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(34, (householdDay.Household.Income >= 600000 && householdDay.Household.Income < 900000).ToFlag()); //GV; 16. april 2013, not significant
//alternative.AddUtilityTerm(35, (householdDay.Household.Income >= 900000).ToFlag()); //GV; 16. april 2013, not significant
alternative.AddUtilityTerm(36, householdDay.PrimaryPriorityTimeFlag);
//alternative.AddUtilityTerm(24, person.IsChildUnder5.ToFlag());
//alternative.AddUtilityTerm(25, person.IsNonworkingAdult.ToFlag());
// PatternType Home (all day)
alternative = choiceProbabilityCalculator.GetAlternative(2, true, choice == 3);
alternative.Choice = 3;
alternative.AddUtilityTerm(41, 1);
alternative.AddUtilityTerm(42, person.WorksAtHome().ToFlag());
//GV: introduced again - 10. june 2016; not sign.
//alternative.AddUtilityTerm(43, person.IsUniversityStudent.ToFlag());
//alternative.AddUtilityTerm(54, (homeCount > 0)? 1 : 0); //GV: can be estimated but the valus is huge
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, PersonDayWrapper personDay, HouseholdDayWrapper householdDay, DayPattern[] dayPatterns, DayPattern choice = null)
{
Framework.DomainModels.Wrappers.IHouseholdWrapper household = personDay.Household;
Framework.DomainModels.Wrappers.IParcelWrapper residenceParcel = household.ResidenceParcel;
Framework.DomainModels.Wrappers.IPersonWrapper person = personDay.Person;
int carsPerDriver = household.GetCarsPerDriver();
double mixedDensity = residenceParcel.MixedUse3Index2();
double intersectionDensity = residenceParcel.IntersectionDensity34Minus1Buffer2();
double[] purposeLogsums = new double[Global.Settings.Purposes.TotalPurposes + 2];
double[] atUsualLogsums = new double[3];
//var carOwnership = person.CarOwnershipSegment; //GV: sat car ownership not to impact logsums
int carOwnership = 0;
int votSegment = person.Household.GetVotALSegment();
int transitAccess = residenceParcel.TransitAccessSegment();
//GV: input 26. july 2013
// household inputs
//var childrenUnder5 = householdTotals.ChildrenUnder5;
//var childrenAge5Through15 = householdTotals.ChildrenAge5Through15;
//var nonworkingAdults = householdTotals.NonworkingAdults;
//var retiredAdults = householdTotals.RetiredAdults;
int onePersonHouseholdFlag = household.IsOnePersonHousehold.ToFlag();
int twoPersonHouseholdFlag = household.IsTwoPersonHousehold.ToFlag();
int householdCars = household.VehiclesAvailable;
//var noCarsInHouseholdFlag = HouseholdWrapper.GetNoCarsInHouseholdFlag(householdCars);
//var carsLessThanDriversFlag = household.GetCarsLessThanDriversFlag(householdCars);
//var carsLessThanWorkersFlag = household.GetCarsLessThanWorkersFlag(householdCars);
int HHwithChildrenFlag = household.HasChildren.ToFlag();
int HHwithSmallChildrenFlag = household.HasChildrenUnder5.ToFlag();
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 primaryFamilyTimeFlag = householdDay.PrimaryPriorityTimeFlag;
//GV: input 26. july 2013
// person inputs
int partTimeWorkerFlag = person.IsPartTimeWorker.ToFlag();
int nonworkingAdultFlag = person.IsNonworkingAdult.ToFlag();
int universityStudentFlag = person.IsUniversityStudent.ToFlag();
int retiredAdultFlag = person.IsRetiredAdult.ToFlag();
int fullTimeWorkerFlag = person.IsFulltimeWorker.ToFlag();
int childAge5Through15Flag = person.IsChildAge5Through15.ToFlag();
int childUnder5Flag = person.IsChildUnder5.ToFlag();
int adultFlag = person.IsAdult.ToFlag();
int maleFlag = person.IsMale.ToFlag();
int femaleFlag = person.IsFemale.ToFlag();
if (person.UsualWorkParcel == null || person.UsualWorkParcelId == household.ResidenceParcelId)
{
purposeLogsums[Global.Settings.Purposes.Work] = 0;
}
else
{
int destinationArrivalTime = ChoiceModelUtility.GetDestinationArrivalTime(Global.Settings.Models.WorkTourModeModel);
int destinationDepartureTime = ChoiceModelUtility.GetDestinationDepartureTime(Global.Settings.Models.WorkTourModeModel);
//JLB 201406
//var nestedAlternative = Global.ChoiceModelSession.Get<WorkTourModeModel>().RunNested(personDay, residenceParcel, person.UsualWorkParcel, destinationArrivalTime, destinationDepartureTime, household.VehiclesAvailable);
//JLB 201602
//var nestedAlternative = Global.ChoiceModelSession.Get<WorkTourModeTimeModel>().RunNested(personDay, residenceParcel, person.UsualWorkParcel, destinationArrivalTime, destinationDepartureTime, household.VehiclesAvailable);
ChoiceProbabilityCalculator.Alternative nestedAlternative = Global.ChoiceModelSession.Get <TourModeTimeModel>().RunNested(personDay, residenceParcel, person.UsualWorkParcel, destinationArrivalTime, destinationDepartureTime, household.VehiclesAvailable, Global.Settings.Purposes.Work);
purposeLogsums[Global.Settings.Purposes.Work] = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
atUsualLogsums[Global.Settings.Purposes.Work] = Global.AggregateLogsums[person.UsualWorkParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votSegment][person.UsualWorkParcel.TransitAccessSegment()];
}
if (person.UsualSchoolParcel == null || person.UsualSchoolParcelId == household.ResidenceParcelId)
{
purposeLogsums[Global.Settings.Purposes.School] = 0;
}
else
{
int destinationArrivalTime = ChoiceModelUtility.GetDestinationArrivalTime(Global.Settings.Models.SchoolTourModeModel);
int destinationDepartureTime = ChoiceModelUtility.GetDestinationDepartureTime(Global.Settings.Models.SchoolTourModeModel);
//JLB 201406
//var nestedAlternative = Global.ChoiceModelSession.Get<SchoolTourModeModel>().RunNested(personDay, residenceParcel, person.UsualSchoolParcel, destinationArrivalTime, destinationDepartureTime, household.VehiclesAvailable);
//JLB 201602
//var nestedAlternative = Global.ChoiceModelSession.Get<SchoolTourModeTimeModel>().RunNested(personDay, residenceParcel, person.UsualSchoolParcel, destinationArrivalTime, destinationDepartureTime, household.VehiclesAvailable);
ChoiceProbabilityCalculator.Alternative nestedAlternative = Global.ChoiceModelSession.Get <TourModeTimeModel>().RunNested(personDay, residenceParcel, person.UsualSchoolParcel, destinationArrivalTime, destinationDepartureTime, household.VehiclesAvailable, Global.Settings.Purposes.School);
purposeLogsums[Global.Settings.Purposes.School] = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
atUsualLogsums[Global.Settings.Purposes.School] = Global.AggregateLogsums[person.UsualSchoolParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votSegment][person.UsualSchoolParcel.TransitAccessSegment()];
}
double compositeLogsum = Global.AggregateLogsums[household.ResidenceZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votSegment][transitAccess];
purposeLogsums[Global.Settings.Purposes.Escort] = Global.AggregateLogsums[household.ResidenceZoneId][Global.Settings.Purposes.Escort][carOwnership][votSegment][transitAccess];
purposeLogsums[Global.Settings.Purposes.PersonalBusiness] = Global.AggregateLogsums[household.ResidenceZoneId][Global.Settings.Purposes.PersonalBusiness][carOwnership][votSegment][transitAccess];
purposeLogsums[Global.Settings.Purposes.Shopping] = Global.AggregateLogsums[household.ResidenceZoneId][Global.Settings.Purposes.Shopping][carOwnership][votSegment][transitAccess];
purposeLogsums[Global.Settings.Purposes.Social] = Global.AggregateLogsums[household.ResidenceZoneId][Global.Settings.Purposes.Social][carOwnership][votSegment][transitAccess];
for (int xPurpose = Global.Settings.Purposes.Escort; xPurpose <= Global.Settings.Purposes.Social + 10; xPurpose++)
{
// extra components 1-5 are for 2,3,4,5,6 tour purposes
// extra components 6-10 are for 2,3,4,5,6 stop puroposes
// recode purpose to match coefficients
int purpose = xPurpose <= Global.Settings.Purposes.Social ? xPurpose :
xPurpose <= Global.Settings.Purposes.Social + 5 ? Global.Settings.Purposes.Social + 1 :
Global.Settings.Purposes.Social + 2;
// get correct multiplier on coefficients.
double xMultiplier = xPurpose <= Global.Settings.Purposes.Social ? 1.0 :
xPurpose <= Global.Settings.Purposes.Social + 5 ? Math.Log(xPurpose - Global.Settings.Purposes.Social + 1) :
Math.Log(xPurpose - Global.Settings.Purposes.Social - 5 + 1);
choiceProbabilityCalculator.CreateUtilityComponent(xPurpose);
ChoiceProbabilityCalculator.Component component = choiceProbabilityCalculator.GetUtilityComponent(xPurpose);
component.AddUtilityTerm(100 * purpose + 1, xMultiplier * person.IsFulltimeWorker.ToFlag());
component.AddUtilityTerm(100 * purpose + 2, xMultiplier * person.IsPartTimeWorker.ToFlag());
component.AddUtilityTerm(100 * purpose + 3, xMultiplier * person.IsRetiredAdult.ToFlag());
component.AddUtilityTerm(100 * purpose + 4, xMultiplier * person.IsNonworkingAdult.ToFlag());
component.AddUtilityTerm(100 * purpose + 5, xMultiplier * person.IsUniversityStudent.ToFlag());
component.AddUtilityTerm(100 * purpose + 6, xMultiplier * person.IsDrivingAgeStudent.ToFlag());
component.AddUtilityTerm(100 * purpose + 7, xMultiplier * person.IsChildAge5Through15.ToFlag());
component.AddUtilityTerm(100 * purpose + 8, xMultiplier * person.IsChildUnder5.ToFlag());
component.AddUtilityTerm(100 * purpose + 9, xMultiplier * HHwithLowIncomeFlag);
component.AddUtilityTerm(100 * purpose + 10, xMultiplier * HHwithMidleIncomeFlag);
component.AddUtilityTerm(100 * purpose + 11, xMultiplier * HHwithHighIncomeFlag);
//component.AddUtilityTerm(100 * purpose + 12, xMultiplier * carsPerDriver);
component.AddUtilityTerm(100 * purpose + 12, xMultiplier * householdCars);
component.AddUtilityTerm(100 * purpose + 13, xMultiplier * person.IsOnlyAdult().ToFlag());
component.AddUtilityTerm(100 * purpose + 14, xMultiplier * person.IsOnlyFullOrPartTimeWorker().ToFlag());
component.AddUtilityTerm(100 * purpose + 15, xMultiplier * 0);
component.AddUtilityTerm(100 * purpose + 16, xMultiplier * person.IsFemale.ToFlag() * person.IsAdult.ToFlag() * (!household.HasChildrenUnder16).ToFlag());
component.AddUtilityTerm(100 * purpose + 17, xMultiplier * person.IsFemale.ToFlag() * person.IsAdult.ToFlag() * household.HasChildrenUnder5.ToFlag());
component.AddUtilityTerm(100 * purpose + 18, xMultiplier * person.IsFemale.ToFlag() * person.IsAdult.ToFlag() * household.HasChildrenAge5Through15.ToFlag());
component.AddUtilityTerm(100 * purpose + 19, xMultiplier * person.IsMale.ToFlag() * person.IsAdult.ToFlag() * household.HasChildrenUnder5.ToFlag());
component.AddUtilityTerm(100 * purpose + 20, xMultiplier * person.IsMale.ToFlag() * person.IsAdult.ToFlag() * household.HasChildrenAge5Through15.ToFlag());
//component.AddUtilityTerm(100 * purpose + 21, xMultiplier * primaryFamilyTimeFlag); //GV: wrong sign
//component.AddUtilityTerm(100 * purpose + 21, xMultiplier * person.AgeIsBetween18And25.ToFlag());
//component.AddUtilityTerm(100 * purpose + 22, xMultiplier * person.AgeIsBetween26And35.ToFlag());
//component.AddUtilityTerm(100 * purpose + 23, xMultiplier * person.AgeIsBetween51And65.ToFlag());
component.AddUtilityTerm(100 * purpose + 24, xMultiplier * person.WorksAtHome().ToFlag());
component.AddUtilityTerm(100 * purpose + 25, xMultiplier * mixedDensity);
component.AddUtilityTerm(100 * purpose + 26, xMultiplier * intersectionDensity);
//component.AddUtilityTerm(100 * purpose + 27, xMultiplier * purposeLogsums[purpose]); //GV: 17.08.2013, the logsums are wrong
//component.AddUtilityTerm(100 * purpose + 28, xMultiplier * person.TransitPassOwnershipFlag);
}
// tour utility
const int tourComponentIndex = 18;
choiceProbabilityCalculator.CreateUtilityComponent(tourComponentIndex);
ChoiceProbabilityCalculator.Component tourComponent = choiceProbabilityCalculator.GetUtilityComponent(tourComponentIndex);
//tourComponent.AddUtilityTerm(1701, carsPerDriver);
tourComponent.AddUtilityTerm(1701, householdCars);
tourComponent.AddUtilityTerm(1702, person.WorksAtHome().ToFlag());
tourComponent.AddUtilityTerm(1703, mixedDensity);
tourComponent.AddUtilityTerm(1704, mixedDensity * person.IsChildAge5Through15.ToFlag());
tourComponent.AddUtilityTerm(1705, compositeLogsum);
//tourComponent.AddUtilityTerm(1706, person.TransitPassOwnershipFlag);
tourComponent.AddUtilityTerm(1706, primaryFamilyTimeFlag);
// stop utility
const int stopComponentIndex = 19;
choiceProbabilityCalculator.CreateUtilityComponent(stopComponentIndex);
ChoiceProbabilityCalculator.Component stopComponent = choiceProbabilityCalculator.GetUtilityComponent(stopComponentIndex);
//stopComponent.AddUtilityTerm(1711, carsPerDriver);
stopComponent.AddUtilityTerm(1711, householdCars);
stopComponent.AddUtilityTerm(1712, person.WorksAtHome().ToFlag());
stopComponent.AddUtilityTerm(1713, mixedDensity);
stopComponent.AddUtilityTerm(1714, mixedDensity * person.IsChildAge5Through15.ToFlag());
stopComponent.AddUtilityTerm(1715, compositeLogsum);
//stopComponent.AddUtilityTerm(1716, person.TransitPassOwnershipFlag);
//stopComponent.AddUtilityTerm(1716, primaryFamilyTimeFlag); //GV: 17.08.2013, the logsums are wrong
for (int alternativeIndex = 0; alternativeIndex < TOTAL_ALTERNATIVES; alternativeIndex++)
{
DayPattern dayPattern = dayPatterns[alternativeIndex];
bool available = dayPattern.Available;
ChoiceProbabilityCalculator.Alternative alternative = choiceProbabilityCalculator.GetAlternative(alternativeIndex, available, choice != null && choice.Equals(dayPattern));
if (!Global.Configuration.IsInEstimationMode && !alternative.Available)
{
continue;
}
alternative.Choice = dayPattern;
// components for the purposes
for (int purpose = Global.Settings.Purposes.Escort; purpose <= Global.Settings.Purposes.Social; purpose++)
{
if (dayPattern.Tours[purpose] > 0 || dayPattern.Stops[purpose] > 0)
{
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(purpose));
if (dayPattern.Tours[purpose] > 0)
{
alternative.AddUtilityTerm(100 * purpose + 50, 1); // tour purpose ASC
//alternative.AddUtilityTerm(100 * purpose + 51, purposeLogsums[purpose]); // tour purpose logsum GV: 17.08.2013, the logsums are wrong
}
if (dayPattern.Stops[purpose] > 0)
{
alternative.AddUtilityTerm(100 * purpose + 60, 1); // stop purpose ASC
//alternative.AddUtilityTerm(100 * purpose + 61, purposeLogsums[purpose]); // stop purpose logsum GV: 17.08.2013, the logsums are wrong
}
if (Global.Configuration.IsInEstimationMode)
{
//GV commented out
//alternative.AddUtilityTerm(100 * purpose + 70, 1 - person.PaperDiary);
//GV commented out
//alternative.AddUtilityTerm(100 * purpose + 71, person.ProxyResponse);
}
}
}
// multiple tour purposes component
if (dayPattern.TotalTourPurposes > 1)
{
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(Global.Settings.Purposes.Social + (dayPattern.TotalTourPurposes - 1)));
}
// multiple stop purposes component
if (dayPattern.TotalStopPurposes > 1)
{
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(Global.Settings.Purposes.Social + 5 + (dayPattern.TotalStopPurposes - 1)));
}
for (int tourPurpose = Global.Settings.Purposes.Work; tourPurpose <= Global.Settings.Purposes.Social; tourPurpose++)
{
for (int stopPurpose = Global.Settings.Purposes.Work; stopPurpose <= Global.Settings.Purposes.Social - 1; stopPurpose++)
{
if (tourPurpose > Global.Settings.Purposes.School && stopPurpose <= Global.Settings.Purposes.School)
{
continue;
}
if (dayPattern.Tours[tourPurpose] > 0 && dayPattern.Stops[stopPurpose] > 0)
{
alternative.AddUtilityTerm(1200 + 10 * tourPurpose + stopPurpose, 1); // tour-stop comb. utility
}
}
}
for (int tourPurpose = Global.Settings.Purposes.Work; tourPurpose <= Global.Settings.Purposes.School; tourPurpose++)
{
if (dayPattern.Tours[tourPurpose] == 1 && dayPattern.TotalStopPurposes >= 1)
{
alternative.AddUtilityTerm(1300 + 10 * tourPurpose, purposeLogsums[tourPurpose]); // usual location logsum x presence of stops in work or school pattern
alternative.AddUtilityTerm(1300 + 10 * tourPurpose + 1, compositeLogsum); // home aggregate logsum x presence of stops in work or school pattern
//alternative.AddUtilityTerm(1300 + 10 * tourPurpose + 2, atUsualLogsums[tourPurpose]); // at usual location aggregate logsum x presence of stops in work or school pattern GV: commented out as the sign is wrong
}
}
for (int tourPurpose = Global.Settings.Purposes.Work; tourPurpose <= Global.Settings.Purposes.Social - 2; tourPurpose++)
{
for (int tourPurpose2 = tourPurpose + 1; tourPurpose2 <= Global.Settings.Purposes.Social; tourPurpose2++)
{
if (dayPattern.Tours[tourPurpose] > 0 && dayPattern.Tours[tourPurpose2] > 0)
{
alternative.AddUtilityTerm(1400 + 10 * tourPurpose + tourPurpose2, 1); // tour-tour comb. utility
}
}
}
for (int tourPurpose = Global.Settings.Purposes.Business; tourPurpose <= Global.Settings.Purposes.Business; tourPurpose++)
{
for (int tourPurpose2 = Global.Settings.Purposes.Escort; tourPurpose2 <= Global.Settings.Purposes.Social; tourPurpose2++)
{
if (dayPattern.Tours[tourPurpose] > 0 && dayPattern.Tours[tourPurpose2] > 0)
{
alternative.AddUtilityTerm(1461, 1); // tour-tour comb. utility for business combos
}
}
}
for (int stopPurpose = Global.Settings.Purposes.Work; stopPurpose <= Global.Settings.Purposes.Social - 2; stopPurpose++)
{
for (int stopPurpose2 = stopPurpose + 1; stopPurpose2 <= Global.Settings.Purposes.Social; stopPurpose2++)
{
if (dayPattern.Stops[stopPurpose] > 0 && dayPattern.Stops[stopPurpose2] > 0)
{
alternative.AddUtilityTerm(1500 + 10 * stopPurpose + stopPurpose2, 1); // stop-stop comb. utility
}
}
}
if (dayPattern.TotalTourPurposes > 0 && dayPattern.TotalStopPurposes > 0)
{
int totalStopPurposes = dayPattern.TotalStopPurposes;
if (totalStopPurposes > 3)
{
totalStopPurposes = 3;
}
alternative.AddUtilityTerm(1600 + 10 * dayPattern.TotalTourPurposes + totalStopPurposes, 1); // nttour-ntstop utility
}
if (dayPattern.TotalTourPurposes - dayPattern.Tours[Global.Settings.Purposes.Work]
- dayPattern.Tours[Global.Settings.Purposes.Business]
- dayPattern.Tours[Global.Settings.Purposes.School] > 0)
{
// tour utility
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(tourComponentIndex));
}
if (dayPattern.TotalStopPurposes - dayPattern.Stops[Global.Settings.Purposes.Work]
- dayPattern.Stops[Global.Settings.Purposes.Business]
- dayPattern.Stops[Global.Settings.Purposes.School] > 0)
{
// stop utility
alternative.AddUtilityComponent(choiceProbabilityCalculator.GetUtilityComponent(stopComponentIndex));
}
}
}