using SIIPExamples
using PowerSystems
using PowerSimulations
using PowerSystemCaseBuilder

sys = build_system(PSYTestSystems, "tamu_ACTIVSg2000_sys")
transform_single_time_series!(sys, 2, Hour(1))

using Ipopt
solver = optimizer_with_attributes(Ipopt.Optimizer)

print_tree(PowerSimulations.PM.AbstractPowerModel)

ed_template = OperationsProblemTemplate(QCLSPowerModel)
set_device_model!(ed_template, ThermalStandard, ThermalDispatch)
set_device_model!(ed_template, PowerLoad, StaticPowerLoad)
#set_device_model!(ed_template, FixedAdmittance, StaticPowerLoad) #TODO add constructor for shunts in PSI

problem = OperationsProblem(
    ed_template,
    sys,
    horizon = 1,
    optimizer = solver,
    balance_slack_variables = true,
)
build!(problem, output_dir = mktempdir())

solve!(problem)