protected virtual void setCols()
{
try
{//figure out config dates and get comp columns
HeightConfigCollection configcollection = new HeightConfigCollection(_collection);
List <IConfig> configs = configcollection.GetConfigs();
Console.WriteLine(" doing " + _summarytype);
if (configs.Count > 1)
{
Console.WriteLine("configs found " + configs.Count);
foreach (IConfig c in configs)
{
Console.WriteLine("start date " + c.StartDate);
}
throw new ApplicationException("Can not process grid. System not yet set up to handle multpile configurations");
}
_configstart = configs[0].StartDate;
Console.WriteLine("start date in summary grid " + _configstart);
//wscomps
//SortedDictionary<double, ISessionColumn> wscomps
// = _collection.GetColumnsByType(SessionColumnType.WSAvg, _configstart, true);
_lowerwscompcol = _collection[_collection.LowerWSComp(_configstart)];
_upperwscompcol = _collection[_collection.UpperWSComp(_configstart)];
_lowerwsht = _lowerwscompcol.getConfigAtDate(_configstart).Height;
_upperwsht = _upperwscompcol.getConfigAtDate(_configstart).Height;
//Console.WriteLine("Summary grid found " + wscomps.Values.Count + " ws comps");
//wdcomp
List <ISessionColumn> wdcomp =
_collection.GetColumnsByType(SessionColumnType.WDAvg);
var compcol = from c in wdcomp.AsEnumerable()
where c.IsComposite
select c;
_wdcompcol = compcol.ToList()[0];
Console.WriteLine("Summary grid found " + compcol.ToList().Count + " wd comps");
//shear
_shearwscol = _collection.GetColumnsByType(_sheartype)[0];
_shearht = _shearwscol.getConfigAtDate(_configstart).Height;
Console.WriteLine("Summary grid found " + _collection.GetColumnsByType(_sheartype).Count + " shear");
}
catch
{
throw;
}
}
//methods
public bool CalculateComposites()
{
//get configs
//add columns to datatable and column collection
//check deadranges
//compare sensor a aand b acorsss all configs
//add columns and update data through the dataview that is each config
SelectValue UseColumn = default(SelectValue);
double missing = Settings.Default.MissingValue;
IWindSpeedConfig A = new WindSpeedConfig();
int AIndex;
double AVal;
IWindSpeedConfig B = new WindSpeedConfig();
int BIndex;
double BVal;
ISessionColumn ACol;
ISessionColumn BCol;
int compositeIndex = 0;
IConfigCollection configCollection = new HeightConfigCollection(_collection);
string datecol = _data.Columns [_collection.DateIndex].ColumnName;
string wdcolname = Enum.GetName(typeof(SessionColumnType), SessionColumnType.WDAvg)
+ Settings.Default.CompColName;
int WDIndex = _collection[wdcolname].ColIndex;
double WDVal;
double workVal;
//find each config across height changes..
List <HeightConfig> workConfigs
= configCollection.GetConfigs().ConvertAll <HeightConfig>(t => (HeightConfig)t).ToList();
//Console.WriteLine("workconfigs count " + workConfigs.Count);
//loop through each config
foreach (HeightConfig htConfig in workConfigs)
{ //loop each ht in config
DataView configData = _data.AsDataView();
configData.Sort = Settings.Default.TimeStampName + " asc ";
//filter the data based on start and end of the current ht config
string filter = datecol + ">= '" + htConfig.StartDate.ToString() + "' and " + datecol + " <='"
+ htConfig.EndDate.ToString() + "'";
Console.WriteLine(filter);
configData.RowFilter = filter;
if (configData.Count > _data.Rows.Count)
{
throw new ApplicationException("Error filtering dataset by configuration dates");
}
int RowNum = configData.Count;
Console.WriteLine("datatable and filtered dataviewcounts: " + _data.Rows.Count + " " + RowNum);
List <double> resultArray = new List <double>(RowNum);
Console.WriteLine(" Heights found " + htConfig.Columns.Count);
//loop each height's columns in config
//each htconfig stores a collection of collumns with those hts
foreach (KeyValuePair <double, IList <ISessionColumn> > kv in htConfig.Columns)
{
var a = from f in kv.Value.AsEnumerable()
where f.ColumnType == SessionColumnType.WSAvg
select f;
List <ISessionColumn> thisresult = a.ToList();
//get sensor pairs that are ws avg types
if (a.Count() == 2)
{
Console.WriteLine(thisresult[0].ColumnType);
Console.WriteLine(thisresult[1].ColumnType);
ACol = thisresult[0];
BCol = thisresult[1];
AIndex = ACol.ColIndex;
BIndex = BCol.ColIndex;
//get the specific config for the date of the ht config
A = (IWindSpeedConfig)ACol.getConfigAtDate(htConfig.StartDate);
B = (IWindSpeedConfig)BCol.getConfigAtDate(htConfig.StartDate);
//change the config in the sensor config data if there are dead zones
CheckDeadRange(A, B);
Console.WriteLine("A Good Range: " + A.GoodSector.SectorStart + " to " + A.GoodSector.SectorEnd);
Console.WriteLine("A Bad Range: " + A.ShadowSector.SectorStart + " to " + A.ShadowSector.SectorEnd);
Console.WriteLine("B Good Range: " + B.GoodSector.SectorStart + " to " + B.GoodSector.SectorEnd);
Console.WriteLine("B Bad Range: " + B.ShadowSector.SectorStart + " to " + B.ShadowSector.SectorEnd);
}
else
{
continue;
}
double height = A.Height;
IList <ISessionColumn> thisColumnSet = kv.Value;
//add column to datatable and column collection if necessary
WindSpeedCompColumns wscompcols = new WindSpeedCompColumns(_collection, _data);
wscompcols.Add(height, thisColumnSet, htConfig.StartDate, htConfig.EndDate);
if (NewCompositeColumnAdded != null)
{
NewCompositeColumnAdded("Wind Speed Composite Columns Added");
}
#region decide on value
//begin calculating comps for each row
if (DeterminingWindSpeedCompositeValues != null)
{
DeterminingWindSpeedCompositeValues("Assigning Wind Speed Composite Values");
}
foreach (DataRowView rowview in configData)
{
if (!Double.TryParse(rowview[AIndex].ToString(), out AVal))
{
AVal = missing;
}
else
{
if (AVal < 0)
{
AVal = missing;
}
}
if (!Double.TryParse(rowview[BIndex].ToString(), out BVal))
{
BVal = missing;
}
else
{
if (BVal < 0)
{
BVal = missing;
}
}
if (!double.TryParse(rowview[WDIndex].ToString(), out WDVal) || WDVal < 0)
{
WDVal = missing;
UseColumn = SelectValue.missing;
goto assignvalue;
}
InSector SensorA = A.BelongsToSector(WDVal);
InSector SensorB = B.BelongsToSector(WDVal);
//in A not in B
if (SensorA.Equals(InSector.Not_Shadowed) &&
(SensorB.Equals(InSector.Shadowed) |
SensorB.Equals(InSector.Neither)))
{
workVal = AVal;
if (workVal.Equals(missing))
{
UseColumn = SelectValue.missing;
}
else
{
UseColumn = SelectValue.A;
}
goto assignvalue;
}
//in B not in A take b
if (SensorB.Equals(InSector.Not_Shadowed) &&
(SensorA.Equals(InSector.Shadowed) |
SensorA.Equals(InSector.Neither)))
{
workVal = BVal;
if (workVal < 0)
{
UseColumn = SelectValue.missing;
}
else
{
UseColumn = SelectValue.B;
}
goto assignvalue;
}
//in both A and B avg the two
if (SensorB.Equals(InSector.Not_Shadowed) &&
SensorA.Equals(InSector.Not_Shadowed))
{
if (AVal < 0 && BVal >= 0)
{
UseColumn = SelectValue.B;
goto assignvalue;
}
if (BVal < 0 && AVal >= 0)
{
UseColumn = SelectValue.A;
goto assignvalue;
}
//if wdval is on a boundary use the other sensor
if (WDVal == A.GoodSector.SectorStart | WDVal == A.GoodSector.SectorEnd)
{
UseColumn = SelectValue.B;
goto assignvalue;
}
if (WDVal == B.GoodSector.SectorStart | WDVal == B.GoodSector.SectorEnd)
{
UseColumn = SelectValue.A;
goto assignvalue;
}
List <double> avg = new List <double>()
{
BVal, AVal
};
UseColumn = SelectValue.both;
workVal = avg.Average();
goto assignvalue;
}
//otherwise set to missing
workVal = missing;
goto assignvalue;
#endregion
#region assign value
assignvalue:
{
double useval = missing;
double usevalB = missing;
string heightstring = height.ToString().Replace(".", "_");
string parentColName = heightstring + Enum.GetName(typeof(SessionColumnType), SessionColumnType.WSAvg) + Settings.Default.CompColName;
if (_collection[parentColName] != null)
{
compositeIndex = _collection[parentColName].ColIndex;
}
else
{
throw new ApplicationException(" Parent composite column " + parentColName + " not found");
}
switch (UseColumn)
{
case SelectValue.A:
rowview[compositeIndex] = AVal;
foreach (ISessionColumn child in ACol.ChildColumns)
{
string thischildcompname = heightstring
+ Enum.GetName(typeof(SessionColumnType), child.ColumnType)
+ Settings.Default.CompColName;
if (_collection[thischildcompname] != null)
{
if (!double.TryParse(rowview[child.ColIndex].ToString(), out useval))
{
rowview[thischildcompname] = missing;
}
else
{
rowview[thischildcompname] = useval;
}
}
else
{
throw new ApplicationException(" Child composite column " + thischildcompname + " not found");
}
}
break;
case SelectValue.B:
rowview[compositeIndex] = BVal;
foreach (ISessionColumn child in BCol.ChildColumns)
{
string thischildcompname = heightstring
+ child.ColumnType
+ Settings.Default.CompColName;
if (_collection[thischildcompname] != null)
{
if (!double.TryParse(rowview[child.ColIndex].ToString(), out useval))
{
rowview[thischildcompname] = missing;
}
else
{
rowview[thischildcompname] = useval;
}
}
else
{
throw new ApplicationException(" Child composite column " + thischildcompname + " not found");
}
}
break;
case SelectValue.missing:
rowview[compositeIndex] = missing;
//loop through comp cols for both a and b in case one column has more children than the other
foreach (ISessionColumn child in ACol.ChildColumns)
{
string thischildcompname = heightstring
+ Enum.GetName(typeof(SessionColumnType), child.ColumnType)
+ Settings.Default.CompColName;
if (_collection[thischildcompname] != null)
{
rowview[thischildcompname] = missing;
}
else
{
throw new ApplicationException(" Child composite column " + thischildcompname + " not found");
}
}
foreach (ISessionColumn child in BCol.ChildColumns)
{
string thischildcompname = heightstring
+ Enum.GetName(typeof(SessionColumnType), child.ColumnType)
+ Settings.Default.CompColName;
if (_collection[thischildcompname] != null)
{
rowview[thischildcompname] = missing;
}
else
{
throw new ApplicationException(" Child composite column " + thischildcompname + " not found");
}
}
break;
case SelectValue.both:
double averageWs = new double[2] {
AVal, BVal
}.Average();
rowview[compositeIndex] = averageWs;
//loop through comp cols for both a and b in case one column has more children than the other
foreach (ISessionColumn child in ACol.ChildColumns)
{
string thischildcompname = heightstring
+ Enum.GetName(typeof(SessionColumnType), child.ColumnType)
+ Settings.Default.CompColName;
if (_collection[thischildcompname] != null)
{
var result = from c in BCol.ChildColumns.AsEnumerable()
where c.ColumnType == child.ColumnType
select c;
List <ISessionColumn> columnfound = result.ToList();
if (columnfound.Count == 1)
{
//Console.WriteLine(" ****************** average taken");
if (!double.TryParse(rowview[columnfound[0].ColIndex].ToString(), out useval) &
!double.TryParse(rowview[child.ColIndex].ToString(), out usevalB))
{
rowview[thischildcompname] = missing;
}
else
{
if (child.ColumnType == SessionColumnType.WSStd)
{ //combine ws standard deviations more correctly than straight average
double pooled_var = (Math.Pow(useval, 2.0) + Math.Pow(usevalB, 2.0)
+ Math.Pow((AVal - averageWs), 2.0)
+ Math.Pow((BVal - averageWs), 2.0)) / 2;
double pooled_std = Math.Sqrt(pooled_var);
if (pooled_std < 100.0)
{
rowview[thischildcompname] = pooled_std;
}
else
{
rowview[thischildcompname] = missing;
}
}
else
{
rowview[thischildcompname] = new double[2] {
useval, usevalB
}.Average();
}
}
}
else
{
Console.WriteLine(" missing taken A value =" + AVal + " B Val " + BVal);
rowview[thischildcompname] = missing;
}
}
else
{
throw new ApplicationException(" Child composite column " + thischildcompname + " not found");
}
}
break;
default:
break;
}
}
}
#endregion
}
}
if (CompletedWindSpeedCompositeValues != null)
{
CompletedWindSpeedCompositeValues("Completed Generating Wind Speed Composites");
}
return(true);
}
