comp.lang.idl-pvwave archive: archive » Re: Svdc

Home » Public Forums » archive » Re: Svdc_complex

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Svdc_complex [message #16151]

Fri, 02 July 1999 00:00

Javier Sanchez Gonzal is currently offline

Javier Sanchez Gonzal
Messages: 3
Registered: July 1999

Junior Member

Hi all,
if somebody wants to test this I send the procedure to calculate the
singular value of a complex matrix.

pro svdc_complex, M, Ew2, Eu2, Ev2

xdim=(size(M))(1)
ydim=(size(M))(2)
Ew2=complexarr(xdim)
Eu2=complexarr(xdim,ydim)
Ev2=complexarr(xdim,ydim)
M2=dblarr(2*xdim,2*ydim)
M2(0:xdim-1,0:ydim-1)=double(M)
M2(xdim:2*xdim-1,0:ydim-1)=-imaginary(M)
M2(0:xdim-1,ydim:2*ydim-1)=imaginary(M)
M2(xdim:2*xdim-1,ydim:2*ydim-1)=double(M)
catch, Error_status
if Error_status ne 0 then begin
; print, !err_string
Minv=0.*transpose(M)
end else begin
svdc, m2,w,u,v,/double
for i=0, xdim-1 do begin
Ew2(i)=w(2*i)
Eu2(i,*)=u(2*i,0:ydim-1)+dcomplex(0,1)*u(2*i,ydim:2*ydim-1)
Ev2(i,*)=v(2*i,0:ydim-1)+dcomplex(0,1)*v(2*i,ydim:2*ydim-1)
endfor
end

end

Thanks all

Javier Sanchez Gonzalez

Report message to a moderator

[Message index]

Re: Svdc_complex

By: Javier Sanchez Gonzal on Fri, 02 July 1999 00:00

Previous Topic:	Can an object's structure be redifined without restarting IDL
Next Topic:	Iso -Surface in Object Graphics

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Oct 12 09:47:27 PDT 2025

Total time taken to generate the page: 1.19906 seconds