sf_double_complex sf_dcmul(sf_double_complex a, sf_double_complex b) /*< complex multiplication >*/ { ... }