WIP : point spread + cross talk functions

WIP needs more work
use linalg from scipy (not numpy)