Skip to content

simplenet.io.matpower

Parser for MATPOWER .m case files.

Handles the inline assignment style used in test_9bus_case.m and case9.m files - mpc.bus = [ ... ]; matrix literals plus the single-line mpc.baseMVA = 100; scalar assignment. Lines whose RHS calls readmatrix / readcell are skipped silently so case files such as case_ACTIVSg10kCopy2.m (which sources data from matlab2.xlsx) can still be parsed - the caller is then expected to populate those fields via the xlsx loader.

load_m 

load_m(path: str | Path) -> PowerCase

Load a MATPOWER .m case file into a :class:PowerCase.

Lines that source data from external files (readmatrix, readcell, etc.) are ignored and leave the corresponding field empty so the caller can fill it in (typically via the xlsx loader).

Source code in src/simplenet/io/matpower.py 
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
def load_m(path: str | Path) -> PowerCase:
    """Load a MATPOWER ``.m`` case file into a :class:`PowerCase`.

    Lines that source data from external files (``readmatrix``,
    ``readcell``, etc.) are ignored and leave the corresponding field
    empty so the caller can fill it in (typically via the xlsx
    loader).
    """

    text = Path(path).read_text(encoding="utf-8", errors="replace")
    case = PowerCase()

    pos = 0
    n = len(text)
    while pos < n:
        next_nl = text.find("\n", pos)
        if next_nl == -1:
            next_nl = n
        raw_line = text[pos:next_nl]
        pos = next_nl + 1
        line = _strip_comment(raw_line)
        m = _ASSIGN_RE.match(line)
        if not m:
            continue
        field_name = m.group("field")
        rhs = m.group("rhs").strip()
        if not rhs:
            continue
        if rhs.endswith(";"):
            rhs = rhs[:-1].rstrip()

        if rhs.startswith("[") and not rhs.endswith("]"):
            buf = StringIO()
            buf.write(rhs + "\n")
            depth = rhs.count("[") - rhs.count("]")
            while depth > 0 and pos < n:
                nxt_nl = text.find("\n", pos)
                if nxt_nl == -1:
                    nxt_nl = n
                segment = _strip_comment(text[pos:nxt_nl])
                pos = nxt_nl + 1
                buf.write(segment + "\n")
                depth += segment.count("[") - segment.count("]")
            rhs = buf.getvalue().rstrip().rstrip(";").rstrip()
        elif rhs.startswith("{") and not rhs.endswith("}"):
            buf = StringIO()
            buf.write(rhs + "\n")
            depth = rhs.count("{") - rhs.count("}")
            while depth > 0 and pos < n:
                nxt_nl = text.find("\n", pos)
                if nxt_nl == -1:
                    nxt_nl = n
                segment = _strip_comment(text[pos:nxt_nl])
                pos = nxt_nl + 1
                buf.write(segment + "\n")
                depth += segment.count("{") - segment.count("}")
            rhs = buf.getvalue().rstrip().rstrip(";").rstrip()

        if "readmatrix" in rhs or "readcell" in rhs or "readtable" in rhs:
            continue

        if field_name == "version":
            case.version = rhs.strip().strip("'").strip('"')
        elif field_name == "baseMVA":
            try:
                case.base_mva = float(rhs.strip())
            except ValueError:
                continue
        elif field_name == "bus":
            case.bus = _parse_matrix(rhs)
        elif field_name == "gen":
            case.gen = _parse_matrix(rhs)
        elif field_name == "branch":
            case.branch = _parse_matrix(rhs)
        elif field_name == "gencost":
            case.gencost = _parse_matrix(rhs)
        elif field_name == "dcline":
            case.dcline = _parse_matrix(rhs)
        elif field_name in ("gentype",):
            case.gentype = _parse_cell(rhs)
        elif field_name in ("genfuel",):
            case.genfuel = _parse_cell(rhs)
        elif field_name in ("bus_name",):
            case.bus_name = _parse_cell(rhs)

    return case