void sf_banded_define (sf_bands slv, float* diag /* diagonal [n] */, float** offd /* off-diagonal [band][n] */) /*< define the matrix >*/ { ... }