float sf_cell_update1 (int dim, float s, float v, float *p, const float *g) /*< symplectic second-order: step 1 >*/ { ... }