activsg200.jl

using PowerModelsGMDLib
const _PMGLib = PowerModelsGMDLib

using PowerModelsGMD
const _PMG = PowerModelsGMD

import InfrastructureModels
const _IM = InfrastructureModels

import PowerModels
const _PM = PowerModels

import JSON
import JuMP
import Ipopt
import Juniper
import LinearAlgebra
import SparseArrays
using Test
import GZip
import Memento

Memento.setlevel!(Memento.getlogger(_PMG), "error")
Memento.setlevel!(Memento.getlogger(_IM), "error")
Memento.setlevel!(Memento.getlogger(_PM), "error")

_PMG.logger_config!("error")
const TESTLOG = Memento.getlogger(_PMG)
Memento.setlevel!(TESTLOG, "error")

ipopt_solver = JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol" => 1e-4, "print_level" => 0, "sb" => "yes")
juniper_solver = JuMP.optimizer_with_attributes(Juniper.Optimizer, "nl_solver" => _PM.optimizer_with_attributes(Ipopt.Optimizer, "tol" => 1e-4, "print_level" => 0, "sb" => "yes"), "log_levels" => [])
setting = Dict{String,Any}("output" => Dict{String,Any}("branch_flows" => true))

file = "../data/ACTIVSg200/activsg200.m.gz"
case = _PMGLib.parse_file(file)
_PMG.add_gmd_3w_branch!(case)

sol= _PMG.solve_gmd(case) # linear solver
# sol=  _PMG.solve_gmd(case, ipopt_solver; setting=setting) # for opt solver

high_error = 1e-2 # abs(value) >= .0001
low_error = 1 # abs(value) < .0001

@testset "solve of gmd" begin
    @testset "dc bus voltage" begin
        @test isapprox(sol["solution"]["gmd_bus"]["1"]["gmd_vdc"], -0.00038335, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["2"]["gmd_vdc"], -0.00042969, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["3"]["gmd_vdc"], 1.582e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["4"]["gmd_vdc"], -1.997e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["5"]["gmd_vdc"], -0.00244804, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["6"]["gmd_vdc"], 0.00014362, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["7"]["gmd_vdc"], -0.91744745, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["8"]["gmd_vdc"], -0.00350741, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["9"]["gmd_vdc"], -0.00015459, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["10"]["gmd_vdc"], 0.00240021, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["11"]["gmd_vdc"], 0.00219002, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["12"]["gmd_vdc"], -1.62e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["13"]["gmd_vdc"], 5.891e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["14"]["gmd_vdc"], -0.00016201, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["15"]["gmd_vdc"], -0.00107522, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["16"]["gmd_vdc"], -0.00063821, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["17"]["gmd_vdc"], 7.94e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["18"]["gmd_vdc"], -0.00032006, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["19"]["gmd_vdc"], 6.227e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["20"]["gmd_vdc"], -3.141e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["21"]["gmd_vdc"], 3.73771811, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["22"]["gmd_vdc"], -0.00072417, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["23"]["gmd_vdc"], 0.54608577, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["24"]["gmd_vdc"], -0.0001252, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["25"]["gmd_vdc"], 2.98244643, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["26"]["gmd_vdc"], -0.00197493, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["27"]["gmd_vdc"], -1.15e-06, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["28"]["gmd_vdc"], -2.754e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["29"]["gmd_vdc"], -4.07e-06, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["30"]["gmd_vdc"], -4.92e-06, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["31"]["gmd_vdc"], -1.198e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["32"]["gmd_vdc"], -8.091e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["33"]["gmd_vdc"], 4.50473213, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["34"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["35"]["gmd_vdc"], -2.644e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["36"]["gmd_vdc"], -3.65473819, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["37"]["gmd_vdc"], -0.00023495, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["38"]["gmd_vdc"], 1.44499445, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["39"]["gmd_vdc"], -1.92e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["40"]["gmd_vdc"], 3.887e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["41"]["gmd_vdc"], -2.175e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["42"]["gmd_vdc"], -8.4e-06, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["43"]["gmd_vdc"], -15.05074883, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["44"]["gmd_vdc"], -3.43058801, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["45"]["gmd_vdc"], 5.9466691, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["46"]["gmd_vdc"], -0.00013385, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["47"]["gmd_vdc"], 2.059e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["48"]["gmd_vdc"], -0.00010929, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["49"]["gmd_vdc"], -4.78124809, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["50"]["gmd_vdc"], 0.00212545, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["51"]["gmd_vdc"], -0.00016913, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["52"]["gmd_vdc"], 2.84264874, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["53"]["gmd_vdc"], 0.00017216, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["54"]["gmd_vdc"], 1.094e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["55"]["gmd_vdc"], -0.00043891, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["56"]["gmd_vdc"], -0.00100942, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["57"]["gmd_vdc"], -2.425e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["58"]["gmd_vdc"], 0.00014033, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["59"]["gmd_vdc"], -1.40395248, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["60"]["gmd_vdc"], -2.07247853, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["61"]["gmd_vdc"], -0.00031819, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["62"]["gmd_vdc"], -4.07e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["63"]["gmd_vdc"], -0.00024564, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["64"]["gmd_vdc"], -2.97019434, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["65"]["gmd_vdc"], -4.06109142, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["66"]["gmd_vdc"], 1.24523878, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["67"]["gmd_vdc"], 0.0001572, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["68"]["gmd_vdc"], -0.00013163, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["69"]["gmd_vdc"], 2.55e-06, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["70"]["gmd_vdc"], -8.71586514, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["71"]["gmd_vdc"], -0.00013669, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["72"]["gmd_vdc"], 8.531e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["73"]["gmd_vdc"], -0.0001826, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["74"]["gmd_vdc"], 5.394e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["75"]["gmd_vdc"], -2.067e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["76"]["gmd_vdc"], -1.831e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["77"]["gmd_vdc"], 0.00015058, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["78"]["gmd_vdc"], -2.054e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["79"]["gmd_vdc"], -0.00053575, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["80"]["gmd_vdc"], -15.89966679, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["81"]["gmd_vdc"], 0.00015061, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["82"]["gmd_vdc"], -9.33065987, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["83"]["gmd_vdc"], -2.71218681, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["84"]["gmd_vdc"], -3.692e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["85"]["gmd_vdc"], -6.059e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["86"]["gmd_vdc"], -4.77413368, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["87"]["gmd_vdc"], 9.969e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["88"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["89"]["gmd_vdc"], -0.00017627, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["90"]["gmd_vdc"], 0.0002242, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["91"]["gmd_vdc"], 0.00012665, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["92"]["gmd_vdc"], 0.00102077, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["93"]["gmd_vdc"], 0.00013173, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["94"]["gmd_vdc"], -0.00245528, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["95"]["gmd_vdc"], -2.081e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["96"]["gmd_vdc"], 16.7992363, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["97"]["gmd_vdc"], 0.00011209, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["98"]["gmd_vdc"], -3.62125754, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["99"]["gmd_vdc"], 7.166e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["100"]["gmd_vdc"], -0.00058082, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["101"]["gmd_vdc"], -2.948e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["102"]["gmd_vdc"], -3.70574188, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["103"]["gmd_vdc"], -4.166e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["104"]["gmd_vdc"], 0.00045644, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["105"]["gmd_vdc"], 1.519e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["106"]["gmd_vdc"], -0.00010309, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["107"]["gmd_vdc"], -8.778e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["108"]["gmd_vdc"], 0.60580927, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["109"]["gmd_vdc"], -4.193e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["110"]["gmd_vdc"], 0.00396924, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["111"]["gmd_vdc"], -7.069e-05, rtol=low_error)
        @test isapprox(sol["solution"]["gmd_bus"]["112"]["gmd_vdc"], -4.49861383, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["113"]["gmd_vdc"], -4.49860859, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["114"]["gmd_vdc"], -10.84522057, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["115"]["gmd_vdc"], -10.84521008, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["116"]["gmd_vdc"], 2.488832, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["117"]["gmd_vdc"], 2.4888289, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["118"]["gmd_vdc"], -1.3866359, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["119"]["gmd_vdc"], -1.38663435, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["120"]["gmd_vdc"], -6.52635336, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["121"]["gmd_vdc"], -6.52634668, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["122"]["gmd_vdc"], 3.62638617, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["123"]["gmd_vdc"], 3.62638259, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["124"]["gmd_vdc"], 3.62638211, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["125"]["gmd_vdc"], -0.18062465, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["126"]["gmd_vdc"], -1.27687013, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["127"]["gmd_vdc"], -1.27686977, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["128"]["gmd_vdc"], -17.83521652, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["129"]["gmd_vdc"], -17.83520126, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["130"]["gmd_vdc"], -17.83520126, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["131"]["gmd_vdc"], -17.83519745, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["132"]["gmd_vdc"], -17.83519745, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["133"]["gmd_vdc"], -5.53569412, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["134"]["gmd_vdc"], -5.53568792, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["135"]["gmd_vdc"], -5.53568745, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["136"]["gmd_vdc"], 23.00612831, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["137"]["gmd_vdc"], 23.00610352, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["138"]["gmd_vdc"], 7.1170063, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["139"]["gmd_vdc"], 7.11699867, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["140"]["gmd_vdc"], -1.78770936, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["141"]["gmd_vdc"], -1.78770709, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["142"]["gmd_vdc"], 0.73229289, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["143"]["gmd_vdc"], 0.73229206, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["144"]["gmd_vdc"], 0.73229223, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["145"]["gmd_vdc"], -7.51812363, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["146"]["gmd_vdc"], -7.518116, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["147"]["gmd_vdc"], -10.73323727, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["148"]["gmd_vdc"], -10.73322678, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["149"]["gmd_vdc"], -10.73322487, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["150"]["gmd_vdc"], -15.49197388, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["151"]["gmd_vdc"], -15.49195957, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["152"]["gmd_vdc"], 14.90454006, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["153"]["gmd_vdc"], -10.2237978, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["154"]["gmd_vdc"], 6.94479656, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["155"]["gmd_vdc"], -30.23489571, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["156"]["gmd_vdc"], 4.81281519, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["157"]["gmd_vdc"], 4.29150248, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["158"]["gmd_vdc"], -8.68870163, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["159"]["gmd_vdc"], 0.72702789, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["160"]["gmd_vdc"], 0.54608577, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["161"]["gmd_vdc"], 0.54608577, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["162"]["gmd_vdc"], 0.54608577, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["163"]["gmd_vdc"], 0.54608577, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["164"]["gmd_vdc"], 0.54608577, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["165"]["gmd_vdc"], -9.40817451, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["166"]["gmd_vdc"], 7.46221161, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["167"]["gmd_vdc"], 6.95648432, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["168"]["gmd_vdc"], -4.242486, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["169"]["gmd_vdc"], -0.17857276, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["170"]["gmd_vdc"], -3.65038824, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["171"]["gmd_vdc"], -8.06392288, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["172"]["gmd_vdc"], -1.3050983, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["173"]["gmd_vdc"], -1.62395358, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["174"]["gmd_vdc"], -7.77426434, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["175"]["gmd_vdc"], 5.99116039, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["176"]["gmd_vdc"], 4.50473213, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["177"]["gmd_vdc"], -4.32405663, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["178"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["179"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["180"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["181"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["182"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["183"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["184"]["gmd_vdc"], -4.19008255, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["185"]["gmd_vdc"], -6.99832535, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["186"]["gmd_vdc"], -5.89647579, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["187"]["gmd_vdc"], -3.65473819, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["188"]["gmd_vdc"], -3.65473819, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["189"]["gmd_vdc"], -3.65473819, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["190"]["gmd_vdc"], -3.65473819, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["191"]["gmd_vdc"], -13.42639637, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["192"]["gmd_vdc"], 5.94141722, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["193"]["gmd_vdc"], 3.06791472, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["194"]["gmd_vdc"], -6.47550631, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["195"]["gmd_vdc"], 9.87188053, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["196"]["gmd_vdc"], -27.86432076, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["197"]["gmd_vdc"], -0.47087705, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["198"]["gmd_vdc"], -15.64640427, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["199"]["gmd_vdc"], -15.30696297, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["200"]["gmd_vdc"], -5.08353567, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["201"]["gmd_vdc"], -3.43058801, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["202"]["gmd_vdc"], -3.43058801, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["203"]["gmd_vdc"], -3.43058801, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["204"]["gmd_vdc"], 14.04657364, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["205"]["gmd_vdc"], 5.9466691, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["206"]["gmd_vdc"], -10.78509331, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["207"]["gmd_vdc"], 10.55762386, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["208"]["gmd_vdc"], -27.74376678, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["209"]["gmd_vdc"], -5.18367338, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["210"]["gmd_vdc"], -5.187603, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["211"]["gmd_vdc"], 14.52148914, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["212"]["gmd_vdc"], -4.95997429, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["213"]["gmd_vdc"], 3.72418976, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["214"]["gmd_vdc"], 3.85200214, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["215"]["gmd_vdc"], 2.84264874, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["216"]["gmd_vdc"], 2.84264874, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["217"]["gmd_vdc"], 4.55678511, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["218"]["gmd_vdc"], 5.58110619, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["219"]["gmd_vdc"], -17.62033653, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["220"]["gmd_vdc"], -9.95981407, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["221"]["gmd_vdc"], -4.53042126, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["222"]["gmd_vdc"], 10.89848328, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["223"]["gmd_vdc"], -2.46164012, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["224"]["gmd_vdc"], 0.2449422, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["225"]["gmd_vdc"], -1.40395248, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["226"]["gmd_vdc"], -1.40395248, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["227"]["gmd_vdc"], -2.52953553, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["228"]["gmd_vdc"], -2.90283608, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["229"]["gmd_vdc"], -16.67985153, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["230"]["gmd_vdc"], -7.61108589, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["231"]["gmd_vdc"], -32.44390488, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["232"]["gmd_vdc"], -3.47038746, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["233"]["gmd_vdc"], -3.19161725, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["234"]["gmd_vdc"], -4.62116861, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["235"]["gmd_vdc"], -4.92776871, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["236"]["gmd_vdc"], -4.06109142, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["237"]["gmd_vdc"], -4.06109142, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["238"]["gmd_vdc"], -4.06109142, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["239"]["gmd_vdc"], 1.19278777, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["240"]["gmd_vdc"], 1.61868978, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["241"]["gmd_vdc"], 18.29226685, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["242"]["gmd_vdc"], -8.70999622, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["243"]["gmd_vdc"], 0.89899296, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["244"]["gmd_vdc"], -8.99282265, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["245"]["gmd_vdc"], -9.04585457, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["246"]["gmd_vdc"], -8.71586514, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["247"]["gmd_vdc"], -8.71586514, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["248"]["gmd_vdc"], -9.42301369, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["249"]["gmd_vdc"], 2.25957966, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["250"]["gmd_vdc"], -7.22554398, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["251"]["gmd_vdc"], 8.81908226, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["252"]["gmd_vdc"], -3.73242521, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["253"]["gmd_vdc"], -2.13307357, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["254"]["gmd_vdc"], 30.97347641, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["255"]["gmd_vdc"], -0.8164289, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["256"]["gmd_vdc"], -27.94487762, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["257"]["gmd_vdc"], -16.08408165, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["258"]["gmd_vdc"], -15.89966679, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["259"]["gmd_vdc"], 12.19706917, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["260"]["gmd_vdc"], -9.82914734, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["261"]["gmd_vdc"], -9.20613289, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["262"]["gmd_vdc"], -9.33065987, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["263"]["gmd_vdc"], -9.33065987, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["264"]["gmd_vdc"], -9.33065987, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["265"]["gmd_vdc"], -9.33065987, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["266"]["gmd_vdc"], -9.33065987, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["267"]["gmd_vdc"], -4.42287302, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["268"]["gmd_vdc"], -3.10276675, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["269"]["gmd_vdc"], -5.28528786, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["270"]["gmd_vdc"], -5.51869869, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["271"]["gmd_vdc"], -4.99704885, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["272"]["gmd_vdc"], -4.77413368, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["273"]["gmd_vdc"], 8.8538599, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["274"]["gmd_vdc"], 10.26072121, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["275"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["276"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["277"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["278"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["279"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["280"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["281"]["gmd_vdc"], 9.21918011, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["282"]["gmd_vdc"], -16.60688019, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["283"]["gmd_vdc"], 33.06455612, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["284"]["gmd_vdc"], 16.03970146, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["285"]["gmd_vdc"], 5.32272434, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["286"]["gmd_vdc"], 23.2712822, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["287"]["gmd_vdc"], -19.31587219, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["288"]["gmd_vdc"], -6.89282417, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["289"]["gmd_vdc"], 19.5370388, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["290"]["gmd_vdc"], 18.50569534, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["291"]["gmd_vdc"], 25.26343155, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["292"]["gmd_vdc"], -6.35538769, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["293"]["gmd_vdc"], -3.62125754, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["294"]["gmd_vdc"], -3.62125754, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["295"]["gmd_vdc"], 10.4335947, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["296"]["gmd_vdc"], -15.21825695, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["297"]["gmd_vdc"], -6.4817028, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["298"]["gmd_vdc"], -3.84731245, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["299"]["gmd_vdc"], -4.06449556, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["300"]["gmd_vdc"], -3.70574188, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["301"]["gmd_vdc"], -12.53932381, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["302"]["gmd_vdc"], 45.75274277, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["303"]["gmd_vdc"], 2.99225783, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["304"]["gmd_vdc"], -5.6210413, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["305"]["gmd_vdc"], -9.55320358, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["306"]["gmd_vdc"], 0.8253988, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["307"]["gmd_vdc"], 0.60580927, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["308"]["gmd_vdc"], 0.60580927, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["309"]["gmd_vdc"], -10.852911, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["310"]["gmd_vdc"], 24.91715431, rtol=high_error)
        @test isapprox(sol["solution"]["gmd_bus"]["311"]["gmd_vdc"], -35.51470566, rtol=high_error)
    end
    @testset "auto transformers" begin
        @test isapprox(sol["solution"]["qloss"]["25"], 0.9564708285775501, rtol=high_error) || isapprox(sol["solution"]["qloss"]["25"], 0.9585803326291501, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["25"], 2.1854641400, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["70"], 0.9429352963637999, rtol=high_error) || isapprox(sol["solution"]["qloss"]["70"], 0.9406617915937199, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["70"], 2.1704027700, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["85"], 0.7122036585233702, rtol=high_error) || isapprox(sol["solution"]["qloss"]["85"], 0.710904929765355, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["85"], 1.6438099100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["122"], 3.7923462314880445, rtol=high_error) || isapprox(sol["solution"]["qloss"]["122"], 3.7941088093373847, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["122"], 8.7627449000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["132"], 0.58922885940552, rtol=high_error) || isapprox(sol["solution"]["qloss"]["132"], 0.5869732110714, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["132"], 1.3535671200, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["148"], 0.5058792309748501, rtol=high_error) || isapprox(sol["solution"]["qloss"]["148"], 0.50242476280455, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["148"], 1.1832468500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["157"], 0.423088972210215, rtol=high_error) || isapprox(sol["solution"]["qloss"]["157"], 0.4216124139668551, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["157"], 0.9752893400, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["171"], 0.4153595083644601, rtol=high_error) || isapprox(sol["solution"]["qloss"]["171"], 0.4071476586553201, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["171"], 0.9626939900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["176"], 0.7135876425795751, rtol=high_error) || isapprox(sol["solution"]["qloss"]["176"], 0.705918492721665, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["176"], 1.6709897500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["183"], 0.5252850277571399, rtol=high_error) || isapprox(sol["solution"]["qloss"]["183"], 0.5221383878424, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["183"], 1.2029095900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["186"], 1.3909799261651998, rtol=high_error) || isapprox(sol["solution"]["qloss"]["186"], 1.37275042117296, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["186"], 3.2194948200, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["193"], 0.056753975763509996, rtol=high_error) || isapprox(sol["solution"]["qloss"]["193"], 0.055585525076265005, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["193"], 0.1316089200, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["196"], 1.13319111293388, rtol=high_error) || isapprox(sol["solution"]["qloss"]["196"], 1.1285472215186398, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["196"], 2.6177458800, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["211"], 0.998280090117825, rtol=high_error) || isapprox(sol["solution"]["qloss"]["211"], 0.996380278761075, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["211"], 2.2870965000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["219"], 2.3623327708289095, rtol=high_error) || isapprox(sol["solution"]["qloss"]["219"], 2.3544057308871, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["219"], 5.4670534100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["237"], 4.517205801359999, rtol=high_error) || isapprox(sol["solution"]["qloss"]["237"], 4.511456741508001, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["237"], 10.4554147700, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["243"], 0.08397893497365, rtol=high_error) || isapprox(sol["solution"]["qloss"]["243"], 0.08328937531733999, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["243"], 0.1919332300, rtol=high_error)
    end
    @testset "d-y transformers" begin
    end
    @testset "y-d transformers" begin
        @test isapprox(sol["solution"]["qloss"]["78"], 0.021839478815700003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["78"], 0.02208981022794, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["78"], 0.1007913600, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["79"], 0.027299348519625007, rtol=high_error) || isapprox(sol["solution"]["qloss"]["79"], 0.02754018339345, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["79"], 0.1259892100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["80"], 0.021839478815700003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["80"], 0.022122560825820002, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["80"], 0.1007913600, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["81"], 0.027299348519625007, rtol=high_error) || isapprox(sol["solution"]["qloss"]["81"], 0.027527532095699997, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["81"], 0.1259892100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["82"], 0.027954717495075002, rtol=high_error) || isapprox(sol["solution"]["qloss"]["82"], 0.028346758909740002, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["82"], 0.1290129400, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["103"], 3.8186149523215045, rtol=high_error) || isapprox(sol["solution"]["qloss"]["103"], 3.85865829159006, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["103"], 17.5617713900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["105"], 0.0211803677496, rtol=high_error) || isapprox(sol["solution"]["qloss"]["105"], 0.021247199182800002, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["105"], 0.0967051300, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["106"], 0.06440729009436, rtol=high_error) || isapprox(sol["solution"]["qloss"]["106"], 0.06461051751498001, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["106"], 0.2940599300, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["107"], 0.10063629357911999, rtol=high_error) || isapprox(sol["solution"]["qloss"]["107"], 0.10095383611716001, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["107"], 0.4594686300, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["108"], 0.06440729009436, rtol=high_error) || isapprox(sol["solution"]["qloss"]["108"], 0.06461051751498001, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["108"], 0.2940599300, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["109"], 0.12579614452044, rtol=high_error) || isapprox(sol["solution"]["qloss"]["109"], 0.12619307514642, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["109"], 0.5743358100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["110"], 0.12579614452044, rtol=high_error) || isapprox(sol["solution"]["qloss"]["110"], 0.12619307514642, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["110"], 0.5743358100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["111"], 0.06440729009436, rtol=high_error) || isapprox(sol["solution"]["qloss"]["111"], 0.06461051751498001, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["111"], 0.2940599300, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["114"], 0.18902030235046502, rtol=high_error) || isapprox(sol["solution"]["qloss"]["114"], 0.191986957551615, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["114"], 0.8813830000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["115"], 0.09073321677121499, rtol=high_error) || isapprox(sol["solution"]["qloss"]["115"], 0.09295883642465999, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["115"], 0.4230814000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["116"], 0.8509453768032299, rtol=high_error) || isapprox(sol["solution"]["qloss"]["116"], 0.8509453768032299, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["116"], 3.9678738100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["117"], 0.68075569238241, rtol=high_error) || isapprox(sol["solution"]["qloss"]["117"], 0.68075569238241, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["117"], 3.1742992400, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["137"], 0.08927471934045, rtol=high_error) || isapprox(sol["solution"]["qloss"]["137"], 0.09138791648508002, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["137"], 0.4159287500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["138"], 0.21797801424157504, rtol=high_error) || isapprox(sol["solution"]["qloss"]["138"], 0.2207275322463, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["138"], 1.0155566900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["139"], 0.27459649172025, rtol=high_error) || isapprox(sol["solution"]["qloss"]["139"], 0.27459649172025, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["139"], 1.2793419400, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["142"], 2.50765006471392, rtol=high_error) || isapprox(sol["solution"]["qloss"]["142"], 2.5200613430745604, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["142"], 11.4694604900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["158"], 0.688485359719725, rtol=high_error) || isapprox(sol["solution"]["qloss"]["158"], 0.6974213731760999, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["158"], 1.5870753500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["159"], 1.4905373470376102, rtol=high_error) || isapprox(sol["solution"]["qloss"]["159"], 1.50988338192756, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["159"], 3.4359364500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["172"], 0.02875534015788, rtol=high_error) || isapprox(sol["solution"]["qloss"]["172"], 0.02880493365015, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["172"], 0.0666469900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["173"], 1.3478985286981204, rtol=high_error) || isapprox(sol["solution"]["qloss"]["173"], 1.37284050719844, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["173"], 3.1240775600, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["187"], 0.20421041290964995, rtol=high_error) || isapprox(sol["solution"]["qloss"]["187"], 0.20770307201142002, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["187"], 0.4726549100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["188"], 0.61814973437385, rtol=high_error) || isapprox(sol["solution"]["qloss"]["188"], 0.6287220958183801, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["188"], 1.4307392800, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["189"], 0.9148616682239998, rtol=high_error) || isapprox(sol["solution"]["qloss"]["189"], 0.9305087642112001, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["189"], 2.1174941100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["197"], 1.0541316252285, rtol=high_error) || isapprox(sol["solution"]["qloss"]["197"], 1.0700819588430002, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["197"], 2.4351124800, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["198"], 1.0541316252285, rtol=high_error) || isapprox(sol["solution"]["qloss"]["198"], 1.0700819588430002, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["198"], 2.4351124800, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["209"], 0.23943809577613495, rtol=high_error) || isapprox(sol["solution"]["qloss"]["209"], 0.24444872453199, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["209"], 1.1155530200, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["212"], 0.03628644426567001, rtol=high_error) || isapprox(sol["solution"]["qloss"]["212"], 0.036532078594920005, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["212"], 0.0831336900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["213"], 0.33036810184452003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["213"], 0.33260446720752, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["213"], 0.7568838000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["214"], 0.149839245321315, rtol=high_error) || isapprox(sol["solution"]["qloss"]["214"], 0.15085355419794003, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["214"], 0.3432861300, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["215"], 0.33036810184452003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["215"], 0.33260446720752, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["215"], 0.7568838000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["216"], 0.33036810184452003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["216"], 0.33260446720752, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["216"], 0.7568838000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["220"], 0.52572490891908, rtol=high_error) || isapprox(sol["solution"]["qloss"]["220"], 0.5328772595047201, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["220"], 2.4252646000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["222"], 0.41160922632818997, rtol=high_error) || isapprox(sol["solution"]["qloss"]["222"], 0.41160922632818997, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["222"], 1.9191424800, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["223"], 0.4565172284356199, rtol=high_error) || isapprox(sol["solution"]["qloss"]["223"], 0.4565172284356199, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["223"], 2.1285262100, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["224"], 0.20639164089501003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["224"], 0.20639164089501003, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["224"], 0.9623102500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["225"], 0.32248760511020996, rtol=high_error) || isapprox(sol["solution"]["qloss"]["225"], 0.32548741572636003, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["225"], 1.5036097800, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["226"], 0.18767984645276994, rtol=high_error) || isapprox(sol["solution"]["qloss"]["226"], 0.18767984645276994, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["226"], 0.8750608000, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["227"], 0.20639164089501003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["227"], 0.20639164089501003, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["227"], 0.9623102500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["228"], 0.20639164089501003, rtol=high_error) || isapprox(sol["solution"]["qloss"]["228"], 0.21000064287492004, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["228"], 0.9623102500, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["240"], 0.8147962769714251, rtol=high_error) || isapprox(sol["solution"]["qloss"]["240"], 0.82665488820558, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["240"], 3.7623205200, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["241"], 1.9352101585618497, rtol=high_error) || isapprox(sol["solution"]["qloss"]["241"], 1.9624181245655101, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["241"], 8.9358367900, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["244"], 0.780630068555325, rtol=high_error) || isapprox(sol["solution"]["qloss"]["244"], 0.78401384984562, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["244"], 1.7841280700, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["245"], 0.25195826850293995, rtol=high_error) || isapprox(sol["solution"]["qloss"]["245"], 0.25195826850293995, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["245"], 1.1652016600, rtol=high_error)
        @test isapprox(sol["solution"]["qloss"]["246"], 0.25195826850293995, rtol=high_error) || isapprox(sol["solution"]["qloss"]["246"], 0.25601783015316004, rtol=high_error) 
        @test isapprox(sol["solution"]["ieff"]["246"], 1.1652016600, rtol=high_error)
    end
    @testset "y-y transformers" begin
    end
    @testset "d-d transformers" begin
    end
    @testset "lines" begin
        @test isapprox(sol["solution"]["qloss"]["1"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["1"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["2"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["2"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["3"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["3"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["4"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["4"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["5"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["5"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["6"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["6"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["7"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["7"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["8"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["8"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["9"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["9"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["10"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["10"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["11"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["11"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["12"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["12"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["13"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["13"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["14"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["14"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["15"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["15"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["16"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["16"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["17"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["17"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["18"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["18"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["19"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["19"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["20"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["20"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["21"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["21"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["22"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["22"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["23"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["23"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["24"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["24"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["26"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["26"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["27"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["27"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["28"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["28"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["29"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["29"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["30"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["30"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["31"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["31"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["32"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["32"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["33"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["33"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["34"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["34"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["35"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["35"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["36"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["36"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["37"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["37"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["38"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["38"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["39"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["39"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["40"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["40"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["41"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["41"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["42"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["42"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["43"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["43"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["44"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["44"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["45"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["45"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["46"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["46"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["47"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["47"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["48"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["48"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["49"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["49"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["50"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["50"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["51"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["51"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["52"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["52"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["53"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["53"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["54"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["54"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["55"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["55"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["56"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["56"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["57"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["57"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["58"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["58"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["59"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["59"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["60"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["60"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["61"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["61"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["62"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["62"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["63"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["63"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["64"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["64"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["65"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["65"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["66"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["66"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["67"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["67"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["68"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["68"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["69"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["69"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["71"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["71"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["72"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["72"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["73"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["73"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["74"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["74"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["75"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["75"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["76"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["76"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["77"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["77"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["83"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["83"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["84"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["84"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["86"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["86"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["87"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["87"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["88"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["88"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["89"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["89"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["90"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["90"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["91"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["91"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["92"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["92"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["93"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["93"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["94"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["94"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["95"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["95"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["96"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["96"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["97"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["97"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["98"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["98"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["99"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["99"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["100"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["100"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["101"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["101"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["102"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["102"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["104"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["104"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["112"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["112"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["113"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["113"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["118"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["118"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["119"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["119"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["120"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["120"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["121"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["121"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["123"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["123"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["124"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["124"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["125"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["125"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["126"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["126"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["127"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["127"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["128"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["128"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["129"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["129"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["130"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["130"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["131"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["131"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["133"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["133"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["134"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["134"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["135"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["135"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["136"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["136"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["140"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["140"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["141"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["141"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["143"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["143"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["144"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["144"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["145"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["145"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["146"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["146"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["147"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["147"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["149"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["149"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["150"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["150"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["151"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["151"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["152"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["152"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["153"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["153"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["154"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["154"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["155"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["155"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["156"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["156"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["160"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["160"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["161"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["161"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["162"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["162"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["163"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["163"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["164"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["164"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["165"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["165"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["166"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["166"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["167"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["167"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["168"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["168"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["169"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["169"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["170"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["170"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["174"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["174"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["175"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["175"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["177"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["177"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["178"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["178"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["179"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["179"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["180"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["180"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["181"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["181"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["182"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["182"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["184"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["184"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["185"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["185"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["190"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["190"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["191"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["191"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["192"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["192"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["194"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["194"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["195"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["195"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["199"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["199"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["200"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["200"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["201"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["201"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["202"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["202"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["203"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["203"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["204"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["204"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["205"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["205"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["206"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["206"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["207"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["207"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["208"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["208"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["210"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["210"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["217"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["217"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["218"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["218"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["221"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["221"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["229"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["229"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["230"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["230"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["231"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["231"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["232"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["232"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["233"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["233"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["234"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["234"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["235"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["235"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["236"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["236"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["238"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["238"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["239"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["239"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["qloss"]["242"], 0.0, rtol=low_error)
        @test isapprox(sol["solution"]["ieff"]["242"], 0.0, rtol=low_error)
    end
end