public static void CalcularPesosEAbscissas(MetodoCoeficientesDosMomentos metodoCoeficientes, double[] momentos,
out double[] pesos, out double[] abscissas)
{
int numeroDePontos = momentos.Length / 2;
switch (metodoCoeficientes)
{
case MetodoCoeficientesDosMomentos.PDA:
AlgoritmoPD algoritmoPd = new AlgoritmoPD(numeroDePontos, momentos);
algoritmoPd.CalcularPesosEAbscissas(out pesos, out abscissas);
break;
case MetodoCoeficientesDosMomentos.ChebyShev:
var chebyChev = new ChebyshevModificado(numeroDePontos,momentos,ChebyshevModificado.TipoDeDominio.Finito);
//var chebyChev = new ChebyshevModificado(momentos, numeroDePontos, ChebyshevModificado.TipoDeDominio.Finito);
chebyChev.CalcularPesosEAbscissas(out pesos, out abscissas);
break;
case MetodoCoeficientesDosMomentos.Wheeler:
Wheeler wheeler = new Wheeler(numeroDePontos, momentos);
wheeler.CalcularPesosEAbscissas(out pesos, out abscissas);
break;
default:
throw new ArgumentOutOfRangeException("metodoCoeficientes", metodoCoeficientes, null);
}
}
private void CalcularPesosEAbscissas()
{
double[] pf;
double[] abf;
var momentos = CalcularMomentosFinais();
Wheeler wheeler = new Wheeler(momentos.Length / 2, momentos);
wheeler.CalcularPesosEAbscissas(out pf, out abf);
double[] pi;
double[] abi;
double[] momentosIniciais = new[]
{
1.0000000000000000,
0.7343097109235880,
0.5502539727012860,
0.4201449510879860,
0.3264374407964930,
0.2577688630354670,
0.2066322315214840,
0.1679756057694100,
0.1383405829521260,
0.1153216484404730,
0.0972205496278511,
0.0828215151094570,
0.0712427053220418,
0.0618365584375131,
0.0541219897678821,
0.0477376477314290
};
wheeler = new Wheeler(momentosIniciais.Length / 2, momentosIniciais);
wheeler.CalcularPesosEAbscissas(out pi, out abi);
pesos = new double[pi.Length];
abscissas = new double[abi.Length];
pesos = pf;
abscissas = abf;
//for (int i = 0; i < pi.Length; i++)
//{
// pesos[i] = (pf[i] + pi[i])/2;
// abscissas[i] = (abf[i] + abi[i]) / 2;
//}
//Até aqui OK
}