Rosenfeld, loop starts with j = 0