voly 0.0.127__py3-none-any.whl → 0.0.128__py3-none-any.whl

voly/client.py

@@ -368,40 +368,3 @@ class VolyClient:
             'x_surface': x_surface,
             'moments': moments
         }
-
-    @staticmethod
-    def get_rv_surface(model_results: pd.DataFrame,
-                       pdf_surface: Dict[str, np.ndarray],
-                       x_surface: Dict[str, np.ndarray],
-                       domain_params: Tuple[float, float, int] = (-1.5, 1.5, 1000),
-                       return_domain: str = 'log_moneyness') -> Dict[str, Any]:
-        """
-        Transform historical density surface into a volatility surface.
-
-        Parameters:
-        - model_results: DataFrame from fit_model() with maturity information
-        - pdf_surface: Dictionary mapping maturity names to historical density arrays
-        - x_surface: Dictionary mapping maturity names to x domain arrays
-        - domain_params: Tuple of (min, max, num_points) for the log-moneyness grid
-        - return_domain: Domain for x-axis values ('log_moneyness', 'moneyness', 'returns', 'strikes', 'delta')
-
-        Returns:
-        - Dictionary with implied volatility surface information
-        """
-
-        logger.info("Calculating realized volatility surface")
-
-        # Generate the surface
-        fit_results, rv_surface, x_surface = get_rv_surface(
-            model_results=model_results,
-            pdf_surface=pdf_surface,
-            x_surface=x_surface,
-            domain_params=domain_params,
-            return_domain=return_domain
-        )
-
-        return {
-            'fit_results': fit_results,
-            'rv_surface': rv_surface,
-            'x_surface': x_surface
-        }

voly/core/hd.py

@@ -205,169 +205,3 @@ def get_hd_surface(model_results: pd.DataFrame,
     moments = pd.DataFrame(all_moments).T
 
     return pdf_surface, cdf_surface, x_surface, moments
-
-
-@catch_exception
-def get_rv_surface(model_results: pd.DataFrame,
-                   pdf_surface: Dict[str, np.ndarray],
-                   x_surface: Dict[str, np.ndarray],
-                   domain_params: Tuple[float, float, int] = (-1.5, 1.5, 1000),
-                   return_domain: str = 'log_moneyness') -> Tuple[
-    Dict[str, np.ndarray], Dict[str, np.ndarray], pd.DataFrame]:
-
-    # Check if required columns are present
-    required_columns = ['s', 't', 'r']
-    missing_columns = [col for col in required_columns if col not in model_results.columns]
-    if missing_columns:
-        raise VolyError(f"Required columns missing in model_results: {missing_columns}")
-
-    rv_surface = {}
-    new_x_surface = {}
-    all_params = {}
-
-    # Check if pdf_surface is empty
-    if not pdf_surface:
-        logger.warning("Historical density surface is empty.")
-        return {}, {}, pd.DataFrame()
-
-    # Process each maturity
-    for i in model_results.index:
-        if i not in pdf_surface:
-            logger.warning(f"No historical density available for maturity {i}, skipping.")
-            continue
-
-        # Get parameters for this maturity
-        s = model_results.loc[i, 's']
-        r = model_results.loc[i, 'r']
-        t = model_results.loc[i, 't']
-
-        # Get historical density for this maturity
-        pdf = pdf_surface[i]
-        x = x_surface[i]
-
-        # Calculate x_domain grids
-        LM = get_domain(domain_params, s, r, None, t, 'log_moneyness')
-        M = get_domain(domain_params, s, r, None, t, 'moneyness')
-        R = get_domain(domain_params, s, r, None, t, 'returns')
-        K = get_domain(domain_params, s, r, None, t, 'log_moneyness')
-
-        # Recover call prices from the PDF
-        c_recovered = np.zeros_like(LM)
-        for j, lm_k in enumerate(LM):
-            mask = LM >= lm_k
-            if np.any(mask):
-                integrand = s * (np.exp(LM[mask]) - np.exp(lm_k)) * pdf[mask]
-                c_recovered[j] = np.exp(-r * t) * np.trapz(integrand, LM[mask])
-
-        # Ensure call prices are at least the intrinsic value
-        intrinsic_values = np.maximum(s - K, 0)
-        c_recovered = np.maximum(c_recovered, intrinsic_values)
-
-        # Determine min_lm and max_lm based on days to expiry (DTE)
-        dte = t * 365.25
-        if dte <= 30:
-            min_lm, max_lm = -0.3, 0.3
-        elif dte <= 90:
-            min_lm, max_lm = -0.6, 0.6
-        else:
-            min_lm, max_lm = -0.9, 0.9
-
-        # Generate key log-moneyness points
-        key_lm_points = generate_lm_points(min_lm, max_lm)
-
-        # Find the indices of the key log-moneyness points
-        key_indices = [np.argmin(np.abs(LM - lm)) for lm in key_lm_points]
-        key_lm_actual = LM[key_indices]  # Actual log moneyness values
-        key_strikes = K[key_indices]  # Corresponding strikes
-
-        # Extract call prices at key log-moneyness points
-        key_call_prices = c_recovered[key_indices]
-
-        # Compute IV at key log-moneyness points using our own iv function
-        key_ivs = []
-        for j, idx in enumerate(key_indices):
-            call_price = key_call_prices[j]
-            strike = key_strikes[j]
-            if call_price <= 0:
-                iv_value = 0.01  # Minimum IV of 1%
-            else:
-                try:
-                    iv_value = iv(option_price=call_price, s=s, K=strike, r=r, t=t, option_type='call')
-                    iv_value = max(0.01, min(iv_value, 3.0))  # Clamp between 1% and 300%
-                except Exception as e:
-                    logger.warning(f"IV calculation failed for strike {strike}: {str(e)}")
-                    iv_value = 0.01  # Fallback to 1%
-            key_ivs.append(iv_value)
-        key_ivs = np.array(key_ivs)
-
-        # Create a synthetic option chain for SVI fitting
-        # Convert to DataFrame columns that fit_model expects
-        synthetic_chain = pd.DataFrame({
-            'maturity_name': [i] * len(key_strikes),
-            'maturity_date': pd.Timestamp.now() + pd.Timedelta(days=int(t * 365.25)),
-            'index_price': s,
-            'underlying_price': s,
-            'strike': key_strikes,
-            'log_moneyness': key_lm_actual,
-            'mark_iv': key_ivs,
-            't': t,
-            'r': r,
-            'option_type': 'call'
-        })
-
-        # Fit the SVI model to the recovered IVs
-        fit_results_rv = fit_model(option_chain=synthetic_chain, model_name='svi')
-
-        # Get the parameters for this maturity
-        a = fit_results_rv.loc[i, 'a']
-        b = fit_results_rv.loc[i, 'b']
-        sigma = fit_results_rv.loc[i, 'sigma']
-        rho = fit_results_rv.loc[i, 'rho']
-        m = fit_results_rv.loc[i, 'm']
-
-        # Store the parameters
-        params = {
-            's': s,
-            'r': r,
-            't': t,
-            'a': a,
-            'b': b,
-            'sigma': sigma,
-            'rho': rho,
-            'm': m
-        }
-
-        nu, psi, p, c, nu_tilde = SVIModel.raw_to_jw_params(params['a'], params['b'], params['sigma'], params['rho'], params['m'], params['t'])
-        params.update({
-            'nu': nu,
-            'psi': psi,
-            'p': p,
-            'c': c,
-            'nu_tilde': nu_tilde,
-        })
-        all_params[i] = params
-
-        # Calculate implied volatility using SVI model
-        w = np.array([SVIModel.svi(lm, a, b, sigma, rho, m) for lm in LM])
-        o_recovered = np.sqrt(w / t)
-
-        # Store results
-        rv_surface[i] = o_recovered
-
-        if return_domain == 'log_moneyness':
-            x = LM
-        elif return_domain == 'moneyness':
-            x = M
-        elif return_domain == 'returns':
-            x = R
-        elif return_domain == 'strikes':
-            x = K
-        elif return_domain == 'delta':
-            x = get_domain(domain_params, s, r, o_recovered, t, 'delta')
-
-        new_x_surface[i] = x
-
-    # Create a DataFrame with parameters
-    fit_results = pd.DataFrame(all_params).T
-    x_surface = new_x_surface
-    return fit_results, rv_surface, x_surface

{voly-0.0.127.dist-info → voly-0.0.128.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: voly
-Version: 0.0.127
+Version: 0.0.128
 Summary: Options & volatility research package
 Author-email: Manu de Cara <manu.de.cara@gmail.com>
 License: MIT

{voly-0.0.127.dist-info → voly-0.0.128.dist-info}/RECORD

@@ -1,5 +1,5 @@
 voly/__init__.py,sha256=8xyDk7rFCn_MOD5hxuv5cxxKZvBVRiSIM7TgaMPpwpw,211
-voly/client.py,sha256=R5QSKSJhwYYoJe9LFPL9dynOVwz0q6LI8TU--vnsILs,14887
+voly/client.py,sha256=CJjQzes3K6ZdzPV1JFBoXQ35_EDz05MyP26YQKJe9dc,13402
 voly/exceptions.py,sha256=PBsbn1vNMvKcCJwwJ4lBO6glD85jo1h2qiEmD7ArAjs,92
 voly/formulas.py,sha256=HejPfgVh6-hmWDhvwbHgAIfQMt8iPbIiaKqcBrB5Spw,10766
 voly/models.py,sha256=BF-O7BjGf0BLMpw4rCtfwW7s8_f4iyUZdUY6q1dVxLs,3363
@@ -7,13 +7,13 @@ voly/core/__init__.py,sha256=bu6fS2I1Pj9fPPnl-zY3L7NqrZSY5Zy6NY2uMUvdhKs,183
 voly/core/charts.py,sha256=E21OZB5lTY4YL2flgaFJ6s5g3_ExtAQT2zryZZxLPyM,12735
 voly/core/data.py,sha256=pDeuYhP0GX4RbtlqByvsE3rfHcIkix0BU5MLW8sKIeI,8935
 voly/core/fit.py,sha256=Tb9eeG7e_2dQTcqt6aqEwFrZdy6jR9rSNqe6tzOdVhQ,9245
-voly/core/hd.py,sha256=nU_sjOn-DVRFNmjUiZfeMf-RJehm-I0HrqfOmupINiw,13459
+voly/core/hd.py,sha256=spxNmqw_L3xym-CcPQ_GJEpIFibsuP9bPGLZj0YhDtw,7423
 voly/core/interpolate.py,sha256=JkK172-FXyhesW3hY4pEeuJWG3Bugq7QZXbeKoRpLuo,5305
 voly/core/rnd.py,sha256=0VE77lxesx_BPAO46QwKpcauZNaHnPTiDhmRbSURn3c,10022
 voly/utils/__init__.py,sha256=E05mWatyC-PDOsCxQV1p5Xi1IgpOomxrNURyCx_gB-w,200
 voly/utils/logger.py,sha256=4-_2bVJmq17Q0d7Rd2mPg1AeR8gxv6EPvcmBDMFWcSM,1744
-voly-0.0.127.dist-info/licenses/LICENSE,sha256=wcHIVbE12jfcBOai_wqBKY6xvNQU5E909xL1zZNq_2Q,1065
-voly-0.0.127.dist-info/METADATA,sha256=8s8R7TsvT7fLIhVumViAOTfJ3ofzCJshQBc2QrFOx3k,4115
-voly-0.0.127.dist-info/WHEEL,sha256=1tXe9gY0PYatrMPMDd6jXqjfpz_B-Wqm32CPfRC58XU,91
-voly-0.0.127.dist-info/top_level.txt,sha256=ZfLw2sSxF-LrKAkgGjOmeTcw6_gD-30zvtdEY5W4B7c,5
-voly-0.0.127.dist-info/RECORD,,
+voly-0.0.128.dist-info/licenses/LICENSE,sha256=wcHIVbE12jfcBOai_wqBKY6xvNQU5E909xL1zZNq_2Q,1065
+voly-0.0.128.dist-info/METADATA,sha256=izR6T5uT6uB0FAEb5OvLB9q1uva2THo3kYt5q-O_7wQ,4115
+voly-0.0.128.dist-info/WHEEL,sha256=DK49LOLCYiurdXXOXwGJm6U4DkHkg4lcxjhqwRa0CP4,91
+voly-0.0.128.dist-info/top_level.txt,sha256=ZfLw2sSxF-LrKAkgGjOmeTcw6_gD-30zvtdEY5W4B7c,5
+voly-0.0.128.dist-info/RECORD,,

{voly-0.0.127.dist-info → voly-0.0.128.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (77.0.3)
+Generator: setuptools (78.0.2)
 Root-Is-Purelib: true
 Tag: py3-none-any