In recent years, coupled physics imaging techniques have been developed to produce clearer images than those produced by electrical impedance tomography. This paper focuses on the inverse problem arising in current density impedance imaging and magneto-acousto-electric tomography. 
 We consider the electrostatic equation $\nabla\cdot(\sigma\nabla w_b) = 0$ in a bounded domain $\Omega\subset\mathbb{R}^3$ with either the Dirichlet or Neumann boundary condition $b$, where $\sigma$ is a scalar conductivity function. The inverse problem is formulated as recovering $\sigma$ from vector fields $J_b = \sigma\nabla w_b$ with different boundary conditions $b$. 
 We provide a local Lipschitz stability, stating that near some known $\sigma_0$ and under some regularity assumptions, we can find $b_1$ and $b_2$ by constructing complex geometrical optics (CGO) solutions such that $\|\ln\sigma^{(1)} - \ln\sigma^{(2)}\|_{C^{m,\alpha}(\overline{\Omega})} \le C\sum_{j=1}^2 \|J_{b_j}^{(1)} - J_{b_j}^{(2)}\|_{C^{m,\alpha}(\overline{\Omega})}$. 
 Furthermore, we modify the CGO solutions using the reflection method to make $b_1$ and $b_2$ vanish on a portion of a plane, and prove a local Lipschitz stability with partial data.