module Factory.Math.Implementations.Pi.Ramanujan.Implementation(
openR
) where
import qualified Control.Parallel.Strategies
import qualified Factory.Math.Implementations.Pi.Ramanujan.Series as Math.Implementations.Pi.Ramanujan.Series
import qualified Factory.Math.Precision as Math.Precision
import qualified Factory.Math.Summation as Math.Summation
openR
:: Math.Implementations.Pi.Ramanujan.Series.Series squareRootAlgorithm factorialAlgorithm
-> squareRootAlgorithm
-> factorialAlgorithm
-> Math.Precision.DecimalDigits
-> Rational
openR :: Series squareRootAlgorithm factorialAlgorithm
-> squareRootAlgorithm
-> factorialAlgorithm
-> DecimalDigits
-> Rational
openR Math.Implementations.Pi.Ramanujan.Series.MkSeries {
terms :: forall squareRootAlgorithm factorialAlgorithm.
Series squareRootAlgorithm factorialAlgorithm
-> factorialAlgorithm -> [Rational]
Math.Implementations.Pi.Ramanujan.Series.terms = factorialAlgorithm -> [Rational]
terms,
getSeriesScalingFactor :: forall squareRootAlgorithm factorialAlgorithm.
Series squareRootAlgorithm factorialAlgorithm
-> squareRootAlgorithm -> DecimalDigits -> Rational
Math.Implementations.Pi.Ramanujan.Series.getSeriesScalingFactor = squareRootAlgorithm -> DecimalDigits -> Rational
getSeriesScalingFactor,
convergenceRate :: forall squareRootAlgorithm factorialAlgorithm.
Series squareRootAlgorithm factorialAlgorithm -> ConvergenceRate
Math.Implementations.Pi.Ramanujan.Series.convergenceRate = ConvergenceRate
convergenceRate
} squareRootAlgorithm
squareRootAlgorithm factorialAlgorithm
factorialAlgorithm DecimalDigits
decimalDigits = (Rational -> Rational -> Rational)
-> (Rational, Rational) -> Rational
forall a b c. (a -> b -> c) -> (a, b) -> c
uncurry Rational -> Rational -> Rational
forall a. Fractional a => a -> a -> a
(/) ((Rational, Rational) -> Rational)
-> (Rational, Rational) -> Rational
forall a b. (a -> b) -> a -> b
$ Strategy (Rational, Rational)
-> (Rational, Rational) -> (Rational, Rational)
forall a. Strategy a -> a -> a
Control.Parallel.Strategies.withStrategy (
Strategy Rational
-> Strategy Rational -> Strategy (Rational, Rational)
forall a b. Strategy a -> Strategy b -> Strategy (a, b)
Control.Parallel.Strategies.parTuple2 Strategy Rational
forall a. NFData a => Strategy a
Control.Parallel.Strategies.rdeepseq Strategy Rational
forall a. NFData a => Strategy a
Control.Parallel.Strategies.rdeepseq
) (
squareRootAlgorithm -> DecimalDigits -> Rational
getSeriesScalingFactor squareRootAlgorithm
squareRootAlgorithm DecimalDigits
decimalDigits,
DecimalDigits -> [Rational] -> Rational
forall i.
(Integral i, NFData i) =>
DecimalDigits -> [Ratio i] -> Ratio i
Math.Summation.sumR DecimalDigits
64 ([Rational] -> Rational)
-> ([Rational] -> [Rational]) -> [Rational] -> Rational
forall b c a. (b -> c) -> (a -> b) -> a -> c
. DecimalDigits -> [Rational] -> [Rational]
forall a. DecimalDigits -> [a] -> [a]
take (
ConvergenceRate -> DecimalDigits -> DecimalDigits
forall i. Integral i => ConvergenceRate -> DecimalDigits -> i
Math.Precision.getTermsRequired ConvergenceRate
convergenceRate DecimalDigits
decimalDigits
) ([Rational] -> Rational) -> [Rational] -> Rational
forall a b. (a -> b) -> a -> b
$ factorialAlgorithm -> [Rational]
terms factorialAlgorithm
factorialAlgorithm
)