Added an operator for simultaneous computation of f_w and df_w/ds to save time