Start of example_full_model, with point spread function
Minor adjustments to gibbs_test_full_model (fixed messed-up inversion, made it run faster)