ImportForm.cpp

/*
 *  Copyright (C) 2018  Thales Lima Oliveira <thales@ufu.br>
 *
 *  This program is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation; either version 2 of the License, or
 *  any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <https://www.gnu.org/licenses/>.
 */
 
#include "ImportForm.h"
 
#include "../elements/powerElement/Bus.h"
#include "../elements/powerElement/Capacitor.h"
#include "../elements/powerElement/IndMotor.h"
#include "../elements/powerElement/Inductor.h"
#include "../elements/powerElement/Line.h"
#include "../elements/powerElement/Load.h"
#include "../elements/powerElement/SyncGenerator.h"
#include "../elements/powerElement/SyncMotor.h"
#include "../elements/powerElement/Transformer.h"
 
#include "../utils/PropertiesData.h"
#include "../utils/GraphAutoLayout.h"
 
#include "../editors/Workspace.h"
 
ImportForm::ImportForm(wxWindow* parent, Workspace* workspace) : ImportFormBase(parent)
{
    SetInitialSize();
 
    m_parent = parent;
    m_workspace = workspace;
}
 
ImportForm::~ImportForm() {}
 
void ImportForm::OnButtonCancelClick(wxCommandEvent& event)
{
    EndModal(wxID_CANCEL);
    //if(m_workspace) delete m_workspace;
}
 
void ImportForm::OnButtonOKClick(wxCommandEvent& event)
{
    if(ImportSelectedFiles())
        EndModal(wxID_OK);
    else {
        // Error message
        wxMessageDialog msgDialog(this, _("It was not possible to import the selected files."), _("Error"),
                                  wxOK | wxCENTRE | wxICON_ERROR);
        msgDialog.ShowModal();
    }
}
 
bool ImportForm::ImportSelectedFiles()
{
    switch(m_notebook->GetSelection()) {
        case 0: {
            return ImportCEPELFiles();
            break;
        }
        case 1: {
            return ImportMatpowerFiles();
            break;
        }
        default:
            break;
    }
    return false;
}
 
Bus* ImportForm::GetBusFromID(std::vector<Bus*> busList, int id)
{
    for(auto it = busList.begin(), itEnd = busList.end(); it != itEnd; ++it) {
        if((*it)->GetID() == id) return *it;
    }
    return nullptr;
}
 
bool ImportForm::ImportCEPELFiles()
{
    ParseAnarede parseAnarede(m_filePickerANAREDELST->GetFileName(), m_filePickerANAREDEPWF->GetFileName());
    if(!parseAnarede.Parse()) return false;
    double scale = 1.25;
 
    double mvaBasePower = parseAnarede.GetMVAPowerBase();
    auto simProp = m_workspace->GetProperties()->GetSimulationPropertiesData();
    simProp.basePower = mvaBasePower;
    m_workspace->GetProperties()->SetSimulationPropertiesData(simProp);
 
    std::vector<Element*> elementList;
 
    std::vector<Bus*> busList;
    std::vector<SyncGenerator*> syncGeneratorList;
    std::vector<SyncMotor*> syncMotorList;
    std::vector<Load*> loadList;
    std::vector<Inductor*> indList;
    std::vector<Capacitor*> capList;
    std::vector<IndMotor*> indMotorList;
    std::vector<Transformer*> transformerList;
    std::vector<Line*> lineList;
 
    auto components = parseAnarede.GetComponents();
    for(auto it = components.begin(), itEnd = components.end(); it != itEnd; ++it) {
        switch((*it)->type) {
            case ANA_BUS: {
                Bus* bus = new Bus(wxPoint2DDouble((*it)->position.m_x * scale, (*it)->position.m_y * scale));
                bus->SetInserted();
                bus->SetID((*it)->id);
                bus->SetWidth((*it)->length * scale);
                bus->StartMove(bus->GetPosition());
                if((*it)->rotationID == 0) {
                    bus->Move(bus->GetPosition() + wxPoint2DDouble(bus->GetWidth() / 2, bus->GetHeight() / 2));
                } else if((*it)->rotationID == 1) {
                    for(int i = 0; i < (*it)->rotationID * 2; ++i) bus->Rotate();
                    bus->Move(bus->GetPosition() + wxPoint2DDouble(-bus->GetHeight() / 2, bus->GetWidth() / 2));
                } else if((*it)->rotationID == 2) {
                    for(int i = 0; i < (*it)->rotationID * 2; ++i) bus->Rotate();
                    bus->Move(bus->GetPosition() + wxPoint2DDouble(-bus->GetWidth() / 2, -bus->GetHeight() / 2));
                } else if((*it)->rotationID == 3) {
                    for(int i = 0; i < (*it)->rotationID * 2; ++i) bus->Rotate();
                    bus->Move(bus->GetPosition() + wxPoint2DDouble(-bus->GetHeight() / 2, -bus->GetWidth() / 2));
                }
 
                // Electrical data
                auto data = bus->GetElectricalData();
                ParseAnarede::BusData* busData = parseAnarede.GetBusDataFromID((*it)->electricalID);
                if(busData) {
                    data.number = busData->id;
                    data.name = busData->busName;
                    switch(busData->type) {
                        case 0:
                        case 3: {
                            data.isVoltageControlled = false;
                            data.slackBus = false;
                        } break;
                        case 1: {
                            data.isVoltageControlled = true;
                            data.slackBus = false;
                        } break;
                        case 2: {
                            data.isVoltageControlled = true;
                            data.slackBus = true;
                        } break;
                        default: {
                            return false;
                        } break;
                    }
                    data.voltage = std::complex<double>(busData->voltage * std::cos(wxDegToRad(busData->angle)),
                                                        busData->voltage * std::sin(wxDegToRad(busData->angle)));
                    data.controlledVoltage = busData->voltage;
                } else
                    return false;
 
                bus->SetElectricalData(data);
                busList.push_back(bus);
            } break;
            case ANA_GENERATOR:
            case ANA_IND_GENERATOR: {
                // Find parent bus
                Bus* parentBus = GetBusFromID(busList, (*it)->busConnectionID[0].first);
                wxPoint2DDouble nodePos =
                    parseAnarede.GetNodePositionFromID(parentBus, scale, (*it)->busConnectionNode[0].second);
 
                ParseAnarede::BusData* busData = parseAnarede.GetBusDataFromID((*it)->busConnectionID[0].second);
                bool isMotor = false;
                if(busData->genPower.real() <= 0.0) isMotor = true;
                Machines* machine = nullptr;
 
                if(!isMotor) {
                    SyncGenerator* syncGenerator = new SyncGenerator();
 
                    auto data = syncGenerator->GetElectricalData();
                    data.name =
                        wxString::Format("%s %u (%s)", _("Generator"), syncGeneratorList.size() + 1, busData->busName);
                    if((*it)->type == ANA_IND_GENERATOR) {
                        ParseAnarede::IndGenData* genData =
                            static_cast<ParseAnarede::IndGenData*>(parseAnarede.GetIndElementDataFromID(
                                (*it)->electricalID, (*it)->busConnectionID[0].second, ANA_IND_GENERATOR));
                        double numUnits = static_cast<double>(genData->numUnits);
                        data.activePower = genData->power.real() * numUnits;
                        data.reactivePower = genData->power.imag() * numUnits;
                        data.maxReactive = genData->maxReactivePower * numUnits;
                        data.minReactive = genData->minReactivePower * numUnits;
                        data.nominalPower = genData->ratedPower * numUnits;
                        double xd = genData->xd / (100.0 * numUnits);
                        data.syncXd = xd > 1e-3 ? xd : 1.0;
                        data.syncXq = genData->xq / (100.0 * numUnits);
                        data.potierReactance = genData->xl / (100.0 * numUnits);
                    } else {
                        data.activePower = busData->genPower.real();
                        data.reactivePower = busData->genPower.imag();
                        data.minReactive = busData->minReactivePower;
                        data.maxReactive = busData->maxReactivePower;
                        data.useMachineBase = false;  //  ANAREDE use system's base
                    }
 
                    syncGenerator->SetElectricalData(data);
 
                    machine = syncGenerator;
                } else {
                    SyncMotor* syncMotor = new SyncMotor();
 
                    auto data = syncMotor->GetElectricalData();
                    data.name = wxString::Format("%s %u (%s)", _("Synchronous compensator"), syncMotorList.size() + 1,
                                                 busData->busName);
                    data.activePower = busData->genPower.real() == 0.0 ? 0.0 : -busData->genPower.real();
                    data.reactivePower = busData->genPower.imag();
                    data.minReactive = busData->minReactivePower;
                    data.maxReactive = busData->maxReactivePower;
                    data.useMachineBase = false;  //  ANAREDE use system's base
 
                    syncMotor->SetElectricalData(data);
 
                    machine = syncMotor;
                }
 
                if(machine) {
                    machine->SetID((*it)->id);
                    machine->AddParent(parentBus, nodePos);
 
                    machine->StartMove(machine->GetPosition());
                    machine->Move(wxPoint2DDouble((*it)->position.m_x * scale, (*it)->position.m_y * scale));
 
                    for(int i = 0; i < 2; ++i) machine->Rotate(false);  // Set to ANAREDE default rotation
                    for(int i = 0; i < (*it)->rotationID * 2; ++i) machine->Rotate();
 
                    if(!isMotor)
                        syncGeneratorList.push_back(static_cast<SyncGenerator*>(machine));
                    else
                        syncMotorList.push_back(static_cast<SyncMotor*>(machine));
                }
            } break;
            case ANA_LOAD:
            case ANA_IND_LOAD: {
                Bus* parentBus = GetBusFromID(busList, (*it)->busConnectionID[0].first);
                wxPoint2DDouble nodePos =
                    parseAnarede.GetNodePositionFromID(parentBus, scale, (*it)->busConnectionNode[0].second);
 
                ParseAnarede::BusData* busData = parseAnarede.GetBusDataFromID((*it)->busConnectionID[0].second);
 
                Load* load = new Load();
                load->SetID((*it)->id);
                load->AddParent(parentBus, nodePos);
 
                load->StartMove(load->GetPosition());
                load->Move(wxPoint2DDouble((*it)->position.m_x * scale, (*it)->position.m_y * scale));
 
                auto data = load->GetElectricalData();
                data.name = wxString::Format("%s %u (%s)", _("Load"), loadList.size() + 1, busData->busName);
                data.activePower = busData->loadPower.real();
                data.reactivePower = busData->loadPower.imag();
                load->SetElectricalData(data);
 
                for(int i = 0; i < (*it)->rotationID * 2; ++i) load->Rotate();
 
                loadList.push_back(load);
            } break;
            case ANA_SHUNT:
            case ANA_IND_SHUNT: {
                Bus* parentBus = GetBusFromID(busList, (*it)->busConnectionID[0].first);
                wxPoint2DDouble nodePos =
                    parseAnarede.GetNodePositionFromID(parentBus, scale, (*it)->busConnectionNode[0].second);
 
                ParseAnarede::BusData* busData = parseAnarede.GetBusDataFromID((*it)->busConnectionID[0].second);
                bool isInductor = false;
                Shunt* shuntElement = nullptr;
 
                double ql = 0.0;
                if((*it)->type == ANA_IND_SHUNT) {
                    ParseAnarede::IndElementData* shuntData = parseAnarede.GetIndElementDataFromID(
                        (*it)->electricalID, (*it)->busConnectionID[0].second, ANA_IND_SHUNT);
                    if(shuntData) {
                        double numUnits = static_cast<double>(shuntData->numUnits);
                        ql = shuntData->power.imag() * numUnits;
                    }
                } else {
                    ql = busData->shuntReactive;
                }
                if(ql <= 0.0) isInductor = true;
 
                if(!isInductor) {
                    Capacitor* cap = new Capacitor();
                    auto data = cap->GetElectricalData();
                    data.name = wxString::Format("%s %u (%s)", _("Capacitor"), indList.size() + 1, busData->busName);
                    data.reactivePower = ql;
                    cap->SetElectricalData(data);
                    shuntElement = cap;
                } else {
                    Inductor* ind = new Inductor();
                    auto data = ind->GetElectricalData();
                    data.name = wxString::Format("%s %u (%s)", _("Inductor"), indList.size() + 1, busData->busName);
                    data.reactivePower = std::abs(ql);
                    ind->SetElectricalData(data);
                    shuntElement = ind;
                }
 
                if(shuntElement) {
                    shuntElement->SetID((*it)->id);
                    shuntElement->AddParent(parentBus, nodePos);
 
                    shuntElement->StartMove(shuntElement->GetPosition());
                    // Offset (ind->GetHeight() / 2 + 10) to adequate the y coordinate
                    shuntElement->Move(wxPoint2DDouble(
                        (*it)->position.m_x * scale, (*it)->position.m_y * scale + shuntElement->GetHeight() / 2 + 10));
                    if((*it)->rotationID != 0) {
                        // Remove offset in position
                        shuntElement->Move(wxPoint2DDouble((*it)->position.m_x * scale, (*it)->position.m_y * scale));
                        // Get the rotated point ralated to the offset
                        wxPoint2DDouble movePt = shuntElement->RotateAtPosition(
                            shuntElement->GetPosition() + wxPoint2DDouble(0, shuntElement->GetHeight() / 2 + 10),
                            (*it)->rotationID * 90.0, true);
                        shuntElement->Move(movePt);
 
                        for(int i = 0; i < (*it)->rotationID * 2; ++i) shuntElement->Rotate();
                    }
 
                    if(!isInductor)
                        capList.push_back(static_cast<Capacitor*>(shuntElement));
                    else
                        indList.push_back(static_cast<Inductor*>(shuntElement));
                }
            } break;
            case ANA_MIT: {
                // Find parent bus
                Bus* parentBus = GetBusFromID(busList, (*it)->busConnectionID[0].first);
                wxPoint2DDouble nodePos =
                    parseAnarede.GetNodePositionFromID(parentBus, scale, (*it)->busConnectionNode[0].second);
 
                IndMotor* indMotor = new IndMotor();
                indMotor->SetID((*it)->id);
                indMotor->AddParent(parentBus, nodePos);
 
                indMotor->StartMove(indMotor->GetPosition());
                indMotor->Move(wxPoint2DDouble((*it)->position.m_x * scale, (*it)->position.m_y * scale));
 
                auto data = indMotor->GetElectricalData();
                data.name = wxString::Format("%s %u (%s)", _("Motor"), indMotorList.size() + 1,
                                             parentBus->GetElectricalData().name);
                data.activePower = 0.0;
                data.reactivePower = 0.0;
                indMotor->SetElectricalData(data);
 
                for(int i = 0; i < 2; ++i) indMotor->Rotate(false);  // Set to ANAREDE default rotation
                for(int i = 0; i < (*it)->rotationID * 2; ++i) indMotor->Rotate();
 
                indMotorList.push_back(indMotor);
            } break;
            case ANA_TRANSFORMER: {
                Bus* parentBus1 = GetBusFromID(busList, (*it)->busConnectionID[0].first);
                Bus* parentBus2 = GetBusFromID(busList, (*it)->busConnectionID[1].first);
                wxPoint2DDouble nodePos1 =
                    parseAnarede.GetNodePositionFromID(parentBus1, scale, (*it)->busConnectionNode[0].second);
                wxPoint2DDouble nodePos2 =
                    parseAnarede.GetNodePositionFromID(parentBus2, scale, (*it)->busConnectionNode[1].second);
 
                ParseAnarede::BranchData* branchData =
                    parseAnarede.GetBranchDataFromID((*it)->electricalID, (*it)->busConnectionID[0].second,
                                                     (*it)->busConnectionID[1].second, ANA_TRANSFORMER);
 
                Transformer* transformer = new Transformer();
                transformer->SetID((*it)->id);
                transformer->AddParent(parentBus1, nodePos1);
                transformer->AddParent(parentBus2, nodePos2);
 
                transformer->StartMove(transformer->GetPosition());
                transformer->Move(wxPoint2DDouble((*it)->position.m_x * scale, (*it)->position.m_y * scale));
 
                for(int i = 0; i < 2; ++i) transformer->Rotate();  // Set to ANAREDE default rotation
                for(int i = 0; i < (*it)->rotationID * 2; ++i) transformer->Rotate();
 
                auto data = transformer->GetElectricalData();
                data.name =
                    wxString::Format("%s %u (%s - %s)", _("Transformer"), transformerList.size() + 1,
                                     parentBus1->GetElectricalData().name, parentBus2->GetElectricalData().name);
                data.resistance = branchData->resistance / 100.0;
                data.indReactance = branchData->indReactance / 100.0;
                data.turnsRatio = branchData->tap;
                data.phaseShift = branchData->phaseShift;
                transformer->SetElectricaData(data);
 
                transformerList.push_back(transformer);
            } break;
            default: {
            } break;
        }
    }
 
    auto powerLines = parseAnarede.GetLines();
    for(auto it = powerLines.begin(), itEnd = powerLines.end(); it != itEnd; ++it) {
        if((*it)->type == ANA_LINE) {
            Bus* parentBus1 = GetBusFromID(busList, (*it)->busConnectionID[0].first);
            Bus* parentBus2 = GetBusFromID(busList, (*it)->busConnectionID[1].first);
            wxPoint2DDouble nodePos1 =
                parseAnarede.GetNodePositionFromID(parentBus1, scale, (*it)->busConnectionNode[0].second);
            wxPoint2DDouble nodePos2 =
                parseAnarede.GetNodePositionFromID(parentBus2, scale, (*it)->busConnectionNode[1].second);
 
            ParseAnarede::BranchData* branchData = parseAnarede.GetBranchDataFromID(
                (*it)->electricalID, (*it)->busConnectionID[0].second, (*it)->busConnectionID[1].second, ANA_LINE);
 
            Line* line = new Line();
            line->SetID((*it)->id);
            line->AddParent(parentBus1, nodePos1);
            for(unsigned int i = 0; i < (*it)->nodesPosition.size(); ++i)
                line->AddPoint(
                    wxPoint2DDouble((*it)->nodesPosition[i].m_x * scale, (*it)->nodesPosition[i].m_y * scale));
            line->AddParent(parentBus2, nodePos2);
 
            auto data = line->GetElectricalData();
            data.name = wxString::Format("%s %u (%s - %s)", _("Line"), lineList.size() + 1,
                                         parentBus1->GetElectricalData().name, parentBus2->GetElectricalData().name);
            data.resistance = branchData->resistance / 100.0;
            data.indReactance = branchData->indReactance / 100.0;
            data.capSusceptance = branchData->capSusceptance / mvaBasePower;
            line->SetElectricalData(data);
 
            lineList.push_back(line);
        }
    }
 
    // Search for bus data without component
    std::vector<ParseAnarede::BusData*> busDataVector = parseAnarede.GetBusData();
    for(auto it = busDataVector.begin(), itEnd = busDataVector.end(); it != itEnd; ++it) {
        ParseAnarede::BusData* busData = *it;
 
        // Search for bus
        Bus* bus = nullptr;
        for(auto itB = busList.begin(), itBEnd = busList.end(); itB != itBEnd; ++itB) {
            if((*itB)->GetElectricalData().number == busData->id) {
                bus = *itB;
                break;
            }
        }
        if(bus) {
            // Check load data
            if(std::abs(busData->loadPower.real()) > 1e-5 ||
               std::abs(busData->loadPower.imag()) > 1e-5) {  // Have loads
                // Find load associated with the bus
                Load* load = nullptr;
                for(auto itL = loadList.begin(), itLEnd = loadList.end(); itL != itLEnd; ++itL) {
                    if((*itL)->GetParentList().size() > 0) {     // Don't search in empty vectors
                        if((*itL)->GetParentList()[0] == bus) {  // Found load
                            load = *itL;
                            break;
                        }
                    }
                }
                if(!load) {  // The load don't exists, create a new one.
                    Load* newLoad =
                        new Load(wxString::Format("%s %u (%s)", _("Load"), loadList.size() + 1, busData->busName));
                    newLoad->AddParent(bus, bus->GetPosition());
                    auto data = newLoad->GetElectricalData();
                    data.activePower = busData->loadPower.real();
                    data.reactivePower = busData->loadPower.imag();
                    newLoad->SetElectricalData(data);
 
                    loadList.push_back(newLoad);
                }
            }
            // Check generation data
            if(std::abs(busData->genPower.real()) > 1e-5 || std::abs(busData->genPower.imag()) > 1e-5) {
                // Find synchornous machine associated with bus
                if(busData->genPower.real() > 0.0) {  // Synchronous generator
                    SyncGenerator* syncGenerator = nullptr;
                    for(auto itSM = syncGeneratorList.begin(), itSMEnd = syncGeneratorList.end(); itSM != itSMEnd;
                        ++itSM) {
                        if((*itSM)->GetParentList().size() > 0) {     // Don't search in empty vectors
                            if((*itSM)->GetParentList()[0] == bus) {  // Found load
                                syncGenerator = *itSM;
                                break;
                            }
                        }
                    }
                    if(!syncGenerator) {
                        SyncGenerator* newSyncGenerator = new SyncGenerator(wxString::Format(
                            "%s %u (%s)", _("Generator"), syncGeneratorList.size() + 1, busData->busName));
                        newSyncGenerator->AddParent(bus, bus->GetPosition());
                        auto data = newSyncGenerator->GetElectricalData();
                        data.activePower = busData->genPower.real();
                        data.reactivePower = busData->genPower.imag();
                        newSyncGenerator->SetElectricalData(data);
 
                        syncGeneratorList.push_back(newSyncGenerator);
                    }
                } else {
                    SyncMotor* syncMotor = nullptr;
                    for(auto itSM = syncMotorList.begin(), itSMEnd = syncMotorList.end(); itSM != itSMEnd; ++itSM) {
                        if((*itSM)->GetParentList().size() > 0) {     // Don't search in empty vectors
                            if((*itSM)->GetParentList()[0] == bus) {  // Found load
                                syncMotor = *itSM;
                                break;
                            }
                        }
                    }
                    if(!syncMotor) {
                        SyncMotor* newSyncMotor = new SyncMotor(wxString::Format(
                            "%s %u (%s)", _("Synchronous compensator"), syncMotorList.size() + 1, busData->busName));
                        newSyncMotor->AddParent(bus, bus->GetPosition());
                        auto data = newSyncMotor->GetElectricalData();
                        data.activePower = std::abs(busData->genPower.real());
                        data.reactivePower = busData->genPower.imag();
                        newSyncMotor->SetElectricalData(data);
 
                        syncMotorList.push_back(newSyncMotor);
                    }
                }
            }
        }
    }
    // Check for branch data without component (?)
 
    for(auto it = busList.begin(), itEnd = busList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = transformerList.begin(), itEnd = transformerList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = lineList.begin(), itEnd = lineList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = syncGeneratorList.begin(), itEnd = syncGeneratorList.end(); it != itEnd; ++it)
        elementList.push_back(*it);
    for(auto it = syncMotorList.begin(), itEnd = syncMotorList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = loadList.begin(), itEnd = loadList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = indList.begin(), itEnd = indList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = capList.begin(), itEnd = capList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = indMotorList.begin(), itEnd = indMotorList.end(); it != itEnd; ++it) elementList.push_back(*it);
 
    m_workspace->SetElementList(elementList);
    m_workspace->SetName(parseAnarede.GetProjectName());
    return true;
}
 
bool ImportForm::ImportMatpowerFiles()
{
    ParseMatpower parseMatpower(m_filePickerMatpowerM->GetFileName());
    if(!parseMatpower.Parse()) return false;
 
    double mvaBasePower = parseMatpower.GetMVAPowerBase();
    auto simProp = m_workspace->GetProperties()->GetSimulationPropertiesData();
    simProp.basePower = mvaBasePower;
    m_workspace->GetProperties()->SetSimulationPropertiesData(simProp);
 
    std::vector<Element*> elementList;
 
    std::vector<Bus*> busList;
    std::vector<SyncGenerator*> syncGeneratorList;
    std::vector<SyncMotor*> syncMotorList;
    std::vector<Load*> loadList;
    std::vector<Inductor*> indList;
    std::vector<Capacitor*> capList;
    std::vector<IndMotor*> indMotorList;
    std::vector<Transformer*> transformerList;
    std::vector<Line*> lineList;
 
    int iterations = wxAtoi(m_textCtrlIterations->GetValue());
    double scale;
    if(!m_textCtrlScale->GetValue().ToDouble(&scale)) return false;
 
    // Automatically calculate buses positions using weighted graph to determine the layout
    GraphAutoLayout gal(parseMatpower.GetBusData(), parseMatpower.GetBranchData());
    gal.CalculatePositions(iterations, scale);
 
    // Fill bus list
    auto busDataList = parseMatpower.GetBusData();
    for(auto it = busDataList.begin(), itEnd = busDataList.end(); it != itEnd; ++it) {
        ParseMatpower::BusData* busData = *it;
 
        Bus* bus = new Bus(busData->busPosition);
        bus->SetID((*it)->id);
 
        // Electrical data
        auto data = bus->GetElectricalData();
        if(busData) {
            data.number = busData->id;
            data.name = busData->busName;
            switch(busData->type) {
                case 3: {
                    data.isVoltageControlled = true;
                    data.slackBus = true;
                } break;
                case 1: {
                    data.isVoltageControlled = false;
                    data.slackBus = false;
                } break;
                case 2: {
                    data.isVoltageControlled = true;
                    data.slackBus = false;
                } break;
                default: {
                    return false;
                } break;
            }
            data.voltage = std::complex<double>(busData->voltage * std::cos(wxDegToRad(busData->angle)),
                                                busData->voltage * std::sin(wxDegToRad(busData->angle)));
            data.controlledVoltage = busData->voltage;
            // TODO: Nominal voltage are getting errors (e.g. 118bus.m)
            // if(busData->baseVoltage > 1e-3) { data.nominalVoltage = busData->baseVoltage; }
        } else
            return false;
 
        bus->SetElectricalData(data);
        busList.push_back(bus);
    }
    // Fill lines list
    auto branchDataList = parseMatpower.GetBranchData();
    int lineID = 0;
    int transformerID = 0;
    for(auto it = branchDataList.begin(), itEnd = branchDataList.end(); it != itEnd; ++it) {
        ParseMatpower::BranchData* branchData = *it;
        Bus* fstParentBus = GetBusFromID(busList, branchData->busConnections.first);
        Bus* sndParentBus = GetBusFromID(busList, branchData->busConnections.second);
        if(branchData->tap > 1e-3) {  // Transformer
            Transformer* transformer = new Transformer();
            transformer->SetID(transformerID);
            transformer->AddParent(fstParentBus, fstParentBus->GetPosition());
            transformer->AddParent(sndParentBus, sndParentBus->GetPosition());
 
            auto data = transformer->GetElectricalData();
            data.resistance = branchData->resistance;
            data.indReactance = branchData->indReactance;
            data.turnsRatio = branchData->tap;
            data.phaseShift = branchData->phaseShift;
            data.name =
                wxString::Format("%s %u (%s - %s)", _("Transfomer"), transformerList.size() + 1,
                                 fstParentBus->GetElectricalData().name, sndParentBus->GetElectricalData().name);
 
            transformer->SetElectricaData(data);
 
            transformer->SetOnline(branchData->isOnline);
 
            transformerList.push_back(transformer);
            transformerID++;
        } else {  // Line
            Line* line = new Line();
            line->SetID(lineID);
            line->AddParent(fstParentBus, fstParentBus->GetPosition());
            line->AddParent(sndParentBus, sndParentBus->GetPosition());
 
            auto data = line->GetElectricalData();
            data.resistance = branchData->resistance;
            data.indReactance = branchData->indReactance;
            data.capSusceptance = branchData->capSusceptance;
            data.name =
                wxString::Format("%s %u (%s - %s)", _("Line"), lineList.size() + 1,
                                 fstParentBus->GetElectricalData().name, sndParentBus->GetElectricalData().name);
 
            line->SetElectricalData(data);
 
            line->SetOnline(branchData->isOnline);
 
            lineList.push_back(line);
            lineID++;
        }
    }
 
    // Connect Generators
    auto genDataList = parseMatpower.GetGenData();
    int genID = 0;
    for(auto it = genDataList.begin(), itEnd = genDataList.end(); it != itEnd; ++it) {
        SyncGenerator* generator = new SyncGenerator();
        ParseMatpower::GenData* genData = *it;
        generator->SetID(genID);
 
        Bus* parentBus = GetBusFromID(busList, (*it)->busID);
        generator->AddParent(parentBus, parentBus->GetPosition());
 
        auto data = generator->GetElectricalData();
        data.name = wxString::Format("%s %u (%s)", _("Machine"), syncGeneratorList.size() + 1,
                                     parentBus->GetElectricalData().name);
        data.activePower = genData->pg;
        data.reactivePower = genData->qg;
        data.maxReactive = genData->maxReactivePower;
        data.minReactive = genData->minReactivePower;
        data.nominalPower = genData->baseMVA;
        generator->SetElectricalData(data);
 
        syncGeneratorList.push_back(generator);
 
        genID++;
    }
 
    // Connect Loads and capacitors
    for(auto it = busList.begin(), itEnd = busList.end(); it != itEnd; ++it) {
        Bus* bus = *it;
        ParseMatpower::BusData* busData = parseMatpower.GetBusDataFromID(bus->GetID());
        if(busData->pd > 1e-3 || busData->qd > 1e-3) {
            // The bus have load
            Load* load = new Load();
 
            load->AddParent(bus, bus->GetPosition());
 
            auto data = load->GetElectricalData();
            data.name = wxString::Format("%s %u (%s)", _("Load"), loadList.size() + 1, busData->busName);
            data.activePower = busData->pd;
            data.reactivePower = busData->qd;
            load->SetElectricalData(data);
 
            loadList.push_back(load);
        }
        if(std::abs(busData->gs) > 1e-3) {
            // The bus have constant impedance load
            Load* load = new Load();
 
            load->AddParent(bus, bus->GetPosition());
 
            auto data = load->GetElectricalData();
            data.name = wxString::Format("%s %u (%s)", _("Load"), loadList.size() + 1, busData->busName);
            data.activePower = busData->gs;
            data.reactivePower = busData->bs;
            data.loadType = CONST_IMPEDANCE;
            load->SetElectricalData(data);
 
            loadList.push_back(load);
        } else if(std::abs(busData->bs) > 1e-3) {
            // The bus have capacitor or inductor connected
            if(busData->bs < 1e-3) {
                // Inductor
                Inductor* inductor = new Inductor();
 
                inductor->AddParent(bus, bus->GetPosition());
 
                auto data = inductor->GetElectricalData();
                data.name = wxString::Format("%s %u (%s)", _("Inductor"), indList.size() + 1, busData->busName);
                data.reactivePower = std::abs(busData->bs);
                inductor->SetElectricalData(data);
 
                indList.push_back(inductor);
            } else {
                // Capacitor
                Capacitor* capacitor = new Capacitor();
 
                capacitor->AddParent(bus, bus->GetPosition());
 
                auto data = capacitor->GetElectricalData();
                data.name = wxString::Format("%s %u (%s)", _("Capacitor"), capList.size() + 1, busData->busName);
                data.reactivePower = std::abs(busData->bs);
                capacitor->SetElectricalData(data);
 
                capList.push_back(capacitor);
            }
        }
    }
 
    // Adjust generators positions
    for(auto it = syncGeneratorList.begin(), itEnd = syncGeneratorList.end(); it != itEnd; ++it) {
        SyncGenerator* generator = *it;
        generator->StartMove(generator->GetPosition());
        // TODO: Check why node position have 100 pts offset in x axis
        generator->MoveNode(generator->GetParentList()[0],
                            generator->GetParentList()[0]->GetPosition() - wxPoint2DDouble(140, 0));
        generator->Move(generator->GetParentList()[0]->GetPosition() - wxPoint2DDouble(40, 100));
        generator->Rotate();
        generator->Rotate();
    }
 
    // Adjust loads positions
    for(auto it = loadList.begin(), itEnd = loadList.end(); it != itEnd; ++it) {
        Load* load = *it;
        // Move load to the left of the bus
        load->StartMove(wxPoint2DDouble(0, 0));
        load->MoveNode(load->GetParentList()[0], wxPoint2DDouble(-load->GetParentList()[0]->GetWidth() / 2 + 10, 0));
        load->Move(wxPoint2DDouble(-load->GetParentList()[0]->GetWidth() / 2 + 10, 0));
    }
 
    // Adjust capacitors positions
    for(auto it = capList.begin(), itEnd = capList.end(); it != itEnd; ++it) {
        Capacitor* capacitor = *it;
        // Move capacitor to the right of the bus capacitor->StartMove(wxPoint2DDouble(0, 0));
        capacitor->StartMove(wxPoint2DDouble(0, 0));
        capacitor->MoveNode(capacitor->GetParentList()[0],
                            wxPoint2DDouble(capacitor->GetParentList()[0]->GetWidth() / 2 - 20, 0));
        capacitor->Move(wxPoint2DDouble(capacitor->GetParentList()[0]->GetWidth() / 2 - 20, 0));
    }
 
    // Adjust inductors positions
    for(auto it = indList.begin(), itEnd = indList.end(); it != itEnd; ++it) {
        Inductor* inductor = *it;
        // Move indutor to the far right of the bus inductor->StartMove(wxPoint2DDouble(0, 0));
        inductor->StartMove(wxPoint2DDouble(0, 0));
        inductor->MoveNode(inductor->GetParentList()[0],
                           wxPoint2DDouble(inductor->GetParentList()[0]->GetWidth() / 2 - 10, 0));
        inductor->Move(wxPoint2DDouble(inductor->GetParentList()[0]->GetWidth() / 2 + 10, 0));
    }
 
    // Adjust branches
    for(auto it = busList.begin(), itEnd = busList.end(); it != itEnd; ++it) {
        Bus* bus = *it;
        int numberOfConnectedBranches = 0;
        std::vector<Line*> linesConnected;
        std::vector<Transformer*> transformersConnected;
 
        std::vector<Element*> childElements = bus->GetChildList();
        for(unsigned int i = 0; i < childElements.size(); ++i) {
            if(Line* line = dynamic_cast<Line*>(childElements[i])) {
                linesConnected.push_back(line);
                numberOfConnectedBranches++;
            } else if(Transformer* transformer = dynamic_cast<Transformer*>(childElements[i])) {
                transformersConnected.push_back(transformer);
                numberOfConnectedBranches++;
            }
        }
        /*for(auto itc = bus->GetChildList().begin(), itEnd = bus->GetChildList().end(); itc != itEnd; ++itc) {
            if(Line* line = dynamic_cast<Line*>(*itc)) {
                linesConnected.push_back(line);
                numberOfConnectedLines++;
            }
        }*/
 
        if(numberOfConnectedBranches > 0) {
            double dx = (bus->GetWidth() - 30) / (static_cast<double>(numberOfConnectedBranches + 1));
            int cont = 0;
            for(unsigned int i = 0; i < linesConnected.size(); ++i) {
                Line* lineToMove = linesConnected[i];
                // Line nove move in x axis
                wxPoint2DDouble newPos(dx * static_cast<double>(i + 1), 0);
                lineToMove->StartMove(bus->GetPosition());
                lineToMove->MoveNode(bus, bus->GetPosition() - wxPoint2DDouble(bus->GetWidth() / 2 - 10, 0) + newPos);
                cont++;
            }
            for(unsigned int i = 0; i < transformersConnected.size(); ++i) {
                Transformer* trafoToMove = transformersConnected[i];
                wxPoint2DDouble newPos(dx * static_cast<double>(i + cont + 1), 0);
                trafoToMove->StartMove(bus->GetPosition());
                trafoToMove->MoveNode(bus, bus->GetPosition() - wxPoint2DDouble(bus->GetWidth() / 2 - 10, 0) + newPos);
                trafoToMove->SetBestPositionAndRotation();
            }
        }
    }
 
    for(auto it = busList.begin(), itEnd = busList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = transformerList.begin(), itEnd = transformerList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = lineList.begin(), itEnd = lineList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = syncGeneratorList.begin(), itEnd = syncGeneratorList.end(); it != itEnd; ++it)
        elementList.push_back(*it);
    for(auto it = syncMotorList.begin(), itEnd = syncMotorList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = loadList.begin(), itEnd = loadList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = indList.begin(), itEnd = indList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = capList.begin(), itEnd = capList.end(); it != itEnd; ++it) elementList.push_back(*it);
    for(auto it = indMotorList.begin(), itEnd = indMotorList.end(); it != itEnd; ++it) elementList.push_back(*it);
 
    m_workspace->SetElementList(elementList);
 
    return true;
}
 
ParseAnarede::ParseAnarede(wxFileName lstFile, wxFileName pwfFile)
{
    m_lstFile = lstFile;
    m_pwfFile = pwfFile;
}
 
bool ParseAnarede::Parse()
{
    wxTextFile lst(m_lstFile.GetFullPath());
    wxTextFile pwf(m_pwfFile.GetFullPath());
    if(!lst.Open()) return false;
    if(!pwf.Open()) return false;
 
    // Parse LST file
    for(wxString line = lst.GetFirstLine(); !lst.Eof(); line = lst.GetNextLine()) {
        // Current line
        switch(static_cast<char>(line[0])) {
            case 'C': {  // Component
                int parsePosition = 1;
 
                Component* component = new Component();
 
                component->id = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                // Check if is component type is valid
                if(StrToElementType(GetLSTLineNextValue(line, parsePosition), component->type)) {
                    if(component->type == ANA_BUS) {
                        if(!GetLenghtAndRotationFromBusCode(GetLSTLineNextValue(line, parsePosition), component->length,
                                                            component->rotationID))
                            return false;
                    } else {
                        component->rotationID = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                    }
                    parsePosition += 2;  // Jump to position
                    if(!GetLSTLineNextValue(line, parsePosition).ToCDouble(&component->position.m_x)) return false;
                    if(!GetLSTLineNextValue(line, parsePosition).ToCDouble(&component->position.m_y)) return false;
                    parsePosition += 8;  // Jump to electrical ID
                    component->electricalID = wxAtoi(GetLSTLineNextValue(line, parsePosition));
 
                    // Bus connection IDs
                    if(component->type != ANA_BUS) {
                        int fromBus = wxAtoi(GetLSTLineNextValue(line, parsePosition));  // Origin bus
                        int toBus = 0;
                        // If the component is a transformer, parse the next value (toBus)
                        if(component->type == ANA_TRANSFORMER) {
                            toBus = wxAtoi(GetLSTLineNextValue(line, parsePosition));  // Destiny bus
                        }
                        // Iterate through the already parsed components (the buses is saved first)
                        for(auto it = m_components.begin(), itEnd = m_components.end(); it != itEnd; ++it) {
                            if((*it)->type == ANA_BUS) {
                                if(fromBus == (*it)->electricalID) {
                                    component->busConnectionID[0] = std::make_pair((*it)->id, fromBus);
                                } else if(toBus == (*it)->electricalID) {
                                    component->busConnectionID[1] = std::make_pair((*it)->id, toBus);
                                }
                            }
                        }
                    }
                    m_components.push_back(component);
                }
            } break;
            case 'L': {  // Line
                int parsePosition = 1;
 
                PowerLine* pLine = new PowerLine();
 
                pLine->id = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                pLine->type = static_cast<ElementTypeAnarede>(wxAtoi(GetLSTLineNextValue(line, parsePosition)));
                parsePosition += 4;  // Jump to from bus
 
                int fromBus = wxAtoi(GetLSTLineNextValue(line, parsePosition));  // Origin bus
                int toBus = wxAtoi(GetLSTLineNextValue(line, parsePosition));    // Destiny bus
 
                // Iterate through the already parsed components
                for(auto it = m_components.begin(), itEnd = m_components.end(); it != itEnd; ++it) {
                    if((*it)->type == ANA_BUS) {
                        if(fromBus == (*it)->id) {
                            pLine->busConnectionID[0] = std::make_pair((*it)->id, (*it)->electricalID);
                        } else if(toBus == (*it)->id) {
                            pLine->busConnectionID[1] = std::make_pair((*it)->id, (*it)->electricalID);
                        }
                    }
                }
                pLine->electricalID = wxAtoi(GetLSTLineNextValue(line, parsePosition));
 
                m_lines.push_back(pLine);
            } break;
            case 'U': {  // Connection
                // The connection data depends on the component type, so this data definition must be more generic.
                int parsePosition = 1;
 
                int id = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                (void)id;            // Unused id.
                parsePosition += 5;  // Jump to next relevant data
                int data1 = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                int data2 = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                int data3 = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                int data4 = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                int data5 = wxAtoi(GetLSTLineNextValue(line, parsePosition));
                int data6 = wxAtoi(GetLSTLineNextValue(line, parsePosition));
 
                // If data5 differs from -1, the connection is power line data with data5's value as its ID, otherwise
                // data1 is the ID of a component.
                if(data5 != -1) {
                    // Iterate through parsed lines
                    for(auto it = m_lines.begin(), itEnd = m_lines.end(); it != itEnd; ++it) {
                        if(data5 == (*it)->id) {
                            (*it)->busConnectionNode[0] = std::make_pair(data1, data2);
                            (*it)->busConnectionNode[1] = std::make_pair(data3, data4);
                            for(int i = 0; i < data6; ++i) {
                                wxPoint2DDouble nodePt;
                                if(!GetLSTLineNextValue(line, parsePosition).ToCDouble(&nodePt.m_x)) return false;
                                if(!GetLSTLineNextValue(line, parsePosition).ToCDouble(&nodePt.m_y)) return false;
                                (*it)->nodesPosition.push_back(nodePt);
                            }
                        }
                    }
                } else {
                    // Iterate through parsed components
                    for(auto it = m_components.begin(), itEnd = m_components.end(); it != itEnd; ++it) {
                        if(data1 == (*it)->id) {
                            if((*it)->type == ANA_BUS) {
                                // If the data1 is a bus ID, the element is a element with different data order.
                                // Find the correct element ID with data3
                                for(auto itDiff = m_components.begin(), itDiffEnd = m_components.end();
                                    itDiff != itDiffEnd; ++itDiff) {
                                    if(data3 == (*itDiff)->id) {
                                        (*itDiff)->busConnectionNode[data4 - 1] = std::make_pair(data1, data2);
                                        break;
                                    }
                                }
                            } else {
                                (*it)->busConnectionNode[data2 - 1] = std::make_pair(data3, data4);
                                break;
                            }
                        }
                    }
                }
            } break;
            default:
                break;
        }
    }
    // Last line
 
    // Parse PWF file
    for(wxString line = pwf.GetFirstLine(); !pwf.Eof(); line = pwf.GetNextLine()) {
        // Current line
        if(line == "TITU") {  // Title
            m_projectName = pwf.GetNextLine().Trim();
        } else if(line != "") {
            if(line == "FIM") break;
            if(line[0] != '(') {  // Ignore commented line
                wxString exeCode = line;
                wxString data = "";
                for(wxString lineData = pwf.GetNextLine(); lineData != "99999" && !pwf.Eof();
                    lineData = pwf.GetNextLine()) {
                    if(lineData == "FIM") break;
                    if(lineData[0] != '(') {  // Ignore commented line
                        data += lineData + "\n";
                    }
                }
                if(data != "" || data != "FIM") {  // Don't parse empty data
                    if(!ParsePWFExeCode(data, exeCode)) return false;
                }
            }
        }
    }
    // Last line
 
    return true;
}
 
wxString ParseAnarede::GetLSTLineNextValue(wxString line, int& currentPos)
{
    // Parse each line column
    wxString strValue = "=";
    for(; currentPos < static_cast<int>(line.length()); ++currentPos) {
        if(line[currentPos] != ' ') {  // Skip the firsts whitespaces
            wxString parsedLine = line.Mid(currentPos);
            strValue = parsedLine.BeforeFirst(' ');  // Get the character before the first whitespace
            currentPos += strValue.length();         // Set the current parse position
            break;
        }
    }
    return strValue;
}
 
bool ParseAnarede::GetLenghtAndRotationFromBusCode(wxString code, double& lenght, int& rotationID)
{
    // Get code in decimal.
    long longCode;
    if(!code.ToLong(&longCode)) return false;
    // Convert code to hex
    std::stringstream hexCode;
    hexCode << std::hex << longCode;
    // Convert code to string
    wxString hexCodeStr = hexCode.str();
    // Get length on the substring
    wxString lenghtStr = hexCodeStr.Left(hexCodeStr.length() - 1);
    // Convert lenght string (hex) to decimal
    hexCode.str("");
    hexCode << lenghtStr.ToStdString();
    int intLenght;
    hexCode >> std::hex >> intLenght;
    lenght = static_cast<double>(intLenght);
    // Get rotation ID from substring
    wxString rotationIDStr = hexCodeStr.Right(1);
    // Convert rotation ID string to int
    long rotationIDLong;
    if(!rotationIDStr.ToLong(&rotationIDLong)) return false;
    rotationID = static_cast<int>(rotationIDLong);
    // Calculate the true value of lenght
    lenght = 3.0 * lenght - 16.0;
    return true;
}
 
wxPoint2DDouble ParseAnarede::GetNodePositionFromID(Bus* bus, double scale, int nodeID)
{
    wxPoint2DDouble nodePt;
 
    double offset = (static_cast<double>(nodeID) - 1.0) * 16.0;
    nodePt = bus->GetPosition() + wxPoint2DDouble(offset * scale, 0) -
             wxPoint2DDouble(bus->GetWidth() / 2, bus->GetHeight() / 2);
    nodePt = bus->RotateAtPosition(nodePt, bus->GetAngle(), true);
 
    return nodePt;
}
 
bool ParseAnarede::StrToElementType(wxString strType, ElementTypeAnarede& type)
{
    // Check if type exists, othewise return false
    int typeInt = wxAtoi(strType);
    switch(typeInt) {
        case ANA_BUS:
        case ANA_GENERATOR:
        case ANA_LOAD:
        case ANA_SHUNT:
        case ANA_MIT:
        case ANA_TRANSFORMER:
        case ANA_LINE:
        case ANA_IND_LOAD:
        case ANA_IND_SHUNT:
        case ANA_IND_GENERATOR: {
            type = static_cast<ElementTypeAnarede>(typeInt);
            return true;
        } break;
        default:
            break;
    }
    return false;
}
 
bool ParseAnarede::ParsePWFExeCode(wxString data, wxString exeCode)
{
    if(exeCode == "DCTE") {
        // Parse data
        wxString code = "";
        wxString value = "";
        for(unsigned int i = 0; i < data.length(); ++i) {
            code = data.Mid(i, 4);
            value = data.Mid(i + 5, 6).Trim(false);
            i += 11;
            if(code == "BASE") {
                if(!value.ToCDouble(&m_mvaBase)) return false;
                break;
            }
        }
    } else if(exeCode == "DBAR") {
        wxString lineData = "";
        while(data != "") {
            wxString rest = "";
            lineData = data.BeforeFirst('\n', &rest);
            data = rest;
            BusData* busData = new BusData();
 
            if(!GetPWFStructuredData(lineData, 0, 5, busData->id)) return false;
            wxString isOnlineCode = lineData.Mid(6, 1);
            if(isOnlineCode == "L" || isOnlineCode == " ")
                busData->isOnline = true;
            else if(isOnlineCode == "D")
                busData->isOnline = false;
            else
                return false;
            if(!GetPWFStructuredData(lineData, 7, 1, busData->type)) return false;
            busData->voltageBase = lineData.Mid(8, 2);
            busData->busName = lineData.Mid(10, 12).Trim();
            if(!GetPWFStructuredData(lineData, 24, 4, busData->voltage, 1)) return false;
            if(!GetPWFStructuredData(lineData, 28, 4, busData->angle)) return false;
            double pg, qg;
            if(!GetPWFStructuredData(lineData, 32, 5, pg)) return false;
            if(!GetPWFStructuredData(lineData, 37, 5, qg)) return false;
            busData->genPower = std::complex<double>(pg, qg);
            if(!GetPWFStructuredData(lineData, 42, 5, busData->minReactivePower)) return false;
            if(!GetPWFStructuredData(lineData, 47, 5, busData->maxReactivePower)) return false;
            if(!GetPWFStructuredData(lineData, 52, 6, busData->ctrlBusID)) return false;
            double pl, ql;
            if(!GetPWFStructuredData(lineData, 58, 5, pl)) return false;
            if(!GetPWFStructuredData(lineData, 63, 5, ql)) return false;
            busData->loadPower = std::complex<double>(pl, ql);
            if(!GetPWFStructuredData(lineData, 68, 5, busData->shuntReactive)) return false;
 
            m_busData.push_back(busData);
        }
    } else if(exeCode == "DLIN") {
        wxString lineData = "";
        while(data != "") {
            wxString rest = "";
            lineData = data.BeforeFirst('\n', &rest);
            data = rest;
            BranchData* branchData = new BranchData();
 
            if(!GetPWFStructuredData(lineData, 15, 2, branchData->id)) return false;
            int from, to;
            if(!GetPWFStructuredData(lineData, 0, 5, from)) return false;
            if(!GetPWFStructuredData(lineData, 10, 5, to)) return false;
            branchData->busConnections = std::make_pair(from, to);
 
            wxString isOnlineCode = lineData.Mid(6, 1) + lineData.Mid(9, 1) + lineData.Mid(17, 1);
            if(isOnlineCode.Find('D') == wxNOT_FOUND)
                branchData->isOnline = true;
            else
                branchData->isOnline = false;
            if(!GetPWFStructuredData(lineData, 20, 6, branchData->resistance, 4)) return false;
            if(!GetPWFStructuredData(lineData, 26, 6, branchData->indReactance, 4)) return false;
            if(!GetPWFStructuredData(lineData, 32, 6, branchData->capSusceptance, 4)) return false;
            if(!GetPWFStructuredData(lineData, 38, 5, branchData->tap, 2)) return false;
            if(!GetPWFStructuredData(lineData, 53, 5, branchData->phaseShift, 3)) return false;
 
            if(branchData->tap < 1e-3)
                branchData->type = ANA_LINE;
            else
                branchData->type = ANA_TRANSFORMER;
 
            m_branchData.push_back(branchData);
        }
    } else if(exeCode == "DCAI") {
        wxString lineData = "";
        while(data != "") {
            wxString rest = "";
            lineData = data.BeforeFirst('\n', &rest);
            data = rest;
            IndElementData* indData = new IndElementData();
            indData->type = ANA_IND_LOAD;
 
            if(!GetPWFStructuredData(lineData, 9, 2, indData->id)) return false;
            if(!GetPWFStructuredData(lineData, 0, 5, indData->busConnection)) return false;
 
            wxString isOnlineCode = lineData.Mid(13, 1);
            if(isOnlineCode == "L" || isOnlineCode == " ")
                indData->isOnline = true;
            else if(isOnlineCode == "D")
                indData->isOnline = false;
            else
                return false;
 
            double pl, ql;
            if(!GetPWFStructuredData(lineData, 22, 5, pl)) return false;
            if(!GetPWFStructuredData(lineData, 28, 5, ql)) return false;
            indData->power = std::complex<double>(pl, ql);
            // Unities in operation
            if(!GetPWFStructuredData(lineData, 18, 3, indData->numUnits)) return false;
 
            m_indElementData.push_back(indData);
        }
    } else if(exeCode == "DGEI") {
        wxString lineData = "";
        while(data != "") {
            wxString rest = "";
            lineData = data.BeforeFirst('\n', &rest);
            data = rest;
            IndGenData* genData = new IndGenData();
            genData->type = ANA_IND_GENERATOR;
 
            if(!GetPWFStructuredData(lineData, 9, 2, genData->id)) return false;
            if(!GetPWFStructuredData(lineData, 0, 5, genData->busConnection)) return false;
 
            wxString isOnlineCode = lineData.Mid(12, 1);
            if(isOnlineCode == "L" || isOnlineCode == " ")
                genData->isOnline = true;
            else if(isOnlineCode == "D")
                genData->isOnline = false;
            else
                return false;
 
            double pg, qg;
            if(!GetPWFStructuredData(lineData, 22, 5, pg)) return false;
            if(!GetPWFStructuredData(lineData, 27, 5, qg)) return false;
            genData->power = std::complex<double>(pg, qg);
            // Unities in operation
            if(!GetPWFStructuredData(lineData, 16, 3, genData->numUnits)) return false;
            if(!GetPWFStructuredData(lineData, 32, 5, genData->minReactivePower)) return false;
            if(!GetPWFStructuredData(lineData, 37, 5, genData->maxReactivePower)) return false;
            if(!GetPWFStructuredData(lineData, 42, 6, genData->xt, 2)) return false;
            if(!GetPWFStructuredData(lineData, 49, 5, genData->xd, 4)) return false;
            if(!GetPWFStructuredData(lineData, 54, 5, genData->xq, 4)) return false;
            if(!GetPWFStructuredData(lineData, 59, 5, genData->xl, 4)) return false;
            if(!GetPWFStructuredData(lineData, 69, 5, genData->ratedPower, 3)) return false;
 
            m_indElementData.push_back(genData);
        }
    } else if(exeCode == "DBSH") {
        int lineNumber = 1;
        wxString lineData = "";
        int busNumber = -1;
 
        while(data != "") {
            wxString rest = "";
            lineData = data.BeforeFirst('\n', &rest);
            data = rest;
 
            if(lineNumber == 1) {
                if(!GetPWFStructuredData(lineData, 0, 5, busNumber)) return false;
            } else if(lineData == "FBAN") {
                lineNumber = 0;
                busNumber = -1;
            } else if(lineNumber >= 2 && busNumber != -1) {
                IndElementData* indShunt = new IndElementData();
                indShunt->type = ANA_IND_SHUNT;
                indShunt->busConnection = busNumber;
 
                if(!GetPWFStructuredData(lineData, 0, 2, indShunt->id)) return false;
 
                wxString isOnlineCode = lineData.Mid(6, 1);
                if(isOnlineCode == "L" || isOnlineCode == " ")
                    indShunt->isOnline = true;
                else if(isOnlineCode == "D")
                    indShunt->isOnline = false;
                else
                    return false;
 
                if(!GetPWFStructuredData(lineData, 12, 3, indShunt->numUnits)) return false;
                double ql;
                if(!GetPWFStructuredData(lineData, 16, 6, ql)) return false;
                indShunt->power = std::complex<double>(0.0, ql);
                m_indElementData.push_back(indShunt);
            }
            ++lineNumber;
        }
    }
 
    return true;
}
 
bool ParseAnarede::GetPWFStructuredData(wxString data,
                                        unsigned int startPos,
                                        unsigned int dataLenght,
                                        int& value,
                                        int decimalPos)
{
    if(data.length() < startPos || data.length() < (startPos + dataLenght)) {
        value = 0;
        return true;
    }
    wxString strValue = data.Mid(startPos, dataLenght);
    strValue.Replace(' ', '0');
    long lValue = 0;
    if(!strValue.ToLong(&lValue)) return false;
    value = static_cast<int>(lValue);
    return true;
}
 
bool ParseAnarede::GetPWFStructuredData(wxString data,
                                        unsigned int startPos,
                                        unsigned int dataLenght,
                                        double& value,
                                        int decimalPos)
{
    if(data.length() < startPos || data.length() < (startPos + dataLenght)) {
        value = 0.0;
        return true;
    }
    wxString strValue = data.Mid(startPos, dataLenght);
    if(strValue.Find('-') == wxNOT_FOUND) strValue.Replace(' ', '0');
    if(decimalPos != -1 && strValue.Find('.') == wxNOT_FOUND) strValue.insert(decimalPos, '.');
    if(!strValue.ToCDouble(&value)) return false;
    return true;
}
 
ParseAnarede::BusData* ParseAnarede::GetBusDataFromID(int id)
{
    for(auto it = m_busData.begin(), itEnd = m_busData.end(); it != itEnd; ++it) {
        if((*it)->id == id) return *it;
    }
    return nullptr;
}
 
ParseAnarede::BranchData* ParseAnarede::GetBranchDataFromID(int id, int fromBus, int toBus, ElementTypeAnarede type)
{
    for(auto it = m_branchData.begin(), itEnd = m_branchData.end(); it != itEnd; ++it) {
        if((*it)->id == id && (*it)->busConnections.first == fromBus && (*it)->busConnections.second == toBus &&
           (*it)->type == type)
            return *it;
    }
    return nullptr;
}
 
ParseAnarede::IndElementData* ParseAnarede::GetIndElementDataFromID(int id, int busID, ElementTypeAnarede type)
{
    for(auto it = m_indElementData.begin(), itEnd = m_indElementData.end(); it != itEnd; ++it) {
        if((*it)->id == id && (*it)->busConnection == busID && (*it)->type == type) return *it;
    }
    return nullptr;
}
 
void ParseAnarede::ClearData()
{
    for(auto it = m_components.begin(), itEnd = m_components.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_components.clear();
    for(auto it = m_lines.begin(), itEnd = m_lines.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_lines.clear();
    for(auto it = m_busData.begin(), itEnd = m_busData.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_busData.clear();
    for(auto it = m_branchData.begin(), itEnd = m_branchData.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_branchData.clear();
    for(auto it = m_indElementData.begin(), itEnd = m_indElementData.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_indElementData.clear();
}
 
ParseMatpower::ParseMatpower(wxFileName mFile) { m_mFile = mFile; }
 
void ParseMatpower::ClearData()
{
    // Clear data to avoid memory leak
    for(auto it = m_branchData.begin(), itEnd = m_branchData.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_branchData.clear();
    for(auto it = m_branchData.begin(), itEnd = m_branchData.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_branchData.clear();
    for(auto it = m_genData.begin(), itEnd = m_genData.end(); it != itEnd; ++it) {
        if(*it) delete *it;
    }
    m_genData.clear();
}
 
bool ParseMatpower::Parse()
{
    wxTextFile mFile(m_mFile.GetFullPath());
    if(!mFile.Open()) return false;
 
    // Parse M file
    wxString fileTxt = "";
    for(wxString line = mFile.GetFirstLine(); !mFile.Eof(); line = mFile.GetNextLine()) {
        // Current line
        if(line.Trim()[0] != '%') {  // Check if the line is commented (with %)
            // Parse data:
            if(line.Find("mpc.baseMVA") != wxNOT_FOUND) {                                       // Found baseMVA
                if(!line.AfterFirst('=').BeforeFirst(';').ToCDouble(&m_mvaBase)) return false;  // No Trim() needed?
            }
            if(line.Find("mpc.bus ") != wxNOT_FOUND) {  // Found bus
                wxStringTokenizer busStrTok = GetMFileTokenData(mFile, line);
                while(busStrTok.HasMoreTokens()) {
                    BusData* busData = new BusData();
                    wxString busDataStr = busStrTok.GetNextToken();
                    char tokenChar = '\t';
                    if(busDataStr.Find('\t') == wxNOT_FOUND) { tokenChar = ' '; }
                    // Tokenize using tabulation or space (each branch data)
                    wxStringTokenizer busDataStrTok(busDataStr, tokenChar);
 
                    busData->id = wxAtoi(busDataStrTok.GetNextToken());                      // Bus number
                    busData->type = wxAtoi(busDataStrTok.GetNextToken());                    // Bus type
                    if(!busDataStrTok.GetNextToken().ToCDouble(&busData->pd)) return false;  // Real power demand
                    if(!busDataStrTok.GetNextToken().ToCDouble(&busData->qd)) return false;  // Reactive power demand
                    if(!busDataStrTok.GetNextToken().ToCDouble(&busData->gs)) return false;  // Shunt condutance
                    if(!busDataStrTok.GetNextToken().ToCDouble(&busData->bs)) return false;  // Shunt susceptance
                    busData->area = wxAtoi(busDataStrTok.GetNextToken());                    // Bus area
                    if(!busDataStrTok.GetNextToken().ToCDouble(&busData->voltage)) return false;  // Bus abs voltage
                    if(!busDataStrTok.GetNextToken().ToCDouble(&busData->angle)) return false;    // Angle of voltage
                    if(!busDataStrTok.GetNextToken().ToCDouble(&busData->baseVoltage)) return false;  // Base Voltage
 
                    m_busData.push_back(busData);  // Save new bus data at vector
                }
            }
            if(line.Find("mpc.gen ") != wxNOT_FOUND) {  // Found generator
                wxStringTokenizer genStrTok = GetMFileTokenData(mFile, line);
                while(genStrTok.HasMoreTokens()) {
                    GenData* genData = new GenData();
                    wxString genDataStr = genStrTok.GetNextToken();
                    char tokenChar = '\t';
                    if(genDataStr.Find('\t') == wxNOT_FOUND) { tokenChar = ' '; }
                    // Tokenize using tabulation or space (each branch data)
                    wxStringTokenizer genDataStrTok(genDataStr, tokenChar);
                    genData->busID = wxAtoi(genDataStrTok.GetNextToken());                   // Bus number
                    if(!genDataStrTok.GetNextToken().ToCDouble(&genData->pg)) return false;  // Real power output
                    if(!genDataStrTok.GetNextToken().ToCDouble(&genData->qg)) return false;  // Reative power output
                    if(!genDataStrTok.GetNextToken().ToCDouble(&genData->maxReactivePower))
                        return false;  // Maximum reactive power (MVAr)
                    if(!genDataStrTok.GetNextToken().ToCDouble(&genData->minReactivePower))
                        return false;              // Minumum reactive power (MVAr)
                    genDataStrTok.GetNextToken();  // Voltage magnitude setpoint (p.u.), skip
                    if(!genDataStrTok.GetNextToken().ToCDouble(&genData->baseMVA)) return false;  // Power base (MVA)
                    double machineStatus = 0;
                    if(!genDataStrTok.GetNextToken().ToCDouble(&machineStatus))
                        return false;  // Machine status (> 0 = machine in-service; <= 0 = machine out-of-service)
                    genData->isOnline = machineStatus > 1e-3 ? true : false;
 
                    m_genData.push_back(genData);
                }
            }
            if(line.Find("mpc.branch") != wxNOT_FOUND) {  // Found branch
                wxStringTokenizer branchStrTok = GetMFileTokenData(mFile, line);
                while(branchStrTok.HasMoreTokens()) {
                    BranchData* branchData = new BranchData();
                    wxString branchDataStr = branchStrTok.GetNextToken();
                    char tokenChar = '\t';
                    if(branchDataStr.Find('\t') == wxNOT_FOUND) { tokenChar = ' '; }
                    // Tokenize using tabulation or space (each branch data)
                    wxStringTokenizer branchDataStrTok(branchDataStr, tokenChar);
 
                    int fromBus = wxAtoi(branchDataStrTok.GetNextToken());  // From bus number
                    int toBus = wxAtoi(branchDataStrTok.GetNextToken());    // To bus number
                    branchData->busConnections = std::make_pair(fromBus, toBus);
                    if(!branchDataStrTok.GetNextToken().ToCDouble(&branchData->resistance)) return false;  // Resistance
                    if(!branchDataStrTok.GetNextToken().ToCDouble(&branchData->indReactance))
                        return false;  // Reactance
                    if(!branchDataStrTok.GetNextToken().ToCDouble(&branchData->capSusceptance))
                        return false;                 // Line charging susceptance
                    branchDataStrTok.GetNextToken();  // MVA rating A, skip
                    branchDataStrTok.GetNextToken();  // MVA rating B, skip
                    branchDataStrTok.GetNextToken();  // MVA rating C, skip
                    if(!branchDataStrTok.GetNextToken().ToCDouble(&branchData->tap)) return false;  // TAP
                    if(!branchDataStrTok.GetNextToken().ToCDouble(&branchData->phaseShift))
                        return false;  // Transformer phase-shift angle
                    double branchStatus = 0;
                    if(!branchDataStrTok.GetNextToken().ToCDouble(&branchStatus))
                        return false;  // Branch status (> 0 = machine in-service; <= 0 = machine out-of-service)
                    branchData->isOnline = branchStatus > 1e-3 ? true : false;
 
                    m_branchData.push_back(branchData);
                }
            }
            if(line.Find("mpc.bus_name") != wxNOT_FOUND) {  // Found buses names
                wxString dataStr = "";
                while(line.Find('}', true) == wxNOT_FOUND) {  // Concatenate all lines until find ']'
                    dataStr += line;
                    line = mFile.GetNextLine();
                }
                dataStr += line;
                dataStr = dataStr.AfterFirst('{').BeforeFirst('}');  // Get only the element data
                wxStringTokenizer dataStrTok(dataStr, ";");          // Tokenize using ';' char (one element data)
                for(auto it = m_busData.begin(), itEnd = m_busData.end(); it != itEnd; ++it) {
                    if(dataStrTok.HasMoreTokens())
                        (*it)->busName = dataStrTok.GetNextToken().AfterFirst('\'').BeforeFirst('\'');
                }
            }
        }
    }
    // Last line
 
    return true;
}
 
wxStringTokenizer ParseMatpower::GetMFileTokenData(wxTextFile& mFile, wxString currentLine)
{
    wxString dataStr = "";
    while(currentLine.Find(']', true) == wxNOT_FOUND) {  // Concatenate all lines until find ']'
        dataStr += currentLine;
        if(currentLine.Find("mpc") == wxNOT_FOUND && currentLine.Find(';') == wxNOT_FOUND) dataStr += ';';
        currentLine = mFile.GetNextLine();
    }
    dataStr += currentLine + ";";
    if(currentLine.Find(';') == wxNOT_FOUND) dataStr += ';';
    dataStr = dataStr.AfterFirst('[').BeforeFirst(']');  // Get only the element data
    wxStringTokenizer dataStrTok(dataStr, ';');          // Tokenize using ';' char (one element data)
    return dataStrTok;
}
 
ParseMatpower::BusData* ParseMatpower::GetBusDataFromID(int id)
{
    for(auto it = m_busData.begin(), itEnd = m_busData.end(); it != itEnd; ++it) {
        if((*it)->id == id) return *it;
    }
    return nullptr;
}