//Birth
private static int[] ChooseMaleFathers(SimParams par, Population pop,
int[] vacant, List <int> notVacant)
{
/*Picks fathers to sire chicks into vacancies,
* Uses only living males.
* It is likely for high quality males to sire multiple
* offspring in one sim step for both local and global scopes*/
//Remove songless birds from notVacant, and mark all not in that list as unavailable
HashSet <int> PotentialFathersTemp = new HashSet <int>(notVacant.Where(x => pop.SyllableRepertoire[x] > 0));
HashSet <int> Unavailable = new HashSet <int>(Enumerable.Range(0, par.NumBirds));
Unavailable.ExceptWith(PotentialFathersTemp);
//Get the probabilities
List <int> PotentialFathers = PotentialFathersTemp.ToList();
float[] Probs = ReproductiveProbability(par, pop, PotentialFathers);
//pick fathers
int[] Fathers = new int[vacant.Length];
if (par.LocalBreeding)
{
float[] Probability = Enumerable.Repeat(0f, par.NumBirds).ToArray();
for (int i = 0; i < PotentialFathers.Count; i++)
{
Probability[PotentialFathers[i]] = Probs[i];
}
for (int i = 0; i < Fathers.Length; i++)
{
Fathers[i] = Locations.GetLocalBirds(par, pop, vacant[i],
Unavailable, probs: Probability)[0];
}
}
else
{
//float[] PotenProbs = new float[PotentialFathers.Count];
int[] FatherIndex;
//for(int i=0;i<PotenProbs.Length;i++){
// PotenProbs[i] = Probability[PotentialFathers[i]];
//}
FatherIndex = par.RandomSampleUnequal(Probs, vacant.Length, true);
for (int i = 0; i < Fathers.Length; i++)
{
Fathers[i] = PotentialFathers[FatherIndex[i]];
}
}
int[] lol = new int[Fathers.Length];
for (int i = 0; i < Fathers.Length; i++)
{
lol[i] = pop.SyllableRepertoire[Fathers[i]];
}
return(Fathers);
}
//Getting tutors, the original runs faster, but cannot get more than one tutor
private static int[] ChooseTutors(SimParams par, Population pop,
List <int> learners, List <int> notVacant)
{
//Picks tutos that are alive, not chicks, and not songless,
//Males can tutor multiple learners
/*remove chicks and songless birds + any Misc, mark all excluded
* birds as unavailable*/
HashSet <int> PotentialTutorsTemp = new HashSet <int>(notVacant.Where(x => pop.Age[x] > 0));
PotentialTutorsTemp.ExceptWith(PotentialTutorsTemp.Where(x => pop.SyllableRepertoire[x] == 0).ToArray());
//pick Tutors
int[] Tutors = new int[learners.Count];
if (par.LocalTutor)
{
//Set up unavailable for local testing
HashSet <int> Unavailable = new HashSet <int>(Enumerable.Range(0, par.NumBirds));
Unavailable.ExceptWith(PotentialTutorsTemp);
for (int i = 0; i < learners.Count; i++)
{
if (par.SocialCues)
{
Tutors[i] = Locations.GetLocalBirds(par, pop, learners[i], Unavailable, 1, pop.Bred)[0];
}
else
{
Tutors[i] = Locations.GetLocalBirds(par, pop, learners[i], Unavailable)[0];
}
}
}
else //drawn individually so no learner is his own tutor
{
List <int> PotentialTutors;
for (int i = 0; i < learners.Count; i++)
{
PotentialTutors = PotentialTutorsTemp.ToList();
PotentialTutors.Remove(learners[i]);
if (par.SocialCues)
{
float[] TutorProbs = new float[PotentialTutors.Count];
for (int j = 0; j < TutorProbs.Length; j++)
{
TutorProbs[j] = pop.Bred[PotentialTutors[j]];
}
Tutors[i] = PotentialTutors[par.RandomSampleUnequal(TutorProbs, 1)[0]];
}
else
{
Tutors[i] = PotentialTutors[par.RandomSampleEqualReplace(PotentialTutors, 1)[0]];
}
}
}
return(Tutors);
}
private static List <int>[] ChooseMultipleTutors(SimParams par, Population pop,
List <int> learners, List <int> notVacant, int numTutors)
{
/*remove chicks and songless birds + any Misc, mark all excluded
* birds as unavailable*/
HashSet <int> PotentialTutorsTemp = new HashSet <int>(notVacant.Where(x => pop.Age[x] > 0));
PotentialTutorsTemp.ExceptWith(PotentialTutorsTemp.Where(x => pop.SyllableRepertoire[x] == 0).ToArray());
List <int>[] Tutors = new List <int> [learners.Count];
if (par.LocalTutor)
{
HashSet <int> Unavailable = new HashSet <int>(Enumerable.Range(0, par.NumBirds));
Unavailable.ExceptWith(PotentialTutorsTemp);
for (int i = 0; i < learners.Count; i++)
{
if (par.SocialCues)
{
Tutors[i] = Locations.GetLocalBirds(par, pop, learners[i], Unavailable, numTutors, pop.Bred).ToList();
}
else
{
Tutors[i] = Locations.GetLocalBirds(par, pop, learners[i], Unavailable, numTutors).ToList();
}
}
}
else
{
List <int> PotentialTutors;
for (int i = 0; i < learners.Count; i++)
{
PotentialTutors = PotentialTutorsTemp.ToList();
PotentialTutors.Remove(learners[i]);
if (par.SocialCues)
{
float[] TutorProbs = new float[PotentialTutors.Count];
for (int j = 0; j < TutorProbs.Length; j++)
{
TutorProbs[j] = pop.Bred[PotentialTutors[j]];
}
int[] ChosenIndex = par.RandomSampleUnequal(TutorProbs, numTutors);
for (int j = 0; j < numTutors; j++)
{
ChosenIndex[j] = PotentialTutors[ChosenIndex[j]];
}
Tutors[i] = ChosenIndex.ToList();
}
else
{
Tutors[i] = par.RandomSampleEqualNoReplace(PotentialTutors, numTutors).ToList();
}
}
}
return(Tutors);
}
