WIP: Simulation.

Now, adding to the simulation a "structural" component of perfusion ('func.nii.gz') that comes from surface curvature.  The idea is to show that the patch-based decomposition can separate the structural component of the functional image from the ``purely functional'' part.