void sf_tridiagonal_define (sf_tris slv /* solver object */, float* diag /* diagonal */, float* offd /* off-diagonal */) /*< fill the matrix >*/ { ... }