Subject: [ublas] [bindings] Workspaces
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2008-12-10 13:27:56
Could it be that the option for 3 workspace arguments is missing? E.g.,
AFAICT, cgelsd could be called with three workspace arguments (WORK, RWORK,
IWORK). At the moment, this is not supported through workspace.hpp.
I also see is that dgelsd.f defines two workspace arrays, which seems to be
contradicting the assumption that the number of workspaces is a 1:1
connection to the value_type of the underlying.
Are these correct observations?
Perhaps it's safer to select the number of workspaces by the number of
workspaces defined in the Fortran code :-).
Another question: are there cases where the formulae for the
minimum_workspace requirements differ for each value_type? (I.e., 4 cases)?
Or is two formulae the maximum (real vs. complex)? I remember I saw
something along the lines of template specialisations based on
real/complex/mixed scenarios, which would make it three, or is this based
on something else?