@Cache_this(limit=<NUM_LIT:3>, ignore_args=())<EOL><INDENT>def dK2_drdr_diag(self):<DEDENT>
raise NotImplementedError("<STR_LIT>")<EOL>
Second order derivative of K in r_{i,i}. The diagonal entries are always zero, so we do not give it here.