better indentation invPMatrix
insert TD consist method to compute UW reconstruction