A Quadratically Convergent Proximal Algorithm For Nonnegative Tensor Decomposition