Dear @arousta,
According to About the functions LinearCG a symmetric elimination is possible with tgv=-2, although I never tried myself, and it is not available for the rhs, as far as I remember. On the other hand, the implementation is not that hard; at Exact Dirichlet BC I have presented a naive implementation at the level of FreeFEM itself.
Regards,
/Fotios Kasolis