public void ObterPesosEAbascissasCheb()
        {
            const int numeroDePontos = 4;

            double[] momentosAnaliticos = CalcularMomentosAnaliticos(2 * numeroDePontos);

            var chebyChev = new ChebyshevModificado(numeroDePontos, momentosAnaliticos, ChebyshevModificado.TipoDeDominio.SemiInfinito);

            double[] pesos, abscissas;

            chebyChev.CalcularPesosEAbscissas(out pesos, out abscissas);

            AlgoritmoPD algoritmoPd = new AlgoritmoPD(numeroDePontos, momentosAnaliticos);

            double[] pesosEsperados;
            double[] abscissasEsperadas;

            algoritmoPd.CalcularPesosEAbscissas(out pesosEsperados, out abscissasEsperadas);

            Wheeler wheeler = new Wheeler(numeroDePontos,momentosAnaliticos);

            double[] pesosWheeler;
            double[] abscissasWheeler;

            wheeler.CalcularPesosEAbscissas(out pesosWheeler, out abscissasWheeler);
        }
Esempio n. 2
0
        public void CalcularPesosEAbscissasAlgoritmoPDTeste()
        {
            int numeroDePontos = 4;
            double[] momentos = new double[] { 1, 5, 26, 140, 778, 4450, 26140, 157400 };

            double[] pesosEsperados = new[] { 0.0459, 0.4541, 0.4541, 0.0459 };
            double[] abscissasEsperadas = new[] { 2.6656, 4.2580, 5.7420, 7.3344 };

            AlgoritmoPD algoritmo = new AlgoritmoPD(numeroDePontos, momentos);

            double[] pesos;
            double[] abscissas;

            algoritmo.CalcularPesosEAbscissas(out pesos,out abscissas);

            for (int i = 0; i < numeroDePontos; i++)
            {
                Assert.AreEqual(pesosEsperados[i], pesos[i], 0.0001);
                Assert.AreEqual(abscissasEsperadas[i], abscissas[i], 0.0001);
            }
        }
        public void ObterPesosEAbascissasMcCoyCheb()
        {
            const int numeroDePontos = 4;

            double[] momentosAnaliticos = CalcularMomentosIniciaisMcCoy(2 * numeroDePontos, FuncaoMcCoy);

            var chebyChev = new ChebyshevModificado(numeroDePontos, momentosAnaliticos, ChebyshevModificado.TipoDeDominio.SemiInfinito);

            double[] pesos, abscissas;

            chebyChev.CalcularPesosEAbscissas(out pesos, out abscissas);

            AlgoritmoPD algoritmoPd = new AlgoritmoPD(numeroDePontos, momentosAnaliticos);

            double[] pesosEsperados;
            double[] abscissasEsperadas;

            algoritmoPd.CalcularPesosEAbscissas(out pesosEsperados, out abscissasEsperadas);
        }