octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [OctDev] linear-algebra: new QR updating functions


From: David Bateman
Subject: Re: [OctDev] linear-algebra: new QR updating functions
Date: Tue, 19 Feb 2008 15:53:38 +0100
User-agent: Thunderbird 2.0.0.6 (X11/20070914)

Jaroslav Hajek wrote:
> Hello,
> to anyone who might be interested:
> I have just committed an initial revision of functions
> "qrupdate","qrinsert" and "qrdelete" into the main/linear-algebra
> package. These functions should provide equivalent behaviour to their
> Matlab counterparts, i.e. update a QR decomposition after a rank-1
> update, adding/deleting row/column. qrinsert and qrdelete currently
> support adding/deleting columns, rows are on the way.
>
> enjoy,
>
>   
These are core matlab functions (ie not in a toolbox) and there should
also be considered for inclusion in Octave 3.1.x. Is there any reason
why you wrote the code itself in F77? Is it just that is encapsulates
several blas/lapack functions and so to do that in C++ would require
additional wrapping.

I'd suggest reworking this as a patch against octave, with the F77 code
in a directory libcruft/qrupdate (or something similar) and the C++
interfaces in src/DLD-FUNCTIONS.. The Makefile.in files will need to be
modified appropriately.

John, we should include these functions.. Are you happy with the F77
wrapper to blas/lapack function to have only a single F77 function in
the Octave C++ wrapper or should that be modified to use a pure C++ wrapper?

D.

-- 
David Bateman                                address@hidden
Motorola Labs - Paris                        +33 1 69 35 48 04 (Ph) 
Parc Les Algorithmes, Commune de St Aubin    +33 6 72 01 06 33 (Mob) 
91193 Gif-Sur-Yvette FRANCE                  +33 1 69 35 77 01 (Fax) 

The information contained in this communication has been classified as: 

[x] General Business Information 
[ ] Motorola Internal Use Only 
[ ] Motorola Confidential Proprietary



reply via email to

[Prev in Thread] Current Thread [Next in Thread]