is it possible to vectorise the following. I want to construct all pairs of integers (I,j) with i^2 + j^2 <= r^2, This includes negative i and j. Ideally I want a 2 by k matrix listing my pairs in lexicographic order. At the moment I use for loops.

Question

David Epstein 2018-7-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411348-is-it-possible-to-vectorise-the-following-i-want-to-construct-all-pairs-of-integers-i-j-with-i-2

评论： David Epstein 2018-7-20

Given r, I want to construct all pairs of integers (i,j) with i^2 + j^2 <= r^2, This includes negative i and j. Ideally I want a k by 2 matrix listing my pairs in lexicographic order (either i_1 < i_2 or i_1==i_2 and j_1 < j_2). Can this be vectorised in Matlab? At the moment I use for loops, which doesn't matter, because I don't have to do it often. But I'm thinking a might need the fastest possible method for doing something similar in another program.

2 个评论
显示无隐藏无

David Epstein 2018-7-20

I should have said that I am using Matlab R2018a on a MacBook Pro running MacOS 10.13.6

David Epstein 2018-7-20

mostly quite small at the moment (up to around 20), but I might contemplate going up to r=200. However, in that case someone would probably do the programming for me in python or C++ on a very large and powerful machine.

请先登录，再进行评论。

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

Answer 1

James Tursa 2018-7-20

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411348-is-it-possible-to-vectorise-the-following-i-want-to-construct-all-pairs-of-integers-i-j-with-i-2#answer_329662

在 MATLAB Online 中打开

One way:

[x,y] = meshgrid(-r:r,-r:r);
z = x.*x + y.*y <= r*r;
result = [x(z) y(z)];

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

David Epstein 2018-7-20

Very elegant.

请先登录，再进行评论。

is it possible to vectorise the following. I want to construct all pairs of integers (I,j) with i^2 + j^2 <= r^2, This includes negative i and j. Ideally I want a 2 by k matrix listing my pairs in lexicographic order. At the moment I use for loops.

2 个评论
显示无隐藏无

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

is it possible to vectorise the following. I want to construct all pairs of integers (I,j) with i^2 + j^2 <= r^2, This includes negative i and j. Ideally I want a 2 by k matrix listing my pairs in lexicographic order. At the moment I use for loops.

2 个评论 显示 无隐藏 无

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

2 个评论
显示无隐藏无

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