public static TreeEvent FindTreeEvent(this TreeEvent treeEvent, Func <TreeEvent, bool> findAction)
{
if (findAction(treeEvent))
{
return(treeEvent);
}
if (treeEvent.FailingEvent != null)
{
var possibleEvent = FindTreeEvent(treeEvent.FailingEvent, findAction);
if (possibleEvent != null)
{
return(possibleEvent);
}
}
if (treeEvent.PassingEvent != null)
{
var possibleEvent = FindTreeEvent(treeEvent.PassingEvent, findAction);
if (possibleEvent != null)
{
return(possibleEvent);
}
}
return(null);
}
public static FragilityCurve GetFragilityCurve(this TreeEvent treeEvent, IEnumerable <double> waterLevels)
{
switch (treeEvent.ProbabilitySpecificationType)
{
case ProbabilitySpecificationType.Classes:
var classCurve = new FragilityCurve();
foreach (var waterLevel in waterLevels)
{
classCurve.Add(new FragilityCurveElement(waterLevel, GetClassesBasedProbabilityForWaterLevel(treeEvent, waterLevel)));
}
return(classCurve);
case ProbabilitySpecificationType.FixedFreqeuncy:
// TODO: Interpolate if necessary
return(treeEvent.FixedFragilityCurve);
case ProbabilitySpecificationType.FixedValue:
var curve = new FragilityCurve();
foreach (var waterLevel in waterLevels)
{
curve.Add(new FragilityCurveElement(waterLevel, treeEvent.FixedProbability));
}
return(curve);
default:
throw new NotImplementedException();
}
}
public static FragilityCurve GetLowerFragilityCurve(this TreeEvent treeEvent, IEnumerable <double> orderedWaterLevels)
{
if (treeEvent.ProbabilitySpecificationType == ProbabilitySpecificationType.Classes)
{
return(treeEvent.GetClassBasedLowerFragilityCurve(orderedWaterLevels));
}
return(treeEvent.GetFragilityCurve(orderedWaterLevels));
}
public static FragilityCurve GetClassBasedLowerFragilityCurve(this TreeEvent treeEvent, IEnumerable <double> waterLevels)
{
var curve = new FragilityCurve();
foreach (var waterLevel in waterLevels)
{
curve.Add(new FragilityCurveElement(waterLevel, GetClassesBasedProbabilityForWaterLevel(treeEvent, waterLevel, e => e.MinEstimation)));
}
return(curve);
}
public static Probability GetClassesBasedProbabilityForWaterLevel(this TreeEvent treeEvent, double waterLevel, Func <ExpertClassEstimation, ProbabilityClass> getProbabilityClassFunc = null)
{
if (getProbabilityClassFunc == null)
{
getProbabilityClassFunc = (e) => e.AverageEstimation;
}
var relevantEstimations = treeEvent.ClassesProbabilitySpecification.Where(e => Math.Abs(e.HydraulicCondition.WaterLevel - waterLevel) < 1e-8).ToArray();
if (relevantEstimations.Length == 0)
{
return(Probability.NaN);
}
var allRelevantEstimations = relevantEstimations.Select(e => ClassToProbabilityDouble(getProbabilityClassFunc(e))).ToArray();
var validEstimation = allRelevantEstimations.Where(e => !Double.IsNaN(e)).ToArray();
return(validEstimation.Any() ? (Probability)validEstimation.Average() : Probability.NaN);
}
public static IEnumerable <TreeEvent> GetAllEventsRecursive(this TreeEvent treeEvent)
{
var list = new[] { treeEvent };
if (treeEvent == null)
{
return(list);
}
if (treeEvent.FailingEvent != null)
{
list = list.Concat(GetAllEventsRecursive(treeEvent.FailingEvent)).ToArray();
}
if (treeEvent.PassingEvent != null)
{
list = list.Concat(GetAllEventsRecursive(treeEvent.PassingEvent)).ToArray();
}
return(list);
}