@tomleslie  Thank you for kind explanation. The matrix is actually almost upper triangular, but I see that it takes very huge memories to store it. Maybe I should figure out how to deal with this problem without generating the whole matrix as you mentioned. Thanks!