fixed_point_div_snippet.c

op1 = (int32_t)(f_op1*(double)scaler); // op1 in Qn format
op2 = (int32_t)(f_op2*(double)scaler); // op2 in Qn format

// Convert the dividend to Q(2*n) format
tmp = ((uint64_t)op1 << Qn);
// Round by summing half the divisor
tmp += op2 >> 1;
// Division
tmp /= (uint64_t)(op2);