#!/usr/bin/python -Wall

# ================================================================
# John Kerl
# kerl.john.r@gmail.com
# 2008-02-20
# ================================================================
# This software is released under the terms of the GNU GPL.
# Please see LICENSE.txt in the same directory as this file.
# ================================================================

from __future__ import division # 1/2 = 0.5, not 0.
import sys
from sackmatc_m import *

argc = len(sys.argv)
if (argc != 4):
	print >> sys.stderr, "Usage: %s {A file name} {U file name} {P file name}" % (sys.argv[0])
	print >> sys.stderr, "Use - as file name(s) to print to standard output."
	sys.exit(1)

A_file_name = sys.argv[1]
U_file_name = sys.argv[2]
P_file_name = sys.argv[3]

A = read_matrix(complex, A_file_name)
[U,P] = A.polar_decomp()
write_matrix(U,U_file_name)
if ((U_file_name == "-") and (P_file_name == "-")):
	print
write_matrix(P,P_file_name)