void sf_igrad1_lop(bool adj, bool add, int nx, int ny, float *xx, float *yy) /*< linear operator >*/ { ... }