public void Plot_PTB_Results()
{
if (Direct_Data == null && IS_Data == null && Receiver == null && Parameter_Choice.Text != "Sabine RT" && Parameter_Choice.Text != "Eyring RT")
{
return;
}
double[] Schroeder;
string[] paramtype = new string [] { "T30/s", "EDT/s", "D/%", "C/dB", "TS/ms", "G/dB", "LF%", "LFC%", "IACC" };//LF/% LFC/% IACC
string ReceiverLine = "Receiver{0};";
double[, , ,] ParamValues = new double[SourceList.Items.Count, Recs.Length, 8, paramtype.Length];
for (int s = 0; s < Direct_Data.Length; s++)
{
for (int r = 0; r < Recs.Length; r++)
{
ReceiverLine += r.ToString() + ";";
for (int oct = 0; oct < 8; oct++)
{
double[] ETC = IR_Construction.ETCurve(Direct_Data, IS_Data, Receiver, (int)(CO_TIME.Value / 1000), SampleRate, oct, r, s, false);
Schroeder = AcousticalMath.Schroeder_Integral(ETC);
ParamValues[s, r, oct, 0] = AcousticalMath.T_X(Schroeder, 30, SampleRate);
ParamValues[s, r, oct, 1] = AcousticalMath.EarlyDecayTime(Schroeder, SampleRate);
ParamValues[s, r, oct, 2] = AcousticalMath.Definition(ETC, SampleRate, 0.05, Direct_Data[s].Min_Time(r), false);
ParamValues[s, r, oct, 3] = AcousticalMath.Clarity(ETC, SampleRate, 0.08, Direct_Data[s].Min_Time(r), false);
ParamValues[s, r, oct, 4] = AcousticalMath.Center_Time(ETC, SampleRate, Direct_Data[s].Min_Time(r)) * 1000;
ParamValues[s, r, oct, 5] = AcousticalMath.Strength(ETC, Direct_Data[s].SWL[oct], false);
double azi, alt;
PachTools.World_Angles(Direct_Data[s].Src.Origin(), Recs[r], true, out alt, out azi);
double[][] Lateral_ETC = IR_Construction.ETCurve_1d(Direct_Data, IS_Data, Receiver, CutoffTime, SampleRate, oct, r, new System.Collections.Generic.List <int> {
s
}, false, alt, azi, true);
ParamValues[s, r, oct, 6] = AcousticalMath.Lateral_Fraction(ETC, Lateral_ETC, SampleRate, Direct_Data[s].Min_Time(r), false);
}
}
}
System.Windows.Forms.SaveFileDialog sf = new System.Windows.Forms.SaveFileDialog();
sf.DefaultExt = ".txt";
sf.AddExtension = true;
sf.Filter = "Text File (*.txt)|*.txt|" + "All Files|";
if (sf.ShowDialog() == System.Windows.Forms.DialogResult.OK)
{
try
{
System.IO.StreamWriter SW = new System.IO.StreamWriter(System.IO.File.Open(sf.FileName, System.IO.FileMode.Create));
SW.WriteLine("Pachyderm Acoustic Simulation Results");
SW.WriteLine("Saved {0}", System.DateTime.Now.ToString());
SW.WriteLine("Filename:{0}", Rhino.RhinoDoc.ActiveDoc.Name);
for (int oct = 0; oct < 8; oct++)
{
SW.WriteLine(string.Format(ReceiverLine, oct));
for (int param = 0; param < paramtype.Length; param++)
{
SW.Write(paramtype[param] + ";");
for (int i = 0; i < Direct_Data.Length; i++)
{
for (int q = 0; q < Recs.Length; q++)
{
SW.Write(ParamValues[i, q, oct, param].ToString() + ";");
}
}
SW.WriteLine("");
}
}
SW.Close();
}
catch (System.Exception)
{
Rhino.RhinoApp.WriteLine("File is open, and cannot be written over.");
return;
}
}
}
public void Plot_Results()
{
if (Direct_Data == null && IS_Data == null && Receiver == null && Parameter_Choice.Text != "Sabine RT" && Parameter_Choice.Text != "Eyring RT")
{
return;
}
System.Windows.Forms.SaveFileDialog sf = new System.Windows.Forms.SaveFileDialog();
sf.DefaultExt = ".txt";
sf.AddExtension = true;
sf.Filter = "Text File (*.txt)|*.txt|" + "All Files|";
if (sf.ShowDialog() == System.Windows.Forms.DialogResult.OK)
{
System.IO.StreamWriter SW;
fileread:
try
{
SW = new System.IO.StreamWriter(System.IO.File.Open(sf.FileName, System.IO.FileMode.Create));
}
catch
{
System.Windows.Forms.DialogResult dr = System.Windows.Forms.MessageBox.Show("File is in use. Close any programs using the file, and try again.", "File In Use", System.Windows.Forms.MessageBoxButtons.RetryCancel);
if (dr == System.Windows.Forms.DialogResult.Cancel)
{
return;
}
goto fileread;
}
double[] Schroeder;
double[, ,] EDT = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] T10 = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] T15 = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] T20 = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] T30 = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] G = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] C80 = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] C50 = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] D50 = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] TS = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] LF = new double[this.SourceList.Items.Count, Recs.Length, 8];
double[, ,] LE = new double[this.SourceList.Items.Count, Recs.Length, 8];
for (int s = 0; s < SourceList.Items.Count; s++)
{
for (int r = 0; r < Recs.Length; r++)
{
for (int oct = 0; oct < 8; oct++)
{
double[] ETC = IR_Construction.ETCurve(Direct_Data, IS_Data, Receiver, (int)(CO_TIME.Value / 1000), SampleRate, oct, r, s, false);
Schroeder = AcousticalMath.Schroeder_Integral(ETC);
EDT[s, r, oct] = AcousticalMath.EarlyDecayTime(Schroeder, SampleRate);
T10[s, r, oct] = AcousticalMath.T_X(Schroeder, 10, SampleRate);
T15[s, r, oct] = AcousticalMath.T_X(Schroeder, 15, SampleRate);
T20[s, r, oct] = AcousticalMath.T_X(Schroeder, 20, SampleRate);
T30[s, r, oct] = AcousticalMath.T_X(Schroeder, 30, SampleRate);
G[s, r, oct] = AcousticalMath.Strength(ETC, Direct_Data[s].SWL[oct], false);
C50[s, r, oct] = AcousticalMath.Clarity(ETC, SampleRate, 0.05, Direct_Data[s].Min_Time(r), false);
C80[s, r, oct] = AcousticalMath.Clarity(ETC, SampleRate, 0.08, Direct_Data[s].Min_Time(r), false);
D50[s, r, oct] = AcousticalMath.Definition(ETC, SampleRate, 0.05, Direct_Data[s].Min_Time(r), false);
TS[s, r, oct] = AcousticalMath.Center_Time(ETC, SampleRate, Direct_Data[s].Min_Time(r)) * 1000;
double[] L_ETC = IR_Construction.ETCurve_1d(Direct_Data, IS_Data, Receiver, (int)(CO_TIME.Value), SampleRate, oct, r, new System.Collections.Generic.List <int>()
{
s
}, false, (double)this.Alt_Choice.Value, (double)this.Azi_Choice.Value, true)[1];
LF[s, r, oct] = AcousticalMath.Lateral_Fraction(ETC, L_ETC, SampleRate, Direct_Data[s].Min_Time(r), false) * 1000;
LE[s, r, oct] = AcousticalMath.Lateral_Efficiency(ETC, L_ETC, SampleRate, Direct_Data[s].Min_Time(r), false) * 1000;
}
}
}
SW.WriteLine("Pachyderm Acoustic Simulation Results");
SW.WriteLine("Saved {0}", System.DateTime.Now.ToString());
SW.WriteLine("Filename:{0}", Rhino.RhinoDoc.ActiveDoc.Name);
for (int s = 0; s < SourceList.Items.Count; s++)
{
SW.WriteLine("Source {0};", s);
for (int r = 0; r < Recs.Length; r++)
{
SW.WriteLine("Receiver {0};63 Hz.;125 Hz.;250 Hz.;500 Hz.;1000 Hz.; 2000 Hz.;4000 Hz.; 8000 Hz.", r);
SW.WriteLine("Early Decay Time (EDT);{0};{1};{2};{3};{4};{5};{6};{7}", EDT[s, r, 0], EDT[s, r, 1], EDT[s, r, 2], EDT[s, r, 3], EDT[s, r, 4], EDT[s, r, 5], EDT[s, r, 6], EDT[s, r, 7]);
SW.WriteLine("Reverberation Time (T-10);{0};{1};{2};{3};{4};{5};{6};{7}", T10[s, r, 0], T10[s, r, 1], T10[s, r, 2], T10[s, r, 3], T10[s, r, 4], T10[s, r, 5], T10[s, r, 6], T10[s, r, 7]);
SW.WriteLine("Reverberation Time (T-15);{0};{1};{2};{3};{4};{5};{6};{7}", T15[s, r, 0], T15[s, r, 1], T15[s, r, 2], T15[s, r, 3], T15[s, r, 4], T15[s, r, 5], T15[s, r, 6], T15[s, r, 7]);
SW.WriteLine("Reverberation Time (T-20);{0};{1};{2};{3};{4};{5};{6};{7}", T20[s, r, 0], T20[s, r, 1], T20[s, r, 2], T20[s, r, 3], T20[s, r, 4], T20[s, r, 5], T20[s, r, 6], T20[s, r, 7]);
SW.WriteLine("Reverberation Time (T-30);{0};{1};{2};{3};{4};{5};{6};{7}", T30[s, r, 0], T30[s, r, 1], T30[s, r, 2], T30[s, r, 3], T30[s, r, 4], T30[s, r, 5], T30[s, r, 6], T30[s, r, 7]);
SW.WriteLine("Clarity - 50 ms (C-50);{0};{1};{2};{3};{4};{5};{6};{7}", C50[s, r, 0], C50[s, r, 1], C50[s, r, 2], C50[s, r, 3], C50[s, r, 4], C50[s, r, 5], C50[s, r, 6], C50[s, r, 7]);
SW.WriteLine("Definition/Deutlichkeit (D);{0};{1};{2};{3};{4};{5};{6};{7}", D50[s, r, 0], D50[s, r, 1], D50[s, r, 2], D50[s, r, 3], D50[s, r, 4], D50[s, r, 5], D50[s, r, 6], D50[s, r, 7]);
SW.WriteLine("Clarity - 80 ms (C-80);{0};{1};{2};{3};{4};{5};{6};{7}", C80[s, r, 0], C80[s, r, 1], C80[s, r, 2], C80[s, r, 3], C80[s, r, 4], C80[s, r, 5], C80[s, r, 6], C80[s, r, 7]);
SW.WriteLine("Center Time (TS);{0};{1};{2};{3};{4};{5};{6};{7}", TS[s, r, 0], TS[s, r, 1], TS[s, r, 2], TS[s, r, 3], TS[s, r, 4], TS[s, r, 5], TS[s, r, 6], TS[s, r, 7]);
SW.WriteLine("Strength(G);{0};{1};{2};{3};{4};{5};{6};{7}", G[s, r, 0], G[s, r, 1], G[s, r, 2], G[s, r, 3], G[s, r, 4], G[s, r, 5], G[s, r, 6], G[s, r, 7]);
SW.WriteLine("Lateral Fraction(LF);{0};{1};{2};{3};{4};{5};{6};{7}", LF[s, r, 0], LF[s, r, 1], LF[s, r, 2], LF[s, r, 3], LF[s, r, 4], LF[s, r, 5], LF[s, r, 6], LF[s, r, 7]);
SW.WriteLine("Lateral Efficiency(LE);{0};{1};{2};{3};{4};{5};{6};{7}", LE[s, r, 0], LE[s, r, 1], LE[s, r, 2], LE[s, r, 3], LE[s, r, 4], LE[s, r, 5], LE[s, r, 6], LE[s, r, 7]);
}
}
SW.Close();
}
}
