Re: Vectorizing interpolation?

[Juan Pablo Carbajal]

On Tue, Aug 9, 2016 at 9:03 PM, Nir Krakauer <address@hidden> wrote:
> Is it possible to vectorize something like the following loop?
>
> m = 1000; n = 30;
>
> A = sort (randn(m, n), 2);
>
> b = randn (m, 1);
>
> c = nan (m, 1);
>
> for i = 1:m
>
>   c(i) = interp1 (A(i, :), 1:n, b(i));
>
> endfor
>
>
> _______________________________________________
> Help-octave mailing list
> address@hidden
> https://lists.gnu.org/mailman/listinfo/help-octave
>

It seems the difficulty lies on interp1 using only a single x-vector
(multiple y-vlaues would work). Can you invert the function to
interpolate? i.e. interp1 (1:n, A(i, :), zb(i))?
Otherwise I guess array/cellfun arte you only hope (or they parallel
implementations in the parallel pkg)