add test subDCA
use JAMA matrix transpose
break when reach max loop