void sf_banded_const_define_reflect (sf_bands slv, float diag /* diagonal */, const float* offd /* off-diagonal [band] */) /*< define matrix with constant diagonal coefficients and reflecting b.c. >*/ { ... }