found another error in the dCons/dPrim matrices.  thus concludes our errors, the system (in full magnetic field, in absense of degenerate cases) multiplies to identity now.