s1 = grdfun(*args, r=1)