Exemplo n.º 1
0
        private void mnuScaleToFundamental_Click(object sender, EventArgs e)
        {
            if (Bore == null)
            {
                ShowError("You must calculate the impedance first.");
                return;
            }

            ScaleBoreToFundamental dlg = new ScaleBoreToFundamental((decimal)ImpedanceData.ImpedancePeakFrequencies[0]);
            DialogResult dr = dlg.ShowDialog(this);
            decimal targetFundamental = dlg.SelectedFundamental;

            if (dr == DialogResult.OK)
            {
                Bore previousBore = Bore;
                decimal currentFundamental = (decimal)ImpedanceData.ImpedancePeakFrequencies[0];
                while (currentFundamental != targetFundamental)
                {
                    decimal scaleAmount = currentFundamental / targetFundamental;
                    List<BoreSection> newBoreSections = new List<BoreSection>();
                    foreach (BoreSection boreSection in previousBore.BoreSections)
                        newBoreSections.Add(new BoreSection(boreSection.OpeningRadius, boreSection.ClosingRadius, boreSection.Length * scaleAmount));
                    Bore newBore = new Bore(newBoreSections);
                    if (targetFundamental > currentFundamental)
                    {
                        Complex lastImpedance = newBore.GetImpedance((double)currentFundamental);
                        for (decimal freq = currentFundamental + 1; true; freq++)
                        {
                            Complex currentImpedance = newBore.GetImpedance((double)freq);
                            if (currentImpedance.Magnitude < lastImpedance.Magnitude)
                            {
                                if (newBore.InputImpedance.Count == 2)
                                {
                                    newBore.GetImpedance((double)currentFundamental - 1);
                                }
                                newBore.FindResonances(2);
                                ImpedanceData tempImpedanceData = new ImpedanceData(newBore.InputImpedance);
                                currentFundamental = (decimal)tempImpedanceData.ImpedancePeakFrequencies[0];
                                previousBore = newBore;
                                break;
                            }
                            lastImpedance = currentImpedance;
                        }
                    }
                    else
                    {
                        Complex lastImpedance = newBore.GetImpedance(1);
                        for (decimal freq = 2; true; freq++)
                        {
                            Complex currentImpedance = newBore.GetImpedance((double)freq);
                            if (currentImpedance.Magnitude < lastImpedance.Magnitude)
                            {
                                newBore.FindResonances(2);
                                ImpedanceData tempImpedanceData = new ImpedanceData(newBore.InputImpedance);
                                currentFundamental =(decimal)tempImpedanceData.ImpedancePeakFrequencies[0];
                                previousBore = newBore;
                                break;
                            }
                            lastImpedance = currentImpedance;
                        }
                    }
                }

                didgePropertyEditor_DimensionsChanged();

                List<BoreSection> boreSections = new List<BoreSection>();
                foreach (BoreSection section in previousBore.BoreSections)
                    boreSections.Add(new BoreSection(section.OpeningRadius, section.ClosingRadius, Math.Round(section.Length, 5)));
                didgePropertyEditor.BoreSections = boreSections;
            }
        }
Exemplo n.º 2
0
        private void btnCalculate_Click(object sender, EventArgs e)
        {
            if (!didgeValid)
            {
                ShowError("There is an error with the entered dimensions.");
                return;
            }

            tabPlots.Enabled = true;
            verticalLines.Clear();

            calculatedFrequencyCount = 0;

            DidgeDesignProperties didgeProperties = didgePropertyEditor.DidgeDesignProperties;

            Bore bore = new Bore(didgeProperties.BoreSections);
            if (bore.Length == 0)
            {
                ShowError("The bore has a length of 0.");
                return;
            }
            bore.CalculatedFrequency += new Bore.CalculatedFrequencyDelegate(bore_CalculatedFrequency);

            DidgeData data = (DidgeData)treeDidgeHistory.SelectedNode.Tag;
            data.didge = didgePropertyEditor.DidgeDesignProperties;
            data.bore = bore;

            finishedThreads = 0;

            ParameterizedThreadStart threadStart = (threadNum) =>
            {
                //split the work up among however many threads we have
                int startFrequency = (int)Math.Round((2000.0 / settings.NumberOfThreads) * ((int)threadNum) + 1.0, 0);
                int endFrequency = (int)Math.Round((2000.0 / settings.NumberOfThreads) * ((int)threadNum + 1.0), 0);

                bore.CalculateInputImpedance(startFrequency, endFrequency, 1);
                SetImpedanceData();
            };

            for (int i=0; i<settings.NumberOfThreads; i++)
            {
                Thread thread = new Thread(threadStart);
                thread.Start(i);
            }

            progressDialog = new Progress();
            progressDialog.ShowDialog(this);
        }