probability-polynomial 1.0.1.0 → 1.0.2.0

raw patch · 4 files changed

+138/−5 lines, 4 filesdep ~QuickCheckdep ~deepseq

Dependency ranges changed: QuickCheck, deepseq

Files

• CHANGELOG.md +7/−0
• probability-polynomial.cabal +7/−5
• src/Numeric/Function/Piecewise/Polynomial.hs +43/−0
• test/Numeric/Function/Piecewise/PolynomialSpec.hs +81/−0

@@ -1,5 +1,12 @@ # Revision history for `probability-polynomial` +## 1.0.2.0 - 2026-11-05++### Added++* New module `Numeric.Function.Piecewise.Polynomial` containing:+    * Function `integral` to compute the continuous antiderivative of a piecewise polynomial+ ## 1.0.1.0 — 2025-09-26  ### Added

@@ -5,7 +5,7 @@ -- PVP summary:  +-+------- breaking API changes --               | | +----- non-breaking API additions --               | | | +--- code changes with no API change-version:         1.0.1.0+version:         1.0.2.0 synopsis:        Probability distributions via piecewise polynomials description:   Package for manipulating finite probability distributions.@@ -20,14 +20,14 @@ homepage:        https://github.com/DeltaQ-SD/deltaq license:         BSD-3-Clause license-file:    LICENSE-copyright:       Predictable Network Solutions Ltd., 2020-2025+copyright:       Predictable Network Solutions Ltd., 2020-2026 author:          Peter W. Thompson, Heinrich Apfelmus maintainer:      peter.thompson@pnsol.com extra-doc-files:   CHANGELOG.md   README.md -tested-with:     GHC ==8.10.7 || ==9.6.7 || ==9.10.3+tested-with:     GHC ==8.10.7 || ==9.6.7 || ==9.10.3 || ==9.12.4  common warnings   ghc-options: -Wall@@ -44,12 +44,13 @@   build-depends:     , base                 >=4.14.3.0 && <5     , containers           >=0.6      && <0.9-    , deepseq              >=1.4.4.0  && <1.6+    , deepseq              >=1.4.4.0  && <1.7     , exact-combinatorics  >=0.2      && <0.3    exposed-modules:     Data.Function.Class     Numeric.Function.Piecewise+    Numeric.Function.Piecewise.Polynomial     Numeric.Measure.Discrete     Numeric.Measure.Finite.Mixed     Numeric.Measure.Probability@@ -67,10 +68,11 @@     , containers     , hspec                   >=2.11.0 && <2.12     , probability-polynomial-    , QuickCheck              >=2.14   && <2.17+    , QuickCheck              >=2.14   && <2.19    main-is:            Spec.hs   other-modules:+    Numeric.Function.Piecewise.PolynomialSpec     Numeric.Function.PiecewiseSpec     Numeric.Measure.DiscreteSpec     Numeric.Measure.Finite.MixedSpec

@@ -0,0 +1,43 @@+{-|+Copyright   : Predictable Network Solutions Ltd., 2020-2026+License     : BSD-3-Clause+Description : Functions specific to piecewise polynomials.+-}+module Numeric.Function.Piecewise.Polynomial+    ( -- * Operations+      integral+    ) where++import Data.List+    ( mapAccumL+    )+import Numeric.Function.Piecewise+    ( Piecewise+    , fromAscPieces+    , toAscPieces+    )+import Numeric.Polynomial.Simple+    ( Poly+    )++import qualified Numeric.Polynomial.Simple as Poly++{------------------------------------------------------------------------------+    Operations+    Integration+------------------------------------------------------------------------------}+{- | The continuous antiderivative of a piecewise polynomial,+     taking the value 0 at the start of the first piece.+-}+integral+    :: (Ord a, Fractional a)+    => Piecewise (Poly a) -> Piecewise (Poly a)+integral = fromAscPieces . snd . mapAccumL f Poly.zero . toAscPieces+  where+    -- The accumulator is the antiderivative on the previous piece.+    -- Shift the antiderivative on this piece so that it starts+    -- where the previous one ends.+    f before (x, p) = (q, (x, q))+      where+        pp = Poly.integrate p+        q = pp + Poly.constant (Poly.eval before x - Poly.eval pp x)

@@ -0,0 +1,81 @@+{-# LANGUAGE ScopedTypeVariables #-}++{-|+Copyright   : Predictable Network Solutions Ltd., 2020-2026+License     : BSD-3-Clause+-}+module Numeric.Function.Piecewise.PolynomialSpec+    ( spec+    , genPiecewisePoly+    ) where++import Prelude++import Data.Function.Class+    ( eval+    )+import Numeric.Function.Piecewise+    ( Piecewise+    , mapPieces+    , toAscPieces+    )+import Numeric.Function.Piecewise.Polynomial+    ( integral+    )+import Numeric.Function.PiecewiseSpec+    ( genPiecewise+    )+import Numeric.Polynomial.Simple+    ( Poly+    , differentiate+    , scale+    )+import Numeric.Polynomial.SimpleSpec+    ( genPoly+    )+import Test.Hspec+    ( Spec+    , describe+    , it+    )+import Test.QuickCheck+    ( Gen+    , (===)+    , property+    )++{-----------------------------------------------------------------------------+    Tests+------------------------------------------------------------------------------}+spec :: Spec+spec = do+    describe "integral" $ do+        it "is continuous" $ property $+            \(p :: Piecewise (Poly Rational)) ->+                let pieces = toAscPieces (integral p)+                    befores = 0 : map snd pieces+                in  zipWith (\b (x, _) -> eval b x) befores pieces+                        === map (\(x, q) -> eval q x) pieces++        it "is an antiderivative" $ property $+            \(p :: Piecewise (Poly Rational)) x ->+                eval (mapPieces differentiate (integral p)) x+                    === eval p x++        it "is additive" $ property $+            \p (q :: Piecewise (Poly Rational)) x ->+                eval (integral (p + q)) x+                    === eval (integral p + integral q) x++        it "is homogeneous" $ property $+            \(a :: Rational) (p :: Piecewise (Poly Rational)) x ->+                eval (integral (mapPieces (scale a) p)) x+                    === a * eval (integral p) x++{-----------------------------------------------------------------------------+    Random generators+------------------------------------------------------------------------------}+-- | Generate an arbitrary piecewise polynomial,+-- with sorted breakpoints and an arbitrary polynomial on each piece.+genPiecewisePoly :: Gen (Piecewise (Poly Rational))+genPiecewisePoly = genPiecewise genPoly