More confusion.  I'm looking at the Kaiser & Squires paper, and their eqn. 2.1.12 has a minus sign in front that isn't in this line of code.