How can I store dates as arrays for further process when using Matlab GPU computing?

Question

0 个投票

In the following attached code, I would like to record the variable Z in the while loop as array for further process, for example, finding the max value of Z in the array. But the dynamic array is not supported in the arrayfun for GPU computing. How to solve this problem?

function test_GPU_computing
maxIterations=500;
gridSize=1500;
xlim=[-0.75, -0.73];
ylim=[ 0.12,  0.14];
t=tic();
x=gpuArray.linspace( xlim(1), xlim(2), gridSize );
y=gpuArray.linspace( ylim(1), ylim(2), gridSize );
[xGrid,yGrid] = meshgrid( x, y );
count=arrayfun( @tar_fun,xGrid, yGrid, maxIterations );
count=gather(count); 
gpuArrayfunTime=toc(t)
figure(3)
imagesc(x,y,count)
reset(gpuDevice(1))
function count=tar_fun(x0,y0,maxIterations)
z0=complex(x0,y0);
z=z0;
count=1;
while (count <= maxIterations) && (abs(z) <= 2)
    count=count+1;
    z=z*z+z0;
end
count=log(count);

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

Follow Question

Answer 1

Walter Roberson 2018-12-15

0 个投票

It does not need to be dynamic. You have maxIterations, so you can create it that size. When you need the "dynamic" version of it, index to the number actually used.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Walter Roberson 2018-12-17

Unfortunately my OS is too old to support the CUDA driver needed for R2018b :(

Shan Yin 2018-12-17

Thank you, Walter Roberson.

I have changed my code to try the nested function, but the situation did not improve. Though the array can be indexed if it is not assigned, it does not reach to my purpose.

So, can you tell me whether I can achieve my goal by using the arrayfun funtion?

clear all
clc
maxIterations=100;
gridSize=300;
xlim=[-0.75, -0.73];
ylim=[ 0.12,  0.14];
t=tic();
x=gpuArray.linspace(xlim(1), xlim(2), gridSize);
y=gpuArray.linspace(ylim(1), ylim(2), gridSize);
[xGrid,yGrid]=meshgrid(x,y);
Pos=gpuArray.zeros(maxIterations,1);
count=parent_fun(xGrid,yGrid,maxIterations,Pos);
count=gather(count);
gpuArrayfunTime=toc(t)
figure(1)
imagesc(x,y,count)
reset(gpuDevice(1))
function result=parent_fun(xGrid,yGrid,maxIterations,Pos)
         function count=tar_fun(x0,y0)
                  z0=complex(x0,y0);
                  z=z0;
                  count=1;
                  while (count<=maxIterations) && (abs(z)<=2)
                         count=count+1;
                         z=z*z+z0;
                         % z=z*z+z0+Pos(count,1); % This line can be executed, but the following line can not!
                         Pos(count,1)=z;  % where error arises! 
                  end 
                  count=max(log(count),log(abs(z)));
         end
         result=arrayfun(@tar_fun, xGrid, yGrid);
end

Answer 2

Edric Ellis 2018-12-17

在 MATLAB Online 中打开

0 个投票

If you know the reduction operation you wish to perform on z, then you might be able to update a "running total" as the while loop proceeds, like so:

function [count,maxZ]=tar_fun(x0,y0,maxIterations)
z0=complex(x0,y0);
z=z0;
count=1;
% Initialize maxZ
maxZ = complex(0);
while (count <= maxIterations) && (abs(z) <= 2)
    z=z*z+z0;
    % Update the value of maxZ
    if count == 1 || abs(z) > abs(maxZ)
        maxZ = z;
    end
    count=count+1;
end
count=log(count);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Shan Yin 2018-12-17

Dear Edric Ellis,

Thanks for your kind reply and the code you provide.

In my question, finding the max value of variable z was taken as example to explain my purpose of using an array. However, it is not all of my goal. Indeed, I plan to transfer this simple case to some more realistic and complicated cases, e.g., dynamic systems described by ordinary differential equation. In that situation, an array would be required to record the evolution of every system state variable, to further determine the long-term behavior of the dynamic systems.

So, in this posted question, I do need to know directions for using an array in the arrayfun function, rather than finding the max value of variable z simplely. Hope I have explained my purpose clearly.

Thanks again for your comment. And please tell me if you have more ideas on how to use array in the arrayfun function.

Shan Yin

How can I store dates as arrays for further process when using Matlab GPU computing?

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

回答（2 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

产品

版本

标签

Community Treasure Hunt

How can I store dates as arrays for further process when using Matlab GPU computing?

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

回答（2 个）

3 个评论 显示 1更早的评论 隐藏 1更早的评论

1 个评论 显示 -1更早的评论 隐藏 -1更早的评论

产品

版本

标签

另请参阅

Community Treasure Hunt

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

3 个评论
显示 1更早的评论隐藏 1更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论