; IDL Version 7.1.1 (linux x86_64 m64)
; Journal File for mats@sally
; Working directory: /home/mats/idl/develop/matrixalgebra
; Date: Thu Apr  5 16:19:19 2012
 
;The LA_SVD routine works with the ## conventions.
; 
;------ |A| (op:##) [2x3] ------
;    0     1     2
;    3     4     5
;-------------------------------
;Now call "LA_SVD, A, W, U, V"
; 
; % Program caused arithmetic error: Floating underflow
;------ |U| (op:##) [2x2] ------
;  -0.275   -0.962
;  -0.962    0.275
;-------------------------------
;------ |diag_matrix(W)| (op:##) [2x2] ------
;   7.348    0.000
;   0.000    1.000
;--------------------------------------------
;------ |VT=transpose(V)| (op:##) [2x3] ------
;  -0.393   -0.561   -0.729
;   0.824    0.137   -0.549
;---------------------------------------------
; 
;Now try to multiply U, W, and VT and make sure we
;get the original matrix back.
; 
; 
;------ |U ## diag_matrix(W) ## transpose(V)| (op:##) [2x3] ------
;  -0.000    1.000    2.000
;   3.000    4.000    5.000
;-----------------------------------------------------------------
; 
;------ |A| (op:##) [2x3] ------
;   0.000    1.000    2.000
;   3.000    4.000    5.000
;-------------------------------
; 
;------ |A - (U ## diag_matrix(W) ## transpose(V))| (op:##) [2x3] ------
;   0.000    0.000    0.000
;   0.000    0.000    0.000
;-----------------------------------------------------------------------