void numdif(
	void *instance,
	conjugrad_evaluate_t evaluate,
	conjugrad_float_t *x,
	int b,
	int k,
	int a,
	int j,
	conjugrad_float_t epsilon,
	int n
);