The fibo_log routine is doing exponentiation by squaring as you suggest, but it isn't explained at all in the article. It just says it "uses interesting algorithmic strategies and is supposed to be O(log n)."