In [4]:

    

;cd test


    

    




/Users/prezaei/Dropbox/JuliaPower/test

In [3]:

    

include("data/case6ww_ps_2.jl")


    

In [12]:

    

reload("ACPowerFlow")


    

    




Warning: replacing module ACPowerFlow
Warning: Method definition real(SparseMatrixCSC{Tv,Ti<:Integer},) in module ACPowerFlow at /Users/prezaei/Dropbox/JuliaPower/ACPowerFlow/src/ACPowerFlow.jl:16 overwritten in module ACPowerFlow at /Users/prezaei/Dropbox/JuliaPower/ACPowerFlow/src/ACPowerFlow.jl:16.
Warning: Method definition imag(SparseMatrixCSC{Tv,Ti<:Integer},) in module ACPowerFlow at /Users/prezaei/Dropbox/JuliaPower/ACPowerFlow/src/ACPowerFlow.jl:17 overwritten in module ACPowerFlow at /Users/prezaei/Dropbox/JuliaPower/ACPowerFlow/src/ACPowerFlow.jl:17.

In [9]:

    

alak = sin


    

    
Out[9]:





sin (generic function with 17 methods)

In [10]:

    

methods(alak)


    

    
Out[10]:




17 methods for generic function sin:
 sin(a::Complex{Float16}) at float16.jl:141
 sin(z::Complex{T<:Real}) at complex.jl:518
 sin(x::Float64) at math.jl:122
 sin(x::Float32) at math.jl:123
 sin(a::Float16) at float16.jl:140
 sin(x::BigFloat) at mpfr.jl:488
 sin(x::Real) at math.jl:124
 sin{Tv,Ti}(A::SparseMatrixCSC{Tv,Ti}) at sparse/sparsematrix.jl:447
 sin{T<:FloatingPoint}(d::DataArray{T<:FloatingPoint,N}) at /Users/prezaei/.julia/v0.3/DataArrays/src/operators.jl:256
 sin{T<:Real}(d::DataArray{T<:Real,N}) at /Users/prezaei/.julia/v0.3/DataArrays/src/operators.jl:256
 sin{T<:FloatingPoint}(adv::AbstractDataArray{T<:FloatingPoint,N}) at /Users/prezaei/.julia/v0.3/DataArrays/src/operators.jl:266
 sin{T<:Real}(adv::AbstractDataArray{T<:Real,N}) at /Users/prezaei/.julia/v0.3/DataArrays/src/operators.jl:266
 sin{T<:Number}(::AbstractArray{T<:Number,1}) at operators.jl:359
 sin{T<:Number}(::AbstractArray{T<:Number,2}) at operators.jl:360
 sin{T<:Number}(::AbstractArray{T<:Number,N}) at operators.jl:362
 sin(::NAtype) at /Users/prezaei/.julia/v0.3/DataArrays/src/operators.jl:481
 sin(z::Dual{T<:Real}) at /Users/prezaei/.julia/v0.3/DualNumbers/src/dual.jl:170

In [13]:

    

using DataFrames, HDF5, JLD, Makeps, ACPowerFlow, Optim, ACSIMSEP

# user input to the code
verbose = true
PartFactFlag = false
PolarFlag = true
# ps = load("data/case2383wp_ps.jld")["ps"]
# ps = load("../data/case6ww_ps.jld")["ps"]
# ps = load("../data/case2_ps.jld")["ps"]; 
ps = load("../data/case6ww_ps_2.jld")["ps"]
ps.shunt[[4:6],[:P,:Q]] = 210.
# ps.shunt[[1:3],[:P,:Q]] = 210.
ps.branch[23,:status] = false

LoadFactor = 1.
ps.shunt[:,:P] = ps.shunt[:,:P] * LoadFactor
ps.shunt[:,:Q] = ps.shunt[:,:Q] * LoadFactor

ACPowerFlow.acpf!(ps, PartFactFlag = PartFactFlag)


    

    




Iter   Max(|g|)   |g|_2      max(Jac'*g)      alpha

at /Users/prezaei/Dropbox/JuliaPower/ACPowerFlow/src/nrsolve.jl:35

      34           g2 = sqrt(sum(g.^2))
 -->  35           @bp
      36   

debug:35> STDIN> det(Jac)
1.5752209568095036e9

debug:35> STDIN> q






    




interrupted
while loading In[13], in expression starting on line 19

 in trap at /Users/prezaei/.julia/v0.3/Debug/src/UI.jl:68
 in __nrsolve#270__ at no file:35
 in acpf! at /Users/prezaei/Dropbox/JuliaPower/ACPowerFlow/src/acpf.jl:34

In [39]:

    

ps.bus


    

    
Out[39]:




ID
type
Pd
Qd
Gs
Bs
area
Vmag
Vang
basekV
zone
Vmax
Vmin
lamP
lamQ
muVx
muVn
locX
locY
1
1
3.0
0.0
0.0
0.0
0.0
1.0
1.05
0.0
230.0
1.0
1.05
1.05
0.0
0.0
0.0
0.0
1.5
1.0
2
2
2.0
0.0
0.0
0.0
0.0
1.0
1.05
-0.06407378348652876
230.0
1.0
1.05
1.05
0.0
0.0
0.0
0.0
2.0
1.0
3
3
2.0
0.0
0.0
0.0
0.0
1.0
1.07
-0.0745825871344223
230.0
1.0
1.07
1.07
0.0
0.0
0.0
0.0
2.5
1.0
4
4
1.0
0.0
0.0
0.0
0.0
1.0
0.9893728693486216
-0.07323090241990156
230.0
1.0
1.05
0.95
0.0
0.0
0.0
0.0
1.0
2.0
5
5
1.0
0.0
0.0
0.0
0.0
1.0
0.9854448350522254
-0.09209034811199586
230.0
1.0
1.05
0.95
0.0
0.0
0.0
0.0
2.0
2.0
6
6
1.0
0.0
0.0
0.0
0.0
1.0
1.0044252703572052
-0.10380265472313888
230.0
1.0
1.05
0.95
0.0
0.0
0.0
0.0
3.0
2.0
7
7
2.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
230.0
1.0
1.05
1.05
0.0
0.0
0.0
0.0
1.5
1.0
8
8
2.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
230.0
1.0
1.05
1.05
0.0
0.0
0.0
0.0
2.0
1.0
9
9
2.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
230.0
1.0
1.07
1.07
0.0
0.0
0.0
0.0
2.5
1.0
10
10
1.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
230.0
1.0
1.05
0.95
0.0
0.0
0.0
0.0
1.0
2.0
11
11
1.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
230.0
1.0
1.05
0.95
0.0
0.0
0.0
0.0
2.0
2.0
12
12
1.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
230.0
1.0
1.05
0.95
0.0
0.0
0.0
0.0
3.0
2.0

In [40]:

    

ps.gen


    

    
Out[40]:




bus
P
Q
Qmax
Qmin
Vsp
mBase
status
Pmax
Pmin
muPx
muPn
muQx
muQn
type
cost
1
1
107.87549691563802
15.956206299830427
100.0
-100.0
1.05
100.0
1.0
200.0
50.0
0.0
0.0
0.0
0.0
3.0
11.669
2
2
50.0000000002266
74.35647509772897
100.0
-100.0
1.05
100.0
1.0
150.0
37.5
0.0
0.0
0.0
0.0
2.0
10.333
3
3
60.000000000125596
89.62677449890072
60.0
-100.0
1.07
100.0
1.0
180.0
45.0
0.0
0.0
0.0
0.0
2.0
10.833
4
7
0.0
0.0
100.0
-100.0
1.05
100.0
1.0
200.0
50.0
0.0
0.0
0.0
0.0
3.0
11.669
5
8
0.0
0.0
100.0
-100.0
1.05
100.0
1.0
150.0
37.5
0.0
0.0
0.0
0.0
2.0
10.333
6
9
0.0
0.0
60.0
-100.0
1.07
100.0
1.0
180.0
45.0
0.0
0.0
0.0
0.0
2.0
10.833

In [41]:

    

ps.branch


    

    
Out[41]:




from
to
R
X
B
rateA
rateB
rateC
tap
shift
status
Pf
Qf
Pt
Qt
muF
muT
ImagF
ImagT
switchable
fail_rate
type
contg_st
1
1
2
0.1
0.2
0.04
40.0
40.0
40.0
1.0
0.0
true
28.68967918746788
-15.418662090643664
-27.784737184995734
12.818546095587976
0.0
0.0
0.3101945600481609
0.29142028379812457
0.0
0.0
0.0
0.0
2
1
4
0.05
0.2
0.04
60.0
60.0
60.0
1.0
0.0
true
43.58494942844689
20.12014966750796
-42.497393816493656
-19.93264456890108
0.0
0.0
0.45718930709679295
0.4744392866830986
0.0
0.0
0.0
0.0
3
1
5
0.08
0.3
0.06
40.0
40.0
40.0
1.0
0.0
true
35.60086829972325
11.254718722966182
-34.52732269404179
-13.449727270454023
0.0
0.0
0.35559544355404976
0.3760173476075476
0.0
0.0
0.0
0.0
4
2
3
0.05
0.25
0.06
40.0
40.0
40.0
1.0
0.0
true
2.930319554440286
-12.268749405066632
-2.8900062829704787
5.728115762415653
0.0
0.0
0.12013181900313918
0.05996144458327816
0.0
0.0
0.0
0.0
5
2
4
0.05
0.1
0.02
60.0
60.0
60.0
1.0
0.0
true
33.09094658383731
46.05406848658985
-31.585843319780682
-45.125220633079564
0.0
0.0
0.540092204270718
0.5567295672680539
0.0
0.0
0.0
0.0
6
2
5
0.1
0.3
0.04
30.0
30.0
30.0
1.0
0.0
true
15.51452625587755
15.353153175036585
-15.01657680656928
-18.006507872974048
0.0
0.0
0.2078766068147443
0.23792668507513764
0.0
0.0
0.0
0.0
7
2
6
0.07
0.2
0.05
90.0
90.0
90.0
1.0
0.0
true
26.248944791067107
12.399456745581194
-25.665641584970246
-16.01130146606351
0.0
0.0
0.27647827865786684
0.30117127410844996
0.0
0.0
0.0
0.0
8
3
5
0.12
0.26
0.05
70.0
70.0
70.0
1.0
0.0
true
19.116823909487966
23.174456112664586
-18.023247504648793
-26.09504437617412
0.0
0.0
0.2807642721458668
0.3218259635364828
0.0
0.0
0.0
0.0
9
3
6
0.02
0.1
0.02
80.0
80.0
80.0
1.0
0.0
true
43.773182373608044
60.72420262382019
-42.76979738721798
-57.861047815602106
0.0
0.0
0.6995949990317755
0.7163545321597278
0.0
0.0
0.0
0.0
10
4
5
0.2
0.4
0.08
20.0
20.0
20.0
1.0
0.0
true
4.0832371362880595
-4.942134797899984
-4.047017528564535
-2.785266776789868
0.0
0.0
0.06479594389413164
0.04985410084354745
0.0
0.0
0.0
0.0
11
5
6
0.1
0.3
0.06
40.0
40.0
40.0
1.0
0.0
true
1.6141645340893462
-9.663453703215161
-1.5645610278208564
3.8723492820308523
0.0
0.0
0.09942047816289923
0.04158074008806207
0.0
0.0
0.0
0.0
12
7
8
0.1
0.2
0.04
40.0
40.0
40.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
13
7
10
0.05
0.2
0.04
60.0
60.0
60.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
14
7
11
0.08
0.3
0.06
40.0
40.0
40.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
15
8
9
0.05
0.25
0.06
40.0
40.0
40.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
16
8
10
0.05
0.1
0.02
60.0
60.0
60.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
17
8
11
0.1
0.3
0.04
30.0
30.0
30.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
18
8
12
0.07
0.2
0.05
90.0
90.0
90.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
19
9
11
0.12
0.26
0.05
70.0
70.0
70.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
20
9
12
0.02
0.1
0.02
80.0
80.0
80.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
21
10
11
0.2
0.4
0.08
20.0
20.0
20.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
22
11
12
0.1
0.3
0.06
40.0
40.0
40.0
1.0
0.0
true
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
23
6
11
0.2
0.4
0.08
20.0
20.0
20.0
1.0
0.0
false
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip
&vellip

In [4]:

    

ps2=Caseps[]


    

    
Out[4]:





0-element Array{Caseps,1}

In [6]:

    

typeof(ps1)


    

    
Out[6]:





Caseps (constructor with 1 method)

In [7]:

    

ps2.bus = DataFrame(A = [1,2,3])


    

    




type Array has no field bus
while loading In[7], in expression starting on line 1

In [8]:

    

DataFrame(A = [1,2,3])


    

    
Out[8]:




A
1
1
2
2
3
3

In [ ]:

    

ps2[1]