new procedure
new procedure
added diagnistic printf
use known-before divisor
very last bug fixed
test_hermitian_decomp passed