private IEnumerable<SparseTwinIndex<float>> GPUDoublyConstrained(IZone[] zones, IEnumerator<SparseArray<float>> ep, IEnumerator<SparseArray<float>> ea, IEnumerator<IDemographicCategory> ec)
{
using ( MultiRunGPUGravityModel gm =
new MultiRunGPUGravityModel( zones.Length, ( p => this.Progress = p ), this.Epsilon, this.MaxIterations ) )
{
float[] friction = null;
while ( ep.MoveNext() && ea.MoveNext() && ec.MoveNext() )
{
var production = ep.Current;
var attraction = ea.Current;
var cat = ec.Current;
var ret = production.CreateSquareTwinArray<float>();
friction = this.ComputeFriction( zones, cat, production.GetFlatData(), attraction.GetFlatData(), friction );
yield return gm.ProcessFlow( friction, production, attraction );
}
}
}
public IEnumerable<SparseTwinIndex<float>> Distribute(IEnumerable<SparseArray<float>> productions, IEnumerable<SparseArray<float>> attractions, IEnumerable<IDemographicCategory> cat)
{
var productionEnum = productions.GetEnumerator();
var attractionEnum = attractions.GetEnumerator();
var catEnum = cat.GetEnumerator();
var numberOfZones = this.Root.ZoneSystem.ZoneArray.GetFlatData().Length;
using ( MultiRunGPUGravityModel multiDist = new MultiRunGPUGravityModel( numberOfZones, ( (p) => this.Progress = p ), this.Epsilon, this.MaxIterations ) )
{
var zoneArray = this.Root.ZoneSystem.ZoneArray;
float[] friction = null;
while ( productionEnum.MoveNext() && attractionEnum.MoveNext() && catEnum.MoveNext() )
{
yield return multiDist.ProcessFlow( ( friction = this.ComputeFriction( zoneArray.GetFlatData(), catEnum.Current, friction ) ), productionEnum.Current, attractionEnum.Current );
}
friction = null;
}
}
