Actual source code: davidson.c