The Value of Value Semantics
I am pleased to welcome back Dave Foti to examine some of the ways value semantics work in MATLAB objects and how data copies can be minimized.
As I know some of you may have children who are learning remotely at home now, I thought I would share a MATLAB object I built to help my son practice multiplication. While it is possible to build this into an app, I was looking to make paper tests so this example is designed to make a test that can be printed. Some of you are probably familiar with handle objects in MATLAB or perhaps have used other languages where objects are always references. I want to use my example classes to point out some of the reasons handle or reference semantics are often unnecessary in MATLAB.
Contents
Handles and Values
First, I think of a handle as an object that can access some data that is not part of the object itself but reachable from the object. When I pass around this object, I don’t pass around the actual data, but just an ID that provides the means to reach that data. This means when I change the data owned by the object, I may have passed that ID to many functions. Many variables in various places may now be able to access the same data and see those changes. It often makes sense to represent objects in the real world as handles because you can’t just make copies of real-world objects. For example, if I have a Train object that is designed to track the locations of trains, it might look like this:
classdef Train < handle properties Number Location end methods function p = Train(Number, Location) p.Number = Number; p.Location = Location; end end end
No matter how many variables I use to store handles for the same train, I always see the same state.
X = Train(5402, "Boston"); Y = X; % We can copy the handle or ID, but not the train itself Y.Location = "New York"; X.Location
ans = "New York"
What I mean by a value or value semantics, is the notion that when I pass an object to a function the function receives a copy of all the data contained in that object. Similarly, when I assign an object to a new variable, that variable gets a copy of all the data and while two variables may hold identical data, changing data in one variable’s object never changes the data for a different variable. While numbers are the most common forms of values, many kinds of objects can benefit from the same value semantics. Take my arithmetic test for example, it is just a collection of parameters and arithmetic problems.
test = MultTest("NumProblems", 50, "Max", 10)
test = MultTest with properties: NumProblems: 50 Max: 10 X: [50×1 double] Y: [50×1 double] Z: [50×1 double] OpSymbol: '×'
Now if I want to make a new test that is the same except for going from 1 to 12 instead of 1 to 10, I can assign a new variable to test and then modify it without changing test.
test12 = test;
test12.Max = 12;
test12 = drawProblems(test12); % recreate and draw the problems
My original test is not changed:
test
test = MultTest with properties: NumProblems: 50 Max: 10 X: [50×1 double] Y: [50×1 double] Z: [50×1 double] OpSymbol: '×'
Semantic Copies are Efficient in MATLAB
Sometimes people worry that making all these copies of objects will lead to using lots of memory. It is important to realize that MATLAB optimizes copies in many ways to make them efficient. When an object is copied from one variable to another, MATLAB simply records this fact, but no data is actually copied. Both variables are actually referencing the same object, but it knows it is linked to multiple variables. This doesn't matter until one of those variables is used to modify the object. Even then, MATLAB makes a shallow copy of the object which means that any matrices inside the object that are not being changed remain shared. When one property is modified, only that property's data is copied, not the data for properties that weren't modified. If we look at the above example again:
test12 = test; test12.Max = 12;
when test12.Max is changed, no copy of the problems is made because those arrays can still be shared between the objects. When we redraw the problems, only then will new arrays be created for the new problems. With these simple examples, you won't notice the difference, but when objects store very large arrays, efficient copies make a real difference.
Value Semantics Are Clear about Input and Output
One of the features of MATLAB that I've always appreciated is the fact that you can look at a function call and easily tell what the inputs and outputs are, and even if a variable is both an input and an output. Some languages and systems use inputs for outputs and you have to read the function documentation or input argument annotations like "in", "out" and "in/out" to understand how the arguments are used. In MATLAB, inputs are on the right, outputs are on the left and parameters that are in/out appear on both the left and the right. When a variable is used on the left and right side of a function call, MATLAB can often avoid making a copy even if the variable is modified by the function. Let's look at the first few lines of the drawProblems function for ArithmeticTest:
dbtype ArithmeticTest.m 45:50
45 function t = drawProblems(t, withAnswers) 46 arguments 47 t 48 withAnswers (1,1) logical = false; 49 end 50 t = genProblems(t);
On line 45, we can see that t is both an input and output. This means that MATLAB can sometimes modify t in place without making a copy but only if the call to the function uses the same variable for the input and output. We can see an example of such a call to a different method on line 50. Even through these multiple levels of function calls, it is possible that t is never actually copied because it is possible that MATLAB can tell that it never needs more than one copy to preserve value semantics.
Below you will find the full program listings for reference:
I'd like to hear about any examples of new value types you've created in MATLAB. Please post any descriptions or links you would like to share in here.
Appendix: Code Files
classdef ArithmeticTest % ArithmeticTest Abstract base class for arithmetic tests properties % Number of problems in the test NumProblems (1,1) {mustBeReal, mustBeInteger} = 35; % Limits the values of operands Max (1,1) {mustBeReal, mustBeInteger} = 20; end properties (SetAccess = protected) X Y Z OpSymbol end methods(Abstract) t = genProblems(t) end methods function t = ArithmeticTest(params) % Assign all the name value pairs to properties of the % object being constructed. for s = string(fieldnames(params)') t.(s) = params.(s); end end end methods(Access = protected) function op = getSymbol(t, k) if isscalar(t.OpSymbol) op = t.OpSymbol; else op = t.OpSymbol(k); end end end methods % drawProblems Draw the problems in a figure that can be printed % drawProblems(t, true) will print the problems with the answers % filled in. function t = drawProblems(t, withAnswers) arguments t withAnswers (1,1) logical = false; end t = genProblems(t); if (t.NumProblems <= 48) jN = 6; kN = 8; fontSize = 20; axPos = [0 0 7 10]; elseif (t.NumProblems <= 80) jN = 8; kN = 10; fontSize = 12; axPos = [0.55 0 8.5 10.4]; else jN = 10; kN = 10; fontSize = 12; axPos = [0.55 0 8.5 10.4]; end f = figure('Units', 'inches', 'Position', [3 0 8.5 11]); a = axes('Parent', f, 'Units', 'Inches', 'Position', axPos); a.Visible = 'off'; a.XLim = [0 jN+1]; a.YLim = [0 kN]; n = 1; for j = 1:jN for k = 1:kN if (n > t.NumProblems) break; end x = t.X((j-1)*kN+k); y = t.Y((j-1)*kN+k); if (x < 99) && (y < 99) dx = 2; dy = 2; else dx = 3; dy = 3; end tx1 = text(j-.8, k-.3, ... sprintf(" %" + dx + "d\n%s%" + dy + "d", ... x, t.getSymbol(n), y)); tx1.FontName = 'Lucida Console'; tx1.FontSize = fontSize; tx1.Interpreter = 'none'; if (x < 99) && (y<99) underlineTxt = '___'; else underlineTxt = '____'; end tx2 = text(j-.8, k-.44, underlineTxt); tx2.FontName = 'Courier New'; tx2.FontSize = fontSize; tx2.Interpreter = 'none'; if (withAnswers) tx3 = text(j-.8, k-.7, sprintf(" %"+dx+"d", ... t.Z((j-1)*jN+k))); tx3.FontName = 'Courier New'; tx3.FontSize = fontSize; tx3.Interpreter = 'none'; end n = n + 1; end if (n > t.NumProblems) break; end end end end end classdef MultTest < ArithmeticTest % MultTest Test for multiplication problems methods function t = MultTest(params) arguments params.?MultTest end t@ArithmeticTest(params); t.OpSymbol = char(215); % Unicode code point 00D7 t = drawProblems(t); end function t = genProblems(t) rng('shuffle'); % Start with all combinations of multiplying % numbers from 1 ... t.Max. A1 = 1:t.Max; A2 = 1:t.Max; X = zeros(t.Max^2, 1); Y = zeros(t.Max^2, 1); for k = 1:t.Max for j = 1:t.Max X((k-1)*t.Max+j,1) = A1(k); Y((k-1)*t.Max+j,1) = A2(j); end end % Reduce the number of easy questions involving 10 or 1. idx = find(X == 10); X(idx(1:5)) = []; Y(idx(1:5)) = []; idx = find(Y == 10); X(idx(1:5)) = []; Y(idx(1:5)) = []; idx = find(X == 1); X(idx(1:5)) = []; Y(idx(1:5)) = []; idx = find(Y == 1); X(idx(1:5)) = []; Y(idx(1:5)) = []; % Randomly permute the problems and make sure we have % enough problems. numProbs = numel(X); numSets = ceil(t.NumProblems/numProbs); t.X = zeros(t.NumProblems, 1); t.Y = zeros(t.NumProblems, 1); if numSets > 1 for k = 1:numSets-1 idx = randperm(numProbs, numProbs); t.X((k-1)*numProbs+1:k*numProbs) = X(idx); t.Y((k-1)*numProbs+1:k*numProbs) = Y(idx); end else k = 0; end idx = randperm(numProbs, numProbs); t.X(k*numProbs+1:end) = X(idx(1:t.NumProblems-k*numProbs)); t.Y(k*numProbs+1:end) = Y(idx(1:t.NumProblems-k*numProbs)); % Compute the products t.Z = t.X .* t.Y; end end end
- Category:
- Memory,
- Object-oriented
Comments
To leave a comment, please click here to sign in to your MathWorks Account or create a new one.