Compute F-matrices taking into account per contact source diameters and sources thicknesses

[espenhgn]

In the current implementation inverse methods will fail if contacts are not equally spaced. Furthermore, constant radii (and thickness) of sources are implicitly assumed