Interpolation kernel
The next line constructs the choice of interpolation kernel for this simulation:
79#-------------------------------------------------------------------------------
80# Create our interpolation kernel
81#-------------------------------------------------------------------------------
82WT = TableKernel(WendlandC4Kernel())
83output("WT")
The interpolation kernel is a fundamental piece of the numerics for an SPH type simulation, and the choice of this kernel is intimately tied with the choice of how many points per smoothing scale (the parameter nperh in this example script) should be aimed for during smoothing scale evolution (see interpolation_kernels and smoothing_scales for further discussion). That said, the combination of a Wendland kernel with an nperh=4.01 as used here is a good starting point for most SPH calculations in Spheral.