_assert_not_implemented(i != 0, 'link(x,logsf,phat,i) where i=0 is '