case study 2 error

[chenhu822]

Hi, i try to. run the dcpowergrid case, and i modified partition_and_solve.jl to the following due to some changes in kahypar.jl and schawarzopt.jl. However, when i included the dual and primal links in optimizer, #it shows error:

Optimizing with user provided overlap
Initializing SchwarzOpt...
ERROR: AssertionError: !(incident_node in all_nodes(subproblem_graph))
Stacktrace:
[1] macro expansion at /Users/822qiuqiu/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:216 [inlined]
[2] macro expansion at ./timing.jl:233 [inlined]
[3] _initialize_optimizer!(::SchwarzOpt.Optimizer) at /Users/822qiuqiu/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:136
[4] optimize!(::SchwarzOpt.Optimizer) at /Users/822qiuqiu/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:390
[5] optimize!(::OptiGraph; subgraphs::Array{OptiGraph,1}, sub_optimizer::MathOptInterface.OptimizerWithAttributes, overlap::Int64, max_iterations::Int64, primal_links::Array{LinkConstraintRef,1}, dual_links::Array{LinkConstraintRef,1}) at /Users/822qiuqiu/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:526
[6] top-level scope at /Users/822qiuqiu/Desktop/PlasmoExamples/case_studies/dcpowergrid/partition_and_solve.jl:63

Then i drop the primal and dual links in optimizer, it shows below, and it seemed that the objective function has some problem.

Number of variables: 100203
Number of constraints: 94483
Number of subproblems: 4
Overlap:
Subproblem variables: [27865, 29044, 32810, 34732]
Subproblem constraints: [10493, 10653, 12407, 13112]
┌ Warning: Suboptimal solution detected for subproblem with status INVALID_OPTION
└ @ SchwarzOpt ~/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:309
┌ Warning: Suboptimal solution detected for subproblem with status INVALID_OPTION
└ @ SchwarzOpt ~/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:309
┌ Warning: Suboptimal solution detected for subproblem with status INVALID_OPTION
└ @ SchwarzOpt ~/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:309
┌ Warning: Suboptimal solution detected for subproblem with status INVALID_OPTION
└ @ SchwarzOpt ~/.julia/packages/SchwarzOpt/rfQsH/src/optimizer.jl:309
Iter | Obj | Prf | Duf
1 | 0.00e+00 | 0.00e+00 | 0.00e+00
Number of Iterations: 1
Time spent in subproblems: 239.197291583
Solution Time: 298.77887892723083
EXIT: SchwarzOpt Finished with status: INVALID_OPTION

[chenhu822]

Here is my code:
#NOTE: This file uses n_parts, imbalance, and overlap defined in the calling file

using KaHyPar
using SchwarzOpt
using LightGraphs

include((@dir)*"/powergrid.jl")

max_imbalance = 0.1
n_parts = 4
overlap = 10

hypergraph,hyper_map = gethypergraph(powergrid) #create hypergraph object based on graph
n_vertices = length(vertices(hypergraph))
node_sizes = [num_variables(node) for node in all_nodes(powergrid)]
edge_weights = [num_linkconstraints(edge) for edge in all_edges(powergrid)]

A = LightGraphs.incidence_matrix(hypergraph)
n_vertices = length(vertices(hypergraph))
kgraph = KaHyPar.HyperGraph(A, node_sizes, edge_weights)
node_vector = KaHyPar.partition(kgraph, n_parts; imbalance=max_imbalance, configuration=:edge_cut)

#Create a model partition
#partition = Partition(powergrid,node_vector,hyper_map)
partition = Partition(hypergraph,node_vector,hyper_map)

#Setup subgraphs based on partition
make_subgraphs!(powergrid,partition)

subgraphs = getsubgraphs(powergrid)
expanded_subs = expand.(Ref(powergrid),subgraphs,Ref(overlap))

#Solve with overlapping schwarz algorithm
sub_optimizer = optimizer_with_attributes(Ipopt.Optimizer,"tol" => 1e-12,"print_level" => 0,"linear_solver" => "ma27")

#optimize using schwarz overlapping decomposition
SchwarzOpt.optimize!(powergrid;
subgraphs = expanded_subs,
sub_optimizer = sub_optimizer,
max_iterations = 100,
primal_links = primal_links,dual_links = dual_links)