FIX: for derivatives of force coefficients with array inputs (scalar inputs worked previously, but the shapes were not broadcast properly for array inputs if output was in coefficient form)