Yes, judging from a quick look at the code (and not being familiar with the GiNaC library used) this seems to square the wave function instead of squaring the absolute value.

[1] https://github.com/c0nrad/hydrogen/blob/master/hydrogen.cpp#...