Subject: Initial- and boundary-value problems in Maxima.
From: Richard Fateman
Date: Thu, 15 Sep 2011 11:42:50 -0700
On 9/15/2011 11:24 AM, Raymond Toy wrote:
...
>
> Yes, but it probably won't be fast in colnew. :-( From a quick grep
> of the code, daxpy is called in 3 different places, and each place
> calls daxpy with a slice of a larger array. F2cl tries to be smart
> about this and figures out the underlying array. That could be
> expensive, depending on how many intervening slices have been done.
> It might be cheaper to just use the normal displaced-array access.
The only other hack that comes to mind is to expand daxpy inline,
perhaps as a macro, and that might not help either. Or call fortran.
RJF