octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: GSL in octave


From: Julien Bect
Subject: Re: GSL in octave
Date: Thu, 21 Jul 2016 22:39:22 +0200
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:45.0) Gecko/20100101 Icedove/45.1.0

Le 21/07/2016 à 09:36, Susi Lehtola a écrit :
On 07/21/2016 12:08 AM, Julien Bect wrote:
I think we have a misunderstanding about what "broadcast" means.

Even if you have a function with four arguments, broadcasting (as I
understand it) doesn't mean that you will end up with a 4-dimensional
(or 8-dimensional ???) array.

You'd get a 4-dimensional tensor if your definition was to use a direct tensor product, i.e. every dimension in the input arguments constitutes a dimension of the output.

The templates should be pretty straightforward to modify for broadcasting. What you'd need to do is first figure out the common dimensions, then pad everything to the common dimensions, and just change the calculation and return sections. This only needs to be done in a single place, namely, the generate_template.sh script.

I just realized that all the *.cc.template files can be regenerated by the Makefile. So I pushed a changeset that removes all the templates from the repo and fixes the Makefile. It seems to work ok, but you can check.

I don't have time to work on implementing the broadcasting myself, sorry.

@++
Julien



reply via email to

[Prev in Thread] Current Thread [Next in Thread]