Additionally, I want to enter the components of A matrix and the rhs vector as double complex. So how can build the libraries in double complex. As far as i can see there is only --with-scalar-type=complex option during compiling petsc. Thanks EVRIM