# python testes/teste13.py
#
# 6- Observar a saída exibida no terminal.
from inewave.newave.pmo import PMO # type: ignore
from inewave.newave.parp import PARp # type: ignore
from inewave.config import REES # type: ignore
from typing import Dict
import numpy as np
from parpa.yulewalker import YuleWalkerPARA # type: ignore
# Variáveis auxiliares no processo
diretorio_parpa = "/home/rogerio/ONS/testes_ft/parpa"
# Lê o arquivo pmo.dat
pmo = PMO.le_arquivo(diretorio_parpa)
# Lê o arquivo parp.dat
parp = PARp.le_arquivo(diretorio_parpa)
# Realiza a verificação de igualdade para todas as
# configurações de todas as REEs.
# Variáveis auxiliares para armazenar valores
IDS_REES = range(1, len(REES) + 1)
# Máxima diferença absoluta por REE
max_dif_ree: Dict[int, float] = {ree: -1e4 for ree in IDS_REES}
ano_max_dif_ree: Dict[int, int] = {ree: 0 for ree in IDS_REES}
periodo_max_dif_ree: Dict[int, int] = {ree: 0 for ree in IDS_REES}
ordem_max_dif_ree: Dict[int, int] = {ree: 0 for ree in IDS_REES}
coef_o_max_dif_ree: Dict[int, float] = {ree: 0 for ree in IDS_REES}
#
# 6- Observar a saída exibida no terminal.
# TODO - Atualizar
from inewave.newave.pmo import PMO
from inewave.newave.parp import PARp
from inewave.nwlistop.mediassin import MediasSIN
from inewave.nwlistop.mediasmerc import MediasMerc
# Variáveis auxiliares no processo
diretorio_v270000 = ""
diretorio_v270403 = ""
# Lê os arquivos de cada diretório
pmo_v270000 = PMO.le_arquivo(diretorio_v270000)
pmo_v270403 = PMO.le_arquivo(diretorio_v270403)
parp_v270000 = PARp.le_arquivo(diretorio_v270000)
parp_v270403 = PARp.le_arquivo(diretorio_v270403)
mediassin_v270000 = MediasSIN.le_arquivo(diretorio_v270000)
mediassin_v270403 = MediasSIN.le_arquivo(diretorio_v270403)
mediasmerc_v270000 = MediasMerc.le_arquivo(diretorio_v270000)
mediasmerc_v270403 = MediasMerc.le_arquivo(diretorio_v270403)
# Compara os dados lidos
pmo_iguais = (pmo_v270000.custo_operacao_series_simuladas ==
pmo_v270403.custo_operacao_series_simuladas)
parp_iguais = (
parp_v270000.coeficientes_ree(1) == parp_v270403.coeficientes_ree(1))
mediassin_iguais = mediassin_v270000.medias == mediassin_v270403.medias
mediasmerc_iguais = mediasmerc_v270000.medias == mediasmerc_v270403.medias
from inewave.newave.dger import DGer # type: ignore
from inewave.newave.pmo import PMO # type: ignore
# OBSERVAÇÃO -
# ESTE TESTE NÃO VAI FUNCIONAR POIS ALGUNS RECURSOS DA
# INEWAVE FORAM TEMPORARIAMENTE DESABILITADOS.
# TODO - RESTAURAR TESTE
# Variáveis auxiliares no processo
diretorio_oficial = ""
diretorio_parpa = ""
# Lê os arquivos
dger_oficial = DGer.le_arquivo(diretorio_oficial)
pmo_oficial = PMO.le_arquivo(diretorio_oficial)
dger_parpa = DGer.le_arquivo(diretorio_parpa)
pmo_parpa = PMO.le_arquivo(diretorio_parpa)
# Compara os dados lidos
eco_igual_oficial = dger_oficial.eq_eco_saida(pmo_oficial.dados_gerais)
flags_dger_oficial = dger_oficial.afluencia_anual_parp
flags_pmo_oficial = pmo_oficial.dados_gerais.afluencia_anual_parp
eco_igual_parpa = dger_parpa.eq_eco_saida(pmo_parpa.dados_gerais)
flags_dger_parpa = dger_parpa.afluencia_anual_parp
flags_pmo_parpa = pmo_parpa.dados_gerais.afluencia_anual_parp
print("PMO Oficial")
print(f"Eco de saída igual: {eco_igual_oficial}")
print(f"Flags no dger: AFL. ANUAL = {flags_dger_oficial[0]}" +
f" RED ORDEM = {flags_dger_oficial[1]}")
import os
import numpy as np
import matplotlib.pyplot as plt # type: ignore
from itertools import combinations
# Variáveis auxiliares no processo
# Consideradas três execuções do NEWAVE, por exemplo:
# A - Máquina 1, Nº Processadores N
# B - Máquina 1, Nº Processadores M
# C - Máquina 2, Nº Processadores N ou M
diretorio_execA = ""
diretorio_execB = ""
diretorio_execC = ""
# Lê os arquivos de cada diretório
pmo_execA = PMO.le_arquivo(diretorio_execA)
pmo_execB = PMO.le_arquivo(diretorio_execB)
pmo_execC = PMO.le_arquivo(diretorio_execC)
parp_execA = PARp.le_arquivo(diretorio_execA)
parp_execB = PARp.le_arquivo(diretorio_execB)
parp_execC = PARp.le_arquivo(diretorio_execC)
mediassin_execA = MediasSIN.le_arquivo(diretorio_execA)
mediassin_execB = MediasSIN.le_arquivo(diretorio_execB)
mediassin_execC = MediasSIN.le_arquivo(diretorio_execC)
mediasmerc_execA = MediasMerc.le_arquivo(diretorio_execA)
mediasmerc_execB = MediasMerc.le_arquivo(diretorio_execB)
mediasmerc_execC = MediasMerc.le_arquivo(diretorio_execC)
# Compara os dados lidos
pmo_iguais = all([
d1 == d2 for d1, d2 in combinations([