comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » multiplication by a diagonal matrix
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
multiplication by a diagonal matrix [message #40812] Fri, 03 September 2004 03:53 Go to next message
lbusett is currently offline  lbusett
Messages: 9
Registered: March 2004
Junior Member
Hi all,

I have the following problem: given a matrix A(n,m) and a vector of
weighting factors w(n), i need to multiply each row of the matrix
A(i,*)by the corresponding weighting factor w(i).

I know that I can simply "transform" the w vector into a diagonal
matrix with diag_matrix and then multiply it with A (e.g.: result =
A##diag_matrix(w)), but for large values of n this solution is very
slow.

Can anybody suggest me a faster approach to solve this problem ?

Thanks in advance for the help,

Lorenzo Busetto

Remote Sensing Lab.
University of Milano-Bicocca.
Re: multiplication by a diagonal matrix [message #40927 is a reply to message #40812] Wed, 08 September 2004 06:33 Go to previous message
David Fanning is currently offline  David Fanning
Messages: 11724
Registered: August 2001
Senior Member
Lorenzo Busetto writes:

> Well, when I posted my message I'd never expected to increase the
> speed 800000 times ! Can it be possible ?

You haven't been hanging out very long on the IDL
newsgroup. It's not only possible, it is the humdrum
around here. :-)

http://www.dfanning.com/code_tips/slowloops.html

Cheers,

David

--
David W. Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http:/www.dfanning.com/
Phone: 970-221-0438, IDL Book Orders: 1-888-461-0155
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Re: multiplication by a diagonal matrix
Next Topic: 6.1 IDLDE project amnesia

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Sat Oct 11 15:31:29 PDT 2025

Total time taken to generate the page: 0.72079 seconds