add xassert, read_matrix
add transpose_matrix
add matrix.h