Re: [ublas] [bindings] Routine with 4 workarrays

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