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