这种新镇 - 拉夫森实施对黑人学位的效率暗示着波动？

Question

我已经实施了一种用于计算隐含波动率的算法（如下所示），并正在寻找有关如何改进代码的反馈，无论是速度，组织等。走到任何不良习惯的面前。

请假设最佳编写“ BS_OPTION_PRICE_ANALYTICE”和“ VEGA_BS”的实现。这些用于第一个导数计算。

double newton_raphson_implied_vol(char type, double mkt_price, double S, double K, double T, double sigma0, double r, double q, double tol)
{
    // Initialize variables
    double bs_price_old{bs_option_price_analytical(type, S, K, T, sigma0, r, q)};
    double bs_mkt_diff{bs_price_old - mkt_price}; // To avoid repeated calculations in While condition
    double sigma_new{-bs_price_old/vega_bs(S, K, T, sigma0, r, q) + sigma0}; // First new value (to be updated iteratively later) 

    // Iterate using Newton-Raphson method
    while (abs(bs_mkt_diff/mkt_price) > tol)
    {
        bs_price_old = bs_option_price_analytical(type, S, K, T, sigma_new, r, q);
        bs_mkt_diff = bs_price_old - mkt_price;
        sigma_new = -bs_mkt_diff/vega_bs(S, K, T, sigma_new, r, q) + sigma_new;
    }

    return sigma_new;
}