Combinator Constructors

Combinator constructors use transformations or measurements to produce new transformations or measurements. Combinators are an area of OpenDP that are still heavily in development, but the chainers in particular are foundational and well-tested.

Chaining

Two of the most essential constructors are the “chainers” that chain transformations with transformations, and transformations with measurements. Chainers are used to incrementally piece Transformations or Measurements together that represent longer computational pipelines.

The opendp.comb.make_chain_tt() constructor creates a new Transformation by combining an inner and an outer Transformation. The resulting Transformation contains a function that sequentially executes the function of the constituent Transformations. It also contains a privacy relation that relates an input distance bound on the inner Transformation with an output distance bound on the outer transformation.

The opendp.comb.make_chain_mt() constructor similarly creates a new Measurement by combining an inner Transformation with an outer Measurement. Notice that there is no make_chain_mm for chaining measurements together! Any computation beyond a measurement is postprocessing and need not be governed by relations.

In the following example we chain opendp.meas.make_base_geometric() with opendp.trans.make_bounded_sum().

>>> from opendp.trans import make_bounded_sum
>>> from opendp.meas import make_base_geometric
>>> from opendp.comb import make_chain_mt
...
>>> # call a constructor to produce a transformation
>>> bounded_sum = make_bounded_sum(bounds=(0, 1))
>>> # call a constructor to produce a measurement
>>> base_geometric = make_base_geometric(scale=1.0)
>>> noisy_sum = make_chain_mt(base_geometric, bounded_sum)
...
>>> # investigate the privacy relation
>>> symmetric_distance = 1
>>> epsilon = 1.0
>>> assert noisy_sum.check(d_in=symmetric_distance, d_out=epsilon)
...
>>> # invoke the chained measurement's function
>>> dataset = [0, 0, 1, 1, 0, 1, 1, 1]
>>> release = noisy_sum(dataset)

In practice, these chainers are used so frequently that we’ve written a shorthand (>>). The syntax automatically chooses between make_chain_mt and make_chain_tt.

>>> noisy_sum = bounded_sum >> base_geometric

In this example the chaining was successful because:

  • bounded_sum’s output domain is equivalent to base_geometric’s input domain

  • bounded_sum’s output metric is equivalent to base_geometric’s input metric

Chaining fails if we were to adjust the domains such that they won’t match. In the below example, the adjustment is subtle, but the bounds were adjusted to floats. make_bounded_sum is equally capable of summing floats, but the chaining fails because the sum emits floats and the geometric mechanism expects integers.

>>> from opendp.mod import OpenDPException
>>> try:
...     make_bounded_sum(bounds=(0., 1.)) >> base_geometric
... except OpenDPException as err:
...     print(err.message[:-1])
Intermediate domains don't match. See https://github.com/opendp/opendp/discussions/297
    The structure of the intermediate domains are the same, but the types or parameters differ.
    shared_domain: AllDomain()

Note that noisy_sum’s input domain and input metric come from bounded_sum’s input domain and input metric. This is intended to enable further chaining with preprocessors like make_cast, make_impute_constant, make_clamp and make_bounded_resize. See the section on Transformation Constructors for more information on how to preprocess data in OpenDP.

Composition

OpenDP has a simple composition combinator for building a tuple of measurements out of two measurements. You can nest this composition combinator arbitrarily.

>>> from opendp.comb import make_basic_composition
>>> noisy_sum_pair = make_basic_composition(noisy_sum, noisy_sum)

Progress on more general composition constructors can be found in the following PR’s:

#242

sequential compositor w/ static distances

#264

f-DP

Amplification

Population amplification based on sampling has an implementation but is not yet finalized. Progress on this can be found in the following PR:

#233

population amplification