sf_double_complex sf_dcsub(sf_double_complex a, sf_double_complex b) /*< complex subtraction >*/ { ... }