2

C# で chudnovsky メソッドを使用して Pi を計算するための現在のコードを次に示します。

using System;
using System.Diagnostics;
using System.IO;
using java.math;

namespace pi.chudnovsky
{
    public class Program
    {
        static Double Factorial(Double fact)
        {
            //begin factorial function
            if (fact <= 1)
                return 1;
            else
                return fact * Factorial(fact - 1); //loops multiplication until the factorial is reached
        }
        static Double doSummation(Double maxPower)
        {
            //begin chudnovsky summation function
            Double sum = 0;
            for (int i = 0; i <= maxPower; i++) //starts at i=0
            {
                sum += ((Math.Pow(-1, i)) * Factorial(6 * i) * (13591409 + 5451401 * i)) / (Factorial(3 * i) * Factorial(i) * Factorial(i) * Factorial(i) * Math.Pow(640320, (3 * i + 1.5))); //chudnovsky algorithm
            }
            return sum;
        }
        static void Main(string[] args)
        {
            int num;
            Console.WriteLine("Enter how many terms to compute Chudnovsky summation: ");
            //begin stopwatch
            Stopwatch stopwatch = new Stopwatch();
            stopwatch.Start();
            //parse user input
            num = Convert.ToInt32(Console.ReadLine());
            //perform calculation
            Double inv = 1 / (12 * doSummation(num));
            //stop stopwatch
            stopwatch.Stop();
            //display info
            Console.WriteLine(inv);
            Console.WriteLine("3.14159265358979323846264338327950288419716939937510");
            Console.WriteLine("Time elapsed: {0}", stopwatch.Elapsed.TotalMilliseconds);
            //write to pi.txt
            TextWriter pi = new StreamWriter("pi.txt");
            pi.WriteLine(inv);
            pi.Close();
            //write to stats.txt
            TextWriter stats = new StreamWriter("stats.txt");
            stats.WriteLine(stopwatch.Elapsed.TotalMilliseconds);
            stats.Close();
        }
    }
}

そのため、J# ライブラリを含め、java.math を含めました。すべての「double」を「BigDecimal」に置き換えると、次のコンパイル エラーが発生します。

http://f.cl.ly/items/1r2X26470d0d0n260p0p/Image%202011-11-14%20at%206.16.19%20PM.png

Doubles で完全に機能したため、これがループに Int を使用することの問題ではないことはわかっています。私の質問は、int と BigDecimal に関連するこれらのエラーをどのように解決しますか、または別の任意精度ライブラリを推奨できますか?

XMPIR を使用してみましたが、コンパイルするすべてのものを取得しましたが、次のようになります。

http://f.cl.ly/items/1l3C371j2u3z3n2g3a0j/Image%202011-11-14%20at%206.20.24%20PM.png

p/invoke を使用して xmpir をインクルードし、bigdecimal クラスが何であれ使用できるようにできますか?

お時間をいただきありがとうございます！

4

1 に答える 1

1

それらを比較する前に、int を BigDecimal に変換できませんか?

ここでの問題は、int を受け入れる BigDecimal クラスに、より大きい、より小さいなどの記号に対する演算子のオーバーロードがないことを理解していると思います。

于 2011-11-15T01:35:51.177 に答える