You'll find examples of using all the above variants. Computing TDEs for observation point/source element pairs is really simple: from cutde.fullspace import disp, strain disp = disp(obs_pts, src_tris ...