void sf_cdstep_diag(int nx, float *res /* [nx] */) /*< compute diagonal of the model resolution matrix >*/ { ... }