public HM4Model(
ICalculationsAbstractFactory calculationsAbstractFactory,
IConstraintElementsAbstractFactory constraintElementsAbstractFactory,
IConstraintsAbstractFactory constraintsAbstractFactory,
ICrossJoinElementsAbstractFactory crossJoinElementsAbstractFactory,
ICrossJoinsAbstractFactory crossJoinsAbstractFactory,
IDependenciesAbstractFactory dependenciesAbstractFactory,
IIndexElementsAbstractFactory indexElementsAbstractFactory,
IIndicesAbstractFactory indicesAbstractFactory,
IObjectiveFunctionsAbstractFactory objectiveFunctionsAbstractFactory,
IParameterElementsAbstractFactory parameterElementsAbstractFactory,
IParametersAbstractFactory parametersAbstractFactory,
IVariablesAbstractFactory variablesAbstractFactory,
IHM4InputContext HM4InputContext)
{
this.Context = HM4InputContext;
this.Model = dependenciesAbstractFactory.CreateModelFactory().Create();
// Indices
// d
this.d = indicesAbstractFactory.CreatedFactory().Create(
this.Context.Weekdays
.Select(x => indexElementsAbstractFactory.CreatedIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// j
this.j = indicesAbstractFactory.CreatejFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => x.Item1)
.Select(x => indexElementsAbstractFactory.CreatejIndexElementFactory().Create(x))
.ToImmutableList());
// l
this.l = indicesAbstractFactory.CreatelFactory().Create(
this.Context.LengthOfStayDays
.Select(x => indexElementsAbstractFactory.CreatelIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// r
this.r = indicesAbstractFactory.CreaterFactory().Create(
this.Context.OperatingRooms
.Entry
.Where(x => x.Resource is Location)
.Select(x => indexElementsAbstractFactory.CreaterIndexElementFactory().Create((Location)x.Resource))
.ToImmutableList());
// s
this.s = indicesAbstractFactory.CreatesFactory().Create(
this.Context.Surgeons
.Entry
.Where(x => x.Resource is Organization)
.Select(x => indexElementsAbstractFactory.CreatesIndexElementFactory().Create((Organization)x.Resource))
.ToImmutableList());
// t
this.t = indicesAbstractFactory.CreatetFactory().Create(
this.Context.PlanningHorizon
.Select(x => indexElementsAbstractFactory.CreatetIndexElementFactory().Create(
x.Key.Value.Value,
x.Value))
.ToImmutableList());
// Λ
this.Λ = indicesAbstractFactory.CreateΛFactory().Create(
this.Context.Scenarios
.Select(x => indexElementsAbstractFactory.CreateΛIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// Cross joins
// dt
this.dt = crossJoinsAbstractFactory.CreatedtFactory().Create(
this.d.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatedtCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// rt
this.rt = crossJoinsAbstractFactory.CreatertFactory().Create(
this.r.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatertCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// slΛ
this.slΛ = crossJoinsAbstractFactory.CreateslΛFactory().Create(
this.s.Value
.SelectMany(b => this.l.Value, (a, b) => crossJoinElementsAbstractFactory.CreateslCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreateslΛCrossJoinElementFactory().Create(a.sIndexElement, a.lIndexElement, b))
.ToImmutableList());
// sr
this.sr = crossJoinsAbstractFactory.CreatesrFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// srd
this.srd = crossJoinsAbstractFactory.CreatesrdFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.d.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrdCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// srt
this.srt = crossJoinsAbstractFactory.CreatesrtFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrtCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// st
this.st = crossJoinsAbstractFactory.CreatestFactory().Create(
this.s.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatestCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// sΛ
this.sΛ = crossJoinsAbstractFactory.CreatesΛFactory().Create(
this.s.Value
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesΛCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// tΛ
this.tΛ = crossJoinsAbstractFactory.CreatetΛFactory().Create(
this.t.Value
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreatetΛCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// Parameters
// H
this.H = parametersAbstractFactory.CreateHFactory().Create(
this.Context.TimeBlockLength);
// h(s, Λ)
this.h = parametersAbstractFactory.CreatehFactory().Create(
this.Context.SurgeonScenarioWeightedAverageSurgicalDurations
.Select(x => parameterElementsAbstractFactory.CreatehParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
x.Item3))
.ToImmutableList());
// L(s)
this.L = parametersAbstractFactory.CreateLFactory().Create(
this.Context.SurgeonLengthOfStayMaximums
.Select(x => parameterElementsAbstractFactory.CreateLParameterElementFactory().Create(
this.s.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// n(s, Λ)
this.n = parametersAbstractFactory.CreatenFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatients
.Select(x => parameterElementsAbstractFactory.CreatenParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(PositiveInt)x.Item3))
.ToImmutableList());
// p(s, l, Λ)
this.p = parametersAbstractFactory.CreatepFactory().Create(
this.Context.SurgeonDayScenarioLengthOfStayProbabilities
.Select(x => parameterElementsAbstractFactory.CreatepParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.l.GetElementAt((PositiveInt)x.Item2),
this.Λ.GetElementAt((PositiveInt)x.Item3),
(FhirDecimal)x.Item4))
.ToImmutableList());
// W
this.W = parametersAbstractFactory.CreateWFactory().Create(
(PositiveInt)this.Context.NumberDaysPerWeek);
// y(s, r)
this.y = parametersAbstractFactory.CreateyFactory().Create(
this.Context.SurgeonOperatingRoomAssignments
.Select(x => parameterElementsAbstractFactory.CreateyParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.r.GetElementAt(x.Item2),
(FhirBoolean)x.Item3))
.ToImmutableList());
// β(s, r, d)
this.β = parametersAbstractFactory.CreateβFactory().Create(
this.Context.SurgeonOperatingRoomDayNumberAssignedTimeBlocks
.Select(x => parameterElementsAbstractFactory.CreateβParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.r.GetElementAt(x.Item2),
this.d.GetElementAt((PositiveInt)x.Item3),
(PositiveInt)x.Item4))
.ToImmutableList());
// γ(r, t)
this.γ = parametersAbstractFactory.CreateγFactory().Create(
this.Context.OperatingRoomDayAssignedAvailabilities
.Select(x => parameterElementsAbstractFactory.CreateγParameterElementFactory().Create(
this.r.GetElementAt(x.Item1),
this.t.GetElementAt(x.Item2),
(FhirBoolean)x.Item3))
.ToImmutableList());
// Δ(j)
this.Δ = parametersAbstractFactory.CreateΔFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => parameterElementsAbstractFactory.CreateΔParameterElementFactory().Create(
this.j.GetElementAt(x.Item1),
x.Item2.Select(i => this.s.GetElementAt(i)).ToImmutableList()))
.ToImmutableList());
// δ1(s, r, d)
this.δ1 = calculationsAbstractFactory.Createδ1CalculationFactory().Create().Calculate(
parameterElementsAbstractFactory.Createδ1ParameterElementFactory(),
parametersAbstractFactory.Createδ1Factory(),
calculationsAbstractFactory.Createδ1ParameterElementCalculationFactory().Create(),
this.srd,
this.t,
this.W,
this.β,
this.γ);
// δ2(s, r, d)
this.δ2 = calculationsAbstractFactory.Createδ2CalculationFactory().Create().Calculate(
parameterElementsAbstractFactory.Createδ2ParameterElementFactory(),
parametersAbstractFactory.Createδ2Factory(),
calculationsAbstractFactory.Createδ2ParameterElementCalculationFactory().Create(),
this.srd,
this.t,
this.W,
this.β,
this.γ);
// δ3(s, r, d)
this.δ3 = calculationsAbstractFactory.Createδ3CalculationFactory().Create().Calculate(
parameterElementsAbstractFactory.Createδ3ParameterElementFactory(),
parametersAbstractFactory.Createδ3Factory(),
calculationsAbstractFactory.Createδ3ParameterElementCalculationFactory().Create(),
this.srd,
this.t,
this.W,
this.β,
this.γ);
// δ4(s, r, d, t)
this.δ4 = calculationsAbstractFactory.Createδ4CalculationFactory().Create().Calculate(
parameterElementsAbstractFactory.Createδ4ParameterElementFactory(),
parametersAbstractFactory.Createδ4Factory(),
calculationsAbstractFactory.Createδ4ParameterElementCalculationFactory().Create(),
this.srd,
this.t,
this.W,
this.δ3,
this.γ);
// μ(s, Λ)
this.μ = parametersAbstractFactory.CreateμFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatientMeans
.Select(x => parameterElementsAbstractFactory.CreateμParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(FhirDecimal)x.Item3))
.ToImmutableList());
// Ρ(Λ)
this.Ρ = parametersAbstractFactory.CreateΡFactory().Create(
this.Context.ScenarioProbabilities
.Select(x => parameterElementsAbstractFactory.CreateΡParameterElementFactory().Create(
this.Λ.GetElementAt((PositiveInt)x.Key),
(FhirDecimal)x.Value))
.ToImmutableList());
// σ(s, Λ)
this.σ = parametersAbstractFactory.CreateσFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatientStandardDeviations
.Select(x => parameterElementsAbstractFactory.CreateσParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(FhirDecimal)x.Item3))
.ToImmutableList());
// ΦHat(s, l, Λ)
this.ΦHat = parametersAbstractFactory.CreateΦHatFactory().Create(
this.Context.SurgeonDayScenarioCumulativeNumberPatients
.Select(x => parameterElementsAbstractFactory.CreateΦHatParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.l.GetElementAt((PositiveInt)x.Item2),
this.Λ.GetElementAt((PositiveInt)x.Item3),
(FhirDecimal)x.Item4))
.ToImmutableList());
// Ω
this.Ω = parametersAbstractFactory.CreateΩFactory().Create(
(PositiveInt)this.Context.MaximumNumberRecoveryWardBeds);
// Variables
// IHat(t, Λ)
this.IHat = variablesAbstractFactory.CreateIHatFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.t.Value,
indexSet2: this.Λ.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => double.MaxValue,
variableTypeGenerator: (a, b) => VariableType.Continuous));
// x(s, r, t)
this.x = variablesAbstractFactory.CreatexFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().CreateBinary(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.r.Value,
indexSet3: this.t.Value));
// z(s, t)
this.z = variablesAbstractFactory.CreatezFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().CreateBinary(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.t.Value));
// Constraints
// Constraints 1
this.Model.AddConstraints(
this.srd.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints1ConstraintElementFactory().Create(
x.dIndexElement,
x.rIndexElement,
x.sIndexElement,
this.t,
this.W,
this.β,
this.x)
.Value));
// Constraints 2
this.Model.AddConstraints(
this.st.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints2ConstraintElementFactory().Create(
x.sIndexElement,
x.tIndexElement,
this.r,
this.x,
this.z)
.Value));
// Constraints 3
this.Model.AddConstraints(
this.rt.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints3ConstraintElementFactory().Create(
x.rIndexElement,
x.tIndexElement,
this.s,
this.γ,
this.x)
.Value));
// Constraints 4
this.Model.AddConstraints(
this.srt.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints4ConstraintElementFactory().Create(
x.rIndexElement,
x.sIndexElement,
x.tIndexElement,
this.y,
this.x)
.Value));
// Constraints 5L
this.Model.AddConstraints(
this.srd.Value
.SelectMany(b =>
this.t.GetNthElementsAt(
b.dIndexElement.Value.Value.Value,
this.t.GetT() - this.W.Value.Value.Value + b.dIndexElement.Value.Value.Value,
this.W.Value.Value.Value),
(a, b) => Tuple.Create(a.sIndexElement, a.rIndexElement, a.dIndexElement, b))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints5LConstraintElementFactory().Create(
x.Item3,
x.Item2,
x.Item1,
x.Item4,
this.t,
this.W,
this.δ1,
this.δ4,
this.x)
.Value));
// Constraints 5U
this.Model.AddConstraints(
this.srd.Value
.SelectMany(b =>
this.t.GetNthElementsAt(
b.dIndexElement.Value.Value.Value,
this.t.GetT() - this.W.Value.Value.Value + b.dIndexElement.Value.Value.Value,
this.W.Value.Value.Value),
(a, b) => Tuple.Create(a.sIndexElement, a.rIndexElement, a.dIndexElement, b))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints5UConstraintElementFactory().Create(
x.Item3,
x.Item2,
x.Item1,
x.Item4,
this.t,
this.W,
this.δ2,
this.δ4,
this.x)
.Value));
// Constraints 6
this.Model.AddConstraints(
this.tΛ.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints6ConstraintElementFactory().Create(
x.tIndexElement,
x.ΛIndexElement,
this.l,
this.t,
this.st,
this.ΦHat,
this.IHat,
this.z)
.Value));
// Objective function
this.Model.AddObjective(
objectiveFunctionsAbstractFactory.CreateObjectiveFunctionFactory().Create(
dependenciesAbstractFactory.CreateObjectiveFactory(),
this.t,
this.Λ,
this.Ρ,
this.IHat)
.Value);
}
public HM5Model(
ICalculationsAbstractFactory calculationsAbstractFactory,
IConstraintElementsAbstractFactory constraintElementsAbstractFactory,
IConstraintsAbstractFactory constraintsAbstractFactory,
ICrossJoinElementsAbstractFactory crossJoinElementsAbstractFactory,
ICrossJoinsAbstractFactory crossJoinsAbstractFactory,
IDependenciesAbstractFactory dependenciesAbstractFactory,
IIndexElementsAbstractFactory indexElementsAbstractFactory,
IIndicesAbstractFactory indicesAbstractFactory,
IObjectiveFunctionsAbstractFactory objectiveFunctionsAbstractFactory,
IParameterElementsAbstractFactory parameterElementsAbstractFactory,
IParametersAbstractFactory parametersAbstractFactory,
IVariablesAbstractFactory variablesAbstractFactory,
IHM5InputContext HM5InputContext)
{
this.Context = HM5InputContext;
this.Model = dependenciesAbstractFactory.CreateModelFactory().Create();
// Indices
// d1
this.d1 = indicesAbstractFactory.Created1Factory().Create(
this.Context.Weekdays
.Select(x => indexElementsAbstractFactory.Created1IndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// d2
this.d2 = indicesAbstractFactory.Created2Factory().Create(
this.Context.Weekdays
.Select(x => indexElementsAbstractFactory.Created2IndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// j
this.j = indicesAbstractFactory.CreatejFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => x.Item1)
.Select(x => indexElementsAbstractFactory.CreatejIndexElementFactory().Create(x))
.ToImmutableList());
// l
this.l = indicesAbstractFactory.CreatelFactory().Create(
this.Context.LengthOfStayDays
.Select(x => indexElementsAbstractFactory.CreatelIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// r
this.r = indicesAbstractFactory.CreaterFactory().Create(
this.Context.OperatingRooms
.Entry
.Where(x => x.Resource is Location)
.Select(x => indexElementsAbstractFactory.CreaterIndexElementFactory().Create((Location)x.Resource))
.ToImmutableList());
// s
this.s = indicesAbstractFactory.CreatesFactory().Create(
this.Context.Surgeons
.Entry
.Where(x => x.Resource is Organization)
.Select(x => indexElementsAbstractFactory.CreatesIndexElementFactory().Create((Organization)x.Resource))
.ToImmutableList());
// t
this.t = indicesAbstractFactory.CreatetFactory().Create(
this.Context.PlanningHorizon
.Select(x => indexElementsAbstractFactory.CreatetIndexElementFactory().Create(
x.Key.Value.Value,
x.Value))
.ToImmutableList());
// Λ
this.Λ = indicesAbstractFactory.CreateΛFactory().Create(
this.Context.Scenarios
.Select(x => indexElementsAbstractFactory.CreateΛIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// Cross joins
// rd1
this.rd1 = crossJoinsAbstractFactory.Createrd1Factory().Create(
this.r.Value
.SelectMany(b => this.d1.Value, (a, b) => crossJoinElementsAbstractFactory.Createrd1CrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// rd1d2
this.rd1d2 = crossJoinsAbstractFactory.Createrd1d2Factory().Create(
this.r.Value
.SelectMany(b => this.d1.Value, (a, b) => crossJoinElementsAbstractFactory.Createrd1CrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.d2.Value, (a, b) => crossJoinElementsAbstractFactory.Createrd1d2CrossJoinElementFactory().Create(a.rIndexElement, a.d1IndexElement, b))
.ToImmutableList());
// rd2
this.rd2 = crossJoinsAbstractFactory.Createrd2Factory().Create(
this.r.Value
.SelectMany(b => this.d2.Value, (a, b) => crossJoinElementsAbstractFactory.Createrd2CrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// rt
this.rt = crossJoinsAbstractFactory.CreatertFactory().Create(
this.r.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatertCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// slΛ
this.slΛ = crossJoinsAbstractFactory.CreateslΛFactory().Create(
this.s.Value
.SelectMany(b => this.l.Value, (a, b) => crossJoinElementsAbstractFactory.CreateslCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreateslΛCrossJoinElementFactory().Create(a.sIndexElement, a.lIndexElement, b))
.ToImmutableList());
// srd2
this.srd2 = crossJoinsAbstractFactory.Createsrd2Factory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.d2.Value, (a, b) => crossJoinElementsAbstractFactory.Createsrd2CrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// srd2t
this.srd2t = crossJoinsAbstractFactory.Createsrd2tFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.d2.Value, (a, b) => crossJoinElementsAbstractFactory.Createsrd2CrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.Createsrd2tCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, a.d2IndexElement, b))
.ToImmutableList());
// srt
this.srt = crossJoinsAbstractFactory.CreatesrtFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrtCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// st
this.st = crossJoinsAbstractFactory.CreatestFactory().Create(
this.s.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatestCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// sΛ
this.sΛ = crossJoinsAbstractFactory.CreatesΛFactory().Create(
this.s.Value
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesΛCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// tΛ
this.tΛ = crossJoinsAbstractFactory.CreatetΛFactory().Create(
this.t.Value
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreatetΛCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// Parameters
// H
this.H = parametersAbstractFactory.CreateHFactory().Create(
this.Context.TimeBlockLength);
// h(s, Λ)
this.h = parametersAbstractFactory.CreatehFactory().Create(
this.Context.SurgeonScenarioWeightedAverageSurgicalDurations
.Select(x => parameterElementsAbstractFactory.CreatehParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
x.Item3))
.ToImmutableList());
// L(s)
this.L = parametersAbstractFactory.CreateLFactory().Create(
this.Context.SurgeonLengthOfStayMaximums
.Select(x => parameterElementsAbstractFactory.CreateLParameterElementFactory().Create(
this.s.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// n(s, Λ)
this.n = parametersAbstractFactory.CreatenFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatients
.Select(x => parameterElementsAbstractFactory.CreatenParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(PositiveInt)x.Item3))
.ToImmutableList());
// p(s, l, Λ)
this.p = parametersAbstractFactory.CreatepFactory().Create(
this.Context.SurgeonDayScenarioLengthOfStayProbabilities
.Select(x => parameterElementsAbstractFactory.CreatepParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.l.GetElementAt((PositiveInt)x.Item2),
this.Λ.GetElementAt((PositiveInt)x.Item3),
(FhirDecimal)x.Item4))
.ToImmutableList());
// Δ(j)
// Must be populated before S1
this.Δ = parametersAbstractFactory.CreateΔFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => parameterElementsAbstractFactory.CreateΔParameterElementFactory().Create(
this.j.GetElementAt(x.Item1),
x.Item2.Select(i => this.s.GetElementAt(i)).ToImmutableList()))
.ToImmutableList());
// x(s, r, t)
// Must be populated before S1
this.x = parametersAbstractFactory.CreatexFactory().Create(
this.Context.SurgeonOperatingRoomDayAssignments
.Select(x => parameterElementsAbstractFactory.CreatexParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.r.GetElementAt(x.Item2),
this.t.GetElementAt(x.Item3),
(FhirBoolean)x.Item4))
.ToImmutableList());
// S1
this.S1 = calculationsAbstractFactory.CreateS1CalculationFactory().Create().Calculate(
parameterElementsAbstractFactory.CreateS1ParameterElementFactory(),
parametersAbstractFactory.CreateS1Factory(),
this.r,
this.x,
this.Δ);
// W
// Must be populated before S2
this.W = parametersAbstractFactory.CreateWFactory().Create(
(PositiveInt)this.Context.NumberDaysPerWeek);
// S2
this.S2 = calculationsAbstractFactory.CreateS2CalculationFactory().Create().Calculate(
parameterElementsAbstractFactory.CreateS2ParameterElementFactory(),
parametersAbstractFactory.CreateS2Factory(),
this.d1,
this.d2,
this.r,
this.t,
this.rd1d2,
this.W,
this.x,
this.Δ);
// y(s, r)
this.y = parametersAbstractFactory.CreateyFactory().Create(
this.Context.SurgeonOperatingRoomAssignments
.Select(x => parameterElementsAbstractFactory.CreateyParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.r.GetElementAt(x.Item2),
(FhirBoolean)x.Item3))
.ToImmutableList());
// γ(r, t)
this.γ = parametersAbstractFactory.CreateγFactory().Create(
this.Context.OperatingRoomDayAssignedAvailabilities
.Select(x => parameterElementsAbstractFactory.CreateγParameterElementFactory().Create(
this.r.GetElementAt(x.Item1),
this.t.GetElementAt(x.Item2),
(FhirBoolean)x.Item3))
.ToImmutableList());
// μ(s, Λ)
this.μ = parametersAbstractFactory.CreateμFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatientMeans
.Select(x => parameterElementsAbstractFactory.CreateμParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(FhirDecimal)x.Item3))
.ToImmutableList());
// Ρ(Λ)
this.Ρ = parametersAbstractFactory.CreateΡFactory().Create(
this.Context.ScenarioProbabilities
.Select(x => parameterElementsAbstractFactory.CreateΡParameterElementFactory().Create(
this.Λ.GetElementAt((PositiveInt)x.Key),
(FhirDecimal)x.Value))
.ToImmutableList());
// σ(s, Λ)
this.σ = parametersAbstractFactory.CreateσFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatientStandardDeviations
.Select(x => parameterElementsAbstractFactory.CreateσParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(FhirDecimal)x.Item3))
.ToImmutableList());
// ΦHat(s, l, Λ)
this.ΦHat = parametersAbstractFactory.CreateΦHatFactory().Create(
this.Context.SurgeonDayScenarioCumulativeNumberPatients
.Select(x => parameterElementsAbstractFactory.CreateΦHatParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.l.GetElementAt((PositiveInt)x.Item2),
this.Λ.GetElementAt((PositiveInt)x.Item3),
(FhirDecimal)x.Item4))
.ToImmutableList());
// Ω
this.Ω = parametersAbstractFactory.CreateΩFactory().Create(
(PositiveInt)this.Context.MaximumNumberRecoveryWardBeds);
// Variables
// IHat(t, Λ)
this.IHat = variablesAbstractFactory.CreateIHatFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.t.Value,
indexSet2: this.Λ.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => double.MaxValue,
variableTypeGenerator: (a, b) => VariableType.Continuous));
// xHat(s, r, t)
this.xHat = variablesAbstractFactory.CreatexHatFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.r.Value,
indexSet3: this.t.Value,
lowerBoundGenerator: (a, b, c) => 0,
upperBoundGenerator: (a, b, c) => 1,
variableTypeGenerator: (a, b, c) => VariableType.Binary));
// z(s, t)
this.z = variablesAbstractFactory.CreatezFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.t.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => 1,
variableTypeGenerator: (a, b) => VariableType.Binary));
// α(r, d1, d2)
this.α = variablesAbstractFactory.CreateαFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.r.Value,
indexSet2: this.d1.Value,
indexSet3: this.d2.Value,
lowerBoundGenerator: (a, b, c) => 0,
upperBoundGenerator: (a, b, c) => 1,
variableTypeGenerator: (a, b, c) => VariableType.Binary));
// Constraints
// Constraints 1
this.Model.AddConstraints(
this.st.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints1ConstraintElementFactory().Create(
x.sIndexElement,
x.tIndexElement,
this.r,
this.xHat,
this.z)
.Value));
// Constraints 2
this.Model.AddConstraints(
this.rt.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints2ConstraintElementFactory().Create(
x.rIndexElement,
x.tIndexElement,
this.s,
this.γ,
this.xHat)
.Value));
// Constraints 3
this.Model.AddConstraints(
this.srt.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints3ConstraintElementFactory().Create(
x.rIndexElement,
x.sIndexElement,
x.tIndexElement,
this.y,
this.xHat)
.Value));
// Constraints 4
this.Model.AddConstraints(
this.tΛ.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints4ConstraintElementFactory().Create(
x.tIndexElement,
x.ΛIndexElement,
this.l,
this.t,
this.st,
this.ΦHat,
this.IHat,
this.z)
.Value));
// Constraints 5
this.Model.AddConstraints(
this.rd1.Value
.Where(x => this.S1.IsThereElementAt(x.rIndexElement))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints5ConstraintElementFactory().Create(
x.d1IndexElement,
x.rIndexElement,
this.d2,
this.α)
.Value));
// Constraints 6
this.Model.AddConstraints(
this.rd2.Value
.Where(x => this.S1.IsThereElementAt(x.rIndexElement))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints6ConstraintElementFactory().Create(
x.d2IndexElement,
x.rIndexElement,
this.d1,
this.α)
.Value));
// Constraints 7
this.Model.AddConstraints(
this.rd1d2.Value
.Where(x => this.S1.IsThereElementAt(x.rIndexElement))
.Where(x => this.S2.IsThereElementAt(
x.rIndexElement,
x.d1IndexElement,
x.d2IndexElement))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints7ConstraintElementFactory().Create(
x.rIndexElement,
x.d1IndexElement,
x.d2IndexElement,
this.α)
.Value));
// Constraints 8
this.Model.AddConstraints(
this.srd2.Value
.Where(x => this.S1.IsThereElementAt(x.rIndexElement))
.SelectMany(b =>
this.t.GetNthElementsAt(
0,
this.W.Value.Value.Value * ((this.t.GetT() / this.W.Value.Value.Value) - 1),
this.W.Value.Value.Value),
(a, b) => Tuple.Create(a.sIndexElement, a.rIndexElement, a.d2IndexElement, b))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints8ConstraintElementFactory().Create(
x.Item3,
x.Item2,
x.Item1,
x.Item4,
this.d1,
this.t,
this.x,
this.xHat,
this.α)
.Value));
// Constraints 9
this.Model.AddConstraints(
this.srt.Value
.Where(x => !this.S1.IsThereElementAt(x.rIndexElement))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints9ConstraintElementFactory().Create(
x.rIndexElement,
x.sIndexElement,
x.tIndexElement,
this.x,
this.xHat)
.Value));
// Objective function
this.Model.AddObjective(
objectiveFunctionsAbstractFactory.CreateObjectiveFunctionFactory().Create(
dependenciesAbstractFactory.CreateObjectiveFactory(),
this.t,
this.Λ,
this.Ρ,
this.IHat)
.Value);
}
public HM3AModel(
IConstraintElementsAbstractFactory constraintElementsAbstractFactory,
IConstraintsAbstractFactory constraintsAbstractFactory,
ICrossJoinElementsAbstractFactory crossJoinElementsAbstractFactory,
ICrossJoinsAbstractFactory crossJoinsAbstractFactory,
IDependenciesAbstractFactory dependenciesAbstractFactory,
IIndexElementsAbstractFactory indexElementsAbstractFactory,
IIndicesAbstractFactory indicesAbstractFactory,
IParameterElementsAbstractFactory parameterElementsAbstractFactory,
IParametersAbstractFactory parametersAbstractFactory,
IVariablesAbstractFactory variablesAbstractFactory,
IHM3AInputContext HM3AInputContext)
{
this.Context = HM3AInputContext;
this.Model = dependenciesAbstractFactory.CreateModelFactory().Create();
// Indices
// j
this.j = indicesAbstractFactory.CreatejFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => x.Item1)
.Select(x => indexElementsAbstractFactory.CreatejIndexElementFactory().Create(x))
.ToImmutableList());
// m
this.m = indicesAbstractFactory.CreatemFactory().Create(
this.Context.Machines
.Entry
.Where(x => x.Resource is Device)
.Select(x => indexElementsAbstractFactory.CreatemIndexElementFactory().Create((Device)x.Resource))
.ToImmutableList());
// r
this.r = indicesAbstractFactory.CreaterFactory().Create(
this.Context.OperatingRooms
.Entry
.Where(x => x.Resource is Location)
.Select(x => indexElementsAbstractFactory.CreaterIndexElementFactory().Create((Location)x.Resource))
.ToImmutableList());
// s
this.s = indicesAbstractFactory.CreatesFactory().Create(
this.Context.Surgeons
.Entry
.Where(x => x.Resource is Organization)
.Select(x => indexElementsAbstractFactory.CreatesIndexElementFactory().Create((Organization)x.Resource))
.ToImmutableList());
// t
this.t = indicesAbstractFactory.CreatetFactory().Create(
this.Context.PlanningHorizon
.Select(x => indexElementsAbstractFactory.CreatetIndexElementFactory().Create(
x.Key.Value.Value,
x.Value))
.ToImmutableList());
// Cross joins
// jr
this.jr = crossJoinsAbstractFactory.CreatejrFactory().Create(
this.j.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatejrCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// mr
this.mr = crossJoinsAbstractFactory.CreatemrFactory().Create(
this.m.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatemrCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// rt
this.rt = crossJoinsAbstractFactory.CreatertFactory().Create(
this.r.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatertCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// sr
this.sr = crossJoinsAbstractFactory.CreatesrFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// srj
this.srj = crossJoinsAbstractFactory.CreatesrjFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.j.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrjCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// Parameters
// BBar(j)
this.BBar = parametersAbstractFactory.CreateBBarFactory().Create(
this.Context.SurgicalSpecialtyNumberAssignedTimeBlocks
.Select(x => parameterElementsAbstractFactory.CreateBBarParameterElementFactory().Create(
this.j.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// B(s)
this.Bs = parametersAbstractFactory.CreateBsFactory().Create(
this.Context.SurgeonNumberAssignedTimeBlocks
.Select(x => parameterElementsAbstractFactory.CreateBsParameterElementFactory().Create(
this.s.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// C(m)
this.C = parametersAbstractFactory.CreateCFactory().Create(
this.Context.MachineCosts
.Select(x => parameterElementsAbstractFactory.CreateCParameterElementFactory().Create(
this.m.GetElementAt(x.Key),
x.Value))
.ToImmutableList());
// Δ(j)
this.Δ = parametersAbstractFactory.CreateΔFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => parameterElementsAbstractFactory.CreateΔParameterElementFactory().Create(
this.j.GetElementAt(x.Item1),
x.Item2.Select(i => this.s.GetElementAt(i)).ToImmutableList()))
.ToImmutableList());
// ζ(s, m)
this.ζ = parametersAbstractFactory.CreateζFactory().Create(
this.Context.SurgeonMachineRequirements
.Select(x => parameterElementsAbstractFactory.CreateζParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.m.GetElementAt(x.Item2),
(FhirBoolean)x.Item3))
.ToImmutableList());
// ψ(t)
this.ψ = parametersAbstractFactory.CreateψFactory().Create(
this.Context.DayAvailabilities
.Select(x => parameterElementsAbstractFactory.CreateψParameterElementFactory().Create(
this.t.GetElementAt(x.Key),
(FhirBoolean)x.Value))
.ToImmutableList());
// Variables
// b(s, r)
this.b = variablesAbstractFactory.CreatebFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.r.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => int.MaxValue,
variableTypeGenerator: (a, b) => VariableType.Integer));
// v(m, r)
this.v = variablesAbstractFactory.CreatevFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.m.Value,
indexSet2: this.r.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => 1,
variableTypeGenerator: (a, b) => VariableType.Binary));
// w(j, r)
this.w = variablesAbstractFactory.CreatewFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.j.Value,
indexSet2: this.r.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => 1,
variableTypeGenerator: (a, b) => VariableType.Binary));
// y(s, r)
this.y = variablesAbstractFactory.CreateyFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.r.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => 1,
variableTypeGenerator: (a, b) => VariableType.Binary));
// γ(r, t)
this.γ = variablesAbstractFactory.CreateγFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.r.Value,
indexSet2: this.t.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => 1,
variableTypeGenerator: (a, b) => VariableType.Binary));
// Constraints
// Constraints 1
this.Model.AddConstraints(
this.mr.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints1ConstraintElementFactory().Create(
x.mIndexElement,
x.rIndexElement,
this.m,
this.s,
this.ζ,
this.v,
this.y)
.Value));
// Constraints 2
this.Model.AddConstraints(
this.s.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints2ConstraintElementFactory().Create(
x,
this.r,
this.Bs,
this.b)
.Value));
// Constraints 3
this.Model.AddConstraints(
this.r.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints3ConstraintElementFactory().Create(
x,
this.s,
this.t,
this.b,
this.γ)
.Value));
// Constraints 4
this.Model.AddConstraints(
this.rt.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints4ConstraintElementFactory().Create(
x.rIndexElement,
x.tIndexElement,
this.ψ,
this.γ)
.Value));
// Constraints 5L
this.Model.AddConstraints(
this.srj.Value
.Where(
x => this.Δ.IsSurgeonMemberOfSurgicalSpecialty(
x.jIndexElement,
x.sIndexElement))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints5LConstraintElementFactory().Create(
x.rIndexElement,
x.sIndexElement,
this.b,
this.y)
.Value));
// Constraints 5M
this.Model.AddConstraints(
this.srj.Value
.Where(
x => this.Δ.IsSurgeonMemberOfSurgicalSpecialty(
x.jIndexElement,
x.sIndexElement))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints5MConstraintElementFactory().Create(
x.rIndexElement,
x.sIndexElement,
this.Bs,
this.b,
this.y)
.Value));
// Constraints 5U
this.Model.AddConstraints(
this.srj.Value
.Where(
x => this.Δ.IsSurgeonMemberOfSurgicalSpecialty(
x.jIndexElement,
x.sIndexElement))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints5UConstraintElementFactory().Create(
x.jIndexElement,
x.rIndexElement,
x.sIndexElement,
this.BBar,
this.Bs,
this.w,
this.y)
.Value));
}
public HM1BModel(
IConstraintElementsAbstractFactory constraintElementsAbstractFactory,
IConstraintsAbstractFactory constraintsAbstractFactory,
ICrossJoinElementsAbstractFactory crossJoinElementsAbstractFactory,
ICrossJoinsAbstractFactory crossJoinsAbstractFactory,
IDependenciesAbstractFactory dependenciesAbstractFactory,
IIndexElementsAbstractFactory indexElementsAbstractFactory,
IIndicesAbstractFactory indicesAbstractFactory,
IParameterElementsAbstractFactory parameterElementsAbstractFactory,
IParametersAbstractFactory parametersAbstractFactory,
IVariablesAbstractFactory variablesAbstractFactory,
IHM1BInputContext HM1BInputContext)
{
this.Context = HM1BInputContext;
this.Model = dependenciesAbstractFactory.CreateModelFactory().Create();
// Indices
// j
this.j = indicesAbstractFactory.CreatejFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => x.Item1)
.Select(x => indexElementsAbstractFactory.CreatejIndexElementFactory().Create(x))
.ToImmutableList());
// r
this.r = indicesAbstractFactory.CreaterFactory().Create(
this.Context.OperatingRooms
.Entry
.Where(x => x.Resource is Location)
.Select(x => indexElementsAbstractFactory.CreaterIndexElementFactory().Create((Location)x.Resource))
.ToImmutableList());
// s
this.s = indicesAbstractFactory.CreatesFactory().Create(
this.Context.Surgeons
.Entry
.Where(x => x.Resource is Organization)
.Select(x => indexElementsAbstractFactory.CreatesIndexElementFactory().Create((Organization)x.Resource))
.ToImmutableList());
// t
this.t = indicesAbstractFactory.CreatetFactory().Create(
this.Context.PlanningHorizon
.Select(x => indexElementsAbstractFactory.CreatetIndexElementFactory().Create(
x.Key.Value.Value,
x.Value))
.ToImmutableList());
// Λ
this.Λ = indicesAbstractFactory.CreateΛFactory().Create(
this.Context.Scenarios
.Select(x => indexElementsAbstractFactory.CreateΛIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// υ1
this.υ1 = indicesAbstractFactory.Createυ1Factory().Create(
this.Context.OperatingRoomServiceLevels
.Select(x => indexElementsAbstractFactory.Createυ1IndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// Cross joins
// sΛ
this.sΛ = crossJoinsAbstractFactory.CreatesΛFactory().Create(
this.s.Value
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesΛCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// Parameters
// A(s, υ1)
this.A = parametersAbstractFactory.CreateAFactory().Create(
this.Context.SurgeonServiceLevelNumberTimeBlocks
.Select(x => parameterElementsAbstractFactory.CreateAParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.υ1.GetElementAt((PositiveInt)x.Item2),
(PositiveInt)x.Item3))
.ToImmutableList());
// N(s)
this.N = parametersAbstractFactory.CreateNFactory().Create(
this.Context.SurgeonStrategicTargets
.Select(x => parameterElementsAbstractFactory.CreateNParameterElementFactory().Create(
this.s.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// n(s, Λ)
this.n = parametersAbstractFactory.CreatenFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatients
.Select(x => parameterElementsAbstractFactory.CreatenParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(PositiveInt)x.Item3))
.ToImmutableList());
// Δ(j)
this.Δ = parametersAbstractFactory.CreateΔFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => parameterElementsAbstractFactory.CreateΔParameterElementFactory().Create(
this.j.GetElementAt(x.Item1),
x.Item2.Select(i => this.s.GetElementAt(i)).ToImmutableList()))
.ToImmutableList());
// Ρ(Λ)
this.Ρ = parametersAbstractFactory.CreateΡFactory().Create(
this.Context.ScenarioProbabilities
.Select(x => parameterElementsAbstractFactory.CreateΡParameterElementFactory().Create(
this.Λ.GetElementAt((PositiveInt)x.Key),
(FhirDecimal)x.Value))
.ToImmutableList());
// υ1Star
this.υ1Star = parametersAbstractFactory.Createυ1StarFactory().Create(
this.υ1.GetElementAt((PositiveInt)this.Context.OptimalOperatingRoomServiceLevel));
// υ1StarPlus1
this.υ1StarPlus1 = parametersAbstractFactory.Createυ1StarPlus1Factory().Create(
this.υ1.GetElementAt((PositiveInt)this.Context.NextOperatingRoomServiceLevel));
// ψ(t)
this.ψ = parametersAbstractFactory.CreateψFactory().Create(
this.Context.DayAvailabilities
.Select(x => parameterElementsAbstractFactory.CreateψParameterElementFactory().Create(
this.t.GetElementAt(x.Key),
(FhirBoolean)x.Value))
.ToImmutableList());
// ω(s)
this.ω = parametersAbstractFactory.CreateωFactory().Create(
this.Context.SurgeonPenaltyWeights
.Select(x => parameterElementsAbstractFactory.CreateωParameterElementFactory().Create(
this.s.GetElementAt(x.Key),
(FhirDecimal)x.Value))
.ToImmutableList());
// Variables
// B(s)
this.B = variablesAbstractFactory.CreateBFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
lowerBoundGenerator: (a) => 0,
upperBoundGenerator: (a) => int.MaxValue,
variableTypeGenerator: (a) => VariableType.Integer));
// e1Minus(s, Λ)
this.e1Minus = variablesAbstractFactory.Createe1MinusFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.Λ.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => int.MaxValue,
variableTypeGenerator: (a, b) => VariableType.Integer));
// e1Plus(s, Λ)
this.e1Plus = variablesAbstractFactory.Createe1PlusFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.Λ.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => int.MaxValue,
variableTypeGenerator: (a, b) => VariableType.Integer));
// Constraints
// Constraints 1L
this.Model.AddConstraints(
this.s.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints1LConstraintElementFactory().Create(
x,
this.υ1Star.Value,
this.A,
this.B)
.Value));
// Constraints 1U
this.Model.AddConstraints(
this.s.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints1UConstraintElementFactory().Create(
x,
this.υ1StarPlus1.Value,
this.A,
this.B)
.Value));
// Constraint 2
this.Model.AddConstraint(
constraintsAbstractFactory.CreateConstraint2Factory().Create(
this.r,
this.s,
this.t,
this.ψ,
this.B)
.Value);
// Constraints 3
this.Model.AddConstraints(
this.sΛ.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints3ConstraintElementFactory().Create(
x.sIndexElement,
x.ΛIndexElement,
this.N,
this.n,
this.B,
this.e1Minus,
this.e1Plus)
.Value));
}
public HM3BModel(
IConstraintElementsAbstractFactory constraintElementsAbstractFactory,
IConstraintsAbstractFactory constraintsAbstractFactory,
ICrossJoinElementsAbstractFactory crossJoinElementsAbstractFactory,
ICrossJoinsAbstractFactory crossJoinsAbstractFactory,
IDependenciesAbstractFactory dependenciesAbstractFactory,
IIndexElementsAbstractFactory indexElementsAbstractFactory,
IIndicesAbstractFactory indicesAbstractFactory,
IObjectiveFunctionsAbstractFactory objectiveFunctionsAbstractFactory,
IParameterElementsAbstractFactory parameterElementsAbstractFactory,
IParametersAbstractFactory parametersAbstractFactory,
IVariablesAbstractFactory variablesAbstractFactory,
IHM3BInputContext HM3BInputContext)
{
this.Context = HM3BInputContext;
this.Model = dependenciesAbstractFactory.CreateModelFactory().Create();
// Indices
// d
this.d = indicesAbstractFactory.CreatedFactory().Create(
this.Context.Weekdays
.Select(x => indexElementsAbstractFactory.CreatedIndexElementFactory().Create(
(PositiveInt)x))
.ToImmutableList());
// j
this.j = indicesAbstractFactory.CreatejFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => x.Item1)
.Select(x => indexElementsAbstractFactory.CreatejIndexElementFactory().Create(x))
.ToImmutableList());
// l
this.l = indicesAbstractFactory.CreatelFactory().Create(
this.Context.LengthOfStayDays
.Select(x => indexElementsAbstractFactory.CreatelIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// m
this.m = indicesAbstractFactory.CreatemFactory().Create(
this.Context.Machines
.Entry
.Where(x => x.Resource is Device)
.Select(x => indexElementsAbstractFactory.CreatemIndexElementFactory().Create((Device)x.Resource))
.ToImmutableList());
// r
this.r = indicesAbstractFactory.CreaterFactory().Create(
this.Context.OperatingRooms
.Entry
.Where(x => x.Resource is Location)
.Select(x => indexElementsAbstractFactory.CreaterIndexElementFactory().Create((Location)x.Resource))
.ToImmutableList());
// s
this.s = indicesAbstractFactory.CreatesFactory().Create(
this.Context.Surgeons
.Entry
.Where(x => x.Resource is Organization)
.Select(x => indexElementsAbstractFactory.CreatesIndexElementFactory().Create((Organization)x.Resource))
.ToImmutableList());
// t
this.t = indicesAbstractFactory.CreatetFactory().Create(
this.Context.PlanningHorizon
.Select(x => indexElementsAbstractFactory.CreatetIndexElementFactory().Create(
x.Key.Value.Value,
x.Value))
.ToImmutableList());
// Λ
this.Λ = indicesAbstractFactory.CreateΛFactory().Create(
this.Context.Scenarios
.Select(x => indexElementsAbstractFactory.CreateΛIndexElementFactory().Create((PositiveInt)x))
.ToImmutableList());
// Cross joins
// dt
this.dt = crossJoinsAbstractFactory.CreatedtFactory().Create(
this.d.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatedtCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// mr
this.mr = crossJoinsAbstractFactory.CreatemrFactory().Create(
this.m.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatemrCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// rt
this.rt = crossJoinsAbstractFactory.CreatertFactory().Create(
this.r.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatertCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// sd
this.sd = crossJoinsAbstractFactory.CreatesdFactory().Create(
this.s.Value
.SelectMany(b => this.d.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesdCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// sdt
this.sdt = crossJoinsAbstractFactory.CreatesdtFactory().Create(
this.s.Value
.SelectMany(b => this.d.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesdCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesdtCrossJoinElementFactory().Create(a.sIndexElement, a.dIndexElement, b))
.ToImmutableList());
// slΛ
this.slΛ = crossJoinsAbstractFactory.CreateslΛFactory().Create(
this.s.Value
.SelectMany(b => this.l.Value, (a, b) => crossJoinElementsAbstractFactory.CreateslCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreateslΛCrossJoinElementFactory().Create(a.sIndexElement, a.lIndexElement, b))
.ToImmutableList());
// sr
this.sr = crossJoinsAbstractFactory.CreatesrFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// srd
this.srd = crossJoinsAbstractFactory.CreatesrdFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.d.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrdCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// srj
this.srj = crossJoinsAbstractFactory.CreatesrjFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.j.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrjCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// srt
this.srt = crossJoinsAbstractFactory.CreatesrtFactory().Create(
this.s.Value
.SelectMany(b => this.r.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrCrossJoinElementFactory().Create(a, b))
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesrtCrossJoinElementFactory().Create(a.sIndexElement, a.rIndexElement, b))
.ToImmutableList());
// st
this.st = crossJoinsAbstractFactory.CreatestFactory().Create(
this.s.Value
.SelectMany(b => this.t.Value, (a, b) => crossJoinElementsAbstractFactory.CreatestCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// sΛ
this.sΛ = crossJoinsAbstractFactory.CreatesΛFactory().Create(
this.s.Value
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreatesΛCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// tΛ
this.tΛ = crossJoinsAbstractFactory.CreatetΛFactory().Create(
this.t.Value
.SelectMany(b => this.Λ.Value, (a, b) => crossJoinElementsAbstractFactory.CreatetΛCrossJoinElementFactory().Create(a, b))
.ToImmutableList());
// Parameters
// BBar(j)
this.BBar = parametersAbstractFactory.CreateBBarFactory().Create(
this.Context.SurgicalSpecialtyNumberAssignedTimeBlocks
.Select(x => parameterElementsAbstractFactory.CreateBBarParameterElementFactory().Create(
this.j.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// B(s)
this.Bs = parametersAbstractFactory.CreateBsFactory().Create(
this.Context.SurgeonNumberAssignedTimeBlocks
.Select(x => parameterElementsAbstractFactory.CreateBsParameterElementFactory().Create(
this.s.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// H
this.H = parametersAbstractFactory.CreateHFactory().Create(
this.Context.TimeBlockLength);
// h(s, Λ)
this.h = parametersAbstractFactory.CreatehFactory().Create(
this.Context.SurgeonScenarioWeightedAverageSurgicalDurations
.Select(x => parameterElementsAbstractFactory.CreatehParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
x.Item3))
.ToImmutableList());
// L(s)
this.L = parametersAbstractFactory.CreateLFactory().Create(
this.Context.SurgeonLengthOfStayMaximums
.Select(x => parameterElementsAbstractFactory.CreateLParameterElementFactory().Create(
this.s.GetElementAt(x.Key),
(PositiveInt)x.Value))
.ToImmutableList());
// n(s, Λ)
this.n = parametersAbstractFactory.CreatenFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatients
.Select(x => parameterElementsAbstractFactory.CreatenParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(PositiveInt)x.Item3))
.ToImmutableList());
// p(s, l, Λ)
this.p = parametersAbstractFactory.CreatepFactory().Create(
this.Context.SurgeonDayScenarioLengthOfStayProbabilities
.Select(x => parameterElementsAbstractFactory.CreatepParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.l.GetElementAt((PositiveInt)x.Item2),
this.Λ.GetElementAt((PositiveInt)x.Item3),
(FhirDecimal)x.Item4))
.ToImmutableList());
// W
this.W = parametersAbstractFactory.CreateWFactory().Create(
(PositiveInt)this.Context.NumberDaysPerWeek);
// Δ(j)
this.Δ = parametersAbstractFactory.CreateΔFactory().Create(
this.Context.SurgicalSpecialties
.Select(x => parameterElementsAbstractFactory.CreateΔParameterElementFactory().Create(
this.j.GetElementAt(x.Item1),
x.Item2.Select(i => this.s.GetElementAt(i)).ToImmutableList()))
.ToImmutableList());
// ζ(s, m)
this.ζ = parametersAbstractFactory.CreateζFactory().Create(
this.Context.SurgeonMachineRequirements
.Select(x => parameterElementsAbstractFactory.CreateζParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.m.GetElementAt(x.Item2),
(FhirBoolean)x.Item3))
.ToImmutableList());
// μ(s, Λ)
this.μ = parametersAbstractFactory.CreateμFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatientMeans
.Select(x => parameterElementsAbstractFactory.CreateμParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(FhirDecimal)x.Item3))
.ToImmutableList());
// Ρ(Λ)
this.Ρ = parametersAbstractFactory.CreateΡFactory().Create(
this.Context.ScenarioProbabilities
.Select(x => parameterElementsAbstractFactory.CreateΡParameterElementFactory().Create(
this.Λ.GetElementAt((PositiveInt)x.Key),
(FhirDecimal)x.Value))
.ToImmutableList());
// σ(s, Λ)
this.σ = parametersAbstractFactory.CreateσFactory().Create(
this.Context.SurgeonScenarioMaximumNumberPatientStandardDeviations
.Select(x => parameterElementsAbstractFactory.CreateσParameterElementFactory().Create(
this.s.GetElementAt(x.Item1),
this.Λ.GetElementAt((PositiveInt)x.Item2),
(FhirDecimal)x.Item3))
.ToImmutableList());
// ψ(t)
this.ψ = parametersAbstractFactory.CreateψFactory().Create(
this.Context.DayAvailabilities
.Select(x => parameterElementsAbstractFactory.CreateψParameterElementFactory().Create(
this.t.GetElementAt(x.Key),
(FhirBoolean)x.Value))
.ToImmutableList());
// Ω
this.Ω = parametersAbstractFactory.CreateΩFactory().Create(
(PositiveInt)this.Context.MaximumNumberRecoveryWardBeds);
// Variables
// b(s, r)
this.b = variablesAbstractFactory.CreatebFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.r.Value,
lowerBoundGenerator: (a, b) => 0,
upperBoundGenerator: (a, b) => int.MaxValue,
variableTypeGenerator: (a, b) => VariableType.Integer));
// u(s, d)
this.u = variablesAbstractFactory.CreateuFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().CreateBinary(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.d.Value));
// x(s, r, t)
this.x = variablesAbstractFactory.CreatexFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().CreateBinary(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.r.Value,
indexSet3: this.t.Value));
// z(s, t)
this.z = variablesAbstractFactory.CreatezFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().CreateBinary(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.t.Value));
// β(s, r, d)
this.β = variablesAbstractFactory.CreateβFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().Create(
model: this.Model,
indexSet1: this.s.Value,
indexSet2: this.r.Value,
indexSet3: this.d.Value,
lowerBoundGenerator: (a, b, c) => 0,
upperBoundGenerator: (a, b, c) => int.MaxValue,
variableTypeGenerator: (a, b, c) => VariableType.Integer));
// γ(r, t)
this.γ = variablesAbstractFactory.CreateγFactory().Create(
dependenciesAbstractFactory.CreateVariableCollectionFactory().CreateBinary(
model: this.Model,
indexSet1: this.r.Value,
indexSet2: this.t.Value));
// Constraints
// 1 (v, y), 5 L/M/U (v, w, y), 10 (y)
// 2, 3, 4, 6, 7, 8, 9
// Constraints 2
this.Model.AddConstraints(
this.s.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints2ConstraintElementFactory().Create(
x,
this.r,
this.Bs,
this.b)
.Value));
// Constraints 3
this.Model.AddConstraints(
this.r.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints3ConstraintElementFactory().Create(
x,
this.s,
this.t,
this.b,
this.γ)
.Value));
// Constraints 4
this.Model.AddConstraints(
this.rt.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints4ConstraintElementFactory().Create(
x.rIndexElement,
x.tIndexElement,
this.ψ,
this.γ)
.Value));
// Constraints 6
this.Model.AddConstraints(
this.sr.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints6ConstraintElementFactory().Create(
x.rIndexElement,
x.sIndexElement,
this.d,
this.b,
this.β)
.Value));
// Constraints 7
this.Model.AddConstraints(
this.srd.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints7ConstraintElementFactory().Create(
x.dIndexElement,
x.rIndexElement,
x.sIndexElement,
this.t,
this.W,
this.x,
this.β)
.Value));
// Constraints 8L
this.Model.AddConstraints(
this.sd.Value
.SelectMany(b =>
this.t.GetNthElementsAt(
b.dIndexElement.Value.Value.Value,
this.t.GetT() - this.W.Value.Value.Value + b.dIndexElement.Value.Value.Value,
this.W.Value.Value.Value),
(a, b) => Tuple.Create(a.sIndexElement, a.dIndexElement, b))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints8LConstraintElementFactory().Create(
x.Item2,
x.Item1,
x.Item3,
this.r,
this.x,
this.z)
.Value));
// Constraints 8U
this.Model.AddConstraints(
this.sd.Value
.SelectMany(b =>
this.t.GetNthElementsAt(
b.dIndexElement.Value.Value.Value,
this.t.GetT() - this.W.Value.Value.Value + b.dIndexElement.Value.Value.Value,
this.W.Value.Value.Value),
(a, b) => Tuple.Create(a.sIndexElement, a.dIndexElement, b))
.Select(
x => constraintElementsAbstractFactory.CreateConstraints8UConstraintElementFactory().Create(
x.Item2,
x.Item1,
x.Item3,
this.u,
this.z)
.Value));
// Constraints 9
this.Model.AddConstraints(
this.rt.Value
.Select(
x => constraintElementsAbstractFactory.CreateConstraints9ConstraintElementFactory().Create(
x.rIndexElement,
x.tIndexElement,
this.s,
this.x,
this.γ)
.Value));
// Objective function
this.Model.AddObjective(
objectiveFunctionsAbstractFactory.CreateObjectiveFunctionFactory().Create(
dependenciesAbstractFactory.CreateObjectiveFactory(),
this.sd,
this.u)
.Value);
}