What is drawnow doing?

Question

Alessandro Masullo 2015-11-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/254641-what-is-drawnow-doing

评论： Alessandro Masullo 2015-11-18

Hello everyone.

I'm working on a code where I need to alternate the display of two complex plots in a very short time. The time for creating each plot is quite high (4-5 seconds) and I want to alternate them with a lag of 0.1 seconds. Creating a new plot every 0.1 seconds is therefore completely useless.

My idea was to create the two plots in two different axes and alternate the visibility of them. The problem is that drawnow keeps taking a very long time just to show one axes instead of the other.

Is there any way I can speed up this process?

Example of code:

figure
x1 = rand(100000,1);
y1 = rand(100000,1);
x2 = rand(100000,1);
y2 = rand(100000,1);
t = tic;
delay = 0.5;
flag = 0;
while 1
    if toc(t) > delay
        if flag
            plot(x1,y1,'.')
            flag = 0;
        else
            plot(x2,y2,'.')
            flag = 1;
        end
        drawnow
    end
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2015-11-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/254641-what-is-drawnow-doing#answer_199671

编辑：Walter Roberson 2015-11-14

在 MATLAB Online 中打开

figure
x1 = rand(100000,1);
y1 = rand(100000,1);
x2 = rand(100000,1);
y2 = rand(100000,1);
p1 = plot(x1,y1,'.', 'visible', 'off');
p2 = plot(x2,y2,'.', 'visible', 'off');
t = tic;
delay = 0.5;
flag = 0;
while 1
    if toc(t) > delay
        if flag
            set(p2, 'visible', 'off');
            set(p1, 'visible', 'on');
            flag = 0;
        else
            set(p1, 'visible', 'off');
            set(p2, 'visible', 'on');
            flag = 1;
        end
        drawnow
    end
end

9 个评论
显示 7更早的评论隐藏 7更早的评论

Alessandro Masullo 2015-11-13

在 MATLAB Online 中打开

Thank you for the answer, Walter, but this is the answer that I suggested too (not in the code, though), and it is still too slow.

If you run the following code, you can see that the actual delay between the two plots is twice the lag (in my case). If the plot is even heavier (rand(100000,4)), the time for drawing the axis reaches up to two seconds!

figure,hold on
x1 = rand(100000,2);
y1 = rand(100000,2);
x2 = rand(100000,2);
y2 = rand(100000,2);
p1 = plot(x1,y1,'.', 'visible', 'off');
p2 = plot(x2,y2,'.', 'visible', 'off');
t = tic;
delay = 0.1;
flag = 0;
while 1
    t2 = tic;
    if toc(t) > delay
        if flag
            set(p1, 'visible', 'off');
            set(p2, 'visible', 'on');
            flag = 0;
        else
            set(p1, 'visible', 'on');
            set(p2, 'visible', 'off');
            flag = 1;
        end
        drawnow
        disp(toc(t2))
    end
end

Another solution that I had in mind was using getframe and alternate the two buffered getframes instead of the real plots. Though I wouldn't like to use this solution because I would need to take care of figure resizing and this kind of stuff...

Mike Garrity 2015-11-13

在 MATLAB Online 中打开

You're on the right track here. The basic idea is that we want to park the two plots on the graphics card, and then just tell the card which one to draw each frame. The code you've got almost does that.

Here's a version I tweaked a tiny bit:

figure
hold on
drawnow
rng default
x1 = rand(100000,2);
y1 = rand(100000,2);
x2 = rand(100000,2);
y2 = rand(100000,2);
p1 = plot(x1,y1,'.', 'visible', 'off');
p2 = plot(x2,y2,'.', 'visible', 'off');
nframes = 500;
ft = zeros(1,nframes);
t = tic;
delay = 0.1;
flag = 0;
for i=1:nframes
  t2 = tic;
  if flag
    set(p1,'Visible','off');
    set(p2,'Visible','on');
    flag = 0;
  else
    set(p1,'Visible','on');
    set(p2,'Visible','off');
    flag = 1;
  end
  drawnow
  ft(i) = toc(t2);
end
sprintf('%g FPS', 1/median(ft))

When I run this on my machine (an nVidia K600) I get 38 frames per second.

But we can do better. The problem with this is that we're changing a property on the line object. When we do this, MATLAB decides something's changed and it resends the data to the card. We don't want to do that.

But how do we change the visibility without modifying the object? The trick is to put the line object in a group and toggle the visibility on that. That would look something like this:

figure
drawnow
rng default
x1 = rand(100000,2);
y1 = rand(100000,2);
x2 = rand(100000,2);
y2 = rand(100000,2);
g1 = hggroup;
p1 = plot(x1,y1,'.','Parent',g1);
g2 = hggroup;
p2 = plot(x2,y2,'.','Parent',g2);
nframes = 500;
ft = zeros(1,nframes);
t = tic;
delay = 0.1;
flag = 0;
for i=1:nframes
  t2 = tic;
  if flag
    g1.Visible = 'off';
    g2.Visible = 'on';
    flag = 0;
  else
    g1.Visible = 'on';
    g2.Visible = 'off';
    flag = 1;
  end
  drawnow
  ft(i) = toc(t2);
end
sprintf('%g FPS', 1/median(ft))

When I run this version, I get 90 frames per second.

Does that help?

Mike Garrity 2015-11-17

It's really not an answer yet. We need to figure out why you are getting softwareopengl. There are a couple of reasons this might happen. Two common reasons.

We have a blacklist for some drivers we've seen a lot of problems with.
The system that reports crashes back to us looks at the call stack. If the crash was in your graphics card driver, we'll set your default to softwareopengl to prevent it happening again.

Could you start MATLAB with the -hardwareopengl switch and then do opengl info again? That's probably the easiest way to get all of the details about your graphics card. I'll see if I can find someone who can help figure out what's going on. Hopefully it'll be as simple as updating your driver to a different version.

Alessandro Masullo 2015-11-18