How can i use nonlinear conjugate gradient (ncg) method from snes solver?

Question

I try to use the nonlinear conjugate gradient method using the snes solver (PETScSNESSolver) but is is not obvious how to properly use it.

I have implemented a code where this method is called from scipy.optimize but i would like a parallelized version.

Thank you.