2D regularisation

• Seems like the default methods to determine the regularisation parameter are not optimal
• Improve regularisation on edges of the phase space
• Bootstrap?