Added an explicit field_p parameter for CG solver and fixed the halo updates in the C version so that it now gives the correct answer (for CG)