Problem in using bsxfun

Question

AP 2014-5-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/131939-problem-in-using-bsxfun

编辑： James Tursa 2014-5-31

I have two arrays in MATLAB:

    A; % size(A) = [NX NY NZ 3 3]
    b; % size(b) = [NX NY NZ 3 1]

These two arrays are very big. In fact, in the three dimensional domain, I have two arrays defined for each (i, j, k) which are obtained from above-mentioned arrays A and b, respectively and their sizes are [3 3] and [3 1], respectively. Let's for the sake of example, call these arrays m and n.

    m; % size(m) = [3 3]
    n; % size(n) = [3 1]

How can I solve m\n for each point of the domain in a vectorize fashion? I used bsxfun but I am not successful.

solution = bsxfun( @(A,b) A\b, A, b );

I think the problem is with the expansion of the singleton elements and I don't know how to fix it.

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

请先登录，再进行评论。

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

Answer 1

James Tursa 2014-5-31

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/131939-problem-in-using-bsxfun#answer_139004

编辑：James Tursa 2014-5-31

For a multi-dimensional paged based linear solver, bsxfun is not the correct approach. You will need to either write loops, or use some 3rd party utility software. You might try this package by Bruno Luong:

http://www.mathworks.com/matlabcentral/fileexchange/27762-small-size-linear-solver

However, you will need to permute your data so the 2D pages are the first two dimensions (gets each page contiguous in memory).