title: C++ Libraries for solving Complex Linear systems Ax=b
body: <p>I am interested in solving a sparse complex linear system Ax=b where A is a square matrix of complex numbers and b is vector of complex numbers. </p>
  
  <p>If possible I would like such a library to be templated (for the ease of installation and use)
  sth in the spirit of <a href="http://eigen.tuxfamily.org/index.php?title=Main_Page" rel="nofollow">Eigen</a></p>
  
  <p>I checked out Eigen but it does not, I think,  look like it supports solving linear equations with complex sparse matrices, (although one can create and do elementary operations on complex matrices.) </p>
  
  <p>Another trick someone suggested to me was one can work around this, by solving an extended system of twice the dimension using the fact that (A1 + iA2)(x1 + ix2) = (b1 + ib2)
  but I would prefer some simple black box which gets the job done. </p>
  
  <p>Any suggestions?</p>