fp16_sub
delta = fp16_div(fp16_sub(fp16_mul(y, y), x), y * 2);
y = fp16_sub(y, delta);
m = fp16_sub(m, fp16_cos(d / 4));