add test_psf() for optimization