Performance: scatter vs. line

October 14, 2009 4 Comments

Following my previous article on the undocumented behavior of the scatter function, one of my readers, Benoit Charles, reported a discovery that in many circumstances the line function generates identical plots much faster than scatter.
Unlike scatter, line does not enable specific data-point marker customization, although the colors could be modified. On the other hand, line only uses a single handle object, saving memory and system resources compared to scatter keeping a separate handle for each data point. So, if you just need to quickly plot a bunch of scattered points then line could be a better choice than scatter.
Here is a simple code snippet, which generates identical plots and shows the performance difference:

>> x=rand(1000,1); y=rand(1000,1);
>> tic, for idx=1:100, cla; h=scatter(x,y); end; toc
Elapsed time is 2.521322 seconds.
>> props = {'LineStyle','none','Marker','o','MarkerEdge','b','MarkerSize',6};
>> tic, for idx=1:100, cla; h=line([x,x],[y,y],props{:}); end; toc
Elapsed time is 0.333369 seconds.

In the past, I have posted about other undocumented performance aspects, comparing the documented ismember function with the undocumented ismembc and about cellfun‘s undocumented options. If you are aware of other similar functions having identical outputs and a significant performance difference, please let me know.

4 Responses

Naor October 14, 2009 at 13:05 Reply

Now this is really interesting. Now I’m curious about generating 1000 line objects with one data point/marker each, and saving into a handle vector? Is this what scatter is doing under the hood?

Benoit October 15, 2009 at 02:54 Reply

In this case, lines are faster than scatter only when colors are used and different for each point (???). Observe this script:

%% Scatter creation
T = 0;
cla
for i=1:10
	x = rand(1,1000);
	y = rand(1,1000);
	cla
	t = cputime;
	h = scatter( x,y );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'Scatter creations: %.2fn', T )
 
%% Lines creation
T = 0;
cla
for i=1:10
	x = rand(1,1000);
	y = rand(1,1000);
	cla
	t = cputime;
	h = line( [x;x], [y;y], props{:} );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'Line creations: %.2fn', T )
 
%% Scatter modification
T = 0;
cla
x = rand(1,1000);
y = rand(1,1000);
h = scatter(x,y);
for i=1:10
	y = rand(1,1000);
	t = cputime;
	set( h, 'YData', y );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'nScatter modifications: %.2fn', T )
 
%% Lines modification
T = 0;
cla
x = rand(1,1000);
y = rand(1,1000);
h = line( [x;x], [y;y], props{:} );
for i=1:10
	y = rand(1,1000);
	t = cputime;
	set( h, {'YData'}, num2cell([y;y]',2) );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'Line modifications: %.2fn', T )
 
%% Scatter color creation
T = 0;
cla
for i=1:10
	x = rand(1,1000);
	y = rand(1,1000);
	c = rand(1000,3);
	cla
	t = cputime;
	h = scatter( x,y,50,c );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'nScatter color creation: %.2fn', T )
 
%% Lines color creation
T = 0;
cla
for i=1:10
	x = rand(1,1000);
	y = rand(1,1000);
	c = rand(1000,3);
	cla
	t = cputime;
	h = line( [x;x], [y;y], props{:} );
	set( h, {'MarkerEdgeColor'}, num2cell(c,2) );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'Line color creation: %.2fn', T )
 
%% Scatter color modification
T = 0;
cla
x = rand(1,1000);
y = rand(1,1000);
c = rand(1000,3);
h = scatter(x,y,50,c);
for i=1:10
	y = rand(1,1000);
	t = cputime;
	set( h, 'YData', y );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'nScatter color modification: %.2fn', T )
 
%% Lines modification
T = 0;
cla
x = rand(1,1000);
y = rand(1,1000);
c = rand(1000,3);
h = line( [x;x], [y;y], props{:} );
set( h, {'MarkerEdgeColor'}, num2cell(c,2) );
for i=1:10
	y = rand(1,1000);
	t = cputime;
	set( h, {'YData'}, num2cell([y;y]',2) );
	drawnow
	T = T + (cputime-t);
end
fprintf( 'Line color modification: %.2fn', T )

When you execute this script, you obtain:

Scatter creations: 0.76
Line creations: 9.53
 
Scatter modifications: 0.61
Line modifications: 1.51
 
Scatter color creation: 4.02
Line color creation: 17.07
 
Scatter color modification: 4.34
Line color modification: 1.53

Magic, isn’t it?
Sometimes, I wonder if Matlab was not developed at Hogwarts…

Naor October 15, 2009 at 15:41 Reply

very strange indeed.
damayi May 22, 2012 at 19:59 Reply

In my application, I will use line as much as possible.
What’s more, line is more simple and easy to control.

HTML tags such as <b> or <i> are accepted.
Wrap code fragments inside <pre lang="matlab"> tags, like this:

<pre lang="matlab">
a = magic(3);
disp(sum(a))
</pre>

I reserve the right to edit/delete comments (read the site policies).
Not all comments will be answered. You can always email me (altmany at gmail) for private consulting.

Click here to cancel reply.

Rolf (18 days 2 hours ago): Disclosure: This comes from the developer/vendor of the Code Checker for MATLAB. Have a look at the Code Checker for MATLAB by MonkeyProof Solutions which has many such (configurable)...
Yair Altman (25 days 19 hours ago): In theory you could use the data: URI scheme (e.g. .load('data:text/html,<html ><body>Hello world</body></html>')), but for some reason it doesn’t work well in...
Sebastian (26 days 1 hour ago): Thanks for your update to the now unsupported `HTMLBrowserPanel`. I tried to replace it with the `LightweightHelpPanel` but the application I am working on is showing generated...
Leuze (48 days 2 hours ago): Hello, I would like to add a new object into SLDD and display columns according to my object properties <>
Eric Delgado (77 days 4 hours ago): Hey Yair, I think your site is the right place to share that I wrote ccTools, a lib that allows a lot of customizations of Matlab built-in web components (such as uifigure,...
Eric Delgado (77 days 4 hours ago): Hey guys, first of all, thanks to @Yair, your site saved me a lot of times! 🙂 That’s why I am sharing that I wrote m2mlapp (*.m => *.mlapp) file conversion, that...
François (94 days 2 hours ago): more than 10 years later it helps me after struggling few hours, so thanks you very much !
Yair Altman (95 days 0 hours ago): I am not familiar with a universal catch-all tab-completion mechanism for all functions, and it also doesn’t make much sense to me, because different functions use...
Arye (95 days 18 hours ago): hi, is there a universal functionSignature.json file that will allow my to get the function hints when opening brackets of my funcs? something that automaticly recognize my function...
Sunham (98 days 13 hours ago): This is an old article, but the issue persists even in 2023. 2023a: z = mat2cell(1:1e6,1,repmat(1,1,1e 6)); f = @() cellfun(‘isempty’, z); g = @() cellfun(@isempty,z);...
Yair Altman (108 days 3 hours ago): Robot only runs when you tell it to run a command such as keyPress. If you don’t tell it to run a command, it uses no CPU, so there’s no need to remove the Robot...
Eric (108 days 14 hours ago): Hey @Kevin, can you share your code about create group of figures in the AppContainer? The container of multiples uifigures could be an amazing improvement over AppDesigner and its...
Elsa Smith (109 days 5 hours ago): I recently used java.awt.Robot to perform GUI testing on MATLAB and found it to be an extremely easy and useful way to control mouse movements.
Elsa Smith (109 days 5 hours ago): I’m suspecting that the slow performance of my GUI may be due to the use of java.awt.Robot. Is there a way to cancel/stop/remove the robot after it has been created, or is...
Michelle Kline (109 days 22 hours ago): *edit* tip about fopen(), not about fwrite(). ‘Wb’ vs. ‘wb’

Related posts:

4 Responses

Leave a Reply