ltbams 1.0.11__py3-none-any.whl → 1.0.13__py3-none-any.whl

ams/routines/pypower.py

@@ -32,6 +32,8 @@ class DCPF1(RoutineBase):
     - This class does not implement the AMS-style DC power flow formulation.
     - For detailed mathematical formulations and algorithmic details, refer to the
       MATPOWER User's Manual, section on Power Flow.
+
+    .. versionadded:: 1.0.10
     """
 
     def __init__(self, system, config):
@@ -131,6 +133,15 @@ class DCPF1(RoutineBase):
                         unit='p.u.',
                         name='vBus', tex_name=r'v_{Bus}',
                         src='v', model='Bus',)
+        # --- load ---
+        self.pd = RParam(info='active demand',
+                         name='pd', tex_name=r'p_{d}',
+                         model='StaticLoad', src='p0',
+                         unit='p.u.',)
+        self.qd = RParam(info='reactive demand',
+                         name='qd', tex_name=r'q_{d}',
+                         model='StaticLoad', src='q0',
+                         unit='p.u.',)
         # --- gen ---
         self.pg = Var(info='Gen active power',
                       unit='p.u.',
@@ -267,13 +278,13 @@ class DCPF1(RoutineBase):
             return False
 
     def _get_off_constrs(self):
-        pass
+        logger.debug(f"{self.class_name} does not implement _get_off_constrs.")
 
     def _data_check(self, info=True, **kwargs):
-        pass
+        logger.debug(f"{self.class_name} does not implement _data_check.")
 
     def update(self, params=None, build_mats=False, **kwargs):
-        pass
+        logger.debug(f"{self.class_name} does not implement update.")
 
     def enable(self, name):
         raise NotImplementedError
@@ -282,7 +293,7 @@ class DCPF1(RoutineBase):
         raise NotImplementedError
 
     def _post_add_check(self):
-        pass
+        logger.debug(f"{self.class_name} does not implement _post_add_check.")
 
     def addRParam(self,
                   name: str,
@@ -353,6 +364,8 @@ class PFlow1(DCPF1):
       MATPOWER User's Manual, section on Power Flow.
     - Fast-Decoupled (XB version) and Fast-Decoupled (BX version) algorithms are
       not fully supported yet.
+
+    .. versionadded:: 1.0.10
     """
 
     def __init__(self, system, config):
@@ -427,6 +440,8 @@ class DCOPF1(DCPF1):
     - For detailed mathematical formulations and algorithmic details, refer to the
       MATPOWER User's Manual, section on Optimal Power Flow.
     - Algorithms 400, 500, 600, and 700 are not fully supported yet.
+
+    .. versionadded:: 1.0.10
     """
 
     def __init__(self, system, config):
@@ -577,6 +592,8 @@ class ACOPF1(DCOPF1):
     - This class does not implement the AMS-style AC optimal power flow formulation.
     - For detailed mathematical formulations and algorithmic details, refer to the
       MATPOWER User's Manual, section on Optimal Power Flow.
+
+    .. versionadded:: 1.0.10
     """
 
     def __init__(self, system, config):

ams/routines/routine.py

@@ -3,9 +3,9 @@ Module for routine data.
 """
 
 import logging
-import os
-from typing import Optional, Union, Type, Iterable, Dict
+import json
 from collections import OrderedDict
+from typing import Optional, Union, Type, Iterable, Dict
 
 import numpy as np
 
@@ -19,6 +19,8 @@ from ams.core.service import RBaseService, ValueService
 from ams.opt import OModel
 from ams.opt import Param, Var, Constraint, Objective, ExpressionCalc, Expression
 
+from ams.utils.paths import get_export_path
+
 from ams.shared import pd
 
 logger = logging.getLogger(__name__)
@@ -430,6 +432,90 @@ class RoutineBase:
             logger.warning(msg)
             return False
 
+    def load_json(self, path):
+        """
+        Load scheduling results from a json file.
+
+        Parameters
+        ----------
+        path : str
+            Path of the json file to load.
+
+        Returns
+        -------
+        bool
+            True if the loading is successful, False otherwise.
+        """
+        try:
+            with open(path, 'r') as f:
+                data = json.load(f)
+        except Exception as e:
+            logger.error(f"Failed to load JSON file: {e}")
+            return False
+
+        if not self.initialized:
+            self.init()
+
+        # Unpack variables and expressions from JSON
+        for group, group_data in data.items():
+            if not isinstance(group_data, dict):
+                continue
+            for key, values in group_data.items():
+                if key == 'idx':
+                    continue
+                # Find the corresponding variable or expression
+                if key in self.vars:
+                    var = self.vars[key]
+                    # Assign values to the variable
+                    try:
+                        var.v = np.array(values)
+                    except Exception as e:
+                        logger.warning(f"Failed to assign values to var '{key}': {e}")
+                elif key in self.exprs:
+                    continue
+                elif key in self.exprcs:
+                    exprc = self.exprcs[key]
+                    # Assign values to the expression calculation
+                    try:
+                        exprc.v = np.array(values)
+                    except Exception as e:
+                        logger.warning(f"Failed to assign values to exprc '{key}': {e}")
+        logger.info(f"Loaded results from {path}")
+        return True
+
+    def export_json(self, path=None):
+        """
+        Export scheduling results to a json file.
+
+        Parameters
+        ----------
+        path : str, optional
+            Path of the json file to export.
+
+        Returns
+        -------
+        str
+            The exported json file name
+        """
+        if not self.converged:
+            logger.warning("Routine did not converge, aborting export.")
+            return None
+
+        path, file_name = get_export_path(self.system, self.class_name,
+                                          path=path, fmt='json')
+
+        data_dict = dict()
+
+        group_data(self, data_dict, self.vars, 'v')
+        group_data(self, data_dict, self.exprs, 'v')
+        group_data(self, data_dict, self.exprcs, 'v')
+
+        with open(path, 'w') as f:
+            json.dump(data_dict, f, indent=4,
+                      default=lambda x: x.tolist() if isinstance(x, np.ndarray) else x)
+
+        return file_name
+
     def export_csv(self, path=None):
         """
         Export scheduling results to a csv file.
@@ -441,26 +527,20 @@ class RoutineBase:
 
         Parameters
         ----------
-        path : str
-            path of the csv file to save
+        path : str, optional
+            Path of the csv file to export.
 
         Returns
         -------
-        export_path
-            The path of the exported csv file
+        str
+            The exported csv file name
         """
         if not self.converged:
             logger.warning("Routine did not converge, aborting export.")
             return None
 
-        if not path:
-            if self.system.files.fullname is None:
-                logger.info("Input file name not detacted. Using `Untitled`.")
-                file_name = f'Untitled_{self.class_name}'
-            else:
-                file_name = os.path.splitext(self.system.files.fullname)[0]
-                file_name += f'_{self.class_name}'
-            path = os.path.join(os.getcwd(), file_name + '.csv')
+        path, file_name = get_export_path(self.system, self.class_name,
+                                          path=path, fmt='csv')
 
         data_dict = initialize_data_dict(self)
 
@@ -473,7 +553,7 @@ class RoutineBase:
 
         pd.DataFrame(data_dict).to_csv(path, index=False)
 
-        return file_name + '.csv'
+        return file_name
 
     def summary(self, **kwargs):
         """
@@ -1027,3 +1107,35 @@ def collect_data(rtn: RoutineBase, data_dict: Dict, items: Dict, attr: str):
             logger.debug(f"Error with collecting data for '{key}': {e}")
             data_v = [np.nan] * len(idx_v)
         data_dict.update(OrderedDict(zip([f'{key} {dev}' for dev in idx_v], data_v)))
+
+
+def group_data(rtn: RoutineBase, data_dict: Dict, items: Dict, attr: str):
+    """
+    Collect data for export from grouped items.
+
+    Parameters
+    ----------
+    rtn : ams.routines.routine.RoutineBase
+        The routine to collect data from.
+    data_dict : Dict
+        The data dictionary to populate.
+    items : dict
+        Dictionary of items to collect data from.
+    attr : str
+        Attribute to collect data for.
+    """
+    for key, item in items.items():
+        if item.owner is None:
+            continue
+        if item.owner.class_name not in data_dict.keys():
+            idx_v = item.get_all_idxes()
+            data_dict[item.owner.class_name] = dict(idx=idx_v)
+        else:
+            idx_v = data_dict[item.owner.class_name]['idx']
+        try:
+            data_v = rtn.get(src=key, attr=attr, idx=idx_v,
+                             horizon=rtn.timeslot.v if hasattr(rtn, 'timeslot') else None)
+        except Exception as e:
+            logger.warning(f"Error with collecting data for '{key}': {e}")
+            data_v = [np.nan] * item.owner.n
+        data_dict[item.owner.class_name][key] = data_v

ams/shared.py

@@ -4,6 +4,7 @@ Shared constants and delayed imports.
 This module is supplementary to the ``andes.shared`` module.
 """
 import logging
+import sys
 import unittest
 from functools import wraps
 from time import strftime
@@ -24,6 +25,7 @@ ppoption = LazyImport('from pypower.ppoption import ppoption')
 runpf = LazyImport('from pypower.runpf import runpf')
 runopf = LazyImport('from pypower.runopf import runopf')
 opf = LazyImport('from gurobi_optimods import opf')
+tqdm = LazyImport('from tqdm.auto import tqdm')
 
 # --- an empty ANDES system ---
 empty_adsys = adSystem(autogen_stale=False)
@@ -118,7 +120,7 @@ def skip_unittest_without_PYPOWER(f):
         try:
             import pypower  # NOQA
         except ImportError:
-            raise unittest.SkipTest("PYPOWER is not installed.")
+            raise unittest.SkipTest("PYPOWER is not available.")
         return f(*args, **kwargs)
     return wrapper
 
@@ -131,6 +133,32 @@ def skip_unittest_without_gurobi_optimods(f):
         try:
             import gurobi_optimods  # NOQA
         except ImportError:
-            raise unittest.SkipTest("Gurobi is not installed.")
+            raise unittest.SkipTest("gurobi_optimods is not available.")
         return f(*args, **kwargs)
     return wrapper
+
+
+def _init_pbar(total, unit, no_tqdm):
+    """Initializes and returns a tqdm progress bar."""
+    pbar = tqdm(total=total, unit=unit, ncols=80, ascii=True,
+                file=sys.stdout, disable=no_tqdm)
+    pbar.update(0)
+    return pbar
+
+
+def _update_pbar(pbar, current, total):
+    """Updates and closes the progress bar."""
+    perc = np.round(min((current / total) * 100, 100), 2)
+    if pbar.total is not None:  # Check if pbar is still valid
+        last_pc = pbar.n / pbar.total * 100  # Get current percentage based on updated value
+    else:
+        last_pc = 0
+
+    perc_diff = perc - last_pc
+    if perc_diff >= 1:
+        pbar.update(perc_diff)
+
+    # Ensure pbar finishes at 100% and closes
+    if pbar.n < pbar.total:  # Check if it's not already at total
+        pbar.update(pbar.total - pbar.n)  # Update remaining
+    pbar.close()

ams/system.py

@@ -29,7 +29,7 @@ from ams.report import Report
 from ams.shared import ad_dyn_models
 
 from ams.io.matpower import system2mpc
-from ams.io.matpower import write as wrtite_m
+from ams.io.matpower import write as write_m
 from ams.io.xlsx import write as write_xlsx
 from ams.io.json import write as write_json
 from ams.io.psse import write_raw
@@ -602,6 +602,7 @@ class System(adSystem):
                  **kwargs):
         """
         Convert the AMS system to an ANDES system.
+        Wrapper method for `ams.interface.to_andes`.
 
         A preferred dynamic system file to be added has following features:
         1. The file contains both power flow and dynamic models.
@@ -718,17 +719,36 @@ class System(adSystem):
     def to_mpc(self):
         """
         Export an AMS system to a MATPOWER dict.
+        Wrapper method for `ams.io.matpower.system2mpc`.
 
         Returns
         -------
         dict
             A dictionary representing the MATPOWER case.
+
+        Notes
+        -----
+        - In the `gen` section, slack generators are listed before PV generators.
+        - For uncontrolled generators (`ctrl.v == 0`), their max and min power
+          limits are set to their initial power (`p0.v`) in the converted MPC.
+        - In the converted MPC, the indices of area (`bus[:, 6]`) and zone (`bus[:, 10]`)
+          may differ from the original MPC. However, the mapping relationship is preserved.
+          For example, if the original MPC numbers areas starting from 1, the converted
+          MPC may number them starting from 0.
+        - The coefficients `c2` and `c1` in the generator cost data are scaled by
+          `baseMVA`.
+        - Unlike the XLSX and JSON converters, this implementation uses value providers
+          (`v`) instead of vin. As a result, any changes made through `model.set` will be
+          reflected in the generated MPC.
+
+        .. versionadded:: 1.0.10
         """
         return system2mpc(self)
 
     def to_m(self, outfile: str, overwrite: bool = None):
         """
         Export an AMS system to a MATPOWER M-file.
+        Wrapper method for `ams.io.matpower.write`.
 
         Parameters
         ----------
@@ -736,12 +756,30 @@ class System(adSystem):
             The output file name.
         overwrite : bool, optional
             If True, overwrite the existing file. Default is None.
-        """
-        return wrtite_m(self, outfile=outfile, overwrite=overwrite)
+
+        Notes
+        -----
+        - In the `gen` section, slack generators are listed before PV generators.
+        - For uncontrolled generators (`ctrl.v == 0`), their max and min power
+          limits are set to their initial power (`p0.v`) in the converted MPC.
+        - In the converted MPC, the indices of area (`bus[:, 6]`) and zone (`bus[:, 10]`)
+          may differ from the original MPC. However, the mapping relationship is preserved.
+          For example, if the original MPC numbers areas starting from 1, the converted
+          MPC may number them starting from 0.
+        - The coefficients `c2` and `c1` in the generator cost data are scaled by
+          `baseMVA`.
+        - Unlike the XLSX and JSON converters, this implementation uses value providers
+          (`v`) instead of vin. As a result, any changes made through `model.set` will be
+          reflected in the generated MPC.
+
+        .. versionadded:: 1.0.10
+        """
+        return write_m(self, outfile=outfile, overwrite=overwrite)
 
     def to_xlsx(self, outfile: str, overwrite: bool = None):
         """
         Export an AMS system to an Excel file.
+        Wrapper method for `ams.io.xlsx.write`.
 
         Parameters
         ----------
@@ -749,12 +787,15 @@ class System(adSystem):
             The output file name.
         overwrite : bool, optional
             If True, overwrite the existing file. Default is None.
+
+        .. versionadded:: 1.0.10
         """
         return write_xlsx(self, outfile=outfile, overwrite=overwrite)
 
     def to_json(self, outfile: str, overwrite: bool = None):
         """
         Export an AMS system to a JSON file.
+        Wrapper method for `ams.io.json.write`.
 
         Parameters
         ----------
@@ -762,12 +803,17 @@ class System(adSystem):
             The output file name.
         overwrite : bool, optional
             If True, overwrite the existing file. Default is None.
+
+        .. versionadded:: 1.0.10
         """
         return write_json(self, outfile=outfile, overwrite=overwrite)
 
     def to_raw(self, outfile: str, overwrite: bool = None):
         """
         Export an AMS system to a v33 PSS/E RAW file.
+        Wrapper method for `ams.io.psse.write_raw`.
+
+        This method has not been fully benchmarked yet!
 
         Parameters
         ----------
@@ -775,6 +821,8 @@ class System(adSystem):
             The output file name.
         overwrite : bool, optional
             If True, overwrite the existing file. Default is None.
+
+        .. versionadded:: 1.0.10
         """
         return write_raw(self, outfile=outfile, overwrite=overwrite)
 

ams/utils/paths.py

@@ -255,3 +255,67 @@ def confirm_overwrite(outfile, overwrite=None):
         pass
 
     return True
+
+
+def get_export_path(system, fname, path=None, fmt='csv'):
+    """
+    Get the absolute export path and the derived file name for
+    an AMS system export.
+
+    This function is not intended to be used directly by users.
+
+    Parameters
+    ----------
+    system : ams.system.System
+        The AMS system to export. (Mocked in example)
+    fname : str
+        The descriptive file name, e.g., 'PTDF', or 'DCOPF'.
+    path : str, optional
+        The desired output path.
+        - If it's a directory, the file name will be generated.
+        - If it's a full file path (with extension), the filename part will be used
+          as the base, and the `fmt` argument will determine the final extension.
+        - If None, the current working directory will be used, and the filename will
+          be generated.
+    fmt : str, optional
+        The file format to export, e.g., 'csv', 'json'. Default is 'csv'.
+
+    Returns
+    -------
+    tuple
+        (str, str): A tuple containing:
+            - The **absolute export path** (e.g., '/home/user/project/data_Routine.csv').
+            - The **export file name** (e.g., 'data_Routine.csv'), including the format extension.
+    """
+    # Determine the base name from system.files.fullname or default to "Untitled"
+    if system.files.fullname is None:
+        logger.info("Input file name not detected. Using `Untitled`.")
+        base_name_prefix = f'Untitled_{fname}'
+    else:
+        base_name_prefix = os.path.splitext(os.path.basename(system.files.fullname))[0]
+        base_name_prefix += f'_{fname}'
+
+    target_extension = fmt.lower()  # Ensure consistent extension casing
+
+    if path:
+        abs_path = os.path.abspath(path)  # Resolve to absolute path early
+
+        # Check if the provided path is likely intended as a directory
+        if not os.path.splitext(abs_path)[1]:  # No extension implies it's a directory
+            dir_path = abs_path
+            final_file_name = f"{base_name_prefix}.{target_extension}"
+            full_export_path = os.path.join(dir_path, final_file_name)
+        else:
+            # Path includes a filename. Use its directory, and its base name.
+            dir_path = os.path.dirname(abs_path)
+            # Use the provided filename's base, but enforce the target_extension
+            provided_base_filename = os.path.splitext(os.path.basename(abs_path))[0]
+            final_file_name = f"{provided_base_filename}.{target_extension}"
+            full_export_path = os.path.join(dir_path, final_file_name)
+    else:
+        # No path provided, use current working directory
+        dir_path = os.getcwd()
+        final_file_name = f"{base_name_prefix}.{target_extension}"
+        full_export_path = os.path.join(dir_path, final_file_name)
+
+    return full_export_path, final_file_name

docs/source/index.rst

@@ -89,9 +89,10 @@ dynamic simulator ANDES.
     publication.
 
 
-.. [Wang2025] J. Wang et al., "Dynamics-incorporated Modeling Framework for Stability
-       Constrained Scheduling Under High-penetration of Renewable Energy," in IEEE
-       Transactions on Sustainable Energy, doi: 10.1109/TSTE.2025.3528027.
+.. [Wang2025] J. Wang et al., "Dynamics-Incorporated Modeling Framework for Stability
+       Constrained Scheduling Under High-Penetration of Renewable Energy," in IEEE
+       Transactions on Sustainable Energy,  vol. 16, no. 3, pp. 1673-1685, July 2025,
+       doi: 10.1109/TSTE.2025.3528027.
 
 
 .. toctree::

docs/source/release-notes.rst

@@ -9,6 +9,21 @@ The APIs before v3.0.0 are in beta and may change without prior notice.
 v1.0
 ==========
 
+v1.0.13 (2025-08-18)
+----------------------
+
+- Add methods ``export_npz``, ``load_npz``, ``load_csv`` in ``MatProcessor``
+  to save and load matrices from files
+- Add methods ``export_json``, ``load_json`` in ``RoutineBase``
+  to save and load routine results from JSON files
+
+v1.0.12 (2025-05-29)
+----------------------
+
+- Add RParam pd and qd in ``DCPF1`` for easy access to load
+- Bug fix in ``RoutineBase.export_csv`` when path is specified
+- Fix bug in ``io.matpower.system2mpc`` with multiple PQ at one bus
+
 v1.0.11 (2025-05-23)
 ----------------------
 
@@ -18,8 +33,8 @@ v1.0.11 (2025-05-23)
 v1.0.10 (2025-05-23)
 ----------------------
 
-- Add bus type correction in ``system.System.setup()``
-- Revise ``ams.io.psse.read`` to complete model Zone when necessary
+- Add bus type correction in ``system.System.setup``
+- Revise ``io.psse.read`` to complete model Zone when necessary
 - Use numerical Area and Zone idx in MATPOWER and PSSE RAW file conversion
 - Support JSON format addfile when converting to ANDES case
 - Add PSS/E v33 RAW file writer
@@ -34,7 +49,7 @@ v1.0.10 (2025-05-23)
 - Revise ``andes.common.config.Config.update`` to ensure configuration parameters
   are consistently updated in both the object and its internal ``_dict``
 - Remove legacy revised PYPOWER module
-- Remove function ``ams.shared.ppc2df``
+- Remove function ``shared.ppc2df``
 
 v1.0.9 (2025-04-23)
 --------------------
@@ -62,7 +77,7 @@ v1.0.7 (2025-04-14)
 - Extend common parameters in groups ``StaticGen`` and ``StaticLoad`` with ``area``
 - Set case ``pjm5bus_demo.xlsx`` as a all-inclusive case
 - Include module ``MatProcessor`` in the API documentation
-- Improve Line parameters correction in ``System.setup()``
+- Improve Line parameters correction in ``System.setup``
 - Make func ``interface._to_andes_pflow`` public
 - Discard ``sync_adsys`` step in func ``to_andes_pflow`` to fix mistake in
   parameters conversion
@@ -91,7 +106,7 @@ v1.0.4 (2025-04-05)
 v1.0.3 (2025-03-17)
 --------------------
 
-- Bug fix in function ``ams.interface.parse_addfile``, released in v1.0.3a1
+- Bug fix in function ``interface.parse_addfile``, released in v1.0.3a1
 
 v1.0.2 (2025-02-01)
 --------------------
@@ -136,10 +151,10 @@ v0.9.12 (2024-11-23)
 
 - Refactor ``OModel.initialized`` as a property method
 - Add a demo to show using ``Constraint.e`` for debugging
-- Fix ``ams.opt.omodel.Param.evaluate`` when its value is a number
-- Improve ``ams.opt.omodel.ExpressionCalc`` for better performance
-- Refactor module ``ams.opt``
-- Add class ``ams.opt.Expression``
+- Fix ``opt.omodel.Param.evaluate`` when its value is a number
+- Improve ``opt.omodel.ExpressionCalc`` for better performance
+- Refactor module ``opt``
+- Add class ``opt.Expression``
 - Switch from PYPOWER to ANDES in routine ``PFlow``
 - Switch from PYPOWER to regular formulation in routine ``DCPF``
 - Refactor routines ``DCPF`` and ``DCOPF``
@@ -382,7 +397,7 @@ v0.6.5 (2023-06-27)
 -------------------
 
 - Update documentation with auto-generated model and routine reference
-- Add interface with ANDES ``ams.interop.andes``
+- Add interface with ANDES ``interop.andes``
 - Add routine RTED and example of RTED-TDS co-simulation
 - Draft development documentation
 

{ltbams-1.0.11.dist-info → ltbams-1.0.13.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ltbams
-Version: 1.0.11
+Version: 1.0.13
 Summary: Python software for scheduling modeling and co-simulation with dynamics.
 Home-page: https://github.com/CURENT/ams
 Author: Jinning Wang
@@ -23,6 +23,7 @@ Requires-Dist: openpyxl
 Requires-Dist: andes>=1.9.3
 Requires-Dist: pybind11
 Requires-Dist: cvxpy
+Requires-Dist: cffi
 Provides-Extra: dev
 Requires-Dist: pytest; extra == "dev"
 Requires-Dist: pytest-cov; extra == "dev"
@@ -60,6 +61,7 @@ Python Software for Power System Scheduling Modeling and Co-Simulation with Dyna
 | Badges |  |  |
 |---|---|---|
 | Repo | ![Project Status: Active](https://www.repostatus.org/badges/latest/active.svg) | ![Repo Size](https://img.shields.io/github/repo-size/CURENT/ams) |
+| Python | [![Python Versions](https://img.shields.io/badge/Python-3.9%20%7C%203.10%20%7C%203.11%20%7C%203.12-blue)](https://www.python.org/) |
 | Version | [![PyPI Version](https://img.shields.io/pypi/v/ltbams.svg)](https://pypi.org/project/ltbams/) | [![Anaconda-Server Badge](https://anaconda.org/conda-forge/ltbams/badges/version.svg)](https://anaconda.org/conda-forge/ltbams) |
 | Tag | [![GitHub Tag](https://img.shields.io/github/v/tag/CURENT/ams)](https://github.com/CURENT/ams/tags) | ![GitHub commits since latest release (branch)](https://img.shields.io/github/commits-since/curent/ams/latest/develop) |
 | Documentation | [![Read the Docs](https://img.shields.io/readthedocs/ams?label=stable)](https://ltb.readthedocs.io/projects/ams/en/stable/?badge=stable) | [![Read the Docs](https://img.shields.io/readthedocs/ams?label=develop)](https://ltb.readthedocs.io/projects/ams/en/develop/?badge=develop) |
@@ -217,7 +219,7 @@ sa
 # Citing AMS
 If you use AMS for research or consulting, please cite the following paper in your publication that uses AMS:
 
-> J. Wang et al., "Dynamics-incorporated Modeling Framework for Stability Constrained Scheduling Under High-penetration of Renewable Energy," in IEEE Transactions on Sustainable Energy, doi: 10.1109/TSTE.2025.3528027.
+> J. Wang et al., "Dynamics-Incorporated Modeling Framework for Stability Constrained Scheduling Under High-Penetration of Renewable Energy," in IEEE Transactions on Sustainable Energy, vol. 16, no. 3, pp. 1673-1685, July 2025, doi: 10.1109/TSTE.2025.3528027.
 
 # Sponsors and Contributors
 AMS is the scheduling simulation engine for the CURENT Largescale Testbed (LTB).