Please open notebook rsepython-s4r8.ipynb
We’ve seen that using object orientation can produce quite complex class structures, with classes owning each other, instantiating each other, and inheriting from each other.
There are lots of different ways to design things, and decisions to make.
Programmers have noticed that there are certain ways of arranging classes that work better than others.
These are called “design patterns”.
They were first collected on one of the world’s first Wikis, as the Portland Pattern Repository.
A description of a pattern in a book such as the Gang Of Four book (UCL Library) usually includes:
There are lots and lots of design patterns, and it’s a great literature to get into to read about design questions in programming and learn from other people’s experience.
We’ll just show a few in this reading:
matplotlib inline
from unittest.mock import Mock
import requests
from IPython.display import Image, HTML
def yuml(model):
result=requests.get("http://yuml.me/diagram/boring/class/" + model)
return Image(result.content)
Here’s what the Gang of Four Book says about Factory Method:
Intent: Define an interface for creating an object, but let subclasses decide which class to instantiate.
Factory Method lets a class defer instantiation to subclasses.
Applicability: Use the Factory method pattern when:
This is pretty hard to understand, so let’s look at an example.
yuml("[Product]^-[ConcreteProduct], "+
"[Creator| (v) FactoryMethod()]^-[ConcreteCreator| FactoryMethod()], "+
"[ConcreteCreator]-.->[ConcreteProduct]")
An “agent based model” is one like the Boids model from last week: agents act and interact under certain rules. Complex phenomena can be described by simple agent behaviours.
class AgentModel:
def simulate(self):
for agent in agents:
for target in agents:
agent.interact(target)
agent.simulate()
This logic is common to many kinds of Agent based model (ABM), so we can imagine a common class for agent based models: the constructor could parse a configuration specifying how many agents of each type to create, their initial conditions and so on.
However, this common constructor doesn’t know what kind of agent to create; as a common base, it could be a model of boids, or the agents could be remote agents on foreign servers, or they could even be physical hardware robots connected to the driving model over Wifi!
We need to defer the construction of the agents. We can do this with polymorphism: each derived class of the ABM can have an appropriate method to create its agents:
class AgentModel:
def __init__(self, config):
self.agents=[]
for agent_config in config:
self.agents.append(self.create(**agent_config))
This is the factory method pattern: a common design solution to the need to defer the construction of daughter objects to a derived class. self.create
is not defined here, but in each of the agents that inherits from AgentModel
. Using polimorphism to get deffered behaviour on what you want to create.
The type that is created is different in the different derived classes:
class BirdModel(AgentModel):
def create(self, agent_config):
return Boid(agent_config)
Agents are the base product, boids or robots are a ConcreteProduct.
class WebAgentFactory(AgentModel):
def __init__(self, url):
self.url=url
self.connection=AmazonCompute.connect(url)
AgentModel.__init__(self)
def create(self, agent_config):
return OnlineAgent(agent_config, self.connection)
There is no need to define an explicit base interface for the “Agent” concept in Python: anything that responds to “simulate” and “interact” methods will do: this is our Agent concept.
I personally have got into a terrible tangle trying to make base classes which somehow “promote” themselves into a derived class based on some code in the base class.
This is an example of an “Antipattern”: like a Smell, this is a recognised Wrong Way of doing things.
What I should have written was a Creator with a Factory Method.
Consider the following code:
class AgentModel:
def simulate(self):
for agent in agents:
for target in agents:
agent.interact(target)
agent.simulate()
class BirdModel(AgentModel):
def __init__(self, config):
self.boids = []
for boid_config in config:
self.boids.append(Boid(**boid_config))
class WebAgentFactory(AgentModel):
def __init__(self, url, config):
self.url = url
connection = AmazonCompute.connect(url)
AgentModel.__init__(self)
self.web_agents = []
for agent_config in config:
self.web_agents.append(OnlineAgent(agent_config, connection))
The agent creation loop is almost identical in the two classes; so we can be sure we need to refactor it away; but the type that is created is different in the two cases, so this is the smell that we need a factory pattern.
from unittest.mock import Mock
Intent: Separate the steps for constructing a complex object from its final representation.
yuml("[Director|Construct()]<>->[Builder| (a) BuildPart()],"+
" [Builder]^-[ConcreteBuilder| BuildPart();GetResult() ],"+
"[ConcreteBuilder]-.->[Product]")
Let’s continue our Agent Based modelling example.
There’s a lot more to defining a model than just adding agents of different kinds: we need to define boundary conditions, specify wind speed or light conditions.
We could define all of this for an imagined advanced Model with a very very long constructor, with lots of optional arguments:
class Model:
def __init__(self, xsize, ysize,
agent_count, wind_speed,
agent_sight_range, eagle_start_location):
pass
However, long constructors easily become very complicated. Instead, it can be cleaner to define a Builder for models.
A builder is like a deferred factory: each step of the construction process is implemented as an individual method call, and the completed object is returned when the model is ready.
Model=Mock() # Create a temporary mock so the example works!
class ModelBuilder:
def start_model(self):
self.model=Model()
self.model.xlim = None
self.model.ylim = None
def set_bounds(self, xlim, ylim):
self.model.xlim=xlim
self.model.ylim=ylim
def add_agent(self, xpos, ypos):
pass # Implementation here
def finish(self):
self.validate()
return self.model
def validate(self):
assert(self.model.xlim is not None)
# Check that the all the
# parameters that need to be set
# have indeed been set.
Inheritance of an Abstract Builder for multiple concrete builders could be used where there might be multiple ways to build models with the same set of calls to the builder: for example a version of the model builder yielding models which can be executed in parallel on a remote cluster.
builder=ModelBuilder()
builder.start_model()
builder.set_bounds(500, 500)
builder.add_agent(40, 40)
builder.add_agent(400, 100)
model=builder.finish()
model.simulate()
<Mock name=’mock().simulate()’ id=’4497379624>
We could, of course, just add all the building methods to the model itself, rather than having the model be yielded from a separate builder.
This is an anti-pattern that is often seen: a class whose __init__
constructor alone is insufficient for it to be ready to use. A series of
methods must be called, in the right order, in order for it to be ready to use.
This results in very fragile code: its hard to keep track of whether an object instance is “ready” or not. Use the builder pattern to keep deferred construction in control.
We might ask why we couldn’t just use a validator in all of the methods that must follow the deferred constructors; to check they have been called.
But we’d need to put these in every method of the class, whereas with a builder, we can validate only in the finish
method.
%matplotlib inline
Define a family of algorithms, encapsulate each one, and make them interchangeable.
Strategy lets the algorithm vary independently from clients that use it.
from numpy import linspace,exp,log,sqrt, array
import math
from scipy.interpolate import UnivariateSpline
from scipy.signal import lombscargle
from scipy.integrate import cumtrapz
from numpy.fft import rfft,fft,fftfreq
import csv
from io import StringIO
from datetime import datetime
import requests
import matplotlib.pyplot as plt
%matplotlib inline
Consider the sequence of sunspot observations:
def load_sunspots():
url_base="http://www.quandl.com/api/v1/datasets/SIDC/SUNSPOTS_A.csv"
x=requests.get(url_base,params={'trim_start':'1700-12-31',
'trim_end':'2017-12-01',
'sort_order':'asc'})
data=csv.reader(StringIO(x.text)) #Convert requests result to look
#like a file buffer before
# reading with CSV
next(data) # Skip header row
return [float(row[1]) for row in data]
spots=load_sunspots()
plt.plot(spots)
[<matplotlib.lines.Line2D at 0x10f8606d8>]
spectrum=rfft(spots)
plt.figure()
plt.plot(abs(spectrum))
plt.savefig('fixed.png')
There’s a potential problem with this analysis however:
The Fast Fourier Transform cannot be applied to uneven time series.
We could:
We also want to find the period of the strongest periodic signal in the data, there are various different methods we could use for this also, such as integrating the fourier series by quadrature to find the mean frequency, or choosing the largest single value.
We could implement a base class for our common code between the different approaches, and define derived classes for each different algorithmic approach. However, this has drawbacks:
class SunspotAnalyzerSplineFFTTrapeziumNearMode
is a bit unweildy.#### Apply the strategy pattern:
First, we’ll define a helper class for our time series.
class Series:
"""Enhance NumPy N-d array with some helper functions for clarity"""
def __init__(self, data):
self.data=array(data)
self.count=self.data.shape[0]
self.start=self.data[0,0]
self.end=self.data[-1,0]
self.range=self.end-self.start
self.step=self.range/self.count
self.times=self.data[:,0]
self.values=self.data[:,1]
self.plot_data=[self.times,self.values]
self.inverse_plot_data=[1.0/self.times[20:], self.values[20:]]
Then, our class which contains the analysis code, except the numerical methods
class AnalyseSunspotData:
def format_date(self, date):
date_format="%Y-%m-%d"
return datetime.strptime(date,date_format)
def load_data(self):
start_date_str='1700-12-31'
end_date_str='2014-01-01'
self.start_date=self.format_date(start_date_str)
end_date=self.format_date(end_date_str)
url_base=("http://www.quandl.com/api/v1/datasets/"+
"SIDC/SUNSPOTS_A.csv")
x=requests.get(url_base,params={'trim_start':start_date_str,
'trim_end':end_date_str,
'sort_order':'asc'})
secs_per_year=(datetime(2014,1,1)-datetime(2013,1,1)
).total_seconds()
data=csv.reader(StringIO(x.text)) #Convert requests
#result to look
#like a file buffer before
#reading with CSV
next(data) # Skip header row
self.series=Series([[
(self.format_date(row[0])-self.start_date
).total_seconds()/secs_per_year
,float(row[1])] for row in data])
def __init__(self, frequency_strategy):
self.load_data()
self.frequency_strategy=frequency_strategy
def frequency_data(self):
return self.frequency_strategy.transform(self.series)
Our existing simple fourier strategy
class FourierNearestFrequencyStrategy:
def transform(self, series):
transformed=fft(series.values)[0:series.count/2]
frequencies=fftfreq(series.count, series.step)[0:series.count/2]
return Series(list(zip(frequencies, abs(transformed)/series.count)))
A strategy based on interpolation to a spline:
class FourierSplineFrequencyStrategy:
def next_power_of_two(self, value):
"Return the next power of 2 above value"
return 2**(1+int(log(value)/log(2)))
def transform(self, series):
spline=UnivariateSpline(series.times, series.values)
# Linspace will give us *evenly* spaced points in the series
fft_count= self.next_power_of_two(series.count)
points=linspace(series.start,series.end,fft_count)
regular_xs=[spline(point) for point in points]
transformed=fft(regular_xs)[0:fft_count//2]
frequencies=fftfreq(fft_count,
series.range/fft_count)[0:fft_count//2]
return Series(list(zip(frequencies, abs(transformed)/fft_count)))
A strategy using the Lomb-Scargle Periodogram:
class LombFrequencyStrategy:
def transform(self,series):
frequencies=array(linspace(1.0/series.range,
0.5/series.step,series.count))
result= lombscargle(series.times,
series.values,2.0*math.pi*frequencies)
return Series(list(zip(frequencies, sqrt(result/series.count))))
Define our concrete solutions with particular strategies:
fourier_model=AnalyseSunspotData(FourierSplineFrequencyStrategy())
lomb_model=AnalyseSunspotData(LombFrequencyStrategy())
nearest_model=AnalyseSunspotData(FourierNearestFrequencyStrategy())
Use these new tools to compare solutions
comparison=fourier_model.frequency_data().inverse_plot_data+['r']
comparison+=lomb_model.frequency_data().inverse_plot_data+['g']
comparison+=nearest_model.frequency_data().inverse_plot_data+['b']
deviation=365*(fourier_model.series.times-linspace(
fourier_model.series.start,
fourier_model.series.end,
fourier_model.series.count))
plt.plot(*comparison)
plt.xlim(0,16)
(0, 16)
plt.plot(deviation)
[<matplotlib.lines.Line2D at 0x10f9ad7f0>]
Whenever we are coding a simulation or model we want to:
We often see scientific programs where the code which is used to display what is happening is mixed up with the mathematics of the analysis. This is hard to understand.
We can do better by separating the Model
from the View
, and using a “Controller
” to manage them.
This is where we describe our internal logic, rules, etc.
import numpy as np
class Model:
def __init__(self):
self.positions=np.random.rand(100,2)
self.speeds=(np.random.rand(100,2) +
np.array([-0.5,-0.5])[np.newaxis,:])
self.deltat=0.01
def simulation_step(self):
self.positions += self.speeds * self.deltat
def agent_locations(self):
return self.positions
This is where we describe what the user sees of our Model, what’s displayed. You may have different type of visualisation (e.g., on one type of projection, a 3D view, a surface view, …) which can be implemented in different view classes.
class View:
def __init__(self, model):
from matplotlib import pyplot as plt
self.figure=plt.figure()
axes=plt.axes()
self.model=model
self.scatter=axes.scatter(
model.agent_locations()[:,0],
model.agent_locations()[:,1])
def update(self):
self.scatter.set_offsets(
self.model.agent_locations())
This is the class that tells the view that the models has changed and updates the model with any change the user has input through the view.
class Controller:
def __init__(self):
self.model = Model() # Or use Builder
self.view = View(self.model)
def animate(frame_number):
self.model.simulation_step()
self.view.update()
self.animator = animate
def go(self):
from matplotlib import animation
anim = animation.FuncAnimation(self.view.figure,
self.animator,
frames=200,
interval=50)
return anim.to_jshtml()
contl=Controller()
HTML(contl.go())
Other resources