Hi Ben,
use the CG-Solver instead of GMRes, then it works
Code:
inverse= CGSolver(a.mat, c.mat, precision=1e-6, maxsteps=200)
CG can solve with complex-symmetric matrices.
Exact= (-5.0265482451562615+6.702064326875025j)
Calculated= (-5.026398383236931+6.701864510664208j)
HalfandHalf= (-5.0264730864946525+6.701964115326205j)
Joachim