The structure of the proton and neutron, parameterized by moments of generalized parton distribution functions (GPDs), can be accessed from first principle through the computation of baryon three-point functions with lattice QCD. The numerical effort involved in such computations is quite high and thus an efficient algorithm that extracts most information at given cost is highly desirable.