Changes in version 1.0.0 (2026-03-10)                  

Added

    o   *Breaking:* Added a form of Algorithm 3 in Bernegger (2026) to
	mommb. Instead of a grid-search, the LS algorithm is a nested
	pair of 1-dimensional line searches along g and b.
	Experimentation shows that it almost always converges to the
	“same” values-close to machine precision-as the original
	expectation-maximization (EM) algorithm, and fails to converge
	in similar scenarios as the EM algorithm. Provided for
	completeness and in gratitude to Dr. Bernegger.

    o   *Breaking:* Added new parameters maxb, minb, maxg, and ming.
	The maxb parameter is used in both fitting algorithms. The
	others are only used in the LS algorithm.

    o   Added a fast C implementation of the dilogarithm function for
	use in the LS algorithm.

    o   Added more informative error messages.

    o   Added and updated unit tests.

Changed

    o   *Breaking:* Given the additional algorithm and parameters, the
	call for mommb has been changed. See the documentation for the
	new order of parameters and the options list. This is the first
	“mature” release. Any breaking changes going forward will be
	announced in advance and will have a deprecation period.

    o   Hardened and optimized underlying C and R code for speed and
	error protection.

Fixed

    o   Fixed errors in the documentation including copy-paste errors
	referencing rnorm instead of rmb, etc.

                 Changes in version 0.1.0 (2026-01-08)                  

Added

    o   *Breaking:* The mommb algorithm now calculates E[x^2] as (E[x])
	^ 2 + Var[x] instead of directly. Using the unbiased estimator
	allows for a slightly larger initial estimate. Given that the
	algorithm fails earlier, this allows for a slightly larger
	legal starting location. The results are functionally
	identical, but not always equal within machine tolerances, so
	the change is breaking.

    o   *Breaking:* The mommb algorithm now allows for passing a
	couplet consisting of the first two raw moments through x when
	new parameter m is set to TRUE. When m is FALSE, the default, x
	is treated as a vector of observations regardless of its
	length. This is breaking as the positions of the parameters has
	now changed.

    o   A new parameter trace has been added that when TRUE will cause
	mommb to print the values of g and b at each iteration i.

Changed

    o   *Breaking:* The mommb algorithm now fails immediately when the
	integration would imply a negative p, instead of restarting
	with a different g. Experimentation indicates that the restart
	never helped.

    o   Made tol in mommb a bit less intrusive by passing NULL and then
	checking for it.

    o   Use R API Memzero instead of memset.

    o   Changed GitHub-related files and workflows.

    o   Updated unit tests.

                 Changes in version 0.0.5 (2025-01-31)                  

Changed

    o   Replaced \() shorthand with full function call in mommb to
	remove dependency on R >= 4.1.0.

    o   Updated package metadata test.

    o   Use log1p instead of log 1.0 + ... where appropriate.

                 Changes in version 0.0.4 (2024-07-24)                  

Fixed

    o   Corrected documentation in ecmb.

    o   Fixed one unit test.

Changed

    o   Use log1p instead of log 1.0 + ... where appropriate.

                 Changes in version 0.0.3 (2024-06-20)                  

Changed

    o   Adjust R code for mommb to be more efficient.

    o   Update documentation.

    o   Use native CRAN DOI.

    o   Expand package metadata tests.

                 Changes in version 0.0.2 (2024-05-19)                  

Fixed

    o   Corrected package link in CITATION.

    o   Removed extra parentheses in description of g in docs.

Added

    o   Expanded documentation of ecmb.

    o   Added DOI.

                 Changes in version 0.0.1 (2024-05-15)                  

    o   Initial release.