From 8b449c76d84780357df67c1d41f644ddd6b206f1 Mon Sep 17 00:00:00 2001 From: diluculo Date: Wed, 25 Sep 2019 10:55:44 +0900 Subject: [PATCH] Clean up codes. --- .paket/Paket.Restore.targets | 2 +- .../GaussRule/GaussKronrodPoint.cs | 20 ++++--------------- 2 files changed, 5 insertions(+), 17 deletions(-) diff --git a/.paket/Paket.Restore.targets b/.paket/Paket.Restore.targets index 952ad42f..ca9842d8 100644 --- a/.paket/Paket.Restore.targets +++ b/.paket/Paket.Restore.targets @@ -203,7 +203,7 @@ $([System.String]::Copy('%(PaketReferencesFileLines.Identity)').Split(',')[0]) $([System.String]::Copy('%(PaketReferencesFileLines.Identity)').Split(',')[1]) $([System.String]::Copy('%(PaketReferencesFileLines.Identity)').Split(',')[4]) - $([System.String]::Copy('%(PaketReferencesFileLines.Identity)').Split(',')[5]) + $([System.String]::Copy('%(PaketReferencesFileLines.Identity)').Split(',')[5]) %(PaketReferencesFileLinesInfo.PackageVersion) diff --git a/src/Numerics/Integration/GaussRule/GaussKronrodPoint.cs b/src/Numerics/Integration/GaussRule/GaussKronrodPoint.cs index 49289e0a..8b682b86 100644 --- a/src/Numerics/Integration/GaussRule/GaussKronrodPoint.cs +++ b/src/Numerics/Integration/GaussRule/GaussKronrodPoint.cs @@ -414,7 +414,7 @@ namespace MathNet.Numerics.Integration.GaussRule kronrodAbscissas[(k - 1) / 2] = x0; } - + // Concatenate two abscissas var abscissas = new double[gaussAbscissas.Length + kronrodAbscissas.Length]; @@ -504,26 +504,14 @@ namespace MathNet.Numerics.Integration.GaussRule a[r] = 1.0; for (int k = 1; k < r; k++) { - double ratio = 1d; - BigInteger num = 1; - BigInteger den = 1; - double rat = 1d; - a[r - k] = 0d; + double ratio = 1.0; + a[r - k] = 0.0; for (int i = r + 1 - k; i <= r; i++) { - num = num * (n - q + 2 * (i + k - 1)) * (n + q + 2 * (k - i + 1)) * (n - 1 - q + 2 * (i - k)) * (2 * (k + i - 1) - 1 - q - n); - den = den * (n - q + 2 * (i - k)) * (2 * (k + i - 1) - q - n) * (n + 1 + q + 2 * (k - i)) * (n - 1 - q + 2 * (i + k)); - - var gcd = Euclid.GreatestCommonDivisor(num, den); - num = num / gcd; - den = den / gcd; - - rat = (double)num / (double)den; double numerator = (n - q + 2 * (i + k - 1)) * (n + q + 2 * (k - i + 1)) * (n - 1 - q + 2 * (i - k)) * (2 * (k + i - 1) - 1 - q - n); double denominator = (n - q + 2 * (i - k)) * (2 * (k + i - 1) - q - n) * (n + 1 + q + 2 * (k - i)) * (n - 1 - q + 2 * (i + k)); ratio = ratio * numerator / denominator; - - a[r - k] -= a[i] * rat; // ratio; + a[r - k] -= a[i] * ratio; } }