void sf_tridiagonal_solve (sf_tris slv /* solver object */, float* b /* in - right-hand side, out - solution */) /*< invert the matrix >*/ { int k; ... }