Re: Matrix decompositions.

Salman Javaid <[email protected]>
Newsgroups gmane.comp.lib.boost.ublas
Message-ID <CAEKdXi1iQJZGmzsdU64Wv--Oww16q6AcxD0b5BeDHeXK4n6Y4Q@mail.gmail.com>
Thank you, Oswin for the detailed response. I am going to update the code.

David, Nasos, any advise on coding conventions? Or anything else that you
can possible suggest? I will stand grateful.





Best Regards,
Salman Javaid


On Tue, May 14, 2013 at 10:53 PM, oswin krause <
[email protected]> wrote:

>  Hi,
>
> in the order I stumbled over the things:
>
> main.cpp
> line 44-54: you don't need a copy, instead you should use a combination of
> row/subrange.
> line 58-60: you should take a look at inner_prod
> line 63: 0.00025 is too big.
> line 66: You should never create H explicitly.
> line 67: because you formed H, this step is O(n^3) which makes the whole
> algorithm O(n^4). This can be done in O(n^2)
> line 73-79: same applies here.
>
> Greetings,
> Oswin
>
> On 14.05.2013 22:12, Salman Javaid wrote:
>
>  Hello uBLAS Contributors:
>
>                                       I have applied to GSoC 2013 and
> pitched implementation of SVD factorization for uBLAS. In order to better
> prepare myself and to get my hands dirty at uBLAS, I ended up implementing
> QR Factorization employing Householder Reflections using uBLAS. This is
> only the first draft and will be needing significant improvement, e.g.,
> computation of QR decomposition of 100 * 100 matrix takes around 30
> seconds. But I guess just to get familiar with code base, it was a good
> exercise. Over the next week or two I will be trying to optimize the code.
>
>
>                                       I will be absolutely grateful if
> contributors can have a quick glance at the code, and point me to any
> improvements they can suggest. Particularly in speeding up matrix
> multiplication.
>
>  I used Visual Studio 2010 to compile the code. I will try to get the
> code running on my Ubuntu machine in a couple of days hopefully.
>
>  Here the header file:
> https://github.com/salmanjavaid/QR_Decomposition/blob/master/QR_Header.hpp
>
>
>  The main file:
> https://github.com/salmanjavaid/QR_Decomposition/blob/master/Main.cpp
>
>
> Best Regards,
>  Salman Javaid
>
>
> _______________________________________________
> ublas mailing [email protected]://lists.boost.org/mailman/listinfo.cgi/ublas
> Sent to: [email protected]
>
>
>
> _______________________________________________
> ublas mailing list
> [email protected]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
> Sent to: [email protected]
>

_______________________________________________
ublas mailing list
[email protected]
http://lists.boost.org/mailman/listinfo.cgi/ublas
Sent to: [email protected]
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.