# J05    
# The set of nodes in the network
set N := 1 2 3 4 5 6 7 8 9;

# The set of modular sizes for structures and couplers
set W := 4 8 16 20 40 80;

# The set of OD pairs and 
# the number wavelengths required for each o-d pair
param: D: r :=
1	3	8
1	6	10
1	7	11
2	3	9
2	5	13
2	6	10
3	5	11
3	6	7
4	8	4
5	7	12
6	7	3
6	8	6
7	8	2
7	9	9
8	9	13
;

# The set of available structures
set S := 1 2 3 4 5;

# The set of available couplers
set C := 1 2 3 4 5 6;

# The set of routing paths
set P := 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17;

# The set of pairs of routing paths that intersect
set H :=
(1,2)     (4,17)    (10,12)   (2,14)    (7,13)    (2,17)    (9,14)    (4,15)
(2,9)     (7,9)     (13,17)   (4,8)     (9,11)    (4,11)    (12,16)   (6,15)
(3,17)    (8,13)    (1,13)    (5,10)    (11,13)   (6,10)    (2,7)     (8,10)
(4,16)    (10,11)   (2,13)    (7,12)    (2,3)     (7,16)    (3,13)    (9,16)
(7,8)     (13,16)   (4,7)     (9,10)    (2,16)    (9,13)    (4,13)    (13,14)
(8,11)    (1,9)     (5,7)     (10,16)   (4,10)    (12,13)   (6,13)
(9,17)    (2,12)    (7,11)    (16,17)   (6,8)     (2,6)     (8,9)
(13,15)   (4,6)     (8,17)    (1,17)    (7,15)    (3,9)     (9,15)
(1,3)     (5,6)     (10,15)   (2,15)    (9,12)    (4,12)    (12,17)
(2,11)    (7,10)    (14,17)   (4,9)     (11,15)   (6,11)    (2,8)
(4,5)     (8,15)    (1,14)    (6,7)     (2,4)     (7,17)    (3,14);

# Sets of routing paths serving each OD pair
set J[1,3] := 1;
set J[1,6] := 2;
set J[1,7] := 3;
set J[2,3] := 4;
set J[2,5] := 5;
set J[2,6] := 6;
set J[3,5] := 7;
set J[3,6] := 8;
set J[4,8] := 9;
set J[5,7] := 10;
set J[6,7] := 11 12;
set J[6,8] := 13;
set J[7,8] := 14;
set J[7,9] := 15 16;
set J[8,9] := 17;

# Sets of routing paths using each structure
set K[1] := 4 5 6 7 10;
set K[2] := 2 4 6 7 8 9 10 11 13 15 17;
set K[3] := 1 2 3 9 13 14 17;
set K[4] := 4 6 7 10;
set K[5] := 2 4 7 9 10 12 13 16 17;

# Sets of routing paths using each coupler
set L[1] := 4 6;
set L[2] := 4 6 7 10;
set L[3] := 2 9 13 17;
set L[4] := 4 7 10;
set L[5] := 4 7 10;
set L[6] := 2 9 13 17;

# The structure costs
param a:
	4	8	16	20	40	80	 :=
1	28	50	84	112	224	420
2	25	45	75	100	200	375
3	23	41	69	92	184	345
4	27	49	81	108	216	405
5	26	47	78	104	208	390
;

# The coupler costs
param f:
	4	8	16	20	40	80	 :=
1	2	3	6	8	16	30
2	9	15	26	35	71	134
3	8	14	23	31	63	119
4	14	24	41	55	111	209
5	13	22	38	51	103	194
6	14	24	41	55	111	209
;