GridCalEngine 5.3.19__tar.gz → 5.3.21__tar.gz

gridcalengine-5.3.21/GridCalEngine/Compilers/circuit_to_gslv.py

@@ -0,0 +1,2436 @@
+# This Source Code Form is subject to the terms of the Mozilla Public
+# License, v. 2.0. If a copy of the MPL was not distributed with this
+# file, You can obtain one at https://mozilla.org/MPL/2.0/.
+# SPDX-License-Identifier: MPL-2.0
+from __future__ import annotations
+import os.path
+import warnings
+import numpy as np
+from typing import List, Dict, Union, Tuple, TYPE_CHECKING
+
+import GridCalEngine
+from GridCalEngine import TapModuleControl, TapPhaseControl
+from GridCalEngine.basic_structures import IntVec, Vec
+from GridCalEngine.Devices.profile import Profile
+from GridCalEngine.Devices.multi_circuit import MultiCircuit
+from GridCalEngine.enumerations import (HvdcControlType, SolverType, TimeGrouping,
+                                        ZonalGrouping, MIPSolvers, ContingencyMethod, ContingencyOperationTypes,
+                                        BuildStatus, BranchGroupTypes)
+import GridCalEngine.Devices as dev
+from GridCalEngine.Simulations.PowerFlow.power_flow_options import PowerFlowOptions
+from GridCalEngine.Simulations.PowerFlow.power_flow_results import PowerFlowResults
+
+from GridCalEngine.DataStructures.numerical_circuit import NumericalCircuit
+
+from GridCalEngine.IO.file_system import get_create_gridcal_folder
+from GridCalEngine.basic_structures import ConvergenceReport
+
+if TYPE_CHECKING:  # Only imports the below statements during type checking
+    from GridCalEngine.Simulations.OPF.opf_results import OptimalPowerFlowResults
+    from GridCalEngine.Simulations.OPF.opf_options import OptimalPowerFlowOptions
+    from GridCalEngine.Simulations.LinearFactors.linear_analysis_options import LinearAnalysisOptions
+    from GridCalEngine.Simulations.ContingencyAnalysis.contingency_analysis_options import ContingencyAnalysisOptions
+    from GridCalEngine.Simulations.ContingencyAnalysis.contingency_analysis_results import ContingencyAnalysisResults
+
+GSLV_RECOMMENDED_VERSION = "0.0.4"
+GSLV_VERSION = ''
+GSLV_AVAILABLE = False
+try:
+    import pygslv as pg
+
+    pg.activate(os.path.join(get_create_gridcal_folder(), "license.gslv"), verbose=True)
+
+    # activate
+    if not pg.isLicensed():
+        # license not found
+        GSLV_AVAILABLE = False
+    else:
+        # already activated
+        GSLV_AVAILABLE = True
+        GSLV_VERSION = pg.get_version()
+
+    if GSLV_AVAILABLE:
+        if GSLV_VERSION < GSLV_RECOMMENDED_VERSION:
+            warnings.warn(f"Recommended version for GSLV is {GSLV_RECOMMENDED_VERSION} "
+                          f"instead of {GSLV_VERSION}")
+
+    build_status_dict = {
+        BuildStatus.Planned: pg.BuildStatus.Planned,
+        BuildStatus.Commissioned: pg.BuildStatus.Commissioned,
+        BuildStatus.Candidate: pg.BuildStatus.Candidate,
+        BuildStatus.Decommissioned: pg.BuildStatus.Decommissioned,
+        BuildStatus.PlannedDecommission: pg.BuildStatus.PlannedDecommission,
+    }
+
+    tap_module_control_mode_dict = {
+        TapModuleControl.fixed: pg.TapModuleControl.fixed,
+        TapModuleControl.Qf: pg.TapModuleControl.Qf,
+        TapModuleControl.Qt: pg.TapModuleControl.Qt,
+        TapModuleControl.Vm: pg.TapModuleControl.Vm,
+    }
+
+    tap_phase_control_mode_dict = {
+        TapPhaseControl.fixed: pg.TapPhaseControl.fixed,
+        TapPhaseControl.Pf: pg.TapPhaseControl.Pf,
+        TapPhaseControl.Pt: pg.TapPhaseControl.Pt,
+    }
+
+    hvdc_control_mode_dict = {
+        HvdcControlType.type_0_free: pg.HvdcControlType.type_0_free,
+        HvdcControlType.type_1_Pset: pg.HvdcControlType.type_1_Pset,
+    }
+
+    group_type_dict = {
+        BranchGroupTypes.GenericGroup: pg.BranchGroupTypes.GenericGroup,
+        BranchGroupTypes.TransformerGroup: pg.BranchGroupTypes.TransformerGroup,
+        BranchGroupTypes.LineSegmentsGroup: pg.BranchGroupTypes.LineSegmentsGroup,
+    }
+
+    contingency_ops_type_dict = {
+        ContingencyOperationTypes.Active: pg.ContingencyOperationTypes.Active,
+        ContingencyOperationTypes.PowerPercentage: pg.ContingencyOperationTypes.PowerPercentage,
+    }
+
+except ImportError as e:
+    pg = None
+    GSLV_AVAILABLE = False
+    GSLV_VERSION = ''
+    build_status_dict = dict()
+    tap_module_control_mode_dict = dict()
+    tap_phase_control_mode_dict = dict()
+    contingency_ops_type_dict = dict()
+
+
+def get_gslv_mip_solvers_list() -> List[str]:
+    """
+    Get list of available MIP solvers
+    :return:
+    """
+    if GSLV_AVAILABLE:
+        return list()
+    else:
+        return list()
+
+
+def convert_tap_module_control_mode_dict(data: Dict[int, TapModuleControl]) -> Dict[int, "pg.TapModuleControl"]:
+    """
+    Function to convert a dictionary of TapModuleControl modes to pg.TapModuleControl modes
+    :param data:
+    :return:
+    """
+    return {i: tap_module_control_mode_dict[val] for i, val in data.items()}
+
+
+def convert_tap_module_control_mode_lst(data: List[TapModuleControl]) -> List["pg.TapModuleControl"]:
+    """
+    Function to convert a list of TapModuleControl modes to pg.TapModuleControl modes
+    :param data:
+    :return:
+    """
+    return [tap_module_control_mode_dict[val] for val in data]
+
+
+def convert_tap_phase_control_mode_dict(data: Dict[int, TapPhaseControl]) -> Dict[int, "pg.TapPhaseControl"]:
+    """
+    Function to convert a dictionary of TapPhaseControl modes to pg.TapPhaseControl modes
+    :param data:
+    :return:
+    """
+    return {i: tap_phase_control_mode_dict[val] for i, val in data.items()}
+
+
+def convert_tap_phase_control_mode_lst(data: List[TapPhaseControl]) -> List["pg.TapPhaseControl"]:
+    """
+    Function to convert a list of TapPhaseControl modes to pg.TapPhaseControl modes
+    :param data:
+    :return:
+    """
+    return [tap_phase_control_mode_dict[val] for val in data]
+
+
+def fill_profile(gslv_profile: "pg.Profiledouble|pg.Profilebool|pg.Profileint|pg.Profileuint",
+                 gc_profile: Profile,
+                 use_time_series: bool,
+                 time_indices: Union[IntVec, None],
+                 n_time: int = 1,
+                 default_val: int | float | bool | TapPhaseControl | TapModuleControl = 0) -> None:
+    """
+    Generates a default time series
+    :param gslv_profile: Profile from gslv to fill in
+    :param gc_profile: Profile from gridcal to convert
+    :param use_time_series: use time series?
+    :param time_indices: time series indices if any (optional)
+    :param n_time: number of time steps
+    :param default_val: Default value
+    """
+
+    if use_time_series:
+        if gc_profile.is_sparse:
+            if time_indices is None:
+
+                if isinstance(default_val, TapPhaseControl):
+                    data = convert_tap_phase_control_mode_dict(data=gc_profile.sparse_array.get_map())
+                elif isinstance(default_val, TapModuleControl):
+                    data = convert_tap_module_control_mode_dict(data=gc_profile.sparse_array.get_map())
+                else:
+                    data = gc_profile.sparse_array.get_map()
+
+                # we pick all the profile
+                if len(data) > 0:
+                    gslv_profile.init_sparse(default_val=gc_profile.default_value, data=data)
+
+            else:
+                assert len(time_indices) == n_time
+
+                # we need a sliced version
+                sp_arr2 = gc_profile.sparse_array.slice(time_indices)
+
+                if isinstance(default_val, TapPhaseControl):
+                    data = convert_tap_phase_control_mode_dict(data=sp_arr2.get_map())
+                elif isinstance(default_val, TapModuleControl):
+                    data = convert_tap_module_control_mode_dict(data=sp_arr2.get_map())
+                else:
+                    data = sp_arr2.get_map()
+
+                gslv_profile.init_sparse(default_val=gc_profile.default_value, data=data)
+
+        else:
+            if time_indices is None:
+                # we pick all the profile
+
+                if isinstance(default_val, TapPhaseControl):
+                    data = convert_tap_phase_control_mode_lst(data=gc_profile.dense_array)
+                elif isinstance(default_val, TapModuleControl):
+                    data = convert_tap_module_control_mode_dict(data=gc_profile.dense_array)
+                else:
+                    data = gc_profile.dense_array
+
+                gslv_profile.init_dense(data)
+
+            else:
+                assert len(time_indices) == n_time
+                # we need a sliced version
+                if isinstance(default_val, TapPhaseControl):
+                    data = convert_tap_phase_control_mode_lst(data=gc_profile.dense_array[time_indices])
+                elif isinstance(default_val, TapModuleControl):
+                    data = convert_tap_module_control_mode_dict(data=gc_profile.dense_array[time_indices])
+                else:
+                    data = gc_profile.dense_array[time_indices]
+
+                gslv_profile.init_dense(data)
+
+    else:
+        if isinstance(default_val, TapPhaseControl):
+            gslv_profile.fill(tap_phase_control_mode_dict[default_val])
+
+        elif isinstance(default_val, TapModuleControl):
+            gslv_profile.fill(tap_module_control_mode_dict[default_val])
+
+        else:
+            gslv_profile.fill(default_val)
+
+
+def fill_profile_with_array(gslv_profile: "pg.Profiledouble",
+                            arr: Vec,
+                            use_time_series: bool,
+                            time_indices: Union[IntVec, None],
+                            n_time=1,
+                            default_val=0) -> None:
+    """
+    Generates a default time series
+    :param gslv_profile: Profile from gslv to fill in
+    :param arr:  array to fill in
+    :param use_time_series: use time series?
+    :param time_indices: time series indices if any (optional)
+    :param n_time: number of time steps
+    :param default_val: Default value
+    """
+
+    if use_time_series:
+        if time_indices is None:
+            # we pick all the profile
+            gslv_profile.init_dense(arr)
+        else:
+            assert len(time_indices) == n_time
+            # we need a sliced version
+            gslv_profile.init_dense(arr[time_indices])
+
+    else:
+        gslv_profile.fill(default_val)
+
+
+def convert_area(area: dev.Area) -> "pg.Area":
+    """
+    
+    :param area:
+    :return:
+    """
+    return pg.Area(idtag=area.idtag, code=str(area.code), name=area.name)
+
+
+def add_areas(circuit: MultiCircuit,
+              gslv_grid: "pg.MultiCircuit") -> Dict[dev.Area, "pg.Area"]:
+    """
+    Add GSLV Areas
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :return: Dictionary [GridCal area] -> GSLV Area
+    """
+    d = dict()
+
+    for i, area in enumerate(circuit.areas):
+        elm = convert_area(area)
+        gslv_grid.add_area(elm)
+        d[area] = elm
+
+    return d
+
+
+def convert_zone(zone: dev.Zone) -> "pg.Zone":
+    """
+
+    :param zone:
+    :return:
+    """
+    return pg.Zone(idtag=zone.idtag, code=str(zone.code), name=zone.name)
+
+
+def add_zones(circuit: MultiCircuit,
+              gslv_grid: "pg.MultiCircuit") -> Dict[dev.Zone, "pg.Zone"]:
+    """
+    Add GSLV Zones
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :return: Dictionary [GridCal Zone] -> GSLV Zone
+    """
+    d = dict()
+
+    for i, zone in enumerate(circuit.zones):
+        elm = convert_zone(zone)
+        gslv_grid.add_zone(elm)
+        d[zone] = elm
+
+    return d
+
+
+def convert_country(country: dev.Country) -> "pg.Country":
+    """
+
+    :param country:
+    :return:
+    """
+    return pg.Country(idtag=country.idtag, code=str(country.code), name=country.name)
+
+
+def add_countries(circuit: MultiCircuit,
+                  gslv_grid: "pg.MultiCircuit") -> Dict[dev.Country, "pg.Country"]:
+    """
+    Add GSLV countries
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :return: Dictionary [GridCal country] -> GSLV country
+    """
+    d = dict()
+
+    for i, country in enumerate(circuit.countries):
+        elm = convert_country(country)
+        gslv_grid.add_country(elm)
+        d[country] = elm
+
+    return d
+
+
+def convert_municipality(country: dev.Municipality) -> "pg.Municipality":
+    """
+
+    :param country:
+    :return:
+    """
+    return pg.Municipality(idtag=country.idtag, code=str(country.code), name=country.name)
+
+
+def add_municipalities(circuit: MultiCircuit,
+                       gslv_grid: "pg.MultiCircuit") -> Dict[dev.Country, "pg.Country"]:
+    """
+    Add GSLV countries
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :return: Dictionary [GridCal country] -> GSLV country
+    """
+    d = dict()
+
+    for i, municipality in enumerate(circuit.municipalities):
+        elm = convert_municipality(municipality)
+        gslv_grid.add_municipality(elm)
+        d[municipality] = elm
+
+    return d
+
+
+def convert_region(country: dev.Municipality) -> "pg.Municipality":
+    """
+
+    :param country:
+    :return:
+    """
+    return pg.Municipality(idtag=country.idtag, code=str(country.code), name=country.name)
+
+
+def add_regions(circuit: MultiCircuit,
+                gslv_grid: "pg.MultiCircuit") -> Dict[dev.Country, "pg.Country"]:
+    """
+    Add GSLV countries
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :return: Dictionary [GridCal country] -> GSLV country
+    """
+    d = dict()
+
+    for i, municipality in enumerate(circuit.regions):
+        elm = convert_region(municipality)
+        gslv_grid.add_region(elm)
+        d[municipality] = elm
+
+    return d
+
+
+def convert_branch_group(country: dev.BranchGroup) -> "pg.BranchGroup":
+    """
+
+    :param country:
+    :return:
+    """
+    return pg.BranchGroup(
+        idtag=country.idtag,
+        code=str(country.code),
+        name=country.name,
+        group_type=group_type_dict[country.group_type]
+    )
+
+
+def add_branch_groups(circuit: MultiCircuit,
+                      gslv_grid: "pg.MultiCircuit") -> Dict[dev.BranchGroup, "pg.BranchGroup"]:
+    """
+    Add GSLV countries
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :return: Dictionary [GridCal country] -> GSLV country
+    """
+    d = dict()
+
+    for i, branch_group in enumerate(circuit.branch_groups):
+        elm = convert_branch_group(branch_group)
+        gslv_grid.add_branch_group(elm)
+        d[branch_group] = elm
+
+    return d
+
+
+def convert_substation(se: dev.Substation, n_time: int) -> "pg.Substation":
+    """
+
+    :param se:
+    :param n_time:
+    :return:
+    """
+    return pg.Substation(
+        nt=n_time,
+        idtag=se.idtag,
+        code=str(se.code),
+        name=se.name
+    )
+
+
+def add_substations(circuit: MultiCircuit,
+                    gslv_grid: "pg.MultiCircuit",
+                    n_time: int) -> Dict[dev.Substation, "pg.Substation"]:
+    """
+    Add GSLV substations
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :param n_time: number of time steps
+    :return: Dictionary [GridCal Zone] -> GSLV Zone
+    """
+    d = dict()
+
+    for i, se in enumerate(circuit.substations):
+        elm = convert_substation(se, n_time=n_time)
+        gslv_grid.add_substation(elm)
+        d[se] = elm
+
+    return d
+
+
+def convert_voltage_level(elm: dev.VoltageLevel,
+                          substations_dict: Dict[dev.Substation, "pg.Substation"]) -> "pg.VoltageLevel":
+    """
+
+    :param elm:
+    :param substations_dict:
+    :return:
+    """
+    return pg.VoltageLevel(
+        idtag=elm.idtag,
+        code=str(elm.code),
+        name=elm.name,
+        Vnom=elm.Vnom,
+        substation=substations_dict.get(elm.substation, None)
+    )
+
+
+def add_voltage_levels(
+        circuit: MultiCircuit,
+        gslv_grid: "pg.MultiCircuit",
+        substations_dict: Dict[dev.Substation, "pg.Substation"]
+) -> Dict[dev.VoltageLevel, "pg.VoltageLevel"]:
+    """
+    Add GSLV substations
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :param substations_dict: substations mapping dictionary
+    :return: Dictionary [GridCal Zone] -> GSLV Zone
+    """
+    d = dict()
+
+    for i, vl in enumerate(circuit.voltage_levels):
+        elm = convert_voltage_level(vl, substations_dict=substations_dict)
+        gslv_grid.add_voltage_level(elm)
+        d[vl] = elm
+
+    return d
+
+
+def convert_contingency_groups(elm: dev.ContingencyGroup) -> "pg.ContingencyGroup":
+    """
+
+    :param elm:
+    :return:
+    """
+    return pg.ContingencyGroup(idtag=elm.idtag,
+                               code=str(elm.code),
+                               name=elm.name,
+                               category=elm.category)
+
+
+def add_contingency_groups(circuit: MultiCircuit,
+                           gslv_grid: "pg.MultiCircuit") -> Dict[dev.ContingencyGroup, "pg.ContingencyGroup"]:
+    """
+    Add GSLV ContingenciesGroup
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :return: Dictionary [GridCal ContingenciesGroup] -> GSLV ContingenciesGroup
+    """
+    d = dict()
+
+    for i, elm in enumerate(circuit.get_contingency_groups()):
+        cg = convert_contingency_groups(elm)
+        gslv_grid.add_contingency_group(cg)
+        d[elm] = cg
+
+    return d
+
+
+def convert_contingencies(elm: dev.Contingency,
+                          n_time: int,
+                          groups_dict: Dict[dev.ContingencyGroup, "pg.ContingencyGroup"]) -> "pg.Contingency":
+    """
+
+    :param elm:
+    :param n_time:
+    :param groups_dict:
+    :return:
+    """
+
+
+
+    return pg.Contingency(idtag=elm.idtag,
+                          device_idtag=elm.device_idtag,
+                          name=elm.name,
+                          code=str(elm.code),
+                          prop=contingency_ops_type_dict[elm.prop],
+                          value=elm.value,
+                          group=groups_dict[elm.group])
+
+
+def add_contingencies(circuit: MultiCircuit,
+                      gslv_grid: "pg.MultiCircuit",
+                      n_time: int,
+                      groups_dict: Dict[dev.ContingencyGroup, "pg.ContingencyGroup"], ):
+    """
+    Add GSLV ContingenciesGroup
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV Circuit
+    :param n_time: number of time steps
+    :param groups_dict: Contingency groups dictionary
+    :return: Dictionary [GridCal ContingenciesGroup] -> GSLV ContingenciesGroup
+    """
+    d = dict()
+
+    for i, elm in enumerate(circuit.contingencies):
+        con = convert_contingencies(elm=elm,
+                                    n_time=n_time,
+                                    groups_dict=groups_dict)
+        gslv_grid.add_contingency(con)
+        d[elm] = con
+
+    return d
+
+
+def convert_investment_group(elm: dev.InvestmentsGroup) -> "pg.InvestmentGroup":
+    """
+
+    :param elm:
+    :return:
+    """
+    return pg.InvestmentGroup(idtag=elm.idtag,
+                              code=str(elm.code),
+                              name=elm.name,
+                              category=elm.category)
+
+
+def add_investment_groups(circuit: MultiCircuit,
+                          gslv_grid: "pg.MultiCircuit") -> Dict[dev.InvestmentsGroup, "pg.InvestmentGroup"]:
+    """
+
+    :param circuit:
+    :param gslv_grid:
+    :return:
+    """
+    d = dict()
+
+    for i, elm in enumerate(circuit.investments_groups):
+        ig = convert_investment_group(elm)
+        gslv_grid.add_investment_group(ig)
+        d[elm] = ig
+
+    return d
+
+
+def convert_investment(
+        elm: dev.Investment,
+        groups_dict: Dict[dev.InvestmentsGroup, "pg.InvestmentGroup"]
+) -> "pg.Investment":
+    """
+
+    :param elm:
+    :param groups_dict:
+    :return:
+    """
+    return pg.Investment(idtag=elm.idtag,
+                         code=str(elm.code),
+                         name=elm.name,
+                         device_idtag=elm.device_idtag,
+                         group=groups_dict[elm.group],
+                         CAPEX=elm.CAPEX,
+                         OPEX=elm.OPEX,
+                         status=elm.status, )
+
+
+def add_investments(circuit: MultiCircuit,
+                    gslv_grid: "pg.MultiCircuit",
+                    groups_dict: Dict[dev.InvestmentsGroup, "pg.InvestmentGroup"]):
+    """
+
+    :param circuit:
+    :param gslv_grid:
+    :param groups_dict:
+    :return:
+    """
+    d = dict()
+
+    for i, elm in enumerate(circuit.investments):
+        investment = convert_investment(elm, groups_dict=groups_dict[elm.group])
+        gslv_grid.add_investment(investment)
+        d[elm] = investment
+
+    return d
+
+
+def convert_facility(elm: dev.Facility) -> "pg.Facility":
+    """
+
+    :param elm:
+    :return:
+    """
+    return pg.Facility(idtag=elm.idtag,
+                       code=str(elm.code),
+                       name=elm.name)
+
+
+def add_facilities(circuit: MultiCircuit,
+                   gslv_grid: "pg.MultiCircuit") -> Dict[dev.Facility, "pg.Facility"]:
+    """
+
+    :param circuit:
+    :param gslv_grid:
+    :return:
+    """
+    d = dict()
+
+    for i, elm in enumerate(circuit.facilities):
+        facility = convert_facility(elm)
+        gslv_grid.add_facility(facility)
+        d[elm] = facility
+
+    return d
+
+
+def convert_modelling_authority(elm: dev.ModellingAuthority) -> "pg.ModellingAuthority":
+    """
+
+    :param elm:
+    :return:
+    """
+    return pg.ModellingAuthority(idtag=elm.idtag,
+                                 code=str(elm.code),
+                                 name=elm.name)
+
+
+def add_modelling_authorities(circuit: MultiCircuit,
+                              gslv_grid: "pg.MultiCircuit") -> Dict[dev.ModellingAuthority, "pg.ModellingAuthority"]:
+    """
+
+    :param circuit:
+    :param gslv_grid:
+    :return:
+    """
+    d = dict()
+
+    for i, elm in enumerate(circuit.modelling_authorities):
+        ma = convert_modelling_authority(elm)
+        gslv_grid.add_modelling_authority(ma)
+        d[elm] = ma
+
+    return d
+
+
+def convert_bus(elm: dev.Bus, n_time: int,
+                area_dict: Dict[dev.Area, "pg.Area"],
+                zone_dict: Dict[dev.Zone, "pg.Zone"],
+                substation_dict: Dict[dev.Substation, "pg.Substation"],
+                voltage_level_dict: Dict[dev.VoltageLevel, "pg.VoltageLevel"],
+                country_dict: Dict[dev.Country, "pg.Country"],
+                time_indices: IntVec,
+                use_time_series: bool) -> "pg.Bus":
+    """
+
+    :param elm:
+    :param n_time:
+    :param area_dict:
+    :param zone_dict:
+    :param substation_dict:
+    :param voltage_level_dict:
+    :param country_dict:
+    :param time_indices:
+    :param use_time_series:
+    :return:
+    """
+    bus = pg.Bus(nt=n_time,
+                 name=elm.name,
+                 idtag=elm.idtag,
+                 code=str(elm.code),
+                 Vnom=elm.Vnom,
+                 vmin=elm.Vmin,
+                 vmax=elm.Vmax,
+                 angle_min=elm.angle_min,
+                 angle_max=elm.angle_max,
+                 r_fault=elm.r_fault,
+                 x_fault=elm.x_fault,
+                 active_default=elm.active,
+
+                 is_slack=elm.is_slack,
+                 is_dc=elm.is_dc,
+                 is_internal=elm.internal,
+
+                 area=area_dict.get(elm.area, None),
+                 zone=zone_dict.get(elm.zone, None),
+                 substation=substation_dict.get(elm.substation, None),
+                 voltage_level=voltage_level_dict.get(elm.substation, None),
+                 country=country_dict.get(elm.country, None),
+                 latitude=elm.latitude,
+                 longitude=elm.longitude,
+                 Vm0=elm.Vm0,
+                 Va0=elm.Va0,
+                 )
+
+    fill_profile(gslv_profile=bus.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    return bus
+
+
+def add_buses(
+        circuit: MultiCircuit,
+        gslv_grid: "pg.MultiCircuit",
+        area_dict: Dict[dev.Area, "pg.Area"],
+        zone_dict: Dict[dev.Zone, "pg.Zone"],
+        substation_dict: Dict[dev.Substation, "pg.Substation"],
+        voltage_level_dict: Dict[dev.VoltageLevel, "pg.VoltageLevel"],
+        country_dict: Dict[dev.Country, "pg.Country"],
+        use_time_series: bool,
+        n_time: int = 1,
+        time_indices: Union[IntVec, None] = None,
+) -> Dict[str, "pg.Bus"]:
+    """
+    Convert the buses to GSLV buses
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param use_time_series: compile the time series from GridCal? otherwise, just the snapshot
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    :param area_dict: Area object translation dictionary
+    :param zone_dict: Zone object translation dictionary
+    :param substation_dict: Substation object translation dictionary
+    :param voltage_level_dict: Voltage level object translation dictionary
+    :param country_dict: Country object translation dictionary
+    :return: bus dictionary buses[uuid] -> Bus
+    """
+
+    if time_indices is not None:
+        assert (len(time_indices) == n_time)
+
+    if area_dict is None:
+        area_dict = {elm: k for k, elm in enumerate(circuit.areas)}
+
+    bus_dict: Dict[str, "pg.Bus"] = dict()
+
+    for i, bus in enumerate(circuit.buses):
+        elm = convert_bus(elm=bus, n_time=n_time,
+                          area_dict=area_dict,
+                          zone_dict=zone_dict,
+                          substation_dict=substation_dict,
+                          voltage_level_dict=voltage_level_dict,
+                          country_dict=country_dict,
+                          use_time_series=use_time_series,
+                          time_indices=time_indices)
+
+        gslv_grid.add_bus(elm)
+        bus_dict[bus.idtag] = elm
+
+    return bus_dict
+
+
+def convert_load(k: int, elm: dev.Load, bus_dict: Dict[str, "pg.Bus"], n_time: int,
+                 use_time_series: bool, time_indices: IntVec | None = None,
+                 opf_results: OptimalPowerFlowResults | None = None) -> "pg.Load":
+    """
+
+    :param k:
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :param opf_results:
+    :return:
+    """
+
+    load = pg.Load(
+        nt=n_time,
+        name=elm.name,
+        idtag=elm.idtag,
+        code=str(elm.code),
+        G=elm.G,
+        B=elm.B,
+        Ir=elm.Ir,
+        Ii=elm.Ii,
+        P=elm.P if opf_results is None else elm.P - opf_results.load_shedding[k],
+        Q=elm.Q,
+        Cost=elm.Cost,
+        active=elm.active,
+        mttf=elm.mttf,
+        mttr=elm.mttr,
+        capex=elm.capex,
+        opex=elm.opex,
+        build_status=build_status_dict[elm.build_status],
+    )
+
+    load.bus = bus_dict[elm.bus.idtag]
+
+    fill_profile(gslv_profile=load.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    if opf_results is None:
+        fill_profile(gslv_profile=load.P,
+                     gc_profile=elm.P_prof,
+                     use_time_series=use_time_series,
+                     time_indices=time_indices,
+                     n_time=n_time,
+                     default_val=elm.P)
+    else:
+        fill_profile_with_array(gslv_profile=load.P,
+                                arr=elm.P_prof.toarray() - opf_results.load_shedding[:, k],
+                                use_time_series=use_time_series,
+                                time_indices=time_indices,
+                                n_time=n_time,
+                                default_val=elm.P)
+
+    fill_profile(gslv_profile=load.Q,
+                 gc_profile=elm.Q_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Q)
+
+    fill_profile(gslv_profile=load.G,
+                 gc_profile=elm.G_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.G)
+
+    fill_profile(gslv_profile=load.B,
+                 gc_profile=elm.B_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.B)
+
+    fill_profile(gslv_profile=load.Ir,
+                 gc_profile=elm.Ir_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Ir)
+
+    fill_profile(gslv_profile=load.Ii,
+                 gc_profile=elm.Ii_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Ii)
+
+    fill_profile(gslv_profile=load.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    return load
+
+
+def add_loads(circuit: MultiCircuit,
+              gslv_grid: "pg.MultiCircuit",
+              bus_dict: Dict[str, "pg.Bus"],
+              use_time_series: bool,
+              n_time=1,
+              time_indices: IntVec | None = None,
+              opf_results: OptimalPowerFlowResults | None = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param use_time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param n_time: number of time steps
+    :param time_indices:
+    :param opf_results:
+    :return:
+    """
+    devices = circuit.get_loads()
+    for k, elm in enumerate(devices):
+        load = convert_load(k=k, elm=elm, bus_dict=bus_dict,
+                            n_time=n_time, use_time_series=use_time_series,
+                            time_indices=time_indices, opf_results=opf_results)
+        gslv_grid.add_load(load)
+
+
+def convert_static_generator(elm: dev.StaticGenerator,
+                             bus_dict: Dict[str, "pg.Bus"],
+                             n_time: int,
+                             use_time_series: bool,
+                             time_indices: IntVec | None = None, ) -> "pg.StaticGenerator":
+    """
+
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :return:
+    """
+
+    pe_inj = pg.StaticGenerator(
+        idtag=elm.idtag,
+        code=str(elm.code),
+        name=elm.name,
+        calc_node=bus_dict[elm.bus.idtag],
+        nt=n_time,
+        P=elm.P,
+        Q=elm.Q,
+        build_status=elm.build_status,
+    )
+
+    fill_profile(gslv_profile=pe_inj.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    fill_profile(gslv_profile=pe_inj.P,
+                 gc_profile=elm.P_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.P)
+
+    fill_profile(gslv_profile=pe_inj.Q,
+                 gc_profile=elm.Q_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Q)
+
+    fill_profile(gslv_profile=pe_inj.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    return pe_inj
+
+
+def add_static_generators(circuit: MultiCircuit, gslv_grid: "pg.MultiCircuit",
+                          bus_dict: Dict[str, "pg.Bus"],
+                          time_series: bool,
+                          n_time=1,
+                          time_indices: Union[IntVec, None] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    """
+    devices = circuit.get_static_generators()
+    for k, elm in enumerate(devices):
+        pe_inj = convert_static_generator(elm=elm, bus_dict=bus_dict, n_time=n_time,
+                                          use_time_series=time_series, time_indices=time_indices)
+        gslv_grid.add_static_generator(pe_inj)
+
+
+def convert_shunt(elm: dev.Shunt, bus_dict: Dict[str, "pg.Bus"], n_time: int,
+                  use_time_series: bool, time_indices: IntVec | None = None, ) -> "pg.Shunt":
+    """
+
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :return:
+    """
+    sh = pg.Shunt(
+        nt=n_time,
+        name=elm.name,
+        idtag=elm.idtag,
+        code=str(elm.code),
+        G=elm.G,
+        B=elm.B,
+        build_status=build_status_dict[elm.build_status],
+    )
+
+    sh.bus = bus_dict[elm.bus.idtag]
+
+    fill_profile(gslv_profile=sh.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    fill_profile(gslv_profile=sh.G,
+                 gc_profile=elm.G_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.G)
+
+    fill_profile(gslv_profile=sh.B,
+                 gc_profile=elm.B_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.B)
+
+    fill_profile(gslv_profile=sh.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    return sh
+
+
+def add_shunts(circuit: MultiCircuit,
+               gslv_grid: "pg.MultiCircuit",
+               bus_dict: Dict[str, "pg.Bus"],
+               time_series: bool,
+               n_time=1,
+               time_indices: Union[IntVec, None] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    """
+    devices = circuit.get_shunts()
+    for k, elm in enumerate(devices):
+        sh = convert_shunt(elm=elm, bus_dict=bus_dict, n_time=n_time,
+                           use_time_series=time_series, time_indices=time_indices)
+        gslv_grid.add_shunt(sh)
+
+
+def convert_generator(k: int, elm: dev.Generator, bus_dict: Dict[str, "pg.Bus"], n_time: int,
+                      use_time_series: bool, time_indices: IntVec | None = None,
+                      opf_results: OptimalPowerFlowResults | None = None) -> "pg.Generator":
+    """
+
+    :param k:
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :param opf_results:
+    :return:
+    """
+    gen = pg.Generator(
+        nt=n_time,
+        name=elm.name,
+        idtag=elm.idtag,
+        active=elm.active,
+        P=elm.P,
+        power_factor=elm.Pf,
+        vset=elm.Vset,
+        Pmin=elm.Pmin,
+        Pmax=elm.Pmax,
+        Qmin=elm.Qmin,
+        Qmax=elm.Qmax,
+        Snom=elm.Snom,
+        is_controlled=elm.is_controlled,
+        enabled_dispatch=elm.enabled_dispatch,
+        q_points=elm.q_curve.get_data().tolist(),
+        use_reactive_power_curve=elm.use_reactive_power_curve
+    )
+
+    gen.bus = bus_dict[elm.bus.idtag]
+
+    fill_profile(gslv_profile=gen.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    if opf_results is None:
+        fill_profile(gslv_profile=gen.P,
+                     gc_profile=elm.P_prof,
+                     use_time_series=use_time_series,
+                     time_indices=time_indices,
+                     n_time=n_time,
+                     default_val=elm.P)
+    else:
+        fill_profile_with_array(gslv_profile=gen.P,
+                                arr=opf_results.generator_power[:, k] - opf_results.generator_shedding[:, k],
+                                use_time_series=use_time_series,
+                                time_indices=time_indices,
+                                n_time=n_time,
+                                default_val=elm.P)
+
+    fill_profile(gslv_profile=gen.Pf,
+                 gc_profile=elm.Pf_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Pf)
+
+    fill_profile(gslv_profile=gen.Vset,
+                 gc_profile=elm.Vset_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Vset)
+
+    fill_profile(gslv_profile=gen.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    fill_profile(gslv_profile=gen.Cost0,
+                 gc_profile=elm.Cost0_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost0)
+
+    fill_profile(gslv_profile=gen.Cost2,
+                 gc_profile=elm.Cost2_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost2)
+
+    return gen
+
+
+def add_generators(circuit: MultiCircuit,
+                   gslv_grid: "pg.MultiCircuit",
+                   bus_dict: Dict[str, "pg.Bus"],
+                   time_series: bool,
+                   n_time=1,
+                   time_indices: Union[IntVec, None] = None,
+                   opf_results: Union[None, OptimalPowerFlowResults] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    :param opf_results: OptimalPowerFlowResults (optional)
+    """
+    devices = circuit.get_generators()
+
+    for k, elm in enumerate(devices):
+        gen = convert_generator(k=k, elm=elm, bus_dict=bus_dict,
+                                n_time=n_time, use_time_series=time_series,
+                                time_indices=time_indices, opf_results=opf_results)
+
+        gslv_grid.add_generator(gen)
+
+
+def convert_battery(k: int, elm: dev.Battery, bus_dict: Dict[str, "pg.Bus"], n_time: int,
+                    use_time_series: bool, time_indices: IntVec | None = None,
+                    opf_results: OptimalPowerFlowResults | None = None) -> "pg.Battery":
+    """
+
+    :param k:
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :param opf_results:
+    :return:
+    """
+    gen = pg.Battery(nt=n_time,
+                     name=elm.name,
+                     idtag=elm.idtag,
+                     P=elm.P,
+                     power_factor=elm.Pf,
+                     vset=elm.Vset,
+                     max_soc=elm.max_soc,
+                     min_soc=elm.min_soc,
+                     Qmin=elm.Qmin,
+                     Qmax=elm.Qmax,
+                     Pmin=elm.Pmin,
+                     Pmax=elm.Pmax,
+                     Snom=elm.Snom,
+                     Enom=elm.Enom,
+                     charge_efficiency=elm.charge_efficiency,
+                     discharge_efficiency=elm.discharge_efficiency,
+                     is_controlled=elm.is_controlled, )
+
+    gen.bus = bus_dict[elm.bus.idtag]
+
+    fill_profile(gslv_profile=gen.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    if opf_results is None:
+        fill_profile(gslv_profile=gen.P,
+                     gc_profile=elm.P_prof,
+                     use_time_series=use_time_series,
+                     time_indices=time_indices,
+                     n_time=n_time,
+                     default_val=elm.P)
+    else:
+        fill_profile_with_array(gslv_profile=gen.P,
+                                arr=opf_results.battery_power[:, k],
+                                use_time_series=use_time_series,
+                                time_indices=time_indices,
+                                n_time=n_time,
+                                default_val=elm.P)
+
+    fill_profile(gslv_profile=gen.Pf,
+                 gc_profile=elm.Pf_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Pf)
+
+    fill_profile(gslv_profile=gen.Vset,
+                 gc_profile=elm.Vset_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Vset)
+
+    fill_profile(gslv_profile=gen.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    fill_profile(gslv_profile=gen.Cost0,
+                 gc_profile=elm.Cost0_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost0)
+
+    fill_profile(gslv_profile=gen.Cost2,
+                 gc_profile=elm.Cost2_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost2)
+
+    return gen
+
+
+def add_battery_data(circuit: MultiCircuit,
+                     gslv_grid: "pg.MultiCircuit",
+                     bus_dict: Dict[str, "pg.Bus"],
+                     time_series: bool,
+                     n_time: int = 1,
+                     time_indices: Union[IntVec, None] = None,
+                     opf_results: Union[None, OptimalPowerFlowResults] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    :param opf_results: OptimelPowerFlowResults (optional)
+    """
+    devices = circuit.get_batteries()
+
+    for k, elm in enumerate(devices):
+        batt = convert_battery(k=k, elm=elm, bus_dict=bus_dict,
+                               n_time=n_time, use_time_series=time_series,
+                               time_indices=time_indices, opf_results=opf_results)
+
+        gslv_grid.add_battery(batt)
+
+
+def convert_line(elm: dev.Line,
+                 n_time: int,
+                 bus_dict: Dict[str, "pg.Bus"],
+                 branch_groups_dict: Dict[dev.BranchGroup, "pg.BranchGroup"],
+                 use_time_series: bool, time_indices: IntVec | None = None, ) -> "pg.Line":
+    """
+
+    :param elm:
+    :param n_time:
+    :param bus_dict:
+    :param branch_groups_dict:
+    :param use_time_series:
+    :param time_indices:
+    :return:
+    """
+    lne = pg.Line(
+        idtag=elm.idtag,
+        code=str(elm.code),
+        name=elm.name,
+        bus_from=bus_dict[elm.bus_from.idtag],
+        bus_to=bus_dict[elm.bus_to.idtag],
+        nt=n_time,
+        length=elm.length,
+        rate=elm.rate if elm.rate > 0 else 9999,
+        active=elm.active,
+        r=elm.R,
+        x=elm.X,
+        b=elm.B,
+        monitor_loading=elm.monitor_loading,
+        contingency_enabled=elm.contingency_enabled,
+    )
+
+    lne.group = branch_groups_dict.get(elm.group, None)
+
+    fill_profile(gslv_profile=lne.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    fill_profile(gslv_profile=lne.rate,
+                 gc_profile=elm.rate_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.rate)
+
+    fill_profile(gslv_profile=lne.contingency_factor,
+                 gc_profile=elm.contingency_factor_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.contingency_factor)
+
+    fill_profile(gslv_profile=lne.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    return lne
+
+
+def add_lines(circuit: MultiCircuit,
+              gslv_grid: "pg.MultiCircuit",
+              bus_dict: Dict[str, "pg.Bus"],
+              branch_groups_dict: Dict[dev.BranchGroup, "pg.BranchGroup"],
+              time_series: bool,
+              n_time: int = 1,
+              time_indices: Union[IntVec, None] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param branch_groups_dict: dictionary of converted branch groups
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    """
+
+    # Compile the lines
+    for i, elm in enumerate(circuit.lines):
+        lne = convert_line(elm=elm,
+                           bus_dict=bus_dict,
+                           branch_groups_dict=branch_groups_dict,
+                           n_time=n_time,
+                           use_time_series=time_series,
+                           time_indices=time_indices)
+        gslv_grid.add_line(lne)
+
+
+def convert_transformer(elm: dev.Transformer2W,
+                        bus_dict: Dict[str, "pg.Bus"],
+                        branch_groups_dict: Dict[dev.BranchGroup, "pg.BranchGroup"],
+                        n_time: int,
+                        use_time_series: bool, time_indices: IntVec | None,
+                        override_controls: bool) -> "pg.Transformer2W":
+    """
+
+    :param elm:
+    :param bus_dict:
+    :param branch_groups_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :param override_controls:
+    :return:
+    """
+    tr2 = pg.Transformer2W(idtag=elm.idtag,
+                           code=str(elm.code),
+                           name=elm.name,
+                           bus_from=bus_dict[elm.bus_from.idtag],
+                           bus_to=bus_dict[elm.bus_to.idtag],
+                           nt=n_time,
+                           HV=elm.HV,
+                           LV=elm.LV,
+                           rate=elm.rate if elm.rate > 0 else 9999,
+                           active=elm.active,
+                           r=elm.R,
+                           x=elm.X,
+                           g=elm.G,
+                           b=elm.B,
+                           monitor_loading=elm.monitor_loading,
+                           contingency_enabled=elm.contingency_enabled,
+                           tap_module=elm.tap_module,
+                           tap_phase=elm.tap_phase)
+
+    tr2.tap_phase_min = elm.tap_phase_min
+    tr2.tap_phase_max = elm.tap_phase_max
+    tr2.tap_module_min = elm.tap_module_min
+    tr2.tap_module_max = elm.tap_module_max
+
+    fill_profile(gslv_profile=tr2.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    fill_profile(gslv_profile=tr2.rate,
+                 gc_profile=elm.rate_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.rate)
+
+    fill_profile(gslv_profile=tr2.contingency_factor,
+                 gc_profile=elm.contingency_factor_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.contingency_factor)
+
+    fill_profile(gslv_profile=tr2.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    fill_profile(gslv_profile=tr2.Pset,
+                 gc_profile=elm.Pset_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Pset)
+
+    fill_profile(gslv_profile=tr2.Qset,
+                 gc_profile=elm.Qset_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Qset)
+
+    fill_profile(gslv_profile=tr2.vset,
+                 gc_profile=elm.vset_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.vset)
+
+    fill_profile(gslv_profile=tr2.tap_phase,
+                 gc_profile=elm.tap_phase_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.tap_phase)
+
+    fill_profile(gslv_profile=tr2.tap_phase_control_mode,
+                 gc_profile=elm.tap_phase_control_mode_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.tap_phase_control_mode)
+
+    fill_profile(gslv_profile=tr2.tap_module,
+                 gc_profile=elm.tap_module_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.tap_module)
+
+    fill_profile(gslv_profile=tr2.tap_module_control_mode,
+                 gc_profile=elm.tap_module_control_mode_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.tap_module_control_mode)
+
+    # control vars
+    if override_controls:
+        tr2.tap_module_control_mode = pg.TapModuleControl.fixed
+        tr2.tap_phase_control_mode = pg.TapPhaseControl.fixed
+    else:
+        # tr2.setAllControlMode(ctrl_dict[elm.control_mode])
+        pass
+
+    return tr2
+
+
+def add_transformers(circuit: MultiCircuit,
+                     gslv_grid: "pg.MultiCircuit",
+                     bus_dict: Dict[str, "pg.Bus"],
+                     branch_groups_dict: Dict[dev.BranchGroup, "pg.BranchGroup"],
+                     time_series: bool,
+                     n_time: int = 1,
+                     time_indices: Union[IntVec, None] = None,
+                     override_controls=False):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param branch_groups_dict: dictionary of branch grous converetd
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    :param override_controls: If true the controls are set to Fix
+    """
+
+    for i, elm in enumerate(circuit.transformers2w):
+        tr2 = convert_transformer(elm=elm,
+                                  bus_dict=bus_dict,
+                                  branch_groups_dict=branch_groups_dict,
+                                  n_time=n_time,
+                                  use_time_series=time_series,
+                                  time_indices=time_indices,
+                                  override_controls=override_controls)
+        gslv_grid.add_transformer(tr2)
+
+
+def convert_transformer3w(elm: dev.Transformer3W,
+                          bus_dict: Dict[str, "pg.Bus"],
+                          n_time: int,
+                          use_time_series: bool,
+                          time_indices: IntVec | None,
+                          override_controls: bool) -> "pg.Transformer3W":
+    """
+
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :param override_controls:
+    :return:
+    """
+
+    tr3 = pg.Transformer3W(idtag=elm.idtag,
+                           code=str(elm.code),
+                           name=elm.name,
+                           nt=n_time,
+                           active=elm.active,
+                           bus1=bus_dict[elm.bus1.idtag] if elm.bus1 else None,
+                           bus2=bus_dict[elm.bus2.idtag] if elm.bus2 else None,
+                           bus3=bus_dict[elm.bus3.idtag] if elm.bus3 else None,
+                           V1=elm.V1,
+                           V2=elm.V2,
+                           V3=elm.V3,
+                           r12=elm.r12, r23=elm.r23, r31=elm.r31,
+                           x12=elm.x12, x23=elm.x23, x31=elm.x31,
+                           rate12=elm.rate1, rate23=elm.rate2, rate31=elm.rate3)
+
+    # this is because the central node is in the buses list already from GridCal
+    tr3.central_node = bus_dict[elm.bus0.idtag]
+
+    if use_time_series:
+        pass
+    else:
+        pass
+
+    return tr3
+
+
+def add_transformers3w(circuit: MultiCircuit,
+                       gslv_grid: "pg.MultiCircuit",
+                       bus_dict: Dict[str, "pg.Bus"],
+                       time_series: bool,
+                       n_time=1,
+                       time_indices: Union[IntVec, None] = None,
+                       override_controls=False):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    :param override_controls: If true the controls are set to Fix
+    """
+    for i, elm in enumerate(circuit.transformers3w):
+        tr3 = convert_transformer3w(elm=elm, bus_dict=bus_dict, n_time=n_time,
+                                    use_time_series=time_series, time_indices=time_indices,
+                                    override_controls=override_controls)
+
+        # because the central bus was added already, do not add it here
+        gslv_grid.add_transformer_3w(tr3)
+
+
+def convert_vsc(elm: dev.VSC, bus_dict: Dict[str, "pg.Bus"], n_time: int,
+                use_time_series: bool, time_indices: IntVec | None) -> "pg.Vsc":
+    """
+
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :return:
+    """
+    vsc = pg.Vsc(idtag=elm.idtag,
+                 code=str(elm.code),
+                 name=elm.name,
+                 calc_node_from=bus_dict[elm.bus_from.idtag],
+                 calc_node_to=bus_dict[elm.bus_to.idtag],
+                 nt=n_time,
+                 active=elm.active, )
+
+    vsc.alpha1 = elm.alpha1
+    vsc.alpha2 = elm.alpha2
+    vsc.alpha3 = elm.alpha3
+
+    vsc.setAllMonitorloading(elm.monitor_loading)
+    vsc.setAllContingencyenabled(elm.contingency_enabled)
+
+    fill_profile(gslv_profile=vsc.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    fill_profile(gslv_profile=vsc.rate,
+                 gc_profile=elm.rate_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.rate)
+
+    fill_profile(gslv_profile=vsc.contingency_factor,
+                 gc_profile=elm.contingency_factor_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.contingency_factor)
+
+    fill_profile(gslv_profile=vsc.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    return vsc
+
+
+def add_vscs(circuit: MultiCircuit,
+             gslv_grid: "pg.MultiCircuit",
+             bus_dict: Dict[str, "pg.Bus"],
+             time_series: bool,
+             n_time: int = 1,
+             time_indices: Union[IntVec, None] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    """
+    for i, elm in enumerate(circuit.vsc_devices):
+        vsc = convert_vsc(elm=elm, bus_dict=bus_dict, n_time=n_time,
+                          use_time_series=time_series, time_indices=time_indices)
+        gslv_grid.add_vsc(vsc)
+
+
+def convert_dc_line(elm: dev.DcLine, bus_dict: Dict[str, "pg.Bus"], n_time: int,
+                    use_time_series: bool, time_indices: IntVec | None) -> "pg.DcLine":
+    """
+
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :return:
+    """
+    lne = pg.DcLine(idtag=elm.idtag,
+                    name=elm.name,
+                    calc_node_from=bus_dict[elm.bus_from.idtag],
+                    calc_node_to=bus_dict[elm.bus_to.idtag],
+                    nt=n_time,
+                    length=elm.length,
+                    rate=elm.rate,
+                    active_default=elm.active,
+                    r=elm.R,
+                    monitor_loading_default=elm.monitor_loading,
+                    monitor_contingency_default=elm.contingency_enabled
+                    )
+
+    fill_profile(gslv_profile=lne.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    fill_profile(gslv_profile=lne.rate,
+                 gc_profile=elm.rate_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.rate)
+
+    fill_profile(gslv_profile=lne.contingency_factor,
+                 gc_profile=elm.contingency_factor_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.contingency_factor)
+
+    fill_profile(gslv_profile=lne.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    return lne
+
+
+def add_dc_lines(circuit: MultiCircuit,
+                 gslv_grid: "pg.MultiCircuit",
+                 bus_dict: Dict[str, "pg.Bus"],
+                 time_series: bool,
+                 n_time: int = 1,
+                 time_indices: Union[IntVec, None] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    """
+    # Compile the lines
+    for i, elm in enumerate(circuit.dc_lines):
+        lne = convert_dc_line(elm=elm, bus_dict=bus_dict, n_time=n_time,
+                              use_time_series=time_series, time_indices=time_indices)
+        gslv_grid.add_dc_line(lne)
+
+
+def convert_hvdc_line(elm: dev.HvdcLine, bus_dict: Dict[str, "pg.Bus"], n_time: int,
+                      use_time_series: bool, time_indices: IntVec | None) -> "pg.HvdcLine":
+    """
+
+    :param elm:
+    :param bus_dict:
+    :param n_time:
+    :param use_time_series:
+    :param time_indices:
+    :return:
+    """
+
+    hvdc = pg.HvdcLine(idtag=elm.idtag,
+                       code=str(elm.code),
+                       name=elm.name,
+                       calc_node_from=bus_dict[elm.bus_from.idtag],
+                       calc_node_to=bus_dict[elm.bus_to.idtag],
+                       cn_from=None,
+                       cn_to=None,
+                       nt=n_time,
+                       active_default=int(elm.active),
+                       rate=elm.rate,
+                       contingency_rate=elm.rate * elm.contingency_factor,
+                       monitor_loading_default=1,
+                       monitor_contingency_default=1,
+                       P=elm.Pset,
+                       Vf=elm.Vset_f,
+                       Vt=elm.Vset_t,
+                       r=elm.r,
+                       angle_droop=elm.angle_droop,
+                       length=elm.length,
+                       min_firing_angle_f=elm.min_firing_angle_f,
+                       max_firing_angle_f=elm.max_firing_angle_f,
+                       min_firing_angle_t=elm.min_firing_angle_t,
+                       max_firing_angle_t=elm.max_firing_angle_t,
+                       control_mode=hvdc_control_mode_dict[elm.control_mode])
+
+    fill_profile(gslv_profile=hvdc.active,
+                 gc_profile=elm.active_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.active)
+
+    fill_profile(gslv_profile=hvdc.Vset_f,
+                 gc_profile=elm.Vset_f_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Vset_f)
+
+    fill_profile(gslv_profile=hvdc.Vset_f,
+                 gc_profile=elm.Vset_f_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Vset_f)
+
+    fill_profile(gslv_profile=hvdc.Vset_t,
+                 gc_profile=elm.Vset_t_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Vset_t)
+
+    fill_profile(gslv_profile=hvdc.angle_droop,
+                 gc_profile=elm.angle_droop_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.angle_droop)
+
+    fill_profile(gslv_profile=hvdc.rate,
+                 gc_profile=elm.rate_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.rate)
+
+    fill_profile(gslv_profile=hvdc.contingency_factor,
+                 gc_profile=elm.contingency_factor_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.contingency_factor)
+
+    fill_profile(gslv_profile=hvdc.cost,
+                 gc_profile=elm.Cost_prof,
+                 use_time_series=use_time_series,
+                 time_indices=time_indices,
+                 n_time=n_time,
+                 default_val=elm.Cost)
+
+    return hvdc
+
+
+def add_hvdcs(circuit: MultiCircuit,
+              gslv_grid: "pg.MultiCircuit",
+              bus_dict: Dict[str, "pg.Bus"],
+              time_series: bool,
+              n_time=1,
+              time_indices: Union[IntVec, None] = None):
+    """
+
+    :param circuit: GridCal circuit
+    :param gslv_grid: GSLV circuit
+    :param time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param bus_dict: dictionary of bus id to GSLV bus object
+    :param n_time: number of time steps
+    :param time_indices: Array of time indices
+    """
+
+    for i, elm in enumerate(circuit.hvdc_lines):
+        hvdc = convert_hvdc_line(elm=elm, bus_dict=bus_dict, n_time=n_time,
+                                 use_time_series=time_series, time_indices=time_indices)
+        gslv_grid.add_hvdc_line(hvdc)
+
+
+class GslvDicts:
+
+    def __init__(self):
+        self.area_dict: Dict[dev.Area, "pg.Area"] = dict()
+
+        self.zone_dict: Dict[dev.Zone, "pg.Zone"] = dict()
+
+        self.substation_dict: Dict[dev.Substation, "pg.Substation"] = dict()
+
+        self.voltage_level_dict: Dict[dev.VoltageLevel, "pg.VoltageLevel"] = dict()
+
+        self.country_dict: Dict[dev.Country, "pg.Country"] = dict()
+
+        self.facility_dict: Dict[dev.Facility, "pg.Facility"] = dict()
+
+        self.regions_dict: Dict[dev.Country, "pg.Country"] = dict()
+
+        self.con_groups_dict: Dict[dev.ContingencyGroup, "pg.ContingencyGroup"] = dict()
+
+        self.inv_groups_dict: Dict[dev.InvestmentsGroup, "pg.InvestmentGroup"] = dict()
+
+
+def to_gslv(circuit: MultiCircuit,
+            use_time_series: bool,
+            time_indices: Union[IntVec, None] = None,
+            override_branch_controls=False,
+            opf_results: Union[None, OptimalPowerFlowResults] = None) -> Tuple["pg.MultiCircuit", GslvDicts]:
+    """
+    Convert GridCal circuit to GSLV
+    :param circuit: MultiCircuit
+    :param use_time_series: compile the time series from GridCal? otherwise just the snapshot
+    :param time_indices: Array of time indices
+    :param override_branch_controls: If true the branch controls are set to Fix
+    :param opf_results:
+    :return: pg.MultiCircuit instance
+    """
+
+    dicts = GslvDicts()
+
+    if time_indices is None:
+        n_time = circuit.get_time_number() if use_time_series else 1
+        if n_time == 0:
+            n_time = 1
+    else:
+        n_time = len(time_indices)
+
+    pg_grid = pg.MultiCircuit(name=circuit.name,
+                              nt=n_time,
+                              Sbase=circuit.Sbase,
+                              fBase=circuit.fBase,
+                              idtag=circuit.idtag)
+
+    dicts.area_dict = add_areas(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.zone_dict = add_zones(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.substation_dict = add_substations(circuit=circuit, gslv_grid=pg_grid, n_time=n_time)
+
+    dicts.voltage_level_dict = add_voltage_levels(circuit=circuit, gslv_grid=pg_grid,
+                                                  substations_dict=dicts.substation_dict)
+
+    dicts.country_dict = add_countries(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.facility_dict = add_facilities(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.modelling_authorities_dict = add_modelling_authorities(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.branch_groups_dict = add_branch_groups(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.municipalities_dict = add_municipalities(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.regions_dict = add_regions(circuit=circuit, gslv_grid=pg_grid)
+
+    dicts.con_groups_dict = add_contingency_groups(circuit=circuit, gslv_grid=pg_grid)
+
+    add_contingencies(circuit=circuit, gslv_grid=pg_grid, n_time=n_time, groups_dict=dicts.con_groups_dict)
+
+    dicts.inv_groups_dict = add_investment_groups(circuit=circuit, gslv_grid=pg_grid)
+
+    add_investments(circuit=circuit, gslv_grid=pg_grid, groups_dict=dicts.inv_groups_dict)
+
+    bus_dict = add_buses(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        use_time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices,
+        area_dict=dicts.area_dict,
+        zone_dict=dicts.zone_dict,
+        substation_dict=dicts.substation_dict,
+        voltage_level_dict=dicts.voltage_level_dict,
+        country_dict=dicts.country_dict,
+    )
+
+    add_loads(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        use_time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices,
+    )
+
+    add_static_generators(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices
+
+    )
+
+    add_shunts(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices
+    )
+
+    add_generators(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices,
+        opf_results=opf_results
+    )
+
+    add_battery_data(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices,
+        opf_results=opf_results
+    )
+
+    add_lines(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        branch_groups_dict=dicts.branch_groups_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices
+    )
+
+    add_transformers(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        branch_groups_dict=dicts.branch_groups_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices,
+        override_controls=override_branch_controls
+    )
+
+    add_transformers3w(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices,
+        override_controls=override_branch_controls
+    )
+
+    add_vscs(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices
+    )
+
+    add_dc_lines(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices
+    )
+
+    add_hvdcs(
+        circuit=circuit,
+        gslv_grid=pg_grid,
+        bus_dict=bus_dict,
+        time_series=use_time_series,
+        n_time=n_time,
+        time_indices=time_indices
+    )
+
+    return pg_grid, dicts
+
+
+class FakeAdmittances:
+    """
+    Fake admittances class needed to make the translation
+    """
+
+    def __init__(self) -> None:
+        self.Ybus = None
+        self.Yf = None
+        self.Yt = None
+
+
+def get_snapshots_from_gslv(circuit: MultiCircuit, override_branch_controls=False) -> List[NumericalCircuit]:
+    """
+
+    :param circuit:
+    :param override_branch_controls:
+    :return:
+    """
+
+    gslv_grid, (bus_dict, area_dict, zone_dict) = to_gslv(circuit,
+                                                          use_time_series=False,
+                                                          override_branch_controls=override_branch_controls)
+
+    logger = pg.Logger()
+    npa_data_lst = pg.compile(gslv_grid, logger, t=0).splitIntoIslands()
+
+    data_lst = list()
+
+    for npa_data in npa_data_lst:
+        data = NumericalCircuit(nbus=0,
+                                nbr=0,
+                                nhvdc=0,
+                                nvsc=0,
+                                nload=0,
+                                ngen=0,
+                                nbatt=0,
+                                nshunt=0,
+                                nfluidnode=0,
+                                nfluidturbine=0,
+                                nfluidpump=0,
+                                nfluidp2x=0,
+                                nfluidpath=0,
+                                sbase=0,
+                                t_idx=0)
+
+        conn = npa_data.getConnectivity()
+        inj = npa_data.getInjections()
+        tpes = npa_data.getSimulationIndices(inj.S0)
+        adm = npa_data.getAdmittances(conn)
+        lin = npa_data.getLinearMatrices(conn)
+        series_adm = npa_data.getSeriesAdmittances(conn)
+        fd_adm = npa_data.getFastDecoupledAdmittances(conn, tpes)
+        qlim = npa_data.getQLimits()
+
+        data.Vbus_ = npa_data.Vbus.reshape(-1, 1)
+        data.Sbus_ = inj.S0.reshape(-1, 1)
+        data.Ibus_ = inj.I0.reshape(-1, 1)
+        data.passive_branch_data.names = np.array(npa_data.passive_branch_data.names)
+        data.passive_branch_data.virtual_tap_f = npa_data.passive_branch_data.vtap_f
+        data.passive_branch_data.virtual_tap_t = npa_data.passive_branch_data.vtap_t
+        data.passive_branch_data.original_idx = npa_data.passive_branch_data.original_indices
+
+        data.bus_data.names = np.array(npa_data.bus_data.names)
+        data.bus_data.original_idx = npa_data.bus_data.original_indices
+
+        data.admittances_ = FakeAdmittances()
+        data.admittances_.Ybus = adm.Ybus
+        data.admittances_.Yf = adm.Yf
+        data.admittances_.Yt = adm.Yt
+
+        data.Bbus_ = lin.Bbus
+        data.Bf_ = lin.Bf
+
+        data.Yseries_ = series_adm.Yseries
+        data.Yshunt_ = series_adm.Yshunt
+
+        data.B1_ = fd_adm.B1
+        data.B2_ = fd_adm.B2
+
+        data.Cf_ = conn.Cf
+        data.Ct_ = conn.Ct
+
+        data.bus_data.bus_types = tpes.types
+        data.pq_ = tpes.pq
+        data.pv_ = tpes.pv
+        data.vd_ = tpes.vd
+        data.pqpv_ = tpes.no_slack
+
+        data.Qmax_bus_ = qlim.qmax_bus
+        data.Qmin_bus_ = qlim.qmin_bus
+
+        control_indices = npa_data.getSimulationIndices(Sbus=data.Sbus_[:, 0])
+
+        data.k_pf_tau = control_indices.k_pf_tau
+        data.k_qf_m = control_indices.k_qf_m
+        data.k_zero_beq = control_indices.k_qf_beq
+        data.k_vf_beq = control_indices.k_vf_beq
+        data.k_vt_m = control_indices.k_v_m
+        data.k_qt_m = control_indices.k_qt_m
+        data.k_pf_dp = control_indices.k_pf_dp
+        data.i_vsc = control_indices.i_vsc
+        data.i_vf_beq = control_indices.i_vf_beq
+        data.i_vt_m = control_indices.i_vt_m
+
+        data_lst.append(data)
+
+    return data_lst
+
+
+def get_gslv_pf_options(opt: PowerFlowOptions) -> "pg.PowerFlowOptions":
+    """
+    Translate GridCal power flow options to GSLV power flow options
+    :param opt:
+    :return:
+    """
+    solver_dict = {SolverType.NR: pg.SolverType.NR,
+                   SolverType.DC: pg.SolverType.DC,
+                   SolverType.HELM: pg.SolverType.HELM,
+                   SolverType.IWAMOTO: pg.SolverType.IWAMOTO,
+                   SolverType.LM: pg.SolverType.LM,
+                   SolverType.LACPF: pg.SolverType.LACPF,
+                   SolverType.FASTDECOUPLED: pg.SolverType.FASTDECOUPLED
+                   }
+
+    if opt.solver_type in solver_dict.keys():
+        solver_type = solver_dict[opt.solver_type]
+    else:
+        solver_type = pg.SolverType.NR
+
+    """
+    solver_type: GSLVpa.SolverType = <SolverType.NR: 0>, 
+    retry_with_other_methods: bool = True, 
+    verbose: bool = False, 
+    initialize_with_existing_solution: bool = False, 
+    tolerance: float = 1e-06, 
+    max_iter: int = 15, 
+    control_q_mode: GSLVpa.ReactivePowerControlMode = <ReactivePowerControlMode.NoControl: 0>, 
+    tap_control_mode: GSLVpa.TapsControlMode = <TapsControlMode.NoControl: 0>, 
+    distributed_slack: bool = False, 
+    ignore_single_node_islands: bool = False, 
+    correction_parameter: float = 0.5, 
+    mu0: float = 1.0
+    """
+
+    """
+    solver_type: pygslv.SolverType = <SolverType.NR: 0>, 
+    retry_with_other_methods: bool = True, 
+    verbose: int = 0, 
+    initialize_with_existing_solution: bool = False, 
+    tolerance: float = 1e-06, 
+    max_iter: int = 25, 
+    max_outer_loop_iter: int = 100, 
+    control_Q: bool = True, 
+    control_taps_modules: bool = True, 
+    control_taps_phase: bool = True, 
+    control_remote_voltage: bool = True, 
+    orthogonalize_controls: bool = True, 
+    apply_temperature_correction: bool = True, 
+    branch_impedance_tolerance_mode: pygslv.BranchImpedanceMode = <BranchImpedanceMode.Specified: 0>, 
+    distributed_slack: bool = False, 
+    ignore_single_node_islands: bool = False, 
+    trust_radius: float = 1.0, 
+    backtracking_parameter: float = 0.05, 
+    use_stored_guess: bool = False, 
+    generate_report: bool = False)
+
+    """
+
+    return pg.PowerFlowOptions(
+        solver_type=solver_type,
+        retry_with_other_methods=opt.retry_with_other_methods,
+        verbose=opt.verbose,
+        initialize_with_existing_solution=opt.use_stored_guess,
+        tolerance=opt.tolerance,
+        max_iter=opt.max_iter,
+        control_Q=opt.control_Q,
+        control_taps_modules=opt.control_taps_modules,
+        control_taps_phase=opt.control_taps_phase,
+        control_remote_voltage=opt.control_remote_voltage,
+        orthogonalize_controls=opt.orthogonalize_controls,
+        apply_temperature_correction=opt.orthogonalize_controls,
+        branch_impedance_tolerance_mode=pg.BranchImpedanceMode.Specified,
+        distributed_slack=opt.distributed_slack,
+        ignore_single_node_islands=opt.ignore_single_node_islands,
+        trust_radius=opt.trust_radius,
+        backtracking_parameter=opt.backtracking_parameter,
+        use_stored_guess=opt.use_stored_guess,
+        generate_report=opt.generate_report
+    )
+
+
+def gslv_pf(circuit: MultiCircuit,
+            pf_opt: PowerFlowOptions,
+            time_series: bool = False,
+            time_indices: Union[IntVec, None] = None,
+            opf_results: Union[None, OptimalPowerFlowResults] = None) -> "pg.PowerFlowResults":
+    """
+    GSLV power flow
+    :param circuit: MultiCircuit instance
+    :param pf_opt: Power Flow Options
+    :param time_series: Compile with GridCal time series?
+    :param time_indices: Array of time indices
+    :param opf_results: Instance of
+    :return: GSLV Power flow results object
+    """
+    override_branch_controls = not (pf_opt.control_taps_modules and pf_opt.control_taps_phase)
+
+    gslv_grid, _ = to_gslv(circuit,
+                           use_time_series=time_series,
+                           time_indices=None,
+                           override_branch_controls=override_branch_controls,
+                           opf_results=opf_results)
+
+    pf_options = get_gslv_pf_options(pf_opt)
+
+    if time_series:
+        # it is already sliced to the relevant time indices
+        if time_indices is None:
+            time_indices = [i for i in range(circuit.get_time_number())]
+        else:
+            time_indices = list(time_indices)
+        n_threads = 0  # max threads
+    else:
+        time_indices = [0]
+        n_threads = 1
+
+    pf_res = pg.multi_island_pf(grid=gslv_grid,
+                                options=pf_options,
+                                time_indices=time_indices,
+                                n_threads=n_threads)
+
+    return pf_res
+
+
+def translate_gslv_pf_results(grid: MultiCircuit, res: "pg.PowerFlowResults") -> PowerFlowResults:
+    """
+    Translate the GSLV Power Analytics results back to GridCal
+    :param grid: MultiCircuit instance
+    :param res: GSLV's PowerFlowResults instance
+    :return: PowerFlowResults instance
+    """
+    results = PowerFlowResults(
+        n=grid.get_bus_number(),
+        m=grid.get_branch_number_wo_hvdc(),
+        n_hvdc=grid.get_hvdc_number(),
+        n_vsc=grid.get_vsc_number(),
+        n_gen=grid.get_generators_number(),
+        n_batt=grid.get_batteries_number(),
+        n_sh=grid.get_shunt_like_device_number(),
+        bus_names=grid.get_bus_names(),
+        branch_names=grid.get_branch_names(add_switch=True),
+        hvdc_names=grid.get_hvdc_names(),
+        vsc_names=grid.get_vsc_names(),
+        gen_names=grid.get_generator_names(),
+        batt_names=grid.get_battery_names(),
+        sh_names=grid.get_shunt_like_devices_names(),
+        bus_types=np.ones(grid.get_bus_number(), dtype=int)
+    )
+
+    results.voltage = res.voltage[0, :]
+    results.Sbus = res.S[0, :]
+    results.Sf = res.Sf[0, :]
+    results.St = res.St[0, :]
+    results.loading = res.loading[0, :]
+    results.losses = res.losses[0, :]
+    # results.Vbranch = res.Vbranch[0, :]
+    # results.If = res.If[0, :]
+    # results.It = res.It[0, :]
+    # results.Beq = res.Beq[0, :]
+    results.m = res.tap_module[0, :]
+    results.tap_angle = res.tap_angle[0, :]
+    # results.F = res.F
+    # results.T = res.T
+    # results.hvdc_F = res.hvdc_F[0, :]
+    # results.hvdc_T = res.hvdc_T[0, :]
+    results.Pf_hvdc = res.Pf_hvdc[0, :]
+    results.Pt_hvdc = res.Pt_hvdc[0, :]
+    results.loading_hvdc = res.loading_hvdc[0, :]
+    results.losses_hvdc = res.losses_hvdc[0, :]
+    # results.bus_area_indices = grid.get_bus_area_indices()
+    # results.area_names = [a.name for a in grid.areas]
+    # results.bus_types = convert_bus_types(res.bus_types[0])  # this is a list of lists
+
+    # for rep in res.stats[0]:
+    #     report = ConvergenceReport()
+    #     for i in range(len(rep.converged)):
+    #         report.add(method=rep.solver[i].name,
+    #                    converged=rep.converged[i],
+    #                    error=rep.norm_f[i],
+    #                    elapsed=rep.elapsed[i],
+    #                    iterations=rep.iterations[i])
+    #         results.convergence_reports.append(report)
+
+    return results