Julia implementation of sfs