Fast Forward Food Labs
Probabilistic Cookies!
Feb 14 2018 · by Ryan Micallef
In the spirit of Valentine’s Day, we at Fast Forward Labs thought it would be fun to bake cookies for our sweethearts. Being DIY nerds, we thought we’d math it up a bit.
We used python to generate probability distributions and matplotlib to check our distributions. Then we wrote a python function to generate a SCAD file defining threedimensional shapes from the distributions. Using OpenSCAD, an opensource CAD program, we checked the 3D models and exported them to STL files for printing. We used a 5th Generation MakerBot Replicator to print our 3D models. And we baked cookies. Here’s one of our office dogs (and my best friend), Dogface, admiring the results.
There were a number of challenges involved in generating 3D models and printing them. Here’s how the basic process went.
We chose a beta distribution for our first prototype because it’s well behaved for purposes of making a 3D printed object.
 A beta distribution only has values from 0 to 1. This gives us a fixedwidth shape to work with. (Compare a Gaussian distribution, which has long tails on both sides and thus may not normalize to a good shape across 0 to 1.)
 The area under the curve of a beta distribution is necessarily 1, which helps keep the shape from getting too eccentric while allowing flexibility in choice of parameters, and thus a wider range of shapes.
 If you take two beta distributions, put them xaxis to xaxis, and squint, they look a bit like a heart.
We created a beta distribution in python.
import numpy as np
from scipy.stats import beta
import matplotlib.pyplot as plt
# choose enough points to have a relatively smooth curve without
# creating so many facets that the 3D printer is slow
#
# choose an odd number so there is a definitive peak to the curve
numpoints = 35
# set up linspace
X = np.linspace(0, 1, numpoints)
# beta distribution parameters
a, b = 2, 1.6
# get beta distribution array
betadist = beta.pdf(X, a, b)
And plotted that distribution with some suitable parameters.
# turn off the axes so we only see the curve
plt.axis('off')
# plot the curve itself
plt.plot(X, betadist)
Once we found a set of parameters that we liked, we used the points in the distribution in python to create an OpenSCADformatted SCAD file. A SCAD file is a readable text file that contains a combination of points in space that define shapes and instructions to manipulate those shapes. We made two 2D copies of the distribution, sized one up a bit, and centered them together. We extruded them both into 3D, one linear, and the other with a slight “cone” projection. This gave us one relatively thin edge for piercing cookie dough. Then we subtracted one shape from the other to make a hollow in the larger shape. Rendered in 3D in OpenSCAD, it looks like this:
Note that the 3D model looks more eccentric than the pythongenerated plot. This is because the plot had axes of differing scales.
With a 3D model in hand, it was time to print a physical object. We exported our SCAD object into a stereolithography (STL) file, which we then imported into MakerBot’s printing software.
We made a final print file and extruded our first prototypes at NYC Resistor our friendly neighborhood hackerspace. 3D printing can take a while for large shapes, so we started out small.
Here are the first results (we had begun tinkering with a Gaussian distribution at that point):
Of course these first test shapes are too small for cookies. After a lot of tinkering and refinement, we ended up with beta, Gaussian, and Poisson distribution shapes scaled up for cookie size (about 100mm high).
We made some cookie dough and got to business.
Note that the Poisson distribution (printed in white) has big solid areas at the top. This makes it more of a cookiedough perturber than a cookie cutter. Those solid areas are an artifact of my own ignorance of OpenSCAD. The “clever” OpenSCAD scale()
approach I had been using was a hack; I later learned that the offset()
function is the correct solution.
OpenSCAD issues aside, the cookies turned out fine. Here are some of the results, some decorated with axes and histograms.
If you feel inclined to give it a shot yourself, here’s a jupyter notebook with some code to get you started.
Happy Valentine’s Day from Fast Forward Labs!
More from the Blog

Older ↓
Post
Jan 22 2018
Exploring Recommendation Systems
by — While we commonly associate recommendation systems with ecommerce, their application extends to any decisionmaking problem which requires pairing two types of things together. To understand why recommenders don’t always work as well as we’d like them to, we set out to build some basic recommendation systems using publicly available data. Data The first ingredient for building a recommendati...
…read more
by Shioulin

Archive ↓
interview
Dec 9 2015
Fashion Goes Deep: Data Science at Lyst
with — On November 16, 2015, Lyst, an online fashion marketplace based in the United Kingdom, launched its first advertising campaign. Featuring a series of ironic headlines (one simply says “Ripoff”) etched over beautiful images, the campaign emphasizes the company’s identity as a “challenger brand,” whose success “has been driven by marrying insights from data science with the emotional nature of ...
…read more
with Eddie Bell

Archive ↓
Interview
Apr 28 2017
Eli Bressert on DataDriven Processes at Netflix
with — Despite their best intentions, companies often struggle to develop processes that provide datadriven decisions through partnership across the business. Netflix stands out as a company that excels at deeply integrating its data science teams into all aspects of the company. We spoke with Eli Bressert, Manager, Data Engineering and Analytics at Netflix, to learn more about how they create the cu...
…read more
with Eli Bressert