run_file.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import Pkg;\n",
    "Pkg.add(\"JSON\"); Pkg.add(\"Ipopt\"); Pkg.add(\"JuMP\")\n",
    "Pkg.add(\"PolyhedralRelaxations\"); Pkg.add(\"SparseArrays\")\n",
    "Pkg.add(\"OpenBLAS32_jll\")\n",
    "using JuMP, Ipopt\n",
    "using JSON\n",
    "using PolyhedralRelaxations\n",
    "using SparseArrays\n",
    "using Plots"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "mk_output_dir (generic function with 1 method)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Create a function to write output in a unique named folder\n",
    "function mk_output_dir()\n",
    "    i = 1\n",
    "    while true\n",
    "        dir_name = joinpath(@__DIR__, \"./output\", \"run_$i\")\n",
    "        if !ispath(dir_name)\n",
    "            mkpath(dir_name)\n",
    "            return dir_name\n",
    "        end\n",
    "        i += 1\n",
    "    end\n",
    "end"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 8 Node Case 1\n",
    "\n",
    "Without carbon offset"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This is Ipopt version 3.14.4, running with linear solver MUMPS 5.4.1.\n",
      "\n",
      "Number of nonzeros in equality constraint Jacobian...:      143\n",
      "Number of nonzeros in inequality constraint Jacobian.:       12\n",
      "Number of nonzeros in Lagrangian Hessian.............:       82\n",
      "\n",
      "Total number of variables............................:       40\n",
      "                     variables with only lower bounds:       19\n",
      "                variables with lower and upper bounds:       21\n",
      "                     variables with only upper bounds:        0\n",
      "Total number of equality constraints.................:       33\n",
      "Total number of inequality constraints...............:        6\n",
      "        inequality constraints with only lower bounds:        0\n",
      "   inequality constraints with lower and upper bounds:        0\n",
      "        inequality constraints with only upper bounds:        6\n",
      "\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "   0  4.2908673e-01 4.28e-01 7.81e-01  -1.0 0.00e+00    -  0.00e+00 0.00e+00   0\n",
      "   1  4.7397430e-01 4.17e-01 1.69e+00  -1.0 4.28e+00    -  2.04e-02 2.47e-02f  1\n",
      "   2  5.9510912e-01 3.73e-01 2.72e+01  -1.0 3.29e+00    -  4.61e-02 1.05e-01f  1\n",
      "   3  6.0589982e-01 3.60e-01 3.80e+01  -1.0 3.73e-01   2.0 1.37e-01 3.36e-02f  1\n",
      "   4  6.1839627e-01 3.26e-01 5.51e+02  -1.0 3.67e-01   2.4 9.91e-01 9.60e-02f  1\n",
      "   5  6.0937293e-01 1.10e-01 3.32e+02  -1.0 4.13e-01    -  8.65e-01 6.63e-01h  1\n",
      "   6  9.2320818e-01 6.35e-02 1.98e+03  -1.0 1.46e+00    -  6.49e-01 4.22e-01h  1\n",
      "   7  1.2997007e+00 5.74e-02 9.62e+04  -1.0 7.80e+00    -  2.43e-01 9.51e-02f  2\n",
      "   8  1.6172563e+00 5.66e-02 1.93e+05  -1.0 1.20e+02    -  4.24e-02 1.44e-02f  3\n",
      "   9  1.5930645e+00 2.30e-04 3.70e+04  -1.0 6.31e-02   1.9 1.00e+00 9.97e-01h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  10  2.8986159e+00 2.18e-04 4.36e+06  -1.0 9.36e+00    -  2.56e-01 4.00e-01f  1\n",
      "  11  1.9721587e+01 1.24e-02 1.82e+09  -1.0 6.85e+01    -  1.10e-01 1.00e+00f  1\n",
      "  12  1.9366842e+01 4.50e-03 5.66e+08  -1.0 1.37e-01   2.4 4.72e-01 6.45e-01h  1\n",
      "  13  1.9143564e+01 1.02e-04 8.86e+06  -1.0 1.73e-02   1.9 1.00e+00 1.00e+00f  1\n",
      "  14  1.9168899e+01 1.02e-04 4.11e+05  -1.0 2.19e-02   1.4 1.00e+00 1.00e+00h  1\n",
      "  15  1.9218349e+01 4.11e-04 1.73e+04  -1.0 4.37e-02   0.9 1.00e+00 1.00e+00f  1\n",
      "  16  1.9722715e+01 5.02e-04 3.15e+06  -1.7 6.99e-02   0.5 8.31e-01 1.00e+00f  1\n",
      "  17  2.0070944e+01 3.79e-04 8.19e+03  -1.7 2.18e-01  -0.0 1.00e+00 1.00e+00f  1\n",
      "  18  2.0426557e+01 1.70e-04 5.57e+03  -1.7 6.27e-01  -0.5 1.00e+00 1.00e+00f  1\n",
      "  19  2.1465510e+01 4.02e-05 9.07e+03  -1.7 1.85e+00  -1.0 1.00e+00 1.00e+00f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  20  3.2262205e+01 5.01e-03 1.25e+06  -1.7 2.91e+05    -  2.54e-04 1.44e-04f  1\n",
      "  21  4.0766220e+01 6.28e-03 1.48e+06  -1.7 2.05e+04    -  6.02e-03 1.59e-03f  1\n",
      "  22  4.1070532e+01 6.21e-03 1.46e+06  -1.7 4.99e+01    -  1.08e-01 1.23e-02f  1\n",
      "  23  6.5102543e+01 3.40e-02 1.79e+06  -1.7 1.27e+03    -  2.90e-03 4.84e-02f  1\n",
      "  24  6.4127421e+01 5.41e-03 1.42e+05  -1.7 1.74e-01   0.4 1.00e+00 1.00e+00h  1\n",
      "  25  8.5594629e+01 2.11e-02 4.95e+05  -1.7 1.49e+04    -  4.18e-03 3.83e-03f  1\n",
      "  26  8.5815551e+01 2.10e-02 4.92e+05  -1.7 1.22e+02    -  1.00e+00 4.70e-03f  1\n",
      "  27  8.6487618e+01 5.30e-04 1.75e+04  -1.7 1.30e-01    -  1.00e+00 1.00e+00h  1\n",
      "  28  8.6547376e+01 5.12e-05 3.36e+02  -1.7 2.92e-02    -  1.00e+00 1.00e+00h  1\n",
      "  29  8.6766027e+01 7.32e-04 1.28e+05  -2.5 1.08e-01    -  1.00e+00 1.00e+00f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  30  8.6805697e+01 6.16e-06 2.93e+02  -2.5 1.17e-02    -  1.00e+00 1.00e+00h  1\n",
      "  31  8.6805796e+01 6.86e-10 4.87e-03  -2.5 9.56e-05    -  1.00e+00 1.00e+00h  1\n",
      "  32  8.6849066e+01 2.73e-05 7.07e+03  -6.0 2.16e-02    -  9.83e-01 9.94e-01f  1\n",
      "  33  8.6850673e+01 8.72e-09 1.78e-01  -6.0 4.45e-04    -  1.00e+00 1.00e+00h  1\n",
      "  34  8.6850674e+01 4.61e-16 2.59e-10  -6.0 5.50e-09    -  1.00e+00 1.00e+00h  1\n",
      "\n",
      "Number of Iterations....: 34\n",
      "\n",
      "                                   (scaled)                 (unscaled)\n",
      "Objective...............:  -8.6850673742187780e+01    8.6850673742187780e+01\n",
      "Dual infeasibility......:   2.5935962351819207e-10    2.5935962351819207e-10\n",
      "Constraint violation....:   4.6064404528767429e-16    4.6064404528767429e-16\n",
      "Variable bound violation:   7.4283870166169663e-09    7.4283870166169663e-09\n",
      "Complementarity.........:   9.0911235594907303e-07    9.0911235594907303e-07\n",
      "Overall NLP error.......:   9.0911235594907303e-07    9.0911235594907303e-07\n",
      "\n",
      "\n",
      "Number of objective function evaluations             = 40\n",
      "Number of objective gradient evaluations             = 35\n",
      "Number of equality constraint evaluations            = 40\n",
      "Number of inequality constraint evaluations          = 40\n",
      "Number of equality constraint Jacobian evaluations   = 35\n",
      "Number of inequality constraint Jacobian evaluations = 35\n",
      "Number of Lagrangian Hessian evaluations             = 34\n",
      "Total seconds in IPOPT                               = 0.033\n",
      "\n",
      "EXIT: Optimal Solution Found.\n",
      "done\n"
     ]
    }
   ],
   "source": [
    "include(\"./src/NGHTwoSteadyOpt.jl\")\n",
    "\n",
    "# Choose which network to run\n",
    "file = \"./data/8-node/\"\n",
    "include(\"./data/8-node/CASE1.jl\")\n",
    "\n",
    "ss = initialize_optimizer(file, initial_guess_filename=\"\");\n",
    "\n",
    "type = \"nlp_eq\"\n",
    "solver_options = Dict()\n",
    "\n",
    "report,shadowPrice = run_optimizer(ss, type, solver_options);\n",
    "println(\"done\")\n",
    "\n",
    "# SAVE OUTPUT\n",
    "open(\"./output/8-node/CASE1_physical.json\",\"w\") do f\n",
    "    JSON.print(f, report.sol)\n",
    "end\n",
    "open(\"./output/8-node/CASE1_shadowprice.json\",\"w\") do f\n",
    "    JSON.print(f, shadowPrice)\n",
    "end"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 8 Node Case 2\n",
    "\n",
    "With 0.055 carbon offset"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This is Ipopt version 3.14.4, running with linear solver MUMPS 5.4.1.\n",
      "\n",
      "Number of nonzeros in equality constraint Jacobian...:      143\n",
      "Number of nonzeros in inequality constraint Jacobian.:       12\n",
      "Number of nonzeros in Lagrangian Hessian.............:       82\n",
      "\n",
      "Total number of variables............................:       40\n",
      "                     variables with only lower bounds:       19\n",
      "                variables with lower and upper bounds:       21\n",
      "                     variables with only upper bounds:        0\n",
      "Total number of equality constraints.................:       33\n",
      "Total number of inequality constraints...............:        6\n",
      "        inequality constraints with only lower bounds:        0\n",
      "   inequality constraints with lower and upper bounds:        0\n",
      "        inequality constraints with only upper bounds:        6\n",
      "\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "   0  4.2951742e-01 4.28e-01 7.81e-01  -1.0 0.00e+00    -  0.00e+00 0.00e+00   0\n",
      "   1  4.7576723e-01 4.17e-01 1.73e+00  -1.0 4.28e+00    -  2.03e-02 2.47e-02f  1\n",
      "   2  6.0689935e-01 3.74e-01 2.76e+01  -1.0 3.37e+00    -  4.60e-02 1.04e-01f  1\n",
      "   3  6.1856741e-01 3.61e-01 4.31e+01  -1.0 3.74e-01   2.0 1.39e-01 3.36e-02f  1\n",
      "   4  6.3207217e-01 3.26e-01 5.53e+02  -1.0 3.67e-01   2.4 9.91e-01 9.73e-02f  1\n",
      "   5  6.3278830e-01 1.10e-01 3.51e+02  -1.0 4.27e-01    -  8.52e-01 6.63e-01h  1\n",
      "   6  9.7838879e-01 6.46e-02 2.17e+03  -1.0 1.57e+00    -  6.25e-01 4.12e-01h  1\n",
      "   7  1.3899472e+00 5.88e-02 9.42e+04  -1.0 8.57e+00    -  2.11e-01 8.88e-02f  2\n",
      "   8  1.5670609e+00 5.86e-02 1.92e+05  -1.0 2.52e+02    -  2.13e-02 4.08e-03f  4\n",
      "   9  1.5488045e+00 2.29e-04 4.02e+04  -1.0 6.49e-02   1.9 1.00e+00 9.99e-01h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  10  2.7961190e+00 2.51e-04 4.27e+06  -1.0 7.28e+00    -  2.66e-01 4.01e-01f  1\n",
      "  11  1.9719688e+01 1.91e-02 1.73e+09  -1.0 5.88e+01    -  1.02e-01 1.00e+00f  1\n",
      "  12  1.9372365e+01 1.77e-02 1.65e+09  -1.0 2.56e+00   1.5 1.08e-02 3.80e-02h  1\n",
      "  13  1.8402146e+01 1.52e-04 5.68e+07  -1.0 1.33e-01   1.0 6.99e-02 1.00e+00f  1\n",
      "  14  1.8428699e+01 1.57e-05 4.79e+05  -1.0 8.52e-03   2.3 1.00e+00 1.00e+00h  1\n",
      "  15  1.8476893e+01 6.02e-05 2.72e+04  -1.0 1.46e-02   1.8 1.00e+00 1.00e+00f  1\n",
      "  16  1.8582314e+01 2.38e-04 2.32e+04  -1.0 3.05e-02   1.4 1.00e+00 1.00e+00f  1\n",
      "  17  2.6277046e+01 8.35e-03 4.26e+06  -1.7 9.58e+04    -  8.72e-05 1.70e-04f  1\n",
      "  18  2.6665964e+01 5.66e-03 3.60e+06  -1.7 2.69e-02   0.9 7.73e-01 3.22e-01h  1\n",
      "  19  3.9546924e+01 8.43e-03 3.57e+06  -1.7 4.98e+03    -  1.43e-02 1.08e-02f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  20  4.6661309e+01 9.38e-03 3.58e+06  -1.7 6.01e+03    -  2.04e-02 4.41e-03f  1\n",
      "  21  5.3419020e+01 3.92e-01 3.47e+06  -1.7 5.72e+02    -  2.03e-02 2.83e-02f  1\n",
      "  22  9.0748991e+01 7.30e-02 4.99e+06  -1.7 2.03e+03    -  1.54e-02 4.11e-02f  1\n",
      "  23  9.3791056e+01 9.54e+01 4.86e+06  -1.7 2.05e+03    -  7.36e-02 3.83e-03f  1\n",
      "  24  8.9140899e+01 5.49e-04 1.96e+04  -1.7 4.02e-01    -  1.00e+00 1.00e+00h  1\n",
      "  25  8.9262956e+01 1.21e-06 3.49e+01  -1.7 3.81e-01    -  1.00e+00 1.00e+00h  1\n",
      "  26  8.9452150e+01 5.52e-05 5.76e+04  -3.8 1.02e-01    -  9.73e-01 1.00e+00f  1\n",
      "  27  8.9457928e+01 1.46e-07 1.06e+01  -3.8 1.16e-03    -  1.00e+00 1.00e+00h  1\n",
      "  28  8.9457925e+01 1.14e-12 1.03e-05  -3.8 3.67e-06    -  1.00e+00 1.00e+00h  1\n",
      "  29  8.9459420e+01 3.79e-09 2.37e+00  -6.0 7.67e-04    -  1.00e+00 1.00e+00f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  30  8.9459420e+01 3.89e-16 3.40e-08  -6.0 5.87e-08    -  1.00e+00 1.00e+00h  1\n",
      "\n",
      "Number of Iterations....: 30\n",
      "\n",
      "                                   (scaled)                 (unscaled)\n",
      "Objective...............:  -8.9459420101059720e+01    8.9459420101059720e+01\n",
      "Dual infeasibility......:   3.4027794825457151e-08    3.4027794825457151e-08\n",
      "Constraint violation....:   3.8857805861880479e-16    3.8857805861880479e-16\n",
      "Variable bound violation:   6.8305646694355460e-09    6.8305646694355460e-09\n",
      "Complementarity.........:   9.0909091312485644e-07    9.0909091312485644e-07\n",
      "Overall NLP error.......:   9.0909091312485644e-07    9.0909091312485644e-07\n",
      "\n",
      "\n",
      "Number of objective function evaluations             = 37\n",
      "Number of objective gradient evaluations             = 31\n",
      "Number of equality constraint evaluations            = 37\n",
      "Number of inequality constraint evaluations          = 37\n",
      "Number of equality constraint Jacobian evaluations   = 31\n",
      "Number of inequality constraint Jacobian evaluations = 31\n",
      "Number of Lagrangian Hessian evaluations             = 30\n",
      "Total seconds in IPOPT                               = 0.028\n",
      "\n",
      "EXIT: Optimal Solution Found.\n",
      "done\n"
     ]
    }
   ],
   "source": [
    "include(\"./src/NGHTwoSteadyOpt.jl\")\n",
    "\n",
    "# Choose which network to run\n",
    "file = \"./data/8-node/\"\n",
    "include(\"./data/8-node/CASE2.jl\")\n",
    "\n",
    "ss = initialize_optimizer(file, initial_guess_filename=\"\");\n",
    "\n",
    "type = \"nlp_eq\"\n",
    "solver_options = Dict()\n",
    "\n",
    "report,shadowPrice = run_optimizer(ss, type, solver_options);\n",
    "println(\"done\")\n",
    "\n",
    "# SAVE OUTPUT\n",
    "open(\"./output/8-node/CASE2_physical.json\",\"w\") do f\n",
    "    JSON.print(f, report.sol)\n",
    "end\n",
    "open(\"./output/8-node/CASE2_shadowprice.json\",\"w\") do f\n",
    "    JSON.print(f, shadowPrice)\n",
    "end"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 40 Node Baseline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This is Ipopt version 3.14.4, running with linear solver MUMPS 5.4.1.\n",
      "\n",
      "Number of nonzeros in equality constraint Jacobian...:      851\n",
      "Number of nonzeros in inequality constraint Jacobian.:       84\n",
      "Number of nonzeros in Lagrangian Hessian.............:      489\n",
      "\n",
      "Total number of variables............................:      208\n",
      "                     variables with only lower bounds:      111\n",
      "                variables with lower and upper bounds:       97\n",
      "                     variables with only upper bounds:        0\n",
      "Total number of equality constraints.................:      171\n",
      "Total number of inequality constraints...............:       32\n",
      "        inequality constraints with only lower bounds:        0\n",
      "   inequality constraints with lower and upper bounds:        0\n",
      "        inequality constraints with only upper bounds:       32\n",
      "\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "   0  6.9306641e+00 6.30e-01 3.21e+01  -1.0 0.00e+00    -  0.00e+00 0.00e+00   0\n",
      "   1  6.9568993e+00 6.25e-01 3.21e+01  -1.0 1.76e+01    -  1.70e-02 8.20e-03f  1\n",
      "   2  7.0964820e+00 5.90e-01 3.26e+01  -1.0 1.41e+01    -  1.17e-02 5.52e-02f  1\n",
      "   3  7.1930198e+00 5.49e-01 4.40e+01  -1.0 1.12e+00   0.0 9.22e-02 7.06e-02f  1\n",
      "   4  8.2524242e+00 5.29e-01 1.28e+03  -1.0 2.30e+01    -  4.47e-02 3.59e-02f  1\n",
      "   5  1.0413110e+01 5.07e-01 6.41e+03  -1.0 3.45e+01    -  8.54e-02 4.11e-02f  1\n",
      "   6  1.4576586e+01 4.89e-01 1.32e+05  -1.0 5.51e+01    -  6.15e-02 3.55e-02f  1\n",
      "   7  2.7765344e+01 4.66e-01 7.62e+06  -1.0 1.17e+02    -  3.09e-02 4.82e-02f  1\n",
      "   8  5.2287138e+01 4.21e-01 2.01e+06  -1.0 6.72e+01    -  5.51e-02 9.56e-02f  1\n",
      "   9  7.0481690e+01 3.97e-01 1.70e+06  -1.0 1.52e+02    -  2.58e-02 5.60e-02f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  10  1.0896597e+02 3.25e-01 1.68e+06  -1.0 4.11e+01    -  1.19e-01 1.83e-01f  1\n",
      "  11  1.5910153e+02 2.41e-01 2.87e+06  -1.0 3.53e+01    -  4.45e-01 2.57e-01f  1\n",
      "  12  2.6289403e+02 1.40e-01 4.99e+07  -1.0 5.53e+01    -  4.15e-01 4.22e-01f  1\n",
      "  13  2.6062156e+02 3.28e-03 4.62e+08  -1.0 2.42e-01   0.4 9.92e-01 1.00e+00f  1\n",
      "  14  4.6288866e+02 5.41e-02 4.46e+08  -1.0 1.51e+03    -  2.71e-02 3.39e-02f  1\n",
      "  15  5.4370548e+02 5.71e-02 4.29e+08  -1.0 8.95e+02    -  7.18e-02 3.53e-02f  1\n",
      "  16  5.4973455e+02 5.50e-02 4.02e+08  -1.0 1.46e+02    -  1.28e-01 3.70e-02f  1\n",
      "  17  5.6492348e+02 2.58e+00 6.99e+08  -1.0 6.99e+01    -  1.40e-01 2.19e-01f  1\n",
      "  18  5.6376800e+02 3.18e-02 5.14e+08  -1.0 6.15e-01  -0.1 2.75e-01 2.22e-01h  1\n",
      "  19  5.5626942e+02 2.33e-02 1.55e+07  -1.0 7.14e-01  -0.5 1.00e+00 1.00e+00h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  20  5.6020417e+02 6.81e-03 3.61e+06  -1.0 7.84e+00    -  7.73e-01 7.67e-01f  1\n",
      "  21  5.6073058e+02 1.16e-02 1.96e+07  -1.0 1.93e+00    -  2.45e-01 7.18e-01f  1\n",
      "  22  5.5978167e+02 6.23e-03 2.76e+08  -1.0 9.27e-01  -1.0 1.00e+00 1.00e+00f  1\n",
      "  23  5.6067754e+02 5.52e-03 2.30e+08  -1.0 2.08e+00  -0.6 1.89e-01 9.97e-02f  1\n",
      "  24  5.6537270e+02 3.15e-02 2.35e+07  -1.0 4.31e+00  -1.1 3.94e-01 3.12e-01f  1\n",
      "  25  5.6623272e+02 3.03e-02 2.24e+07  -1.0 4.11e+01    -  3.97e-01 3.70e-02f  1\n",
      "  26  5.6589682e+02 1.12e-02 5.73e+07  -1.0 3.41e+01    -  3.86e-01 7.71e-01f  1\n",
      "  27  5.6699188e+02 1.83e-03 1.76e+07  -1.0 2.26e+01    -  1.00e+00 8.91e-01f  1\n",
      "  28  5.6649884e+02 2.83e-03 9.25e+05  -1.0 2.70e+01    -  2.56e-01 1.00e+00f  1\n",
      "  29  5.6658471e+02 4.04e-04 1.34e+05  -1.0 1.63e-01  -1.5 1.00e+00 8.57e-01h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  30  5.6682438e+02 2.10e-03 9.86e+02  -1.0 2.32e+01    -  1.00e+00 1.00e+00f  1\n",
      "  31  5.7086420e+02 4.10e-04 1.11e+06  -2.5 7.94e+00    -  8.54e-01 8.93e-01f  1\n",
      "  32  5.7133557e+02 4.36e-03 1.88e+05  -2.5 1.97e+01    -  8.90e-01 1.00e+00f  1\n",
      "  33  5.7158826e+02 1.20e-03 2.45e+03  -2.5 1.31e+01    -  1.00e+00 1.00e+00h  1\n",
      "  34  5.7166859e+02 1.79e-05 5.44e+00  -2.5 2.67e+00    -  1.00e+00 1.00e+00h  1\n",
      "  35  5.7180302e+02 5.09e-04 3.67e+04  -5.7 1.06e+01    -  8.81e-01 9.84e-01f  1\n",
      "  36  5.7184206e+02 4.94e-05 1.14e+03  -5.7 6.86e+00    -  9.69e-01 1.00e+00h  1\n",
      "  37  5.7184586e+02 1.37e-03 1.55e+00  -5.7 3.18e+00    -  1.00e+00 9.86e-01h  1\n",
      "  38  5.7184587e+02 3.40e-07 5.49e-01  -5.7 3.21e-01    -  8.15e-01 1.00e+00h  1\n",
      "  39  5.7184587e+02 6.74e-08 3.05e-04  -5.7 1.31e-01    -  1.00e+00 1.00e+00h  1\n",
      "\n",
      "Number of Iterations....: 39\n",
      "\n",
      "                                   (scaled)                 (unscaled)\n",
      "Objective...............:  -5.7184587415461147e+02    5.7184587415461147e+02\n",
      "Dual infeasibility......:   3.0469478002743359e-04    3.0469478002743359e-04\n",
      "Constraint violation....:   6.7403505542576236e-08    6.7403505542576236e-08\n",
      "Variable bound violation:   8.0130954094315143e-09    8.0130954094315143e-09\n",
      "Complementarity.........:   5.2307706421198755e-06    5.2307706421198755e-06\n",
      "Overall NLP error.......:   5.2307706421198755e-06    3.0469478002743359e-04\n",
      "\n",
      "\n",
      "Number of objective function evaluations             = 40\n",
      "Number of objective gradient evaluations             = 40\n",
      "Number of equality constraint evaluations            = 40\n",
      "Number of inequality constraint evaluations          = 40\n",
      "Number of equality constraint Jacobian evaluations   = 40\n",
      "Number of inequality constraint Jacobian evaluations = 40\n",
      "Number of Lagrangian Hessian evaluations             = 39\n",
      "Total seconds in IPOPT                               = 0.078\n",
      "\n",
      "EXIT: Optimal Solution Found.\n",
      "done\n"
     ]
    }
   ],
   "source": [
    "include(\"./src/NGHTwoSteadyOpt.jl\")\n",
    "\n",
    "# Choose which network to run\n",
    "file = \"./data/40-node-baseline/\"\n",
    "include(\"./data/40-node-baseline/CASE1.jl\")\n",
    "\n",
    "ss = initialize_optimizer(file, initial_guess_filename=\"\");\n",
    "\n",
    "type = \"nlp_eq\"\n",
    "solver_options = Dict()\n",
    "\n",
    "report,shadowPrice = run_optimizer(ss, type, solver_options);\n",
    "println(\"done\")\n",
    "\n",
    "# SAVE OUTPUT\n",
    "open(\"./output/40-node-baseline/CASE1_physical.json\",\"w\") do f\n",
    "    JSON.print(f, report.sol)\n",
    "end\n",
    "open(\"./output/40-node-baseline/CASE1_shadowprice.json\",\"w\") do f\n",
    "    JSON.print(f, shadowPrice)\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total NG : 726.7070850904207\n",
      "Total H2 : 44.92613179053029\n",
      "Total Flow : 771.6332168809511  771.6332168809507\n",
      "Total Energy : 38490.97864889379  38490.978648893804\n",
      "Average CI : 0.05191981482793796\n",
      "CO2 : 1998.44448399668\n",
      "Decarb premium : 21.79959230869799  21.799592308697985   21.797064481334655\n"
     ]
    }
   ],
   "source": [
    "summed_ng = 0\n",
    "summed_h2 = 0\n",
    "summed_total = 0\n",
    "summed_energy = 0\n",
    "summed_ci = 0\n",
    "summed_dcp = 0\n",
    "summed_λd = 0\n",
    "for i = 1:length(report.sol[:total_ng])\n",
    "    summed_ng += report.sol[:total_ng][i]\n",
    "    summed_h2 += report.sol[:total_h2][i]\n",
    "    summed_total += report.sol[:withdrawal_flows][i]\n",
    "    summed_energy += report.sol[:withdrawal_energy][i]\n",
    "    summed_ci += report.sol[:carbon_intensity][i]\n",
    "    summed_dcp += report.sol[:node_decarbpremium][i]*report.sol[:withdrawal_energy][i]\n",
    "    summed_λd += shadowPrice[:λd][i]*report.sol[:withdrawal_energy][i]\n",
    "end\n",
    "println(\"Total NG : \",summed_ng)\n",
    "println(\"Total H2 : \", summed_h2)\n",
    "println(\"Total Flow : \", summed_ng+summed_h2, \"  \",summed_total)\n",
    "println(\"Total Energy : \", summed_ng*44.2+summed_h2*141.8, \"  \",summed_energy)\n",
    "println(\"Average CI : \", report.sol[:total_co2][1]/summed_energy)\n",
    "println(\"CO2 : \", report.sol[:total_co2][1])\n",
    "println(\"Decarb premium : \", report.sol[:obj_fn_CEM ], \"  \", summed_dcp, \"   \", summed_λd)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 40 Node Case 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This is Ipopt version 3.14.4, running with linear solver MUMPS 5.4.1.\n",
      "\n",
      "Number of nonzeros in equality constraint Jacobian...:      866\n",
      "Number of nonzeros in inequality constraint Jacobian.:       93\n",
      "Number of nonzeros in Lagrangian Hessian.............:      507\n",
      "\n",
      "Total number of variables............................:      211\n",
      "                     variables with only lower bounds:      114\n",
      "                variables with lower and upper bounds:       97\n",
      "                     variables with only upper bounds:        0\n",
      "Total number of equality constraints.................:      171\n",
      "Total number of inequality constraints...............:       35\n",
      "        inequality constraints with only lower bounds:        0\n",
      "   inequality constraints with lower and upper bounds:        0\n",
      "        inequality constraints with only upper bounds:       35\n",
      "\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "   0  7.8650132e+00 6.30e-01 3.21e+01  -1.0 0.00e+00    -  0.00e+00 0.00e+00   0\n",
      "   1  7.9867187e+00 6.25e-01 3.19e+01  -1.0 9.49e+00    -  9.66e-03 8.11e-03f  1\n",
      "   2  9.0601187e+00 5.99e-01 3.17e+01  -1.0 1.84e+01    -  1.14e-02 4.11e-02f  1\n",
      "   3  1.0401222e+01 5.93e-01 3.10e+01  -1.0 8.17e+01    -  9.62e-03 9.84e-03f  1\n",
      "   4  1.7450167e+01 5.53e-01 8.13e+01  -1.0 4.85e+01    -  1.58e-02 6.74e-02f  1\n",
      "   5  1.8111609e+01 3.92e-01 1.54e+02  -1.0 1.20e+00   0.0 2.44e-01 2.91e-01f  1\n",
      "   6  2.0182495e+01 2.71e-01 1.81e+04  -1.0 1.22e+00  -0.5 3.66e-01 3.08e-01f  1\n",
      "   7  3.4505951e+01 2.39e-01 1.30e+04  -1.0 3.46e+01    -  9.46e-02 1.18e-01f  1\n",
      "   8  7.8730692e+01 1.93e-01 4.29e+03  -1.0 5.03e+01    -  1.17e-01 1.95e-01f  1\n",
      "   9  1.9348234e+02 1.37e-01 3.88e+03  -1.0 6.77e+01    -  1.66e-01 2.90e-01f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  10  3.1980039e+02 8.71e-02 4.61e+03  -1.0 5.62e+01    -  3.54e-01 3.63e-01f  1\n",
      "  11  3.2652233e+02 8.38e-02 4.43e+03  -1.0 5.53e+01    -  2.37e-01 3.78e-02f  1\n",
      "  12  4.7049337e+02 5.92e-02 3.54e+03  -1.0 5.52e+01    -  6.32e-02 7.67e-01f  1\n",
      "  13  4.9826307e+02 4.74e-02 2.58e+03  -1.0 2.50e+01    -  1.64e-01 2.13e-01f  1\n",
      "  14  5.7016944e+02 1.10e-01 6.09e+02  -1.0 2.09e+01    -  5.21e-01 6.36e-01f  1\n",
      "  15  5.7334180e+02 9.98e-02 4.89e+02  -1.0 1.89e+01    -  3.42e-01 9.30e-02f  1\n",
      "  16  5.7329503e+02 9.57e-02 8.49e+02  -1.0 2.92e+00  -0.1 3.06e-01 4.10e-02f  1\n",
      "  17  6.0370660e+02 3.75e-02 2.34e+04  -1.0 2.47e+01    -  1.43e-01 6.15e-01f  1\n",
      "  18  6.0797276e+02 3.19e-02 7.60e+03  -1.0 7.28e+00    -  3.24e-01 1.60e-01f  1\n",
      "  19  6.0652030e+02 1.51e-02 1.63e+03  -1.0 1.24e+00  -0.5 3.36e-01 5.43e-01h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  20  6.0858921e+02 1.12e-02 1.03e+03  -1.0 2.32e+00  -1.0 5.19e-01 2.67e-01f  1\n",
      "  21  6.1256779e+02 7.13e-03 6.64e+02  -1.0 1.76e+01    -  9.34e-01 3.61e-01f  1\n",
      "  22  6.1322953e+02 3.80e-05 2.64e+01  -1.0 4.21e+00    -  1.00e+00 1.00e+00H  1\n",
      "  23  6.1316221e+02 2.40e-05 2.04e-01  -1.0 1.70e+00    -  1.00e+00 1.00e+00f  1\n",
      "  24  6.1690939e+02 2.29e-03 1.96e+00  -1.7 4.45e+00    -  8.97e-01 1.00e+00f  1\n",
      "  25  6.1957862e+02 3.53e-04 7.89e+00  -1.7 2.25e+01    -  1.00e+00 9.02e-01F  1\n",
      "  26  6.1902160e+02 2.63e-03 3.56e-02  -1.7 4.40e+00    -  1.00e+00 1.00e+00f  1\n",
      "  27  6.2065676e+02 1.77e-03 1.83e+00  -2.5 8.02e+00    -  8.92e-01 1.00e+00f  1\n",
      "  28  6.2105251e+02 1.40e-05 2.48e-01  -2.5 4.01e+00    -  1.00e+00 1.00e+00h  1\n",
      "  29  6.2105449e+02 2.64e-07 3.81e-02  -2.5 9.83e-01    -  1.00e+00 1.00e+00h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  30  6.2105449e+02 9.18e-10 1.46e-04  -2.5 1.51e-02    -  1.00e+00 1.00e+00h  1\n",
      "  31  6.2136219e+02 5.67e-05 6.58e-01  -3.8 4.18e+00    -  8.27e-01 9.93e-01f  1\n",
      "  32  6.2137633e+02 1.52e-06 9.31e-01  -3.8 3.50e+00    -  1.00e+00 1.00e+00h  1\n",
      "  33  6.2137636e+02 1.97e-07 4.76e-01  -3.8 1.26e+00    -  1.00e+00 1.00e+00h  1\n",
      "  34  6.2137637e+02 6.92e-09 4.56e-02  -3.8 2.28e-01    -  1.00e+00 1.00e+00h  1\n",
      "  35  6.2137637e+02 8.20e-12 1.27e-04  -3.8 7.82e-03    -  1.00e+00 1.00e+00h  1\n",
      "  36  6.2139429e+02 4.55e-07 1.06e+00  -5.7 1.63e+00    -  9.50e-01 1.00e+00f  1\n",
      "  37  6.2139433e+02 1.07e-07 1.29e+00  -5.7 8.40e-01    -  1.00e+00 1.00e+00h  1\n",
      "  38  6.2139433e+02 2.17e-08 1.09e+00  -5.7 3.87e-01    -  1.00e+00 1.00e+00h  1\n",
      "  39  6.2139433e+02 2.83e-09 5.03e-01  -5.7 1.41e-01    -  1.00e+00 1.00e+00h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  40  6.2139433e+02 9.04e-11 4.41e-02  -5.7 2.45e-02    -  1.00e+00 1.00e+00h  1\n",
      "  41  6.2139433e+02 9.13e-14 1.05e-04  -5.7 7.66e-04    -  1.00e+00 1.00e+00h  1\n",
      "  42  6.2139433e+02 1.04e-15 2.16e-10  -5.7 7.38e-07    -  1.00e+00 1.00e+00h  1\n",
      "\n",
      "Number of Iterations....: 42\n",
      "\n",
      "                                   (scaled)                 (unscaled)\n",
      "Objective...............:  -6.2139433213228403e+02    6.2139433213228403e+02\n",
      "Dual infeasibility......:   2.1551627149563046e-10    2.1551627149563046e-10\n",
      "Constraint violation....:   1.0408340855860843e-15    1.0408340855860843e-15\n",
      "Variable bound violation:   7.8813285897538776e-09    7.8813285897538776e-09\n",
      "Complementarity.........:   1.8449145094796160e-06    1.8449145094796160e-06\n",
      "Overall NLP error.......:   1.8449145094796160e-06    1.8449145094796160e-06\n",
      "\n",
      "\n",
      "Number of objective function evaluations             = 45\n",
      "Number of objective gradient evaluations             = 43\n",
      "Number of equality constraint evaluations            = 45\n",
      "Number of inequality constraint evaluations          = 45\n",
      "Number of equality constraint Jacobian evaluations   = 43\n",
      "Number of inequality constraint Jacobian evaluations = 43\n",
      "Number of Lagrangian Hessian evaluations             = 42\n",
      "Total seconds in IPOPT                               = 0.082\n",
      "\n",
      "EXIT: Optimal Solution Found.\n",
      "done\n"
     ]
    }
   ],
   "source": [
    "include(\"./src/NGHTwoSteadyOpt.jl\")\n",
    "\n",
    "# Choose which network to run\n",
    "file = \"./data/40-node-counter/\"\n",
    "include(\"./data/40-node-counter/CASE1.jl\")\n",
    "\n",
    "ss = initialize_optimizer(file, initial_guess_filename=\"\");\n",
    "\n",
    "type = \"nlp_eq\"\n",
    "solver_options = Dict()\n",
    "\n",
    "report,shadowPrice = run_optimizer(ss, type, solver_options);\n",
    "println(\"done\")\n",
    "\n",
    "open(\"./output/40-node-counter/CASE1_physical.json\",\"w\") do f\n",
    "    JSON.print(f, report.sol)\n",
    "end\n",
    "open(\"./output/40-node-counter/CASE1_shadowprice.json\",\"w\") do f\n",
    "    JSON.print(f, shadowPrice)\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total NG : 687.563285066714\n",
      "Total H2 : 76.39491875928509\n",
      "Total Flow : 763.9582038259991  763.958203825999\n",
      "Total Energy : 41223.096680015384  41223.09668001538\n",
      "Average CI : 0.04586746717769282\n",
      "CO2 : 1890.7990339334633\n",
      "Decarb premium : 37.069251614481395  37.06925161448139   37.06641812884611\n"
     ]
    }
   ],
   "source": [
    "summed_ng = 0\n",
    "summed_h2 = 0\n",
    "summed_total = 0\n",
    "summed_energy = 0\n",
    "summed_ci = 0\n",
    "summed_dcp = 0\n",
    "summed_λd = 0\n",
    "for i = 1:length(report.sol[:total_ng])\n",
    "    summed_ng += report.sol[:total_ng][i]\n",
    "    summed_h2 += report.sol[:total_h2][i]\n",
    "    summed_total += report.sol[:withdrawal_flows][i]\n",
    "    summed_energy += report.sol[:withdrawal_energy][i]\n",
    "    summed_ci += report.sol[:carbon_intensity][i]\n",
    "    summed_dcp += report.sol[:node_decarbpremium][i]*report.sol[:withdrawal_energy][i]\n",
    "    summed_λd += shadowPrice[:λd][i]*report.sol[:withdrawal_energy][i]\n",
    "end\n",
    "println(\"Total NG : \",summed_ng)\n",
    "println(\"Total H2 : \", summed_h2)\n",
    "println(\"Total Flow : \", summed_ng+summed_h2, \"  \",summed_total)\n",
    "println(\"Total Energy : \", summed_ng*44.2+summed_h2*141.8, \"  \",summed_energy)\n",
    "println(\"Average CI : \", report.sol[:total_co2][1]/summed_energy)\n",
    "println(\"CO2 : \", report.sol[:total_co2][1])\n",
    "println(\"Decarb premium : \", report.sol[:obj_fn_CEM ], \"  \", summed_dcp, \"   \", summed_λd)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 40 Node Case 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This is Ipopt version 3.14.4, running with linear solver MUMPS 5.4.1.\n",
      "\n",
      "Number of nonzeros in equality constraint Jacobian...:      866\n",
      "Number of nonzeros in inequality constraint Jacobian.:       93\n",
      "Number of nonzeros in Lagrangian Hessian.............:      507\n",
      "\n",
      "Total number of variables............................:      211\n",
      "                     variables with only lower bounds:      114\n",
      "                variables with lower and upper bounds:       97\n",
      "                     variables with only upper bounds:        0\n",
      "Total number of equality constraints.................:      171\n",
      "Total number of inequality constraints...............:       35\n",
      "        inequality constraints with only lower bounds:        0\n",
      "   inequality constraints with lower and upper bounds:        0\n",
      "        inequality constraints with only upper bounds:       35\n",
      "\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "   0  7.1769421e+00 6.30e-01 3.21e+01  -1.0 0.00e+00    -  0.00e+00 0.00e+00   0\n",
      "   1  7.3050139e+00 6.25e-01 3.19e+01  -1.0 1.01e+01    -  9.66e-03 8.12e-03f  1\n",
      "   2  8.3731193e+00 6.00e-01 3.17e+01  -1.0 1.89e+01    -  1.14e-02 4.04e-02f  1\n",
      "   3  9.6818427e+00 5.93e-01 3.09e+01  -1.0 6.71e+01    -  1.22e-02 1.17e-02f  1\n",
      "   4  1.8054410e+01 5.47e-01 7.23e+01  -1.0 5.11e+01    -  1.45e-02 7.74e-02f  1\n",
      "   5  1.8295656e+01 5.33e-01 7.24e+01  -1.0 7.34e+00   0.0 6.68e-02 2.44e-02h  1\n",
      "   6  2.0637458e+01 4.71e-01 1.32e+02  -1.0 1.08e+01  -0.5 2.69e-01 1.17e-01h  1\n",
      "   7  3.0811377e+01 4.29e-01 1.20e+02  -1.0 3.14e+01    -  1.40e-02 8.91e-02f  1\n",
      "   8  5.9890126e+01 3.67e-01 3.90e+03  -1.0 1.30e+02    -  3.81e-02 1.43e-01f  1\n",
      "   9  6.9780274e+01 1.03e-01 7.74e+02  -1.0 3.53e+00  -1.0 2.73e-01 7.21e-01f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  10  9.4951872e+01 9.07e-02 6.50e+02  -1.0 4.05e+01    -  1.64e-01 1.15e-01f  1\n",
      "  11  2.0205797e+02 5.49e-02 2.58e+02  -1.0 4.80e+01    -  1.98e-01 3.94e-01f  1\n",
      "  12  2.0217986e+02 4.38e-02 1.96e+02  -1.0 1.07e+00  -0.5 4.20e-01 2.03e-01f  1\n",
      "  13  2.2936244e+02 4.14e-02 1.77e+02  -1.0 1.50e+02    -  1.72e-02 5.39e-02f  1\n",
      "  14  2.3095404e+02 3.64e-02 1.55e+02  -1.0 2.55e+00  -1.0 5.13e-02 1.22e-01f  1\n",
      "  15  2.3095085e+02 3.63e-02 1.55e+02  -1.0 5.85e+00   0.3 4.24e-04 1.18e-03h  1\n",
      "  16  2.3095275e+02 3.63e-02 1.55e+02  -1.0 4.66e+00   0.7 2.02e-04 6.46e-04h  1\n",
      "  17  2.3095191e+02 3.63e-02 2.45e+02  -1.0 9.86e+01   0.3 1.86e-06 2.46e-05f  1\n",
      "  18  2.3093711e+02 3.63e-02 1.54e+02  -1.0 8.15e+00   1.6 9.86e-03 1.28e-03h  1\n",
      "  19  2.2905174e+02 9.65e-03 2.71e+08  -1.0 1.53e-01   1.1 9.40e-01 7.34e-01f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  20  2.3961077e+02 9.62e-03 2.68e+08  -1.0 7.88e+02    -  3.57e-02 3.03e-03f  1\n",
      "  21  2.3975279e+02 1.93e-03 1.91e+07  -1.0 2.27e-01   0.6 1.00e+00 1.00e+00f  1\n",
      "  22  3.6860882e+02 9.98e-03 1.87e+07  -1.0 2.62e+03    -  1.61e-02 2.10e-02f  1\n",
      "  23  4.3416935e+02 1.13e+01 1.85e+07  -1.0 2.79e+03    -  2.10e-02 1.02e-02f  1\n",
      "  24  4.6843349e+02 9.30e+00 1.83e+07  -1.0 1.90e+03    -  3.88e-02 1.51e-02f  1\n",
      "  25  4.9388845e+02 1.01e+01 1.73e+07  -1.0 3.41e+02    -  7.95e-02 5.20e-02f  1\n",
      "  26  4.9253913e+02 1.19e-03 4.22e+06  -1.0 1.34e+00   0.2 1.00e+00 9.61e-01h  1\n",
      "  27  5.1532102e+02 7.18e-03 7.25e+05  -1.0 1.31e+02    -  1.17e-01 1.33e-01f  1\n",
      "  28  5.4136140e+02 8.51e-03 5.35e+05  -1.0 2.66e+02    -  5.10e-02 9.28e-02f  1\n",
      "  29  5.5708312e+02 2.20e-01 4.89e+05  -1.0 4.36e+02    -  8.47e-02 8.63e-02f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  30  5.6085900e+02 1.31e-02 4.74e+05  -1.0 7.84e+02    -  1.11e-01 3.10e-02f  1\n",
      "  31  5.6395455e+02 1.38e+00 4.35e+05  -1.0 5.30e+01    -  1.43e-01 8.16e-02f  1\n",
      "  32  5.7294559e+02 5.79e-02 2.73e+05  -1.0 4.77e+01    -  7.76e-02 3.72e-01f  1\n",
      "  33  5.7512300e+02 3.80e-02 1.24e+05  -1.0 1.29e+01    -  4.78e-01 5.45e-01f  1\n",
      "  34  5.7470041e+02 3.00e-02 9.80e+04  -1.0 8.74e+01    -  1.00e+00 2.11e-01h  1\n",
      "  35  5.7181771e+02 9.07e-04 7.79e+02  -1.0 7.86e+00    -  1.00e+00 1.00e+00h  1\n",
      "  36  5.7179390e+02 2.90e-05 1.07e+01  -1.0 3.96e+00    -  1.00e+00 1.00e+00h  1\n",
      "  37  5.7178899e+02 7.55e-08 2.89e-03  -1.0 3.48e-01    -  1.00e+00 1.00e+00h  1\n",
      "  38  5.7690918e+02 3.91e-03 3.57e+00  -2.5 1.34e+01    -  7.97e-01 9.36e-01f  1\n",
      "  39  5.7946809e+02 5.49e-03 4.35e+01  -2.5 8.25e+00    -  1.00e+00 1.00e+00f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  40  5.8035440e+02 5.73e-05 1.49e+00  -2.5 9.82e+00    -  1.00e+00 1.00e+00h  1\n",
      "  41  5.8035582e+02 3.70e-07 1.03e-01  -2.5 1.72e+00    -  1.00e+00 1.00e+00h  1\n",
      "  42  5.8035586e+02 4.57e-09 1.80e-03  -2.5 1.91e-01    -  1.00e+00 1.00e+00h  1\n",
      "  43  5.8067395e+02 3.20e-05 9.19e+00  -3.8 6.88e+00    -  9.04e-01 1.00e+00f  1\n",
      "  44  5.8068254e+02 2.25e-06 1.82e+00  -3.8 4.26e+00    -  1.00e+00 1.00e+00h  1\n",
      "  45  5.8068255e+02 2.81e-07 5.68e-01  -3.8 1.50e+00    -  1.00e+00 1.00e+00h  1\n",
      "  46  5.8068255e+02 7.55e-09 4.54e-02  -3.8 2.46e-01    -  1.00e+00 1.00e+00h  1\n",
      "  47  5.8068255e+02 5.95e-12 8.71e-05  -3.8 6.91e-03    -  1.00e+00 1.00e+00h  1\n",
      "  48  5.8070079e+02 5.51e-07 1.02e+01  -5.7 2.10e+00    -  9.72e-01 1.00e+00f  1\n",
      "  49  5.8070081e+02 1.37e-07 1.61e+00  -5.7 1.05e+00    -  1.00e+00 1.00e+00h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  50  5.8070081e+02 2.85e-08 1.35e+00  -5.7 4.78e-01    -  1.00e+00 1.00e+00h  1\n",
      "  51  5.8070081e+02 3.68e-09 6.12e-01  -5.7 1.72e-01    -  1.00e+00 1.00e+00h  1\n",
      "  52  5.8070081e+02 1.05e-10 5.11e-02  -5.7 2.90e-02    -  1.00e+00 1.00e+00h  1\n",
      "  53  5.8070081e+02 9.34e-14 1.10e-04  -5.7 8.65e-04    -  1.00e+00 1.00e+00h  1\n",
      "\n",
      "Number of Iterations....: 53\n",
      "\n",
      "                                   (scaled)                 (unscaled)\n",
      "Objective...............:  -5.8070081364377120e+02    5.8070081364377120e+02\n",
      "Dual infeasibility......:   1.0969056207232875e-04    1.0969056207232875e-04\n",
      "Constraint violation....:   9.3406663452680829e-14    9.3406663452680829e-14\n",
      "Variable bound violation:   7.8757137478291384e-09    7.8757137478291384e-09\n",
      "Complementarity.........:   1.8449203796936499e-06    1.8449203796936499e-06\n",
      "Overall NLP error.......:   1.8449203796936499e-06    1.0969056207232875e-04\n",
      "\n",
      "\n",
      "Number of objective function evaluations             = 54\n",
      "Number of objective gradient evaluations             = 54\n",
      "Number of equality constraint evaluations            = 54\n",
      "Number of inequality constraint evaluations          = 54\n",
      "Number of equality constraint Jacobian evaluations   = 54\n",
      "Number of inequality constraint Jacobian evaluations = 54\n",
      "Number of Lagrangian Hessian evaluations             = 53\n",
      "Total seconds in IPOPT                               = 0.106\n",
      "\n",
      "EXIT: Optimal Solution Found.\n",
      "done\n"
     ]
    }
   ],
   "source": [
    "include(\"./src/NGHTwoSteadyOpt.jl\")\n",
    "\n",
    "# Choose which network to run\n",
    "file = \"./data/40-node-counter/\"\n",
    "include(\"./data/40-node-counter/CASE2.jl\")\n",
    "\n",
    "ss = initialize_optimizer(file, initial_guess_filename=\"\");\n",
    "\n",
    "type = \"nlp_eq\"\n",
    "solver_options = Dict()\n",
    "\n",
    "report,shadowPrice = run_optimizer(ss, type, solver_options);\n",
    "println(\"done\")\n",
    "\n",
    "open(\"./output/40-node-counter/CASE2_physical.json\",\"w\") do f\n",
    "    JSON.print(f, report.sol)\n",
    "end\n",
    "open(\"./output/40-node-counter/CASE2_shadowprice.json\",\"w\") do f\n",
    "    JSON.print(f, shadowPrice)\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total NG : 673.8045240133874\n",
      "Total H2 : 74.86635210717243\n",
      "Total Flow : 748.6708761205599  748.67087612056\n",
      "Total Energy : 40398.20869018878  40398.20869018878\n",
      "Average CI : 0.04586744068894398\n",
      "CO2 : 1852.9624410368153\n",
      "Decarb premium : 36.32754231290846  36.32754231290847   36.32471597287387\n"
     ]
    }
   ],
   "source": [
    "summed_ng = 0\n",
    "summed_h2 = 0\n",
    "summed_total = 0\n",
    "summed_energy = 0\n",
    "summed_ci = 0\n",
    "summed_dcp = 0\n",
    "summed_λd = 0\n",
    "for i = 1:length(report.sol[:total_ng])\n",
    "    summed_ng += report.sol[:total_ng][i]\n",
    "    summed_h2 += report.sol[:total_h2][i]\n",
    "    summed_total += report.sol[:withdrawal_flows][i]\n",
    "    summed_energy += report.sol[:withdrawal_energy][i]\n",
    "    summed_ci += report.sol[:carbon_intensity][i]\n",
    "    summed_dcp += report.sol[:node_decarbpremium][i]*report.sol[:withdrawal_energy][i]\n",
    "    summed_λd += shadowPrice[:λd][i]*report.sol[:withdrawal_energy][i]\n",
    "end\n",
    "println(\"Total NG : \",summed_ng)\n",
    "println(\"Total H2 : \", summed_h2)\n",
    "println(\"Total Flow : \", summed_ng+summed_h2, \"  \",summed_total)\n",
    "println(\"Total Energy : \", summed_ng*44.2+summed_h2*141.8, \"  \",summed_energy)\n",
    "println(\"Average CI : \", report.sol[:total_co2][1]/summed_energy)\n",
    "println(\"CO2 : \", report.sol[:total_co2][1])\n",
    "println(\"Decarb premium : \", report.sol[:obj_fn_CEM ], \"  \", summed_dcp, \"   \", summed_λd)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 40 Node Case 3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This is Ipopt version 3.14.4, running with linear solver MUMPS 5.4.1.\n",
      "\n",
      "Number of nonzeros in equality constraint Jacobian...:      866\n",
      "Number of nonzeros in inequality constraint Jacobian.:       93\n",
      "Number of nonzeros in Lagrangian Hessian.............:      507\n",
      "\n",
      "Total number of variables............................:      211\n",
      "                     variables with only lower bounds:      114\n",
      "                variables with lower and upper bounds:       97\n",
      "                     variables with only upper bounds:        0\n",
      "Total number of equality constraints.................:      171\n",
      "Total number of inequality constraints...............:       35\n",
      "        inequality constraints with only lower bounds:        0\n",
      "   inequality constraints with lower and upper bounds:        0\n",
      "        inequality constraints with only upper bounds:       35\n",
      "\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "   0  7.1864041e+00 6.30e-01 3.22e+01  -1.0 0.00e+00    -  0.00e+00 0.00e+00   0\n",
      "   1  7.3187659e+00 6.25e-01 3.19e+01  -1.0 1.03e+01    -  9.63e-03 8.12e-03f  1\n",
      "   2  8.4563604e+00 6.00e-01 3.18e+01  -1.0 1.86e+01    -  1.14e-02 4.05e-02f  1\n",
      "   3  9.9645614e+00 5.93e-01 3.09e+01  -1.0 8.31e+01    -  9.65e-03 1.09e-02f  1\n",
      "   4  1.7615585e+01 5.55e-01 8.63e+01  -1.0 6.23e+01    -  1.50e-02 6.38e-02f  1\n",
      "   5  1.8450308e+01 3.55e-01 1.55e+02  -1.0 1.05e+00   0.0 2.36e-01 3.61e-01f  1\n",
      "   6  2.0206637e+01 2.71e-01 9.59e+03  -1.0 1.34e+00  -0.5 2.84e-01 2.36e-01f  1\n",
      "   7  3.3663261e+01 2.44e-01 6.22e+03  -1.0 3.37e+01    -  1.04e-01 9.87e-02f  1\n",
      "   8  7.6402648e+01 2.10e-01 1.22e+03  -1.0 6.32e+01    -  9.40e-02 1.39e-01f  1\n",
      "   9  1.8854590e+02 1.58e-01 6.15e+03  -1.0 7.21e+01    -  1.30e-01 2.48e-01f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  10  3.1635831e+02 8.63e-02 5.73e+03  -1.0 5.24e+01    -  2.98e-01 4.54e-01f  1\n",
      "  11  3.5273293e+02 7.35e-02 4.89e+03  -1.0 4.19e+01    -  2.83e-01 1.49e-01f  1\n",
      "  12  3.5815541e+02 7.14e-02 4.75e+03  -1.0 3.99e+01    -  2.62e-02 2.84e-02f  1\n",
      "  13  3.9846393e+02 5.80e-02 2.44e+06  -1.0 3.70e+01    -  5.25e-02 1.87e-01f  1\n",
      "  14  4.0067328e+02 4.99e-02 1.96e+06  -1.0 2.27e+00  -1.0 5.89e-01 1.40e-01f  1\n",
      "  15  5.1358919e+02 1.48e-01 7.25e+06  -1.0 7.44e+01    -  1.74e-01 2.68e-01f  1\n",
      "  16  5.2807337e+02 1.35e-01 6.65e+06  -1.0 6.41e+01    -  1.37e-01 8.29e-02f  1\n",
      "  17  5.6949135e+02 1.05e-01 5.03e+06  -1.0 6.19e+01    -  7.84e-02 2.59e-01f  1\n",
      "  18  6.1288693e+02 5.80e-02 2.86e+06  -1.0 4.51e+01    -  5.09e-01 4.86e-01f  1\n",
      "  19  6.2052864e+02 4.44e-02 2.20e+06  -1.0 2.13e+01    -  7.54e-02 2.43e-01f  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  20  6.2333110e+02 3.51e-02 1.74e+06  -1.0 2.79e+01    -  1.48e-01 2.11e-01f  1\n",
      "  21  6.2379565e+02 3.49e-02 1.74e+06  -1.0 7.85e+01  -1.4 3.23e-04 4.92e-03f  1\n",
      "  22  6.2345060e+02 3.44e-02 1.71e+06  -1.0 2.28e+01  -1.0 4.87e-03 1.52e-02h  1\n",
      "  23  6.2371698e+02 3.44e-02 1.71e+06  -1.0 1.02e+03  -1.5 4.57e-06 2.49e-04f  1\n",
      "  24  6.2328328e+02 3.37e-02 1.67e+06  -1.0 1.26e+01  -0.2 1.78e-01 2.14e-02h  1\n",
      "  25  6.3157126e+02 1.67e-02 8.32e+05  -1.0 4.15e+01    -  6.38e-01 5.04e-01f  1\n",
      "  26  6.3391483e+02 8.28e-03 4.20e+05  -1.0 1.99e+01    -  1.00e+00 5.10e-01f  1\n",
      "  27  6.3561920e+02 1.09e-03 1.38e+04  -1.0 6.36e+00    -  1.00e+00 1.00e+00f  1\n",
      "  28  6.3556852e+02 1.47e-06 5.64e+00  -1.0 1.28e+00    -  1.00e+00 1.00e+00h  1\n",
      "  29  6.3556991e+02 2.73e-08 5.45e-05  -1.0 7.96e-02    -  1.00e+00 1.00e+00h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  30  6.4496380e+02 8.24e-04 6.94e+00  -2.5 1.77e+01    -  8.44e-01 9.21e-01f  1\n",
      "  31  6.4647594e+02 6.07e-01 1.21e+01  -2.5 4.63e+00    -  9.67e-01 8.30e-01f  1\n",
      "  32  6.4680903e+02 1.08e-05 5.49e+00  -2.5 9.35e+00    -  1.00e+00 1.00e+00h  1\n",
      "  33  6.4680947e+02 5.59e-08 1.73e-01  -2.5 6.18e-01    -  1.00e+00 1.00e+00h  1\n",
      "  34  6.4680954e+02 1.09e-09 2.01e-03  -2.5 9.36e-02    -  1.00e+00 1.00e+00h  1\n",
      "  35  6.4713989e+02 4.28e-06 9.34e+00  -3.8 5.83e+00    -  8.83e-01 1.00e+00f  1\n",
      "  36  6.4714116e+02 1.85e-06 5.27e+00  -3.8 3.86e+00    -  1.00e+00 1.00e+00h  1\n",
      "  37  6.4714117e+02 2.39e-07 2.61e+00  -3.8 1.39e+00    -  1.00e+00 1.00e+00h  1\n",
      "  38  6.4714117e+02 6.81e-09 2.21e-01  -3.8 2.34e-01    -  1.00e+00 1.00e+00h  1\n",
      "  39  6.4714117e+02 5.88e-12 4.68e-04  -3.8 6.87e-03    -  1.00e+00 1.00e+00h  1\n",
      "iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du alpha_pr  ls\n",
      "  40  6.4715959e+02 4.45e-07 1.02e+01  -5.7 1.89e+00    -  9.68e-01 1.00e+00f  1\n",
      "  41  6.4715959e+02 1.13e-07 7.28e+00  -5.7 9.52e-01    -  1.00e+00 1.00e+00h  1\n",
      "  42  6.4715959e+02 2.36e-08 6.11e+00  -5.7 4.35e-01    -  1.00e+00 1.00e+00h  1\n",
      "  43  6.4715959e+02 3.07e-09 2.79e+00  -5.7 1.57e-01    -  1.00e+00 1.00e+00h  1\n",
      "  44  6.4715959e+02 8.88e-11 2.35e-01  -5.7 2.67e-02    -  1.00e+00 1.00e+00h  1\n",
      "  45  6.4715959e+02 8.01e-14 5.15e-04  -5.7 8.01e-04    -  1.00e+00 1.00e+00h  1\n",
      "\n",
      "Number of Iterations....: 45\n",
      "\n",
      "                                   (scaled)                 (unscaled)\n",
      "Objective...............:  -6.4715959449650325e+02    6.4715959449650325e+02\n",
      "Dual infeasibility......:   5.1485463582778834e-04    5.1485463582778834e-04\n",
      "Constraint violation....:   8.0072548354420641e-14    8.0072548354420641e-14\n",
      "Variable bound violation:   7.8788296997700513e-09    7.8788296997700513e-09\n",
      "Complementarity.........:   1.8449205853651211e-06    1.8449205853651211e-06\n",
      "Overall NLP error.......:   1.8449205853651211e-06    5.1485463582778834e-04\n",
      "\n",
      "\n",
      "Number of objective function evaluations             = 46\n",
      "Number of objective gradient evaluations             = 46\n",
      "Number of equality constraint evaluations            = 46\n",
      "Number of inequality constraint evaluations          = 46\n",
      "Number of equality constraint Jacobian evaluations   = 46\n",
      "Number of inequality constraint Jacobian evaluations = 46\n",
      "Number of Lagrangian Hessian evaluations             = 45\n",
      "Total seconds in IPOPT                               = 0.092\n",
      "\n",
      "EXIT: Optimal Solution Found.\n",
      "done\n"
     ]
    }
   ],
   "source": [
    "include(\"./src/NGHTwoSteadyOpt.jl\")\n",
    "\n",
    "# Choose which network to run\n",
    "file = \"./data/40-node-counter/\"\n",
    "include(\"./data/40-node-counter/CASE3.jl\")\n",
    "\n",
    "ss = initialize_optimizer(file, initial_guess_filename=\"\");\n",
    "\n",
    "type = \"nlp_eq\"\n",
    "solver_options = Dict()\n",
    "\n",
    "report,shadowPrice = run_optimizer(ss, type, solver_options);\n",
    "println(\"done\")\n",
    "\n",
    "open(\"./output/40-node-counter/CASE3_physical.json\",\"w\") do f\n",
    "    JSON.print(f, report.sol)\n",
    "end\n",
    "open(\"./output/40-node-counter/CASE3_shadowprice.json\",\"w\") do f\n",
    "    JSON.print(f, shadowPrice)\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total NG : 679.3330160531552\n",
      "Total H2 : 75.48126596546152\n",
      "Total Flow : 754.8142820186167  754.814282018617\n",
      "Total Energy : 40729.76282345191  40729.76282345191\n",
      "Average CI : 0.045867337903340334\n",
      "CO2 : 1868.1657941461776\n",
      "Decarb premium : 103.21849655658184  103.2184965565818   103.21566992066724\n"
     ]
    }
   ],
   "source": [
    "summed_ng = 0\n",
    "summed_h2 = 0\n",
    "summed_total = 0\n",
    "summed_energy = 0\n",
    "summed_ci = 0\n",
    "summed_dcp = 0\n",
    "summed_λd = 0\n",
    "for i = 1:length(report.sol[:total_ng])\n",
    "    summed_ng += report.sol[:total_ng][i]\n",
    "    summed_h2 += report.sol[:total_h2][i]\n",
    "    summed_total += report.sol[:withdrawal_flows][i]\n",
    "    summed_energy += report.sol[:withdrawal_energy][i]\n",
    "    summed_ci += report.sol[:carbon_intensity][i]\n",
    "    summed_dcp += report.sol[:node_decarbpremium][i]*report.sol[:withdrawal_energy][i]\n",
    "    summed_λd += shadowPrice[:λd][i]*report.sol[:withdrawal_energy][i]\n",
    "end\n",
    "println(\"Total NG : \",summed_ng)\n",
    "println(\"Total H2 : \", summed_h2)\n",
    "println(\"Total Flow : \", summed_ng+summed_h2, \"  \",summed_total)\n",
    "println(\"Total Energy : \", summed_ng*44.2+summed_h2*141.8, \"  \",summed_energy)\n",
    "println(\"Average CI : \", report.sol[:total_co2][1]/summed_energy)\n",
    "println(\"CO2 : \", report.sol[:total_co2][1])\n",
    "println(\"Decarb premium : \", report.sol[:obj_fn_CEM ], \"  \", summed_dcp, \"   \", summed_λd)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Julia 1.7.3",
   "language": "julia",
   "name": "julia-1.7"
  },
  "language_info": {
   "file_extension": ".jl",
   "mimetype": "application/julia",
   "name": "julia",
   "version": "1.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}