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

Newer ↓
Demo
Mar 6 2018
Using three.js for 2D data visualization
by — A code walkthrough of how to use three.js for large 2D data visualizations. Three.js can smoothly render a large number of points in the browser but there are some tricks to making it work for interactive 2D visualizations. This is the guide that I wish I would have had when I built our most recent prototype. It shows how to enable mousedirected pan and zoom, circular points, and hover interac...
…read more
by Grant

Newest ↓
Post
Apr 19 2018
Introducing SciFi
by — Today we are launching a minisite featuring our collection of short stories inspired by new developments in machine learning. Beginning with our fourth report, we started including a sciencefiction story along with the technical and strategic overviews that are the bulk of each report. Using these stories, we can look at the technologies we profile from a different angle and explore their c...
…read more
by Grant