Subject: Re: [ublas] [bindings] Routine with 4 workarrays
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2009-01-18 07:04:54
Karl Meerbergen wrote:
>
> In the case of DGGSVP, I would not use 3 but two workarrays: WORK and
> TAU can be held by a single double array, that you will have to
> partition in two subarrays.
Ok, so suppose we increase the size for array work with the size
requirements for tau, and then calling it with something like
for tau: traits::vector_storage(work.select(real_type()))
for work: traits::vector_storage(work.select(real_type())) + N
?
Thanks,
Rutger