#!/bin/csh
# usage: dospd file.rsa.gz
#
# file.rsa.gz is a compressed Harwell/Boeing file containing
# a symmetric positive definite matrix
#
echo '================================================================='
echo 'Matrix:: ' $1:t:r:r

echo '3'	>  tmp/control.umf4 ; # 1: print level
echo '0.2'	>> tmp/control.umf4 ; # 2: dense row control
echo '0.2'	>> tmp/control.umf4 ; # 3: dense col control
echo '0'	>> tmp/control.umf4 ; # 4: pivot tol (offdiag) NON-DEFAULT 
echo '32'	>> tmp/control.umf4 ; # 5: block size
echo '3'	>> tmp/control.umf4 ; # 6: symmetric strategy, NON-DEFAULT
echo '0.7'	>> tmp/control.umf4 ; # 7: initial alloc
echo '2'	>> tmp/control.umf4 ; # 8: max iter. refinement
echo '1'	>> tmp/control.umf4 ; # 9:
echo '0'	>> tmp/control.umf4 ; # 10:
echo '0'	>> tmp/control.umf4 ; # 11:
echo '0'	>> tmp/control.umf4 ; # 12:
echo '0.01'	>> tmp/control.umf4 ; # 13: 2-by-2 tolerance
echo '0'	>> tmp/control.umf4 ; # 14: Q fixed (auto)
echo '10'	>> tmp/control.umf4 ; # 15: AMD dense row control
echo '0'	>> tmp/control.umf4 ; # 16: diag pivot tolerance, NON-DEFAULT
echo '0'	>> tmp/control.umf4 ; # 17: scaling, NON-DEFAULT
echo '0.5'	>> tmp/control.umf4 ; # 18: frontal matrix alloc.
echo '0'	>> tmp/control.umf4 ; # 19: drop tolerance
echo '1'	>> tmp/control.umf4 ; # 20: AMD/COLAMD aggressive absorption

zcat $1 | readhb_nozeros > tmp/A
zcat $1 | readhb_size    > tmp/Asize
umf4 s

echo '================================================================='