drl = gradfun(phatv[i_notfixed])