Hull.py 文件源码

def necessaryDriftAngle(self, U, Fy, scale=1):
        alpha_test = np.linspace(0, np.max(self.alpha), 5)

        Fx_test, Fy_test, Mz_test = self.driftInducedForces(U, alpha_test, scale=scale)

        alpha_spline = interpolate.splrep(Fy_test, alpha_test)

        return interpolate.splev(Fy, alpha_spline)