Methods and systems for iterative computed tomography add an auxiliary variable to the reconstruction process, while retaining all variables in the original formulation, A weighting operator or filter can be applied that causes the Hessian with respect to an image of the cost function to be well-conditioned. An auxiliary sinogram variable distinct from both a set of actual image measurements and from the set of projections computed based on an image can be applied to iteratively update during the statistical iterative image reconstruction, with applied conditions that causes the Hessian with respect to an image of the cost function to be well-conditioned.