Robust calculation of 3D transonic potential flow based on the non-linear FAS multi-grid method and incomplete LU decomposition