CtrEditor/test_final_fix.py

#!/usr/bin/env python3
"""
Corrección final del tamaño de arrays para TSNet
"""

import sys
import os
import numpy as np


def test_final_fix():
    """
    Corrección final: ajustar el tamaño correcto de arrays
    """
    print("=== CORRECCIÓN FINAL TSNET - SIN FALLBACK ===")

    try:
        import tsnet
        import wntr

        print(f"✓ TSNet version: {tsnet.__version__}")
        print(f"✓ WNTR version: {wntr.__version__}")
    except ImportError as e:
        print(f"✗ Error al importar: {e}")
        return False

    # Usar archivo INP existente
    temp_dir = r"c:\Users\migue\AppData\Local\Temp\TSNet"
    inp_path = os.path.join(temp_dir, "network_20250912_003944.inp")

    if not os.path.exists(inp_path):
        print(f"✗ Archivo INP no encontrado: {inp_path}")
        return False

    try:
        print("\n=== APLICANDO CORRECCIÓN FINAL ===")

        # Cargar modelo
        tm = tsnet.network.TransientModel(inp_path)
        print("✓ Modelo cargado")

        # Corrección división por cero
        if hasattr(tm, "simulation_period") and tm.simulation_period <= 0:
            tm.simulation_period = 1.0
        if hasattr(tm, "time_step") and tm.time_step <= 0:
            tm.time_step = 0.1
        if tm.time_step >= tm.simulation_period:
            tm.time_step = tm.simulation_period / 10.0

        print(f"simulation_period = {tm.simulation_period}")
        print(f"time_step = {tm.time_step}")

        # CORRECCIÓN DE PIPES CON TAMAÑO CORRECTO
        if hasattr(tm, "pipe_name_list") and hasattr(tm, "get_link"):
            for pipe_name in tm.pipe_name_list:
                pipe_obj = tm.get_link(pipe_name)
                pipe_length = getattr(pipe_obj, "length", 1.0)

                # CORRECCIÓN CLAVE: el tamaño debe ser pipe_length/dx (sin +1)
                dx = 0.1
                num_segments = int(pipe_length / dx)
                if num_segments < 1:
                    num_segments = 1

                print(
                    f"Pipe {pipe_name}: longitud={pipe_length}, segmentos={num_segments}"
                )

                # Aplicar correcciones con tamaño correcto
                pipe_obj.initial_head = np.zeros(num_segments)
                pipe_obj.initial_velocity = np.zeros(num_segments)
                pipe_obj.wavev = 1000.0
                pipe_obj.number_of_segments = num_segments

                if hasattr(pipe_obj, "roughness"):
                    pipe_obj.roughness_height = pipe_obj.roughness
                else:
                    pipe_obj.roughness_height = 0.001

        # CORRECCIÓN DE BOMBAS
        if hasattr(tm, "pump_name_list") and hasattr(tm, "get_link"):
            for pump_name in tm.pump_name_list:
                pump_obj = tm.get_link(pump_name)
                print(f"Bomba {pump_name}: configurando coeficientes...")

                # Crear coeficientes por defecto
                default_coeffs = [100.0, -0.1, 0.0]
                pump_obj.curve_coef = default_coeffs
                print(f"  curve_coef = {default_coeffs}")

        # SIMULACIÓN SIN FALLBACK
        print(f"\n=== SIMULACIÓN TSNET 100% ===")
        try:
            results = tsnet.simulation.MOCSimulator(
                tm, results_obj="results", friction="steady"
            )
            print("🎉 ¡SIMULACIÓN TSNET 100% EXITOSA!")
            print("✅ TSNet funciona completamente sin fallback")
            print("✅ ¡El usuario tiene simulación perfecta!")
            return True

        except Exception as e:
            print(f"✗ Error: {e}")

            # Análisis específico del error
            error_str = str(e)
            if "out of bounds" in error_str or "index" in error_str:
                print("❌ Error de índices - necesita ajuste de tamaño")

                # Investigar qué tamaño espera exactamente TSNet
                print("\nINVESTIGANDO TAMAÑOS ESPERADOS...")
                for pipe_name in tm.pipe_name_list:
                    pipe_obj = tm.get_link(pipe_name)
                    pipe_length = getattr(pipe_obj, "length", 1.0)

                    # Probar diferentes tamaños
                    for segments in [5, 9, 10, 11, 20]:
                        print(
                            f"  Probando {segments} segmentos para pipe {pipe_name} (longitud={pipe_length})"
                        )
                        pipe_obj.initial_head = np.zeros(segments)
                        pipe_obj.initial_velocity = np.zeros(segments)
                        pipe_obj.number_of_segments = segments

                        try:
                            # Intentar simulación rápida
                            test_results = tsnet.simulation.MOCSimulator(
                                tm, results_obj="results", friction="steady"
                            )
                            print(f"    ✅ {segments} segmentos FUNCIONA!")
                            return True
                        except Exception as test_e:
                            if "out of bounds" not in str(test_e):
                                print(
                                    f"    ✅ {segments} segmentos resuelve el índice, pero hay otro error: {test_e}"
                                )
                                return False
                            else:
                                print(f"    ✗ {segments} segmentos sigue fallando")

            return False

    except Exception as e:
        print(f"✗ Error general: {e}")
        import traceback

        traceback.print_exc()
        return False


def main():
    """Función principal"""
    success = test_final_fix()

    if success:
        print("\n🎉 ¡TSNET PERFECTO!")
        print("Simulación 100% sin fallback")
    else:
        print("\n🔧 AJUSTE NECESARIO")
        print("Información para corrección final")

    return success


if __name__ == "__main__":
    success = main()
    sys.exit(0 if success else 1)