yet more bugs fixed, reconstruction creates A[i][j] = 2 for all i, j, need to look at Eigen routines