FIx depreciation warnings + more #4

frederikgeth · 2016-09-21T13:10:08Z

See
Updating to be compatible with latest versions of Julia/JuMP/JuMPChance
#3

Updating for Julia 0.5.0, JuMPChance 0.3.0 and JuMP 0.14.1.

fix al the depreciation warnings
fix with_probability

See Updating to be compatible with latest versions of Julia/JuMP/JuMPChance mlubin#3 Updating for Julia 0.5.0, JuMPChance 0.3.0 and JuMP 0.14.1. fix al the depreciation warnings fix with_probability

mlubin · 2016-09-22T02:46:22Z

codejl/ccopfmodel_simulation.jl

@@ -15,7 +15,8 @@ THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLI
 # Define packages to be used:
 using JuMP # Optimization package
 using JuMPChance # Chance constraints package
-using Gurobi # Solver, could be altered if needed
+# using Gurobi # Solver, could be altered if needed
+using Mosek # Solver, could be altered if needed


This change should stay local, no reason to change the solver from what we used in the paper

mlubin · 2016-09-22T02:47:31Z

Everything here appears correct. The problem is infeasible with this code?

frederikgeth · 2016-09-22T14:30:11Z

Sorry for the CPLEX code in the pull request, I just wanted to share some
code to allow you to verify. I'll clean it up if I get it working :)

The console output for that code + my installed packages below:

theron:run fgeth$ julia ../codejl/ccopfmodel_simulation.jl
bpa-season1-gamma1.0.dat


reading config bpa-season1-gamma1.0.dat

>>>> line_probability_threshold = 0.01

>>>> gen_probability_threshold = 0.1666

>>>> logfilename = casebpa.log

>>>> loadscale = 1.0

>>>> mvaBase = 100.0

>>>> thermalLimitscale = 0.9

>>>> hydro_indices = ../data/bpa/hydro_power_plants_indices.csv

>>>> nuclear_indices = ../data/bpa/nuclear_power_plants_indices.csv

>>>> coal_indices = ../data/bpa/coal_power_plants_indices.csv

>>>> gas_indices = ../data/bpa/gas_power_plants_indices.csv

>>>> ramp_file = ../data/bpa/rampbpa.csv

>>>> loaded_lines_file = ../data/bpa/overload_line_ind.csv

>>>> wind_forecast_file = ../data/bpa/season1/bpawind_det_forecast2160.csv

>>>> wind_sigma_file = ../data/bpa/season1/bpawind_det_sigma2160.csv

>>>> load_file = ../data/bpa/season1/load2160.csv

>>>> hydro_limit_file = ../data/bpa/season1/hydro_power_plants2160.csv

>>>> loaded_lines_probability_threshold = 0.01

>>>> robust_mean_lower_file =
../data/bpa/season1/bpawind_rob_mean_dn2160.csv

>>>> robust_mean_upper_file =
../data/bpa/season1/bpawind_rob_mean_up2160.csv

>>>> robust_sigma_lower_file =
../data/bpa/season1/bpawind_rob_sigma_dn2160.csv

>>>> robust_sigma_upper_file =
../data/bpa/season1/bpawind_rob_sigma_up2160.csv

>>>> robust_budget = 1.0

>>>> output_matfile = ccopf_bpa_simulation_results_season1_gamma1.0.mat

ROBUST CCOPF with budget 1.0

10 time steps

node 27 is farm 1

node 18 is farm 2

node 23 is farm 3

node 24 is farm 4

node 18 is farm 5

node 907 is farm 6

node 970 is farm 7

node 1040 is farm 8

node 1384 is farm 9

node 22 is farm 10

node 25 is farm 11

node 2186 is farm 12

node 1121 is farm 13

node 1394 is farm 14

node 1390 is farm 15

node 1400 is farm 16

node 21 is farm 17

node 18 is farm 18

node 31 is farm 19

node 31 is farm 20

node 25 is farm 21

node 24 is farm 22

node 970 is farm 23

node 1413 is farm 24



######## T = 1


-> sumvar 6615.7132128661 stddev 81.33703469432666

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 2


-> sumvar 6472.039335426101 stddev 80.44898591919043

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 3


-> sumvar 6509.1704671961 stddev 80.67943026073064

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 4


-> sumvar 6075.5444353361 stddev 77.94577881666268

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 5


-> sumvar 3093.2861953361 stddev 55.617319203069286

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 6


-> sumvar 2815.3004762561004 stddev 53.05940516304438

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 7


-> sumvar 3223.7749969661 stddev 56.778296883281904

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 8


-> sumvar 4248.993750646101 stddev 65.18430601491514

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 9


-> sumvar 5519.839991086101 stddev 74.29562565243057

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible



######## T = 10


-> sumvar 2843.1364111261 stddev 53.321069110869296

Iteration 0: 6259 constraint violations, 0 objective linearization
violations

WARNING: Not solved to optimality, status: Infeasible


----

Julia packages:

   *_**       _ **_**(_)**_**     |  A fresh approach to technical
computing*

  *(_)**     | **(_)* *(_)**    |  Documentation: http://docs.julialang.org
<http://docs.julialang.org>*

*   _ _   _| |_  __ _   |  Type "?help" for help.*

*  | | | | | | |/ _` |  |*

*  | | |_| | | | (_| |  |  Version 0.5.0 (2016-09-19 18:14 UTC)*

* _/ |\__'_|_|_|\__'_|  |  Official http://julialang.org/
<http://julialang.org/> release*

*|__/                   |  x86_64-apple-darwin13.4.0*


*julia> **Pkg.status()*

6 required packages:

 - CPLEX                         0.1.4

 - Debug                         0.1.6

 - Gurobi                        0.2.4

 - JuMPChance                    0.3.0

 - MAT                           0.3.1

 - Mosek                         0.5.0

25 additional packages:

 - BinDeps                       0.4.5

 - Blosc                         0.1.7

 - BufferedStreams               0.2.0

 - Calculus                      0.1.15

 - Compat                        0.9.2

 - DataStructures                0.4.6

 - Distributions                 0.10.2

 - ECOS                          0.6.7

 - ForwardDiff                   0.2.5

 - HDF5                          0.6.6

 - Homebrew                      0.4.0

 - JSON                          0.7.0

 - JuMP                          0.14.1

 - Lazy                          0.11.1

 - Libz                          0.2.0

 - MacroTools                    0.3.2

 - MathProgBase                  0.5.6

 - NaNMath                       0.2.1

 - PDMats                        0.4.2

 - ReverseDiffSparse             0.5.8

 - Rmath                         0.1.3

 - SHA                           0.2.1

 - StatsBase                     0.9.0

 - StatsFuns                     0.3.1

 - URIParser                     0.1.6



On Thu, Sep 22, 2016 at 4:47 AM, Miles Lubin <[email protected]>
wrote:

> Everything here appears correct. The problem is infeasible with this code?
>
> —
> You are receiving this because you authored the thread.
> Reply to this email directly, view it on GitHub
> <https://github.com/mlubin/RobustCCOPFSupplement/pull/4#issuecomment-248796953>,
> or mute the thread
> <https://github.com/notifications/unsubscribe-auth/AL9YFlCROf_f_X6t40ISPloxAoQAwGwQks5qsexDgaJpZM4KCx77>
> .
>

mlubin · 2016-09-22T14:34:30Z

Were you able to successfully run the original code with the listed versions of the packages?

frederikgeth · 2016-09-22T14:51:44Z

no, I got stuck at the following error in Julia 0.4.6 - JuMP 0.9.0 -
JuMPChance 0.1.1 -
Gurobi.jl 0.1.26:

ERROR: LoadError: Cannot construct an affine expression with a term of type
(JuMP.GenericQuadExpr{Float64,JuMP.Variable})*(Float64)

console dump + pkg list below:

electa-05:run fgeth$
/Applications/Julia-0.4.6.app/Contents/Resources/julia/bin/julia
../codejl/ccopfmodel_simulation.jl bpa-season1-gamma1.0.dat

WARNING: Base.String is deprecated, use AbstractString instead.

likely near
/Users/fgeth/JuliaProjects/RobustCCOPFSupplement-master/codejl/ccopfmodel_simulation.jl:24

WARNING: Base.String is deprecated, use AbstractString instead.

likely near
/Users/fgeth/JuliaProjects/RobustCCOPFSupplement-master/codejl/ccopfmodel_simulation.jl:24

WARNING: module JuMP should explicitly import setindex! from Base