public void Run(PersonWrapper person, int sampleSize)
{
if (person == null)
{
throw new ArgumentNullException("person");
}
person.ResetRandom(0);
if (Global.Configuration.IsInEstimationMode)
{
if (!_helpers[ParallelUtility.GetBatchFromThreadId()].ModelIsInEstimationMode)
{
return;
}
}
var choiceProbabilityCalculator = _helpers[ParallelUtility.GetBatchFromThreadId()].GetChoiceProbabilityCalculator(person.Id);
if (_helpers[ParallelUtility.GetBatchFromThreadId()].ModelIsInEstimationMode)
{
if (person.UsualWorkParcel == null)
{
return;
}
var choseHome = person.UsualWorkParcelId == person.Household.ResidenceParcelId; // JLB 20120329 added these two lines
var chosenParcel = choseHome ? null : person.UsualWorkParcel;
//RunModel(choiceProbabilityCalculator, person, sampleSize, person.UsualWorkParcel);
RunModel(choiceProbabilityCalculator, person, sampleSize, chosenParcel, choseHome); // JLB 20120329 replaced above line
// when chosenParcel is null:
// DestinationSampler doesn't try to assign one of the sampled destinations as chosen
// choseHome is NOT null, and RunModel sets the oddball location as chosen
choiceProbabilityCalculator.WriteObservation();
}
else
{
RunModel(choiceProbabilityCalculator, person, sampleSize);
var chosenAlternative = choiceProbabilityCalculator.SimulateChoice(person.Household.RandomUtility);
var choice = (ParcelWrapper)chosenAlternative.Choice;
person.UsualWorkParcelId = choice.Id;
person.UsualWorkParcel = choice;
person.UsualWorkZoneKey = ChoiceModelFactory.ZoneKeys[choice.ZoneId];
var skimValue = ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Sov, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, 1, person.Household.ResidenceParcel, choice);
person.AutoTimeToUsualWork = skimValue.Variable;
person.AutoDistanceToUsualWork = skimValue.BlendVariable;
person.SetWorkParcelPredictions();
}
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, PersonWrapper person, int sampleSize, IParcelWrapper choice = null, bool choseHome = false)
{
var segment = Global.Kernel.Get <SamplingWeightsSettingsFactory>().SamplingWeightsSettings.GetTourDestinationSegment(Global.Settings.Purposes.Work, Global.Settings.TourPriorities.HomeBasedTour, Global.Settings.Modes.Sov, person.PersonType);
var destinationSampler = new DestinationSampler(choiceProbabilityCalculator, segment, sampleSize, person.Household.ResidenceParcel, choice);
var destinationArrivalTime = ChoiceModelUtility.GetDestinationArrivalTime(Global.Settings.Models.WorkTourModeModel);
var destinationDepartureTime = ChoiceModelUtility.GetDestinationDepartureTime(Global.Settings.Models.WorkTourModeModel);
var workLocationUtilites = new WorkLocationUtilities(person, sampleSize, destinationArrivalTime, destinationDepartureTime);
var sampleItems = destinationSampler.SampleAndReturnTourDestinations(workLocationUtilites);
int index = 0;
foreach (var sampleItem in sampleItems)
{
bool available = sampleItem.Key.Available;
bool isChosen = sampleItem.Key.IsChosen;
double adjustmentFactor = sampleItem.Key.AdjustmentFactor;
var destinationParcel = ChoiceModelFactory.Parcels[sampleItem.Key.ParcelId];
var alternative = choiceProbabilityCalculator.GetAlternative(index++, available, isChosen);
if (!available)
{
continue;
}
alternative.Choice = destinationParcel;
var workTourLogsum = 0D;
var nestedAlternative = Global.ChoiceModelSession.Get <WorkTourModeTimeModel>().RunNested(person, person.Household.ResidenceParcel, destinationParcel, destinationArrivalTime, destinationDepartureTime, person.Household.HouseholdTotals.DrivingAgeMembers, 0.0);
workTourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
var votSegment = person.Household.GetVotALSegment();
var taSegment = destinationParcel.TransitAccessSegment();
var aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][Global.Settings.CarOwnerships.OneOrMoreCarsPerAdult][votSegment][taSegment];
var distanceFromOrigin = person.Household.ResidenceParcel.DistanceFromOrigin(destinationParcel, 1);
var distance1 = Math.Min(distanceFromOrigin, .35);
var distance2 = Math.Max(0, Math.Min(distanceFromOrigin - .35, 1 - .35));
var distance3 = Math.Max(0, distanceFromOrigin - 1);
var distanceLog = Math.Log(1 + distanceFromOrigin);
var distanceFromSchool = person.IsFullOrPartTimeWorker ? 0 : person.UsualSchoolParcel.DistanceFromSchoolLog(destinationParcel, 1);
// parcel buffers
var educationBuffer = Math.Log(destinationParcel.EmploymentEducationBuffer2 + 1);
var governmentBuffer = Math.Log(destinationParcel.EmploymentGovernmentBuffer2 + 1);
var officeBuffer = Math.Log(destinationParcel.EmploymentOfficeBuffer2 + 1);
var serviceBuffer = Math.Log(destinationParcel.EmploymentServiceBuffer2 + 1);
var householdsBuffer = Math.Log(destinationParcel.HouseholdsBuffer2 + 1);
// var retailBuffer = Math.Log(destinationParcel.EmploymentRetailBuffer2 + 1);
var industrialAgricultureConstructionBuffer = Math.Log(destinationParcel.EmploymentIndustrialBuffer2 + destinationParcel.EmploymentAgricultureConstructionBuffer2 + 1);
var foodBuffer = Math.Log(destinationParcel.EmploymentFoodBuffer2 + 1);
var medicalBuffer = Math.Log(destinationParcel.EmploymentMedicalBuffer2 + 1);
var employmentTotalBuffer = Math.Log(destinationParcel.EmploymentTotalBuffer2 + 1);
var studentsUniversityBuffer = Math.Log(destinationParcel.StudentsUniversityBuffer2 + 1);
var studentsK12Buffer = Math.Log(destinationParcel.StudentsK8Buffer2 + destinationParcel.StudentsHighSchoolBuffer2 + 1);
// var mixedUse4Index = destinationParcel.MixedUse4Index2();
//size attributes (derived)
var employmentIndustrialAgricultureConstruction = destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction;
// parking attributes
var parcelParkingDensity = destinationParcel.ParcelParkingPerTotalEmployment();
// connectivity attributes
var c34Ratio = destinationParcel.C34RatioBuffer1();
// Stefan
bool isInCopenhagenMunicipality = true; //destinationParcel.Municipality == 101; Need to change this after Municipality property is added to LD parcel file
var employmentCommercial = destinationParcel.EmploymentRetail + destinationParcel.EmploymentService;
var employmentCommercialBuffer1 = destinationParcel.EmploymentRetailBuffer1 + destinationParcel.EmploymentServiceBuffer1;
double beta00002 = -2.53;
double beta00003 = 2.65;
double beta00004 = 1.57;
double beta00005 = -0.18;
double beta00006 = -0.43;
double beta00007 = -0.19;
double beta00008 = 0.33;
double beta00009 = 0.007;
double beta00010 = 0.55;
var stefanUtility =
beta00002 * destinationParcel.Households / destinationParcel.ThousandsSquareLengthUnits +
beta00003 * (person.Household.Income < 300000).ToFlag() * destinationParcel.Households / destinationParcel.ThousandsSquareLengthUnits +
beta00004 * person.IsFemale.ToFlag() * destinationParcel.Households / destinationParcel.ThousandsSquareLengthUnits +
beta00005 * isInCopenhagenMunicipality.ToFlag() +
beta00006 * (person.Household.HasValidIncome && person.Household.Income < 300000).ToFlag() * isInCopenhagenMunicipality.ToFlag() +
beta00007 * (person.Household.HasValidIncome && person.Household.Income >= 300000 && person.Household.Income < 600000).ToFlag() * isInCopenhagenMunicipality.ToFlag() +
beta00008 * (person.Household.HasValidIncome && person.Household.Income >= 900000).ToFlag() * isInCopenhagenMunicipality.ToFlag() +
beta00009 * person.Age * isInCopenhagenMunicipality.ToFlag() +
0.0; // beta00010 * (person.Household.ResidenceParcel.Municipality == destination.Municipality).ToFlag();
//Stefan non-size terms.
alternative.AddUtilityTerm(1, sampleItem.Key.AdjustmentFactor);
//alternative.AddUtilityTerm(2, destinationParcel.Households / destinationParcel.ThousandsSquareLengthUnits);
//alternative.AddUtilityTerm(3, (person.Household.Income < 300000).ToFlag() * destinationParcel.Households / destinationParcel.ThousandsSquareLengthUnits);
//alternative.AddUtilityTerm(4, person.IsFemale.ToFlag() * destinationParcel.Households / destinationParcel.ThousandsSquareLengthUnits);
//alternative.AddUtilityTerm(5, isInCopenhagenMunicipality.ToFlag());
//alternative.AddUtilityTerm(6, (person.Household.HasValidIncome && person.Household.Income < 300000).ToFlag() * isInCopenhagenMunicipality.ToFlag());
//alternative.AddUtilityTerm(7, (person.Household.HasValidIncome && person.Household.Income >= 300000 && person.Household.Income < 600000).ToFlag() * isInCopenhagenMunicipality.ToFlag());
//alternative.AddUtilityTerm(8, (person.Household.HasValidIncome && person.Household.Income >= 900000).ToFlag() * isInCopenhagenMunicipality.ToFlag());
//alternative.AddUtilityTerm(9, person.Age * isInCopenhagenMunicipality.ToFlag());
////alternative.AddUtilityTerm(10, (person.Household.ResidenceParcel.Municipality == destination.Municipality).ToFlag()); // Acivate this after Municipality property is added to LD parcel file
//following logsums replace Stefan's car and public transport times
alternative.AddUtilityTerm(11, (person.Household.HasValidIncome && person.Household.Income < 300000).ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(12, (person.Household.HasValidIncome && person.Household.Income >= 300000 && person.Household.Income < 600000).ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(13, (person.Household.HasValidIncome && person.Household.Income >= 900000).ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(14, person.Household.HasMissingIncome.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(15, person.IsFemale.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(16, person.Age * workTourLogsum);
alternative.AddUtilityTerm(17, (person.MainOccupation == 50).ToFlag() * workTourLogsum); // self-employed
//Stefan's composite term 18 replaces terms 2-10 above
alternative.AddUtilityTerm(18, stefanUtility); // see above for this composite function of StefanMabitt's utility function
//alternative.AddUtilityTerm(2, person.IsFulltimeWorker.ToFlag() * workTourLogsum);
//alternative.AddUtilityTerm(3, person.IsPartTimeWorker.ToFlag() * workTourLogsum);
//alternative.AddUtilityTerm(4, person.IsNotFullOrPartTimeWorker.ToFlag() * workTourLogsum);
//alternative.AddUtilityTerm(5, distanceLog); // for distance calibration
//alternative.AddUtilityTerm(6, person.IsFulltimeWorker.ToFlag() * distance1);
//alternative.AddUtilityTerm(7, person.IsFulltimeWorker.ToFlag() * distance2);
//alternative.AddUtilityTerm(8, person.IsFulltimeWorker.ToFlag() * distance3);
//alternative.AddUtilityTerm(9, person.IsPartTimeWorker.ToFlag() * distanceLog);
//alternative.AddUtilityTerm(10, person.IsNotFullOrPartTimeWorker.ToFlag() * distanceLog);
//alternative.AddUtilityTerm(11, person.Household.Has0To15KIncome.ToFlag() * distanceLog);
//alternative.AddUtilityTerm(12, person.Household.Has50To75KIncome.ToFlag() * distanceLog);
//alternative.AddUtilityTerm(13, person.Household.Has75To100KIncome.ToFlag() * distanceLog);
//alternative.AddUtilityTerm(14, person.IsFemale.ToFlag() * distanceLog);
//alternative.AddUtilityTerm(15, person.IsStudentAge.ToFlag() * distanceFromSchool);
//alternative.AddUtilityTerm(16, person.IsFulltimeWorker.ToFlag() * aggregateLogsum);
//alternative.AddUtilityTerm(17, person.IsPartTimeWorker.ToFlag() * aggregateLogsum);
//alternative.AddUtilityTerm(18, person.IsNotFullOrPartTimeWorker.ToFlag() * aggregateLogsum);
//alternative.AddUtilityTerm(19, parcelParkingDensity);
//alternative.AddUtilityTerm(20, c34Ratio);
//Neighborhood
//alternative.AddUtilityTerm(21, person.Household.HasValidIncome.ToFlag() * serviceBuffer);
//alternative.AddUtilityTerm(22, person.Household.HasValidIncome.ToFlag() * educationBuffer);
//alternative.AddUtilityTerm(23, person.Household.HasValidIncome.ToFlag() * foodBuffer);
//alternative.AddUtilityTerm(24, person.Household.HasValidIncome.ToFlag() * governmentBuffer);
//alternative.AddUtilityTerm(25, person.Household.HasValidIncome.ToFlag() * officeBuffer);
//alternative.AddUtilityTerm(26, person.Household.HasValidIncome.ToFlag() * medicalBuffer);
//alternative.AddUtilityTerm(27, person.Household.HasValidIncome.ToFlag() * householdsBuffer);
//alternative.AddUtilityTerm(28, person.Household.HasValidIncome.ToFlag() * studentsUniversityBuffer);
//alternative.AddUtilityTerm(29, person.Household.HasValidIncome.ToFlag() * person.IsFulltimeWorker.ToFlag() * studentsK12Buffer);
//alternative.AddUtilityTerm(30, person.Household.HasValidIncome.ToFlag() * person.IsFulltimeWorker.ToFlag() * studentsUniversityBuffer);
//alternative.AddUtilityTerm(31, person.Household.HasValidIncome.ToFlag() * person.IsPartTimeWorker.ToFlag() * industrialAgricultureConstructionBuffer);
//alternative.AddUtilityTerm(32, person.Household.HasValidIncome.ToFlag() * person.IsNotFullOrPartTimeWorker.ToFlag() * foodBuffer);
//alternative.AddUtilityTerm(33, person.Household.HasValidIncome.ToFlag() * person.IsNotFullOrPartTimeWorker.ToFlag() * medicalBuffer);
//alternative.AddUtilityTerm(34, person.IsFulltimeWorker.ToFlag() * person.Household.Has75KPlusIncome.ToFlag() * employmentTotalBuffer);
//alternative.AddUtilityTerm(35, person.IsNotFullOrPartTimeWorker.ToFlag() * person.Household.HasIncomeUnder50K.ToFlag() * governmentBuffer);
//alternative.AddUtilityTerm(36, person.IsNotFullOrPartTimeWorker.ToFlag() * person.Household.HasIncomeUnder50K.ToFlag() * employmentTotalBuffer);
//Size
// Stefan size terms.
// Note: the following assumes: (1) Stefan's size variables enter his utility function as a logsum, a la BAL;
// (2) Jobs--commercial and Jobs--finance apply to hh w unknown incomes in Stefan's spec
// If his size variables enter linearly, then if I want to replicate them I should not use alogit size functions.
// If Jobs--commercial and Jobs--finance don't apply to missing incomes, then I need to change the size functions below
alternative.AddUtilityTerm(51, (person.Household.HasValidIncome && person.Household.Income < 300000).ToFlag() * (destinationParcel.EmploymentTotal - employmentCommercial));
alternative.AddUtilityTerm(52, (person.Household.HasValidIncome && person.Household.Income >= 300000 && person.Household.Income < 600000).ToFlag() * (destinationParcel.EmploymentTotal - employmentCommercial));
alternative.AddUtilityTerm(53, (person.Household.HasValidIncome && person.Household.Income >= 600000 && person.Household.Income < 900000).ToFlag() * (destinationParcel.EmploymentTotal - employmentCommercial));
alternative.AddUtilityTerm(54, (person.Household.HasValidIncome && person.Household.Income >= 900000).ToFlag() * (destinationParcel.EmploymentTotal - employmentCommercial));
alternative.AddUtilityTerm(55, (person.Household.HasValidIncome && person.Household.Income >= 900000).ToFlag() * (destinationParcel.EmploymentTotal - employmentCommercial)); // second term allows first one to have base coef of 0
alternative.AddUtilityTerm(56, (person.Household.HasMissingIncome.ToFlag() * (destinationParcel.EmploymentTotal - employmentCommercial)));
alternative.AddUtilityTerm(57, (person.Household.HasValidIncome && person.Household.Income < 300000).ToFlag() * employmentCommercial);
alternative.AddUtilityTerm(58, (person.Household.HasValidIncome && person.Household.Income >= 300000 && person.Household.Income < 600000).ToFlag() * employmentCommercial);
alternative.AddUtilityTerm(59, (person.Household.HasValidIncome && person.Household.Income >= 600000 && person.Household.Income < 900000).ToFlag() * employmentCommercial);
alternative.AddUtilityTerm(60, (person.Household.HasValidIncome && person.Household.Income >= 900000).ToFlag() * employmentCommercial);
alternative.AddUtilityTerm(61, person.Household.HasMissingIncome.ToFlag() * employmentCommercial);
alternative.AddUtilityTerm(62, destinationParcel.EmploymentOffice);
//The following combine with 51-55, 56-60 and 61 to include size of entire buffer region in main size variables
alternative.AddUtilityTerm(63, (destinationParcel.EmploymentTotalBuffer1 - destinationParcel.EmploymentTotal) - (employmentCommercialBuffer1 - employmentCommercial));
alternative.AddUtilityTerm(64, employmentCommercialBuffer1 - employmentCommercial);
alternative.AddUtilityTerm(65, destinationParcel.EmploymentOfficeBuffer1 - destinationParcel.EmploymentOffice);
//alternative.AddUtilityTerm(51, person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentService);
//alternative.AddUtilityTerm(52, person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentEducation);
//alternative.AddUtilityTerm(53, person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentFood);
//alternative.AddUtilityTerm(54, person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentGovernment);
//alternative.AddUtilityTerm(55, person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentOffice);
//alternative.AddUtilityTerm(56, person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentRetail);
//alternative.AddUtilityTerm(57, person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentMedical);
//alternative.AddUtilityTerm(58, person.Household.HasValidIncome.ToFlag() * employmentIndustrialAgricultureConstruction);
//alternative.AddUtilityTerm(59, person.Household.HasValidIncome.ToFlag() * destinationParcel.StudentsUniversity);
//alternative.AddUtilityTerm(60, person.Household.HasValidIncome.ToFlag() * person.IsFulltimeWorker.ToFlag() * destinationParcel.EmploymentGovernment);
//alternative.AddUtilityTerm(61, person.Household.HasValidIncome.ToFlag() * person.IsFulltimeWorker.ToFlag() * employmentIndustrialAgricultureConstruction);
//alternative.AddUtilityTerm(62, person.Household.HasValidIncome.ToFlag() * person.IsPartTimeWorker.ToFlag() * employmentIndustrialAgricultureConstruction);
//alternative.AddUtilityTerm(63, person.Household.HasValidIncome.ToFlag() * person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentEducation);
//alternative.AddUtilityTerm(64, person.Household.HasValidIncome.ToFlag() * person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentFood);
//alternative.AddUtilityTerm(65, person.Household.HasValidIncome.ToFlag() * person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentRetail);
//alternative.AddUtilityTerm(66, person.Household.HasIncomeUnder50K.ToFlag() * destinationParcel.EmploymentRetail);
//alternative.AddUtilityTerm(67, person.Household.HasIncomeUnder50K.ToFlag() * destinationParcel.EmploymentService);
//alternative.AddUtilityTerm(68, person.Household.Has50To75KIncome.ToFlag() * destinationParcel.EmploymentMedical);
//alternative.AddUtilityTerm(69, person.Household.Has50To75KIncome.ToFlag() * destinationParcel.EmploymentOffice);
//alternative.AddUtilityTerm(70, person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentEducation);
//alternative.AddUtilityTerm(71, person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentGovernment);
//alternative.AddUtilityTerm(72, person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentMedical);
//alternative.AddUtilityTerm(73, person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentOffice);
//alternative.AddUtilityTerm(74, person.IsFulltimeWorker.ToFlag() * person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentGovernment);
//alternative.AddUtilityTerm(75, person.IsFulltimeWorker.ToFlag() * (!person.Household.Has75KPlusIncome).ToFlag() * employmentIndustrialAgricultureConstruction);
//alternative.AddUtilityTerm(76, person.IsPartTimeWorker.ToFlag() * (!person.Household.HasIncomeUnder50K).ToFlag() * destinationParcel.EmploymentMedical);
//alternative.AddUtilityTerm(77, (!person.IsFulltimeWorker).ToFlag() * person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentOffice);
//alternative.AddUtilityTerm(78, person.IsNotFullOrPartTimeWorker.ToFlag() * (!person.Household.HasIncomeUnder50K).ToFlag() * destinationParcel.EmploymentRetail);
//alternative.AddUtilityTerm(79, person.Household.HasMissingIncome.ToFlag() * destinationParcel.EmploymentTotal);
//alternative.AddUtilityTerm(80, person.Household.HasMissingIncome.ToFlag() * destinationParcel.StudentsUniversity);
// set shadow price depending on persontype and add it to utility
// we are using the sampling adjustment factor assuming that it is 1
alternative.AddUtilityTerm(1, destinationParcel.ShadowPriceForEmployment);
//remove nesting for estimation of conditional MNL
alternative.AddNestedAlternative(sampleSize + 2, 0, 98);
}
// JLB 20120329 added third call parameter to idenitfy whether this alt is chosen or not
var homeAlternative = choiceProbabilityCalculator.GetAlternative(sampleSize, true, choseHome);
homeAlternative.Choice = person.Household.ResidenceParcel;
homeAlternative.AddUtilityTerm(41, 1);
homeAlternative.AddUtilityTerm(42, (person.MainOccupation == 50).ToFlag()); // self-employed
//homeAlternative.AddUtilityTerm(42, person.IsPartTimeWorker.ToFlag());
//homeAlternative.AddUtilityTerm(43, person.IsStudentAge.ToFlag());
//homeAlternative.AddUtilityTerm(44, person.IsFemale.ToFlag());
homeAlternative.AddUtilityTerm(90, 1);
//make oddball alt unavailable and remove nesting for estimation of conditional MNL
// alternative.Available = false;
homeAlternative.AddNestedAlternative(sampleSize + 3, 1, 98);
}
private void RunModel(ChoiceProbabilityCalculator choiceProbabilityCalculator, TourWrapper tour, HouseholdDayWrapper householdDay, int sampleSize, IParcelWrapper choice = null)
{
timesStartedRunModel++;
HouseholdWrapper household = (HouseholdWrapper)tour.Household;
PersonWrapper person = (PersonWrapper)tour.Person;
PersonDayWrapper personDay = (PersonDayWrapper)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)
{
TourWrapper tInJoint = (TourWrapper)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 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 primaryFlag = ChoiceModelUtility.GetPrimaryFlag(tourCategory);
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);
if (tour.Household.Id == 80049 && tour.PersonDay.Day == 1 && tour.Person.Sequence == 2 && tour.Sequence == 4)
{
bool testbreak = true;
}
var segment = Global.Kernel.Get <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, tour.OriginParcel, choice);
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 index = 0;
foreach (var sampleItem in sampleItems)
{
bool available = sampleItem.Key.Available;
bool isChosen = sampleItem.Key.IsChosen;
double adjustmentFactor = sampleItem.Key.AdjustmentFactor;
var destinationParcel = ChoiceModelFactory.Parcels[sampleItem.Key.ParcelId];
if (isChosen)
{
Global.PrintFile.WriteLine("Sequence {0}: Chosen parcel {1} Available {2} Sample item {3} of {4}", timesStartedRunModel, destinationParcel.Id, available, index, sampleItems.Count);
}
var alternative = choiceProbabilityCalculator.GetAlternative(index++, available, isChosen);
if (!available)
{
continue;
}
var fastestTravelTime =
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, fastestAvailableTimeOfDay, tour.OriginParcel, destinationParcel).Variable +
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, fastestAvailableTimeOfDay, destinationParcel, tour.OriginParcel).Variable;
if (fastestTravelTime >= maxAvailableMinutes)
{
alternative.Available = false;
continue;
}
alternative.Choice = destinationParcel;
double tourLogsum;
if (tour.IsHomeBasedTour)
{
if (tour.DestinationPurpose == Global.Settings.Purposes.Work)
{
//JLB 201406
//var nestedAlternative = Global.ChoiceModelSession.Get<WorkTourModeModel>().RunNested(tour, destinationParcel);
var nestedAlternative = Global.ChoiceModelSession.Get <WorkTourModeTimeModel>().RunNested(tour, destinationParcel, tour.Household.VehiclesAvailable, tour.Person.GetTransitFareDiscountFraction());
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
// JLB201406
//else if (tour.DestinationPurpose == Global.Settings.Purposes.Escort) {
// var nestedAlternative = Global.ChoiceModelSession.Get<EscortTourModeModel>().RunNested(tour, destinationParcel);
// tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
//}
else
{
// JLB201406
//var nestedAlternative = Global.ChoiceModelSession.Get<OtherHomeBasedTourModeModel>().RunNested(tour, destinationParcel);
var nestedAlternative = Global.ChoiceModelSession.Get <OtherHomeBasedTourModeTimeModel>().RunNested(tour, destinationParcel, tour.Household.VehiclesAvailable, tour.Person.GetTransitFareDiscountFraction());
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
}
else
{
// JLB201406
//var nestedAlternative = Global.ChoiceModelSession.Get<WorkBasedSubtourModeModel>().RunNested(tour, destinationParcel);
var nestedAlternative = Global.ChoiceModelSession.Get <WorkBasedSubtourModeTimeModel>().RunNested(tour, destinationParcel, tour.Household.VehiclesAvailable, tour.Person.GetTransitFareDiscountFraction());
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
//var purpose = tour.TourPurposeSegment;
var carOwnership = person.GetCarOwnershipSegment();
var votSegment = tour.GetVotALSegment();
var transitAccess = destinationParcel.TransitAccessSegment();
//var aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][purpose][carOwnership][votSegment][transitAccess];
var aggregateLogsumHomeBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votSegment][transitAccess];
var aggregateLogsumWorkBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.WorkBased][carOwnership][votSegment][transitAccess];
var distanceFromOrigin = tour.OriginParcel.DistanceFromOrigin(destinationParcel, tour.DestinationArrivalTime);
// 1. new from GV: Cph KM-distances
var piecewiseDistanceFrom0To1Km = Math.Min(distanceFromOrigin, .10);
var piecewiseDistanceFrom0To2Km = Math.Min(distanceFromOrigin, .20); //GV: added July 7th
var piecewiseDistanceFrom0To5Km = Math.Min(distanceFromOrigin, .50); //GV: added July 7th
var piecewiseDistanceFrom1To2Km = Math.Max(0, Math.Min(distanceFromOrigin - .1, .2 - .1));
var piecewiseDistanceFrom2To5Km = Math.Max(0, Math.Min(distanceFromOrigin - .2, .5 - .2));
var piecewiseDistanceFrom5To10Km = Math.Max(0, Math.Min(distanceFromOrigin - .5, 1 - .5));
var piecewiseDistanceFrom10To20Km = Math.Max(0, Math.Min(distanceFromOrigin - 1, 2 - 1));
var piecewiseDistanceFrom20KmToInfinity = Math.Max(0, distanceFromOrigin - 2);
var piecewiseDistanceFrom10KmToInfinity = Math.Max(0, distanceFromOrigin - 1);
// 1. finished
var distanceFromOriginLog = Math.Log(1 + distanceFromOrigin);
var distanceFromWorkLog = person.UsualWorkParcel.DistanceFromWorkLog(destinationParcel, 1);
var distanceFromSchoolLog = person.UsualSchoolParcel.DistanceFromSchoolLog(destinationParcel, 1);
var timePressure = Math.Log(1 - fastestTravelTime / maxAvailableMinutes);
// 2. new from GV: Cph buffers for neighborhood effects
// log transforms of buffers for Neighborhood effects
var logOfOnePlusEducationK8Buffer2 = Math.Log(destinationParcel.StudentsK8Buffer2 + 1.0);
var logOfOnePlusEducationUniStuBuffer2 = Math.Log(destinationParcel.StudentsUniversityBuffer2 + 1.0);
var logOfOnePlusEmploymentEducationBuffer2 = Math.Log(destinationParcel.EmploymentEducationBuffer2 + 1.0);
var logOfOnePlusEmploymentGovernmentBuffer2 = Math.Log(destinationParcel.EmploymentGovernmentBuffer2 + 1.0);
var logOfOnePlusEmploymentIndustrialBuffer2 = Math.Log(destinationParcel.EmploymentIndustrialBuffer2 + 1.0);
var logOfOnePlusEmploymentOfficeBuffer2 = Math.Log(destinationParcel.EmploymentOfficeBuffer2 + 1.0);
var logOfOnePlusEmploymentRetailBuffer2 = Math.Log(destinationParcel.EmploymentRetailBuffer2 + 1.0);
var logOfOnePlusEmploymentServiceBuffer2 = Math.Log(destinationParcel.EmploymentServiceBuffer2 + 1.0);
var logOfOnePlusEmploymentAgrConstrBuffer2 = Math.Log(destinationParcel.EmploymentAgricultureConstructionBuffer2 + 1.0);
var logOfOnePlusEmploymentJobsBuffer2 = Math.Log(destinationParcel.EmploymentTotalBuffer2 + 1.0);
var logOfOnePlusHouseholdsBuffer2 = Math.Log(destinationParcel.HouseholdsBuffer2 + 1.0);
// 2. finished
var logOfOnePlusParkingOffStreetDailySpacesBuffer1 = Math.Log(1 + destinationParcel.ParkingOffStreetPaidDailySpacesBuffer1);
// connectivity attributes
var c34Ratio = destinationParcel.C34RatioBuffer1();
var carCompetitionFlag = FlagUtility.GetCarCompetitionFlag(carOwnership); // exludes no cars
var noCarCompetitionFlag = FlagUtility.GetNoCarCompetitionFlag(carOwnership);
var noCarsFlag = FlagUtility.GetNoCarsFlag(carOwnership);
alternative.AddUtilityTerm(2, household.Id);
alternative.AddUtilityTerm(3, personDay.Day);
alternative.AddUtilityTerm(4, person.Sequence);
alternative.AddUtilityTerm(5, tour.Sequence);
alternative.AddUtilityTerm(8, adjustmentFactor);
alternative.AddUtilityTerm(9, tourLogsum);
// 3. new from GV: definition of Cph variables
//alternative.AddUtilityTerm(260, secondaryFlag * workOrSchoolPatternFlag * piecewiseDistanceFrom0To1Km);
alternative.AddUtilityTerm(260, secondaryFlag * workOrSchoolPatternFlag * piecewiseDistanceFrom0To2Km); //GV: added July 7th
//alternative.AddUtilityTerm(261, secondaryFlag * workOrSchoolPatternFlag * piecewiseDistanceFrom1To2Km);
alternative.AddUtilityTerm(262, secondaryFlag * workOrSchoolPatternFlag * piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(263, secondaryFlag * workOrSchoolPatternFlag * piecewiseDistanceFrom5To10Km);
alternative.AddUtilityTerm(264, secondaryFlag * workOrSchoolPatternFlag * piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(265, secondaryFlag * workOrSchoolPatternFlag * piecewiseDistanceFrom20KmToInfinity);
//alternative.AddUtilityTerm(266, secondaryFlag * otherPatternFlag * piecewiseDistanceFrom0To1Km);
alternative.AddUtilityTerm(266, secondaryFlag * otherPatternFlag * piecewiseDistanceFrom0To2Km); //GV: added July 7th
//alternative.AddUtilityTerm(267, secondaryFlag * otherPatternFlag * piecewiseDistanceFrom1To2Km);
alternative.AddUtilityTerm(268, secondaryFlag * otherPatternFlag * piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(269, secondaryFlag * otherPatternFlag * piecewiseDistanceFrom5To10Km);
alternative.AddUtilityTerm(270, secondaryFlag * otherPatternFlag * piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(271, secondaryFlag * otherPatternFlag * piecewiseDistanceFrom20KmToInfinity);
//alternative.AddUtilityTerm(268, (!_tour.IsHomeBasedTour).ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(269, household.Has0To15KIncome.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(270, household.HasMissingIncome.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(271, person.IsRetiredAdult.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(272, person.IsUniversityStudent.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(273, person.IsChildAge5Through15.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(274, person.IsChildUnder5.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(272, (!tour.IsHomeBasedTour).ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(273, (household.Income >= 300000 && household.Income < 600000).ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(274, (household.Income >= 600000 && household.Income < 900000).ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(275, (household.Income >= 900000).ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(276, person.IsChildUnder5.ToFlag() * distanceFromOriginLog); // commented out by GV, July 7th
//alternative.AddUtilityTerm(277, person.IsChildAge5Through15.ToFlag() * distanceFromOriginLog); // commented out by GV, July 7th
//alternative.AddUtilityTerm(278, person.IsChildUnder16.ToFlag() * distanceFromOriginLog); // commented out by GV, July 7th
alternative.AddUtilityTerm(279, person.IsUniversityStudent.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(280, person.IsAdultMale.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(281, person.IsAdultFemale.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(282, person.IsRetiredAdult.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(283, (tour.IsHomeBasedTour).ToFlag() * timePressure); //commented out by GV: 7th July 2013
alternative.AddUtilityTerm(284, (tour.IsHomeBasedTour).ToFlag() * distanceFromSchoolLog);
//alternative.AddUtilityTerm(14, distanceFromWorkLog);
// GV commented out this - on TO DO list
//alternative.AddUtilityTerm(277, (carCompetitionFlag + noCarCompetitionFlag) * destinationParcel.ParkingHourlyEmploymentCommercialMixInParcel());
//alternative.AddUtilityTerm(278, noCarCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixInParcel());
//alternative.AddUtilityTerm(279, carCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixBuffer1());
//alternative.AddUtilityTerm(280, noCarCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixBuffer1());
//alternative.AddUtilityTerm(281, noCarsFlag * c34Ratio);
//alternative.AddUtilityTerm(282, noCarCompetitionFlag * c34Ratio);
//alternative.AddUtilityTerm(283, (carCompetitionFlag + noCarCompetitionFlag) * logOfOnePlusParkingOffStreetDailySpacesBuffer1);
//alternative.AddUtilityTerm(285, jointTourFlag * piecewiseDistanceFrom0To1Km);
//alternative.AddUtilityTerm(286, jointTourFlag * piecewiseDistanceFrom1To2Km);
alternative.AddUtilityTerm(286, jointTourFlag * piecewiseDistanceFrom0To2Km);
alternative.AddUtilityTerm(287, jointTourFlag * piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(288, jointTourFlag * piecewiseDistanceFrom5To10Km);
alternative.AddUtilityTerm(289, jointTourFlag * piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(290, jointTourFlag * piecewiseDistanceFrom20KmToInfinity);
// 3. finished
//4. new from GV: purpose utilities
// COMPAS puposes are: Work, Education, Escort, Shopping, Leisure, Personal business, business
// You need NO "Work" and "Education", their destinations are known in the synthetic population
if (tour.DestinationPurpose == Global.Settings.Purposes.Business)
{
//alternative.AddUtilityTerm(10, piecewiseDistanceFrom0To1Km);
//alternative.AddUtilityTerm(11, piecewiseDistanceFrom1To2Km);
//alternative.AddUtilityTerm(12, piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(12, piecewiseDistanceFrom0To5Km);
alternative.AddUtilityTerm(13, piecewiseDistanceFrom5To10Km);
alternative.AddUtilityTerm(14, piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(15, piecewiseDistanceFrom20KmToInfinity);
alternative.AddUtilityTerm(16, aggregateLogsumWorkBased);
// Neighborhood
//GV: commented out just temp.
//alternative.AddUtilityTerm(20, logOfOnePlusEducationK8Buffer2);
//alternative.AddUtilityTerm(21, logOfOnePlusEducationUniStuBuffer2);
//alternative.AddUtilityTerm(22, logOfOnePlusEmploymentEducationBuffer2);
alternative.AddUtilityTerm(23, logOfOnePlusEmploymentGovernmentBuffer2);
//alternative.AddUtilityTerm(24, logOfOnePlusEmploymentIndustrialBuffer2);
alternative.AddUtilityTerm(25, logOfOnePlusEmploymentOfficeBuffer2);
alternative.AddUtilityTerm(26, logOfOnePlusEmploymentRetailBuffer2);
alternative.AddUtilityTerm(27, logOfOnePlusEmploymentServiceBuffer2);
alternative.AddUtilityTerm(28, logOfOnePlusEmploymentAgrConstrBuffer2);
//alternative.AddUtilityTerm(29, logOfOnePlusEmploymentJobsBuffer2);
// Size terms
alternative.AddUtilityTerm(30, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(31, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(32, destinationParcel.EmploymentIndustrial);
alternative.AddUtilityTerm(33, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(34, destinationParcel.EmploymentRetail);
// GV: 35 is fixed to zero
alternative.AddUtilityTerm(35, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(36, destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(37, destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(38, destinationParcel.Households);
}
else if (tour.DestinationPurpose == Global.Settings.Purposes.Escort)
{
//alternative.AddUtilityTerm(50, piecewiseDistanceFrom0To1Km);
alternative.AddUtilityTerm(51, piecewiseDistanceFrom0To2Km);
//alternative.AddUtilityTerm(52, piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(52, piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(53, piecewiseDistanceFrom5To10Km);
//alternative.AddUtilityTerm(54, piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(55, piecewiseDistanceFrom10KmToInfinity);
//alternative.AddUtilityTerm(55, piecewiseDistanceFrom20KmToInfinity);
alternative.AddUtilityTerm(56, aggregateLogsumHomeBased);
// Neighborhood
//GV: commented out just temp.
alternative.AddUtilityTerm(60, householdHasNoChildren.ToFlag() * logOfOnePlusEmploymentJobsBuffer2);
//alternative.AddUtilityTerm(61, householdHasNoChildren.ToFlag() * logOfOnePlusHouseholdsBuffer2);
//alternative.AddUtilityTerm(62, householdHasChildren.ToFlag() * logOfOnePlusHouseholdsBuffer2);
//alternative.AddUtilityTerm(64, logOfOnePlusEmploymentJobsBuffer2);
// Size terms
// GV: no observations
alternative.AddUtilityTerm(70, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(71, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(72, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentIndustrial);
alternative.AddUtilityTerm(73, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(74, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentRetail);
// GV: 75 is fixed to zero
alternative.AddUtilityTerm(75, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(76, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(77, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(78, (!householdHasChildren).ToFlag() * destinationParcel.Households);
alternative.AddUtilityTerm(80, (householdHasChildren).ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(81, (householdHasChildren).ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(82, (householdHasChildren).ToFlag() * destinationParcel.EmploymentIndustrial);
alternative.AddUtilityTerm(83, (householdHasChildren).ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(84, (householdHasChildren).ToFlag() * destinationParcel.EmploymentRetail);
// GV 85 is fixed to zero at the moment
alternative.AddUtilityTerm(85, (householdHasChildren).ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(86, (householdHasChildren).ToFlag() * destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(87, (householdHasChildren).ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(88, (householdHasChildren).ToFlag() * destinationParcel.Households);
}
else if (tour.DestinationPurpose == Global.Settings.Purposes.PersonalBusiness)
{
alternative.AddUtilityTerm(90, piecewiseDistanceFrom0To2Km);
//alternative.AddUtilityTerm(91, piecewiseDistanceFrom1To2Km);
alternative.AddUtilityTerm(92, piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(93, piecewiseDistanceFrom5To10Km);
alternative.AddUtilityTerm(94, piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(95, piecewiseDistanceFrom20KmToInfinity);
alternative.AddUtilityTerm(96, aggregateLogsumHomeBased);
// Neighborhood
//GV: commented out just temp.
//alternative.AddUtilityTerm(100, logOfOnePlusEmploymentEducationBuffer2);
//alternative.AddUtilityTerm(101, logOfOnePlusEmploymentGovernmentBuffer2);
//alternative.AddUtilityTerm(102, logOfOnePlusEmploymentIndustrialBuffer2);
alternative.AddUtilityTerm(103, logOfOnePlusEmploymentOfficeBuffer2);
alternative.AddUtilityTerm(104, logOfOnePlusEmploymentRetailBuffer2);
alternative.AddUtilityTerm(105, logOfOnePlusEmploymentServiceBuffer2);
//alternative.AddUtilityTerm(106, logOfOnePlusEmploymentAgrConstrBuffer2);
//alternative.AddUtilityTerm(107, logOfOnePlusEmploymentJobsBuffer2);
// Size terms
alternative.AddUtilityTerm(110, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(111, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(112, destinationParcel.EmploymentIndustrial);
alternative.AddUtilityTerm(113, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(114, destinationParcel.EmploymentRetail);
// GV 115 is fixed to zero
alternative.AddUtilityTerm(115, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(116, destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(117, destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(118, destinationParcel.Households);
}
else if (tour.DestinationPurpose == Global.Settings.Purposes.Shopping)
{
//alternative.AddUtilityTerm(120, piecewiseDistanceFrom0To1Km);
alternative.AddUtilityTerm(121, piecewiseDistanceFrom0To2Km);
alternative.AddUtilityTerm(122, piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(123, piecewiseDistanceFrom5To10Km);
alternative.AddUtilityTerm(124, piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(125, piecewiseDistanceFrom20KmToInfinity);
alternative.AddUtilityTerm(126, aggregateLogsumHomeBased);
// Neighborhood
//GV: commented out just temp.
//alternative.AddUtilityTerm(130, logOfOnePlusEmploymentEducationBuffer2);
alternative.AddUtilityTerm(131, logOfOnePlusEmploymentRetailBuffer2);
//alternative.AddUtilityTerm(132, logOfOnePlusEmploymentJobsBuffer2);
// Size terms
alternative.AddUtilityTerm(140, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(141, destinationParcel.EmploymentRetail);
// GV 142 is fixed to zero
alternative.AddUtilityTerm(142, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(143, destinationParcel.EmploymentTotal);
}
else if (tour.DestinationPurpose == Global.Settings.Purposes.Social)
{
//alternative.AddUtilityTerm(170, piecewiseDistanceFrom0To1Km);
//alternative.AddUtilityTerm(171, piecewiseDistanceFrom1To2Km);
alternative.AddUtilityTerm(170, piecewiseDistanceFrom0To2Km);
alternative.AddUtilityTerm(172, piecewiseDistanceFrom2To5Km);
alternative.AddUtilityTerm(173, piecewiseDistanceFrom5To10Km);
alternative.AddUtilityTerm(174, piecewiseDistanceFrom10To20Km);
alternative.AddUtilityTerm(175, piecewiseDistanceFrom20KmToInfinity);
alternative.AddUtilityTerm(176, aggregateLogsumHomeBased);
// Neighborhood
//GV: commented out just temp.
//alternative.AddUtilityTerm(180, logOfOnePlusEmploymentOfficeBuffer2);
alternative.AddUtilityTerm(181, logOfOnePlusEmploymentRetailBuffer2);
alternative.AddUtilityTerm(182, logOfOnePlusEmploymentServiceBuffer2);
//alternative.AddUtilityTerm(183, logOfOnePlusEmploymentJobsBuffer2);
// Size terms
alternative.AddUtilityTerm(190, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(191, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(192, destinationParcel.EmploymentIndustrial);
alternative.AddUtilityTerm(193, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(194, destinationParcel.EmploymentRetail);
// GV 195 is fixed to zero
alternative.AddUtilityTerm(195, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(196, destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(197, destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(198, destinationParcel.Households);
}
}
}
