C# (CSharp) ScalarValue.Lnの例

コード例 #1

        static ScalarValue DiLog1(ScalarValue e)
        {
            ScalarValue f = new ScalarValue(Math.PI * Math.PI / 6.0);

            if (e == 0.0)
            {
                return(f);
            }

            var L  = e.Ln();
            var ek = ScalarValue.One;

            for (int k = 1; k < 250; k++)
            {
                var f_old = f.Clone();
                ek *= e;
                var df = ek * (L - 1.0 / k) / k;
                f += df;

                if (f == f_old)
                {
                    return(f);
                }
            }

            throw new YAMPNotConvergedException("spence");
        }

コード例 #2

ファイル: LambertFunction.cs プロジェクト: FlorianRappl/YAMP

 // series useful for large x
 static ScalarValue SeriesLarge(ScalarValue x)
 {
     var L1 = x.Ln();
     var L2 = L1.Ln();
     var L2L1 = L2 / L1;
     return L1 - L2 + L2L1 + (L2 - 2.0) * L2L1 / L1 / 2.0;
 }

コード例 #3

ファイル: Gamma.cs プロジェクト: CHMP-HLab/HLab.Erp.Lims.Analysis

        static ScalarValue LogGamma_Stirling(ScalarValue z)
        {
            if (z.Im < 0.0)
            {
                return(LogGamma_Stirling(z.Conjugate()).Conjugate());
            }

            var f      = (z - 0.5) * z.Ln() - z + Math.Log(2.0 * Math.PI) / 2.0;
            var reduce = f.Im / (2.0 * Math.PI);

            reduce = f.Im - (int)(reduce) * 2.0 * Math.PI;
            f      = new ScalarValue(f.Re, reduce);

            var zsqu = z * z;
            var zp   = z.Clone();

            for (var i = 1; i < 10; i++)
            {
                var f_old = f.Clone();
                f += Helpers.BernoulliNumbers[i] / (2 * i) / (2 * i - 1) / zp;

                if (f == f_old)
                {
                    return(f);
                }

                zp = zp * zsqu;
            }

            throw new YAMPNotConvergedException("gamma");
        }

コード例 #4

ファイル: LambertFunction.cs プロジェクト: marcusanth/YAMP

        // series useful for large x
        static ScalarValue SeriesLarge(ScalarValue x)
        {
            var L1   = x.Ln();
            var L2   = L1.Ln();
            var L2L1 = L2 / L1;

            return(L1 - L2 + L2L1 + (L2 - 2.0) * L2L1 / L1 / 2.0);
        }

コード例 #5

ファイル: Gamma.cs プロジェクト: CHMP-HLab/HLab.Erp.Lims.Analysis

        static ScalarValue LanczosLogGamma(ScalarValue z)
        {
            ScalarValue sum = new ScalarValue(Helpers.LanczosD[0], 0.0);

            for (int i = 1; i < Helpers.LanczosD.Length; i++)
            {
                sum += Helpers.LanczosD[i] / (z + i);
            }

            sum = (2.0 / Math.Sqrt(Math.PI)) * (sum / z);
            var zshift = z + 0.5;
            var t      = zshift * (zshift + Helpers.LanczosR).Ln() - z;

            return(t + sum.Ln());
        }

コード例 #6

        static ScalarValue DiLog_Log_Series(ScalarValue z)
        {
            var ln  = z.Ln();
            var ln2 = ln * ln;
            var f   = Math.PI * Math.PI / 6.0 + ln * (1.0 - (-ln).Ln()) - ln2 / 4.0;
            var p   = ln.Clone();

            for (int k = 1; k < 17; k++)
            {
                var f_old = f.Clone();
                p *= ln2 / (2 * k + 1) / (2 * k);
                f += (-Helpers.BernoulliNumbers[k] / (2 * k)) * p;

                if (f == f_old)
                {
                    return(f);
                }
            }

            throw new YAMPNotConvergedException("spence");
        }

コード例 #7

ファイル: LnFunction.cs プロジェクト: FlorianRappl/YAMP

 protected override ScalarValue GetValue(ScalarValue value)
 {
     return value.Ln();
 }

コード例 #8

ファイル: Gamma.cs プロジェクト: FlorianRappl/YAMP

        static ScalarValue LogGamma_Stirling(ScalarValue z)
        {
            if (z.Im < 0.0)
            {
                return LogGamma_Stirling(z.Conjugate()).Conjugate();
            }

            var f = (z - 0.5) * z.Ln() - z + Math.Log(2.0 * Math.PI) / 2.0;
            var reduce = f.Im / (2.0 * Math.PI);
            reduce = f.Im - (int)(reduce) * 2.0 * Math.PI;
            f = new ScalarValue(f.Re, reduce);

            var zsqu = z * z;
            var zp = z.Clone();

            for (var i = 1; i < 10; i++)
            {
                var f_old = f.Clone();
                f += Helpers.BernoulliNumbers[i] / (2 * i) / (2 * i - 1) / zp;

                if (f == f_old)
                {
                    return (f);
                }

                zp = zp * zsqu;
            }

            throw new YAMPNotConvergedException("gamma");
        }

コード例 #9

ファイル: Gamma.cs プロジェクト: FlorianRappl/YAMP

        static ScalarValue LanczosLogGamma(ScalarValue z)
        {
            ScalarValue sum = new ScalarValue(Helpers.LanczosD[0], 0.0);

            for (int i = 1; i < Helpers.LanczosD.Length; i++)
                sum += Helpers.LanczosD[i] / (z + i);

            sum = (2.0 / Math.Sqrt(Math.PI)) * (sum / z);
            var zshift = z + 0.5;
            var t = zshift * (zshift + Helpers.LanczosR).Ln() - z;

            return t + sum.Ln();
        }

コード例 #10

ファイル: SpenceFunction.cs プロジェクト: FlorianRappl/YAMP

        static ScalarValue DiLog1(ScalarValue e)
        {
            var f = new ScalarValue(Math.PI * Math.PI / 6.0);

            if (e == 0.0)
            {
                return f;
            }

            var L = e.Ln();
            var ek = ScalarValue.One;

            for (var k = 1; k < 250; k++)
            {
                var f_old = f.Clone();
                ek *= e;
                var df = ek * (L - 1.0 / k) / k;
                f += df;

                if (f == f_old)
                {
                    return f;
                }
            }

            throw new YAMPNotConvergedException("spence");
        }

コード例 #11

ファイル: SpenceFunction.cs プロジェクト: FlorianRappl/YAMP

        static ScalarValue DiLog_Log_Series(ScalarValue z)
        {
            var ln = z.Ln();
            var ln2 = ln * ln;
            var f = Math.PI * Math.PI / 6.0 + ln * (1.0 - (-ln).Ln()) - ln2 / 4.0;
            var p = ln.Clone();

            for (var k = 1; k < 17; k++)
            {
                var f_old = f.Clone();
                p *= ln2 / (2 * k + 1) / (2 * k);
                f += (-Helpers.BernoulliNumbers[k] / (2 * k)) * p;

                if (f == f_old)
                {
                    return f;
                }
            }

            throw new YAMPNotConvergedException("spence");
        }

コード例 #12

ファイル: LnFunction.cs プロジェクト: CHMP-HLab/HLab.Erp.Lims.Analysis

 protected override ScalarValue GetValue(ScalarValue value)
 {
     return(value.Ln());
 }