How to speed up the process of determining if a point is inside a cylinder

Question

Christos 2012-7-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43534-how-to-speed-up-the-process-of-determining-if-a-point-is-inside-a-cylinder

评论： Bruno Luong 2020-9-4

I would like to speed up the process of determining if a point is inside a cylinder. Specifically:

A 3D cubic mesh is used to divide a problem. In 3 separate 3D arrays the cell center coordinates are stored (i.e cx, cy, cz). A cylinder is upright either in the x, y, or z direction. Then based on the cap center coordinates, the height of the cylinder, and the radius of the cylinder, the following code is used to determine if a point is inside the cylinder:

% nx, ny, nz is the number of cells in cx, cy or cz.
% range = distance between the cap centers of the cylinder
% For all points in the mesh
for lind = 1:nx
    for jnd = 1:ny
        for knd = 1:nz
            % Get cell center coordinates (x, y, z)
            apoint = [cx(lind,jnd,knd) cy(lind,jnd,knd) cz(lind,jnd,knd)];
            % Project point on center line
            tmin(lind, jnd, knd) = (dot(apoint, range) - dot(cap1, range))/ ...
            dot(range, range);
            linept = cap1 + (tmin(lind, jnd, knd) * range);
            % Distance form center line
            diff = linept - apoint;
            dist(lind, jnd, knd) = sqrt(dot(diff, diff));
        end
    end
end
% Get all the indices of the points that are part of the cylinder
l = find((((tmin>=0) & (tmin<=1)) & (dist<=radius)));

How can I speed up the computation of determining whether a point is part of a cylinder specifically avoid using for-loops?

Any help would be much appreciated!

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

请先登录，再进行评论。

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

Answer 1

Doug Hull 2012-7-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43534-how-to-speed-up-the-process-of-determining-if-a-point-is-inside-a-cylinder#answer_53457

在 MATLAB Online 中打开

For loops are not always a problem. Run this through the profiler to find where the slow down is. I would do something like this.

pointXVec, pointYVec, pointZVec are locations TopZ, BotZ, are the locations of top and bottom CenterX, CenterY, are location in X and Y Radius

These can be switched around for other orientations.

flagBelowVec = (pointZVec <= topZ);
flagAboveVec = (pointZVec >= botZ);
radialDistanceSquaredVec = (pointXVec-centerX)^2 + (pointYVec-centerY)^2;
flagInsideVec = (radialDistanceSquaredVec <= radius^2);
flagIsIn = (flagBelowVec & flagAboveVec & flagInsideVec);