void sf_mask_lop(bool adj, bool add, int nx, int ny, float *x, float *y) /*< linear operator >*/ { ... }