Boost logo

Ublas :

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(
for work: traits::vector_storage( + N