[Octave-bug-tracker] [bug #38628] bsxfun slow for complex

[Rik]

Follow-up Comment #14, bug #38628 (project octave):

Well, it might get messy in the same way that cellfun gets messy.  In cellfun
certain functions are detected and accelerated, and the rest are processed
normally.

I think it might be similarly useful if functions that are accelerated via
broadcasting, such as @plus, were simply forwarded to the broadcasting
operator.

But, we also need to be able to handle arbitrary binary functions.  For
example, this code


fun = @(a,b) a - exp(b);
a = 1:7; 
b = pi*[0 1/4 1/3 1/2 2/3 3/4 1].';
C = bsxfun(fun,a,b)


For those, we should do everything possible to speed things up.

It is not, in general, possible to know the real/complex nature of the result
without calculating the function.  For specific operations, say @plus, it is
possible to work out the return type in advance given the inputs.  But
something like


f = @(a,b) = sqrt (a - b);


Can produce real or complex output even though the inputs a,b are both real.


    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?38628>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/