  We present the extension of our previous results for three-loop Yukawa
coupling beta-functions to the case of complex Yukawa matrices describing the
flavour structure of the SM. The calculation is carried out in the context of
unbroken phase of the SM with the help of the MINCER program in a general
linear gauge and cross-checked by means of MATAD/BAMBA codes. In addition,
ambiguities in Yukawa matrix beta-functions are studied.
