Control algorithms for power-shaping calculations and critical rod bank position searches are described. The algorithms allow a flux computation module using a finite difference solution of the few-group neutron diffusion equations with coarse-mesh rebalancing to be directly extended to control calculations. The control algorithms are formulated as an extension of the standard critical buckling search by recasting the procedure for calculating the eigenvalue. A comparison of the power-shaping calculation with an equivalent keff calculation demonstrates that the former requires approximately the same number of outer iterations as the latter.