big question: the sparse matrix provided by Matrix package in R don't support complex numbers, so I write a function to create dense matrix with complex number with sme parameters temporarily named powerDenseMatrix