-
Notifications
You must be signed in to change notification settings - Fork 0
/
3PK.SIF
671 lines (515 loc) · 17 KB
/
3PK.SIF
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
***************************
* SET UP THE INITIAL DATA *
***************************
NAME 3PK
* Problem :
* *********
* A problem arising in the estimation of structured O/D matrix
* Source:
* M. Bierlaire, private communication
* see also
* M. Bierlaire and Ph. L. Toint,
* MEUSE: an origin-destination estimator that exploits structure,
* Transportation Research B, 29, 1, 47--60, 1995.
* SIF input: Ph. Toint, Dec 1989, Corrected July 1993.
* classification SBR2-MN-30-0
* Parameters
* Number of parking columns
IE NPKC 3
IA NPKC-1 NPKC -1
IA NPKC+1 NPKC 1
* Number of centroids
IE NCENT 6
IA NCENT-1 NCENT -1
RI RNCENT-1 NCENT-1
* Gamma : general trust in the countings
RE GAMMA 1.0000e+04
* Objective for fill-up coefficients
RE FT0 0.500000
RE FT1 0.500000
RE FT2 0.500000
* Inverse of the weights for fill-up coefficients
RE WFT0 1.000000
RE WFT1 1.000000
RE WFT2 1.000000
* Number of available countings
IE COUNT 9
IA COUNT-1 COUNT -1
* Inverse of the weight on unknow cells
RE DEFW 999.999953
* Trip ends
* Some interesting values...
IE 0 0
IE 1 1
* Countings
RE COU0 910.000000
RE COU1 175.000000
RE COU2 1915.000000
RE COU3 450.000000
RE COU4 260.000000
RE COU5 80.000000
RE COU6 670.000000
RE COU7 1450.000000
RE COU8 990.000000
* Inverse of the trusts in countings
RE PHI0 1.000000
RE PHI1 1.000000
RE PHI2 1.000000
RE PHI3 1.000000
RE PHI4 1.000000
RE PHI5 1.000000
RE PHI6 1.000000
RE PHI7 1.000000
RE PHI8 1.000000
DO I 0 COUNT-1
A/ PHI(I) PHI(I) GAMMA
OD I
VARIABLES
* Coefficients pour les parkings
X A1,0
X A2,0
X A3,0
X A4,0
X A5,0
X A0,1
X A2,1
X A3,1
X A4,1
X A5,1
X A0,2
X A1,2
X A3,2
X A4,2
X A5,2
* O/D matrix
DO J NPKC NCENT-1
IA J+1 J 1
IA J-1 J -1
DO I 0 J-1
X T(I,J)
OD I
DO I J+1 NCENT-1
X T(I,J)
OD I
OD J
GROUPS
* A priori information, but not parkings
XN G0,3 T0,3 0.010000
XN G0,3 'SCALE' 1.0
XN G1,3 T1,3 0.007143
XN G1,3 'SCALE' 1.0
XN G2,3 T2,3 0.008333
XN G2,3 'SCALE' 1.0
XN G4,3 T4,3 0.050000
XN G4,3 'SCALE' 1.0
XN G5,3 T5,3 0.050000
XN G5,3 'SCALE' 1.0
XN G0,4 T0,4 0.005000
XN G0,4 'SCALE' 1.0
XN G1,4 T1,4 0.005556
XN G1,4 'SCALE' 1.0
XN G2,4 T2,4 0.050000
XN G2,4 'SCALE' 1.0
XN G3,4 T3,4 0.001667
XN G3,4 'SCALE' 1.0
XN G5,4 T5,4 0.025000
XN G5,4 'SCALE' 1.0
XN G0,5 T0,5 0.020000
XN G0,5 'SCALE' 1.0
XN G1,5 T1,5 0.033333
XN G1,5 'SCALE' 1.0
XN G2,5 T2,5 0.014286
XN G2,5 'SCALE' 1.0
XN G3,5 T3,5 0.006667
XN G3,5 'SCALE' 1.0
XN G4,5 T4,5 0.050000
XN G4,5 'SCALE' 1.0
* Parking coefficients
AM TMP FT0 5.000000
RD TMP1 TMP 1.0
ZN H0 'SCALE' WFT0
ZN H0 A1,0 TMP1
ZN H0 A2,0 TMP1
ZN H0 A3,0 TMP1
ZN H0 A4,0 TMP1
ZN H0 A5,0 TMP1
AM TMP FT1 5.000000
RD TMP1 TMP 1.0
ZN H1 'SCALE' WFT1
ZN H1 A0,1 TMP1
ZN H1 A2,1 TMP1
ZN H1 A3,1 TMP1
ZN H1 A4,1 TMP1
ZN H1 A5,1 TMP1
AM TMP FT2 5.000000
RD TMP1 TMP 1.0
ZN H2 'SCALE' WFT2
ZN H2 A0,2 TMP1
ZN H2 A1,2 TMP1
ZN H2 A3,2 TMP1
ZN H2 A4,2 TMP1
ZN H2 A5,2 TMP1
* Unknown parking entries
* Countings
DO I 0 COUNT-1
ZN K(I) 'SCALE' PHI(I)
OD I
* Parking columns
* 1,0
RD TMP1 COU7 200.000000
RM TMP TMP1 1.000000
ZN K7 A1,0 TMP
RD TMP1 COU4 200.000000
RM TMP TMP1 1.000000
ZN K4 A1,0 TMP
RD TMP1 COU2 200.000000
RM TMP TMP1 1.000000
ZN K2 A1,0 TMP
* 2,0
RD TMP1 COU8 480.000000
RM TMP TMP1 1.000000
ZN K8 A2,0 TMP
RD TMP1 COU7 480.000000
RM TMP TMP1 1.000000
ZN K7 A2,0 TMP
RD TMP1 COU6 480.000000
RM TMP TMP1 1.000000
ZN K6 A2,0 TMP
RD TMP1 COU2 480.000000
RM TMP TMP1 1.000000
ZN K2 A2,0 TMP
* 3,0
RD TMP1 COU2 120.000000
RM TMP TMP1 1.000000
ZN K2 A3,0 TMP
* 4,0
RD TMP1 COU7 360.000000
RM TMP TMP1 1.000000
ZN K7 A4,0 TMP
RD TMP1 COU2 360.000000
RM TMP TMP1 1.000000
ZN K2 A4,0 TMP
* 5,0
RD TMP1 COU8 560.000000
RM TMP TMP1 1.000000
ZN K8 A5,0 TMP
RD TMP1 COU7 560.000000
RM TMP TMP1 1.000000
ZN K7 A5,0 TMP
RD TMP1 COU2 560.000000
RM TMP TMP1 1.000000
ZN K2 A5,0 TMP
* 0,1
RD TMP1 COU0 240.000000
RM TMP TMP1 1.000000
ZN K0 A0,1 TMP
* 2,1
RD TMP1 COU8 400.000000
RM TMP TMP1 1.000000
ZN K8 A2,1 TMP
RD TMP1 COU7 400.000000
RM TMP TMP1 1.000000
ZN K7 A2,1 TMP
RD TMP1 COU6 400.000000
RM TMP TMP1 1.000000
ZN K6 A2,1 TMP
RD TMP1 COU2 400.000000
RM TMP TMP1 1.000000
ZN K2 A2,1 TMP
RD TMP1 COU0 400.000000
RM TMP TMP1 1.000000
ZN K0 A2,1 TMP
* 3,1
RD TMP1 COU2 420.000000
RM TMP TMP1 1.000000
ZN K2 A3,1 TMP
RD TMP1 COU0 420.000000
RM TMP TMP1 1.000000
ZN K0 A3,1 TMP
* 4,1
RD TMP1 COU7 180.000000
RM TMP TMP1 1.000000
ZN K7 A4,1 TMP
RD TMP1 COU2 180.000000
RM TMP TMP1 1.000000
ZN K2 A4,1 TMP
RD TMP1 COU0 180.000000
RM TMP TMP1 1.000000
ZN K0 A4,1 TMP
* 5,1
RD TMP1 COU8 320.000000
RM TMP TMP1 1.000000
ZN K8 A5,1 TMP
RD TMP1 COU7 320.000000
RM TMP TMP1 1.000000
ZN K7 A5,1 TMP
RD TMP1 COU2 320.000000
RM TMP TMP1 1.000000
ZN K2 A5,1 TMP
RD TMP1 COU0 320.000000
RM TMP TMP1 1.000000
ZN K0 A5,1 TMP
* 0,2
RD TMP1 COU1 20.000000
RM TMP TMP1 1.000000
ZN K1 A0,2 TMP
RD TMP1 COU0 20.000000
RM TMP TMP1 1.000000
ZN K0 A0,2 TMP
* 1,2
RD TMP1 COU1 60.000000
RM TMP TMP1 1.000000
ZN K1 A1,2 TMP
* 3,2
RD TMP1 COU2 40.000000
RM TMP TMP1 1.000000
ZN K2 A3,2 TMP
RD TMP1 COU1 40.000000
RM TMP TMP1 1.000000
ZN K1 A3,2 TMP
RD TMP1 COU0 40.000000
RM TMP TMP1 1.000000
ZN K0 A3,2 TMP
* 4,2
RD TMP1 COU5 120.000000
RM TMP TMP1 1.000000
ZN K5 A4,2 TMP
* 5,2
RD TMP1 COU8 20.000000
RM TMP TMP1 1.000000
ZN K8 A5,2 TMP
RD TMP1 COU5 20.000000
RM TMP TMP1 1.000000
ZN K5 A5,2 TMP
* Other columns
* 0,3
AD TMP COU7 1.000000
ZN K7 T0,3 TMP
AD TMP COU3 1.000000
ZN K3 T0,3 TMP
* 1,3
AD TMP COU7 1.000000
ZN K7 T1,3 TMP
AD TMP COU4 1.000000
ZN K4 T1,3 TMP
* 2,3
AD TMP COU8 1.000000
ZN K8 T2,3 TMP
AD TMP COU7 1.000000
ZN K7 T2,3 TMP
AD TMP COU6 1.000000
ZN K6 T2,3 TMP
* 4,3
AD TMP COU7 1.000000
ZN K7 T4,3 TMP
* 5,3
AD TMP COU8 1.000000
ZN K8 T5,3 TMP
AD TMP COU7 1.000000
ZN K7 T5,3 TMP
* 0,4
AD TMP COU3 1.000000
ZN K3 T0,4 TMP
* 1,4
AD TMP COU4 1.000000
ZN K4 T1,4 TMP
* 2,4
AD TMP COU8 1.000000
ZN K8 T2,4 TMP
AD TMP COU6 1.000000
ZN K6 T2,4 TMP
* 3,4
AD TMP COU3 1.000000
ZN K3 T3,4 TMP
AD TMP COU2 1.000000
ZN K2 T3,4 TMP
* 5,4
AD TMP COU8 1.000000
ZN K8 T5,4 TMP
* 0,5
AD TMP COU6 1.000000
ZN K6 T0,5 TMP
AD TMP COU1 1.000000
ZN K1 T0,5 TMP
AD TMP COU0 1.000000
ZN K0 T0,5 TMP
* 1,5
AD TMP COU6 1.000000
ZN K6 T1,5 TMP
AD TMP COU1 1.000000
ZN K1 T1,5 TMP
* 2,5
AD TMP COU6 1.000000
ZN K6 T2,5 TMP
* 3,5
AD TMP COU6 1.000000
ZN K6 T3,5 TMP
AD TMP COU2 1.000000
ZN K2 T3,5 TMP
AD TMP COU1 1.000000
ZN K1 T3,5 TMP
AD TMP COU0 1.000000
ZN K0 T3,5 TMP
* 4,5
AD TMP COU6 1.000000
ZN K6 T4,5 TMP
AD TMP COU5 1.000000
ZN K5 T4,5 TMP
* Variance minimization
XN L1,0 A1,0 -0.800000
XN L1,0 'SCALE' 0.500000
XN L1,0 A2,0 0.200000
XN L1,0 A3,0 0.200000
XN L1,0 A4,0 0.200000
XN L1,0 A5,0 0.200000
XN L2,0 A1,0 0.200000
XN L2,0 A2,0 -0.800000
XN L2,0 'SCALE' 0.500000
XN L2,0 A3,0 0.200000
XN L2,0 A4,0 0.200000
XN L2,0 A5,0 0.200000
XN L3,0 A1,0 0.200000
XN L3,0 A2,0 0.200000
XN L3,0 A3,0 -0.800000
XN L3,0 'SCALE' 0.500000
XN L3,0 A4,0 0.200000
XN L3,0 A5,0 0.200000
XN L4,0 A1,0 0.200000
XN L4,0 A2,0 0.200000
XN L4,0 A3,0 0.200000
XN L4,0 A4,0 -0.800000
XN L4,0 'SCALE' 0.500000
XN L4,0 A5,0 0.200000
XN L5,0 A1,0 0.200000
XN L5,0 A2,0 0.200000
XN L5,0 A3,0 0.200000
XN L5,0 A4,0 0.200000
XN L5,0 A5,0 -0.800000
XN L5,0 'SCALE' 0.500000
XN L0,1 A0,1 -0.800000
XN L0,1 'SCALE' 0.500000
XN L0,1 A2,1 0.200000
XN L0,1 A3,1 0.200000
XN L0,1 A4,1 0.200000
XN L0,1 A5,1 0.200000
XN L2,1 A0,1 0.200000
XN L2,1 A2,1 -0.800000
XN L2,1 'SCALE' 0.500000
XN L2,1 A3,1 0.200000
XN L2,1 A4,1 0.200000
XN L2,1 A5,1 0.200000
XN L3,1 A0,1 0.200000
XN L3,1 A2,1 0.200000
XN L3,1 A3,1 -0.800000
XN L3,1 'SCALE' 0.500000
XN L3,1 A4,1 0.200000
XN L3,1 A5,1 0.200000
XN L4,1 A0,1 0.200000
XN L4,1 A2,1 0.200000
XN L4,1 A3,1 0.200000
XN L4,1 A4,1 -0.800000
XN L4,1 'SCALE' 0.500000
XN L4,1 A5,1 0.200000
XN L5,1 A0,1 0.200000
XN L5,1 A2,1 0.200000
XN L5,1 A3,1 0.200000
XN L5,1 A4,1 0.200000
XN L5,1 A5,1 -0.800000
XN L5,1 'SCALE' 0.500000
XN L0,2 A0,2 -0.800000
XN L0,2 'SCALE' 0.500000
XN L0,2 A1,2 0.200000
XN L0,2 A3,2 0.200000
XN L0,2 A4,2 0.200000
XN L0,2 A5,2 0.200000
XN L1,2 A0,2 0.200000
XN L1,2 A1,2 -0.800000
XN L1,2 'SCALE' 0.500000
XN L1,2 A3,2 0.200000
XN L1,2 A4,2 0.200000
XN L1,2 A5,2 0.200000
XN L3,2 A0,2 0.200000
XN L3,2 A1,2 0.200000
XN L3,2 A3,2 -0.800000
XN L3,2 'SCALE' 0.500000
XN L3,2 A4,2 0.200000
XN L3,2 A5,2 0.200000
XN L4,2 A0,2 0.200000
XN L4,2 A1,2 0.200000
XN L4,2 A3,2 0.200000
XN L4,2 A4,2 -0.800000
XN L4,2 'SCALE' 0.500000
XN L4,2 A5,2 0.200000
XN L5,2 A0,2 0.200000
XN L5,2 A1,2 0.200000
XN L5,2 A3,2 0.200000
XN L5,2 A4,2 0.200000
XN L5,2 A5,2 -0.800000
XN L5,2 'SCALE' 0.500000
CONSTANTS
* A priori information, but not parkings
DO J NPKC NCENT-1
IA J+1 J 1
IA J-1 J -1
DO I 0 J-1
XN 3PK G(I,J) 1.0
OD I
DO I J+1 NCENT-1
XN 3PK G(I,J) 1.0
OD I
OD J
* Parking coefficients
DO J 0 NPKC-1
XN 3PK H(J) 1.0
OD J
* Countings
DO J 0 COUNT-1
XN 3PK K(J) 1.0
OD J
START POINT
ZV 3PK A1,0 FT0
ZV 3PK A2,0 FT0
ZV 3PK A3,0 FT0
ZV 3PK A4,0 FT0
ZV 3PK A5,0 FT0
ZV 3PK A0,1 FT1
ZV 3PK A2,1 FT1
ZV 3PK A3,1 FT1
ZV 3PK A4,1 FT1
ZV 3PK A5,1 FT1
ZV 3PK A0,2 FT2
ZV 3PK A1,2 FT2
ZV 3PK A3,2 FT2
ZV 3PK A4,2 FT2
ZV 3PK A5,2 FT2
V 3PK T0,3 100.000000
V 3PK T1,3 140.000000
V 3PK T2,3 120.000000
V 3PK T4,3 20.000000
V 3PK T5,3 20.000000
V 3PK T0,4 200.000000
V 3PK T1,4 180.000000
V 3PK T2,4 20.000000
V 3PK T3,4 600.000000
V 3PK T5,4 40.000000
V 3PK T0,5 50.000000
V 3PK T1,5 30.000000
V 3PK T2,5 70.000000
V 3PK T3,5 150.000000
V 3PK T4,5 20.000000
GROUP TYPE
GV SQUARE X
GROUP USES
'DEFAULT' SQUARE
ENDATA
*********************
* SET UP THE GROUPS *
* ROUTINE *
*********************
GROUPS 3PK
INDIVIDUALS
T SQUARE
F X * X
G X + X
H 2.0
ENDATA