test/vivien.c

   1 /* Generated from vivien.cloog by CLooG v0.14.0 gmp bits in 0.62s. */
   2 /* CLooG asked for 1368 KBytes. */
   3 for (p1=-54*n+4;p1<=min(4,4*n+1);p1++) {
   4   if (p1%2 == 0) {
   5     i = (p1-2)/2 ;
   6     S1 ;
   7   }
   8 }
   9 if (n >= 1) {
  10   S3(i = 1) ;
  11 }
  12 if (n >= 2) {
  13   S4(i = 1,j = 2) ;
  14   S1(i = 2) ;
  15   S6(i = 1,j = 2) ;
  16 }
  17 for (p1=max(-54*n+4,4*n+2);p1<=6;p1++) {
  18   if (p1%2 == 0) {
  19     i = (p1-2)/2 ;
  20     S1 ;
  21   }
  22 }
  23 for (p1=7;p1<=min(2*n+2,9);p1++) {
  24   for (p2=ceild(-p1+2,4);p2<=-1;p2++) {
  25     if (p1%2 == 0) {
  26       j = (p1+2*p2)/2 ;
  27       S4(i = -p2) ;
  28     }
  29   }
  30   if ((p1+3)%4 == 0) {
  31     i = (p1-1)/4 ;
  32     S3 ;
  33   }
  34   for (i=ceild(p1+1,4);i<=floord(p1-3,2);i++) {
  35     if ((p1+1)%2 == 0) {
  36       j = (p1-2*i-1)/2 ;
  37       S2 ;
  38     }
  39   }
  40   if (p1%2 == 0) {
  41     i = (p1-2)/2 ;
  42     S1 ;
  43   }
  44   if (p1%2 == 0) {
  45     j = (p1-2)/2 ;
  46     S6(i = 1) ;
  47   }
  48 }
  49 for (p1=2*n+3;p1<=min(4*n-2,9);p1++) {
  50   for (p2=ceild(-p1+2,4);p2<=floord(-p1+2*n,2);p2++) {
  51     if (p1%2 == 0) {
  52       j = (p1+2*p2)/2 ;
  53       S4(i = -p2) ;
  54     }
  55   }
  56   if ((p1+3)%4 == 0) {
  57     i = (p1-1)/4 ;
  58     S3 ;
  59   }
  60   for (i=ceild(p1+1,4);i<=n;i++) {
  61     if ((p1+1)%2 == 0) {
  62       j = (p1-2*i-1)/2 ;
  63       S2 ;
  64     }
  65   }
  66   if (p1%2 == 0) {
  67     i = (p1-2)/2 ;
  68     S1 ;
  69   }
  70 }
  71 for (p1=10;p1<=min(12,2*n+2);p1++) {
  72   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
  73     if (p1%2 == 0) {
  74       j = (p1+2*p2)/2 ;
  75       S4(i = -p2) ;
  76     }
  77   }
  78   for (p2=ceild(-p1+6,4);p2<=-1;p2++) {
  79     if (p1%2 == 0) {
  80       j = (p1+2*p2)/2 ;
  81       S4(i = -p2) ;
  82     }
  83     for (p3=1;p3<=-p2;p3++) {
  84       i = -p2+1 ;
  85       if (p1%2 == 0) {
  86         j = (p1+2*p2-2)/2 ;
  87         S5(k = p3) ;
  88       }
  89     }
  90   }
  91   if (p1%2 == 0) {
  92     j = (p1-4)/2 ;
  93     S6(i = 2) ;
  94   }
  95   for (i=max(3,ceild(p1+1,4));i<=floord(p1-3,2);i++) {
  96     if ((p1+1)%2 == 0) {
  97       j = (p1-2*i-1)/2 ;
  98       S2 ;
  99     }
 100   }
 101   if (p1%2 == 0) {
 102     i = (p1-2)/2 ;
 103     S1 ;
 104   }
 105   if (p1%2 == 0) {
 106     j = (p1-2)/2 ;
 107     S6(i = 1) ;
 108   }
 109 }
 110 if ((n >= 5) && (n <= 5)) {
 111   S3(i = 3) ;
 112   for (i=4;i<=n;i++) {
 113     j = -i+6 ;
 114     S2 ;
 115   }
 116 }
 117 if (n >= 6) {
 118   S3(i = 3) ;
 119   for (i=4;i<=5;i++) {
 120     j = -i+6 ;
 121     S2 ;
 122   }
 123 }
 124 for (p1=max(max(4*n-5,10),2*n+3);p1<=min(12,2*n+4);p1++) {
 125   for (p2=ceild(-p1+2,4);p2<=floord(-p1+2*n,2);p2++) {
 126     if (p1%2 == 0) {
 127       j = (p1+2*p2)/2 ;
 128       S4(i = -p2) ;
 129     }
 130   }
 131   for (p2=max(ceild(-p1+2*n+1,2),ceild(-p1+6,4));p2<=-1;p2++) {
 132     for (p3=1;p3<=-p2;p3++) {
 133       i = -p2+1 ;
 134       if (p1%2 == 0) {
 135         j = (p1+2*p2-2)/2 ;
 136         S5(k = p3) ;
 137       }
 138     }
 139   }
 140   if (p1%2 == 0) {
 141     j = (p1-4)/2 ;
 142     S6(i = 2) ;
 143   }
 144   for (i=max(3,ceild(p1+1,4));i<=n;i++) {
 145     if ((p1+1)%2 == 0) {
 146       j = (p1-2*i-1)/2 ;
 147       S2 ;
 148     }
 149   }
 150   if (p1%2 == 0) {
 151     i = (p1-2)/2 ;
 152     S1 ;
 153   }
 154 }
 155 for (p1=14;p1<=2*n+2;p1++) {
 156   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 157     if (p1%2 == 0) {
 158       j = (p1+2*p2)/2 ;
 159       S4(i = -p2) ;
 160     }
 161   }
 162   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
 163     if (p1%2 == 0) {
 164       j = (p1+2*p2)/2 ;
 165       S4(i = -p2) ;
 166     }
 167     for (p3=1;p3<=-p2;p3++) {
 168       i = -p2+1 ;
 169       if (p1%2 == 0) {
 170         j = (p1+2*p2-2)/2 ;
 171         S5(k = p3) ;
 172       }
 173     }
 174   }
 175   for (p2=ceild(-p1+10,4);p2<=-1;p2++) {
 176     if (p1%2 == 0) {
 177       j = (p1+2*p2)/2 ;
 178       S4(i = -p2) ;
 179     }
 180     i = -p2+2 ;
 181     if (p1%2 == 0) {
 182       j = (p1+2*p2-4)/2 ;
 183       S6 ;
 184     }
 185     for (p3=1;p3<=-p2;p3++) {
 186       i = -p2+1 ;
 187       if (p1%2 == 0) {
 188         j = (p1+2*p2-2)/2 ;
 189         S5(k = p3) ;
 190       }
 191     }
 192   }
 193   if (p1%2 == 0) {
 194     j = (p1-4)/2 ;
 195     S6(i = 2) ;
 196   }
 197   if ((p1+3)%4 == 0) {
 198     i = (p1-1)/4 ;
 199     S3 ;
 200   }
 201   for (i=ceild(p1+1,4);i<=floord(p1-3,2);i++) {
 202     if ((p1+1)%2 == 0) {
 203       j = (p1-2*i-1)/2 ;
 204       S2 ;
 205     }
 206   }
 207   if (p1%2 == 0) {
 208     i = (p1-2)/2 ;
 209     S1 ;
 210   }
 211   if (p1%2 == 0) {
 212     j = (p1-2)/2 ;
 213     S6(i = 1) ;
 214   }
 215 }
 216 if ((n >= 4) && (n <= 4)) {
 217   S3(i = 3) ;
 218   for (i=4;i<=n;i++) {
 219     j = -i+6 ;
 220     S2 ;
 221   }
 222 }
 223 for (p1=2*n+3;p1<=min(4*n-10,2*n+4);p1++) {
 224   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 225     if (p1%2 == 0) {
 226       j = (p1+2*p2)/2 ;
 227       S4(i = -p2) ;
 228     }
 229   }
 230   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
 231     if (p1%2 == 0) {
 232       j = (p1+2*p2)/2 ;
 233       S4(i = -p2) ;
 234     }
 235     for (p3=1;p3<=-p2;p3++) {
 236       i = -p2+1 ;
 237       if (p1%2 == 0) {
 238         j = (p1+2*p2-2)/2 ;
 239         S5(k = p3) ;
 240       }
 241     }
 242   }
 243   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n,2);p2++) {
 244     if (p1%2 == 0) {
 245       j = (p1+2*p2)/2 ;
 246       S4(i = -p2) ;
 247     }
 248     i = -p2+2 ;
 249     if (p1%2 == 0) {
 250       j = (p1+2*p2-4)/2 ;
 251       S6 ;
 252     }
 253     for (p3=1;p3<=-p2;p3++) {
 254       i = -p2+1 ;
 255       if (p1%2 == 0) {
 256         j = (p1+2*p2-2)/2 ;
 257         S5(k = p3) ;
 258       }
 259     }
 260   }
 261   for (p2=ceild(-p1+2*n+1,2);p2<=-1;p2++) {
 262     i = -p2+2 ;
 263     if (p1%2 == 0) {
 264       j = (p1+2*p2-4)/2 ;
 265       S6 ;
 266     }
 267     for (p3=1;p3<=-p2;p3++) {
 268       i = -p2+1 ;
 269       if (p1%2 == 0) {
 270         j = (p1+2*p2-2)/2 ;
 271         S5(k = p3) ;
 272       }
 273     }
 274   }
 275   if (p1%2 == 0) {
 276     j = (p1-4)/2 ;
 277     S6(i = 2) ;
 278   }
 279   if ((p1+3)%4 == 0) {
 280     i = (p1-1)/4 ;
 281     S3 ;
 282   }
 283   for (i=ceild(p1+1,4);i<=n;i++) {
 284     if ((p1+1)%2 == 0) {
 285       j = (p1-2*i-1)/2 ;
 286       S2 ;
 287     }
 288   }
 289   if (p1%2 == 0) {
 290     i = (p1-2)/2 ;
 291     S1 ;
 292   }
 293 }
 294 for (p1=2*n+5;p1<=min(4*n-10,2*n+58);p1++) {
 295   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 296     if (p1%2 == 0) {
 297       j = (p1+2*p2)/2 ;
 298       S4(i = -p2) ;
 299     }
 300   }
 301   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
 302     if (p1%2 == 0) {
 303       j = (p1+2*p2)/2 ;
 304       S4(i = -p2) ;
 305     }
 306     for (p3=1;p3<=-p2;p3++) {
 307       i = -p2+1 ;
 308       if (p1%2 == 0) {
 309         j = (p1+2*p2-2)/2 ;
 310         S5(k = p3) ;
 311       }
 312     }
 313   }
 314   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n,2);p2++) {
 315     if (p1%2 == 0) {
 316       j = (p1+2*p2)/2 ;
 317       S4(i = -p2) ;
 318     }
 319     i = -p2+2 ;
 320     if (p1%2 == 0) {
 321       j = (p1+2*p2-4)/2 ;
 322       S6 ;
 323     }
 324     for (p3=1;p3<=-p2;p3++) {
 325       i = -p2+1 ;
 326       if (p1%2 == 0) {
 327         j = (p1+2*p2-2)/2 ;
 328         S5(k = p3) ;
 329       }
 330     }
 331   }
 332   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 333     i = -p2+2 ;
 334     if (p1%2 == 0) {
 335       j = (p1+2*p2-4)/2 ;
 336       S6 ;
 337     }
 338     for (p3=1;p3<=-p2;p3++) {
 339       i = -p2+1 ;
 340       if (p1%2 == 0) {
 341         j = (p1+2*p2-2)/2 ;
 342         S5(k = p3) ;
 343       }
 344     }
 345   }
 346   for (p2=ceild(-p1+2*n+3,2);p2<=min(floord(-p1+2*n+4,2),-1);p2++) {
 347     i = -p2+2 ;
 348     if (p1%2 == 0) {
 349       j = (p1+2*p2-4)/2 ;
 350       S6 ;
 351     }
 352   }
 353   if ((p1+3)%4 == 0) {
 354     i = (p1-1)/4 ;
 355     S3 ;
 356   }
 357   for (i=ceild(p1+1,4);i<=n;i++) {
 358     if ((p1+1)%2 == 0) {
 359       j = (p1-2*i-1)/2 ;
 360       S2 ;
 361     }
 362   }
 363   if (p1%2 == 0) {
 364     i = (p1-2)/2 ;
 365     S1 ;
 366   }
 367 }
 368 for (p1=max(4*n-9,2*n+3);p1<=min(4*n-8,2*n+4);p1++) {
 369   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 370     if (p1%2 == 0) {
 371       j = (p1+2*p2)/2 ;
 372       S4(i = -p2) ;
 373     }
 374   }
 375   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 376     if (p1%2 == 0) {
 377       j = (p1+2*p2)/2 ;
 378       S4(i = -p2) ;
 379     }
 380     for (p3=1;p3<=-p2;p3++) {
 381       i = -p2+1 ;
 382       if (p1%2 == 0) {
 383         j = (p1+2*p2-2)/2 ;
 384         S5(k = p3) ;
 385       }
 386     }
 387   }
 388   for (p2=ceild(-p1+2*n+1,2);p2<=-1;p2++) {
 389     i = -p2+2 ;
 390     if (p1%2 == 0) {
 391       j = (p1+2*p2-4)/2 ;
 392       S6 ;
 393     }
 394     for (p3=1;p3<=-p2;p3++) {
 395       i = -p2+1 ;
 396       if (p1%2 == 0) {
 397         j = (p1+2*p2-2)/2 ;
 398         S5(k = p3) ;
 399       }
 400     }
 401   }
 402   if (p1%2 == 0) {
 403     j = (p1-4)/2 ;
 404     S6(i = 2) ;
 405   }
 406   if ((p1+3)%4 == 0) {
 407     i = (p1-1)/4 ;
 408     S3 ;
 409   }
 410   for (i=ceild(p1+1,4);i<=n;i++) {
 411     if ((p1+1)%2 == 0) {
 412       j = (p1-2*i-1)/2 ;
 413       S2 ;
 414     }
 415   }
 416   if (p1%2 == 0) {
 417     i = (p1-2)/2 ;
 418     S1 ;
 419   }
 420 }
 421 for (p1=max(4*n-7,14);p1<=2*n+4;p1++) {
 422   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 423     if (p1%2 == 0) {
 424       j = (p1+2*p2)/2 ;
 425       S4(i = -p2) ;
 426     }
 427   }
 428   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 429     if (p1%2 == 0) {
 430       j = (p1+2*p2)/2 ;
 431       S4(i = -p2) ;
 432     }
 433     for (p3=1;p3<=-p2;p3++) {
 434       i = -p2+1 ;
 435       if (p1%2 == 0) {
 436         j = (p1+2*p2-2)/2 ;
 437         S5(k = p3) ;
 438       }
 439     }
 440   }
 441   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+9,4);p2++) {
 442     for (p3=1;p3<=-p2;p3++) {
 443       i = -p2+1 ;
 444       if (p1%2 == 0) {
 445         j = (p1+2*p2-2)/2 ;
 446         S5(k = p3) ;
 447       }
 448     }
 449   }
 450   for (p2=ceild(-p1+10,4);p2<=-1;p2++) {
 451     i = -p2+2 ;
 452     if (p1%2 == 0) {
 453       j = (p1+2*p2-4)/2 ;
 454       S6 ;
 455     }
 456     for (p3=1;p3<=-p2;p3++) {
 457       i = -p2+1 ;
 458       if (p1%2 == 0) {
 459         j = (p1+2*p2-2)/2 ;
 460         S5(k = p3) ;
 461       }
 462     }
 463   }
 464   if (p1%2 == 0) {
 465     j = (p1-4)/2 ;
 466     S6(i = 2) ;
 467   }
 468   if ((p1+3)%4 == 0) {
 469     i = (p1-1)/4 ;
 470     S3 ;
 471   }
 472   for (i=ceild(p1+1,4);i<=n;i++) {
 473     if ((p1+1)%2 == 0) {
 474       j = (p1-2*i-1)/2 ;
 475       S2 ;
 476     }
 477   }
 478   if (p1%2 == 0) {
 479     i = (p1-2)/2 ;
 480     S1 ;
 481   }
 482 }
 483 for (p1=max(4*n-9,2*n+5);p1<=min(4*n-8,2*n+58);p1++) {
 484   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 485     if (p1%2 == 0) {
 486       j = (p1+2*p2)/2 ;
 487       S4(i = -p2) ;
 488     }
 489   }
 490   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 491     if (p1%2 == 0) {
 492       j = (p1+2*p2)/2 ;
 493       S4(i = -p2) ;
 494     }
 495     for (p3=1;p3<=-p2;p3++) {
 496       i = -p2+1 ;
 497       if (p1%2 == 0) {
 498         j = (p1+2*p2-2)/2 ;
 499         S5(k = p3) ;
 500       }
 501     }
 502   }
 503   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 504     i = -p2+2 ;
 505     if (p1%2 == 0) {
 506       j = (p1+2*p2-4)/2 ;
 507       S6 ;
 508     }
 509     for (p3=1;p3<=-p2;p3++) {
 510       i = -p2+1 ;
 511       if (p1%2 == 0) {
 512         j = (p1+2*p2-2)/2 ;
 513         S5(k = p3) ;
 514       }
 515     }
 516   }
 517   for (p2=ceild(-p1+2*n+3,2);p2<=min(floord(-p1+2*n+4,2),-1);p2++) {
 518     i = -p2+2 ;
 519     if (p1%2 == 0) {
 520       j = (p1+2*p2-4)/2 ;
 521       S6 ;
 522     }
 523   }
 524   if ((p1+3)%4 == 0) {
 525     i = (p1-1)/4 ;
 526     S3 ;
 527   }
 528   for (i=ceild(p1+1,4);i<=n;i++) {
 529     if ((p1+1)%2 == 0) {
 530       j = (p1-2*i-1)/2 ;
 531       S2 ;
 532     }
 533   }
 534   if (p1%2 == 0) {
 535     i = (p1-2)/2 ;
 536     S1 ;
 537   }
 538 }
 539 for (p1=max(4*n-7,2*n+5);p1<=min(4*n-6,2*n+58);p1++) {
 540   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 541     if (p1%2 == 0) {
 542       j = (p1+2*p2)/2 ;
 543       S4(i = -p2) ;
 544     }
 545   }
 546   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 547     if (p1%2 == 0) {
 548       j = (p1+2*p2)/2 ;
 549       S4(i = -p2) ;
 550     }
 551     for (p3=1;p3<=-p2;p3++) {
 552       i = -p2+1 ;
 553       if (p1%2 == 0) {
 554         j = (p1+2*p2-2)/2 ;
 555         S5(k = p3) ;
 556       }
 557     }
 558   }
 559   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+9,4);p2++) {
 560     for (p3=1;p3<=-p2;p3++) {
 561       i = -p2+1 ;
 562       if (p1%2 == 0) {
 563         j = (p1+2*p2-2)/2 ;
 564         S5(k = p3) ;
 565       }
 566     }
 567   }
 568   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n+2,2);p2++) {
 569     i = -p2+2 ;
 570     if (p1%2 == 0) {
 571       j = (p1+2*p2-4)/2 ;
 572       S6 ;
 573     }
 574     for (p3=1;p3<=-p2;p3++) {
 575       i = -p2+1 ;
 576       if (p1%2 == 0) {
 577         j = (p1+2*p2-2)/2 ;
 578         S5(k = p3) ;
 579       }
 580     }
 581   }
 582   for (p2=ceild(-p1+2*n+3,2);p2<=min(floord(-p1+2*n+4,2),-1);p2++) {
 583     i = -p2+2 ;
 584     if (p1%2 == 0) {
 585       j = (p1+2*p2-4)/2 ;
 586       S6 ;
 587     }
 588   }
 589   if ((p1+3)%4 == 0) {
 590     i = (p1-1)/4 ;
 591     S3 ;
 592   }
 593   for (i=ceild(p1+1,4);i<=n;i++) {
 594     if ((p1+1)%2 == 0) {
 595       j = (p1-2*i-1)/2 ;
 596       S2 ;
 597     }
 598   }
 599   if (p1%2 == 0) {
 600     i = (p1-2)/2 ;
 601     S1 ;
 602   }
 603 }
 604 for (p1=max(max(4*n-5,14),2*n+5);p1<=min(4*n-2,2*n+58);p1++) {
 605   for (p2=ceild(-p1+2,4);p2<=floord(-p1+2*n,2);p2++) {
 606     if (p1%2 == 0) {
 607       j = (p1+2*p2)/2 ;
 608       S4(i = -p2) ;
 609     }
 610   }
 611   for (p2=max(ceild(-p1+2*n+1,2),ceild(-p1+6,4));p2<=floord(-p1+2*n+2,2);p2++) {
 612     for (p3=1;p3<=-p2;p3++) {
 613       i = -p2+1 ;
 614       if (p1%2 == 0) {
 615         j = (p1+2*p2-2)/2 ;
 616         S5(k = p3) ;
 617       }
 618     }
 619   }
 620   for (p2=max(ceild(-p1+10,4),ceild(-p1+2*n+3,2));p2<=min(floord(-p1+2*n+4,2),-1);p2++) {
 621     i = -p2+2 ;
 622     if (p1%2 == 0) {
 623       j = (p1+2*p2-4)/2 ;
 624       S6 ;
 625     }
 626   }
 627   if ((p1+3)%4 == 0) {
 628     i = (p1-1)/4 ;
 629     S3 ;
 630   }
 631   for (i=ceild(p1+1,4);i<=n;i++) {
 632     if ((p1+1)%2 == 0) {
 633       j = (p1-2*i-1)/2 ;
 634       S2 ;
 635     }
 636   }
 637   if (p1%2 == 0) {
 638     i = (p1-2)/2 ;
 639     S1 ;
 640   }
 641 }
 642 if ((n <= 29) && (n >= 3)) {
 643   p1 = 4*n-1 ;
 644   j = n-1 ;
 645   S2(i = n) ;
 646 }
 647 if (n == 2) {
 648   S2(i = 2,j = 1) ;
 649 }
 650 for (p1=2*n+59;p1<=4*n-10;p1++) {
 651   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 652     if (p1%2 == 0) {
 653       j = (p1+2*p2)/2 ;
 654       S4(i = -p2) ;
 655     }
 656   }
 657   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
 658     if (p1%2 == 0) {
 659       j = (p1+2*p2)/2 ;
 660       S4(i = -p2) ;
 661     }
 662     for (p3=1;p3<=-p2;p3++) {
 663       i = -p2+1 ;
 664       if (p1%2 == 0) {
 665         j = (p1+2*p2-2)/2 ;
 666         S5(k = p3) ;
 667       }
 668     }
 669   }
 670   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n,2);p2++) {
 671     if (p1%2 == 0) {
 672       j = (p1+2*p2)/2 ;
 673       S4(i = -p2) ;
 674     }
 675     i = -p2+2 ;
 676     if (p1%2 == 0) {
 677       j = (p1+2*p2-4)/2 ;
 678       S6 ;
 679     }
 680     for (p3=1;p3<=-p2;p3++) {
 681       i = -p2+1 ;
 682       if (p1%2 == 0) {
 683         j = (p1+2*p2-2)/2 ;
 684         S5(k = p3) ;
 685       }
 686     }
 687   }
 688   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 689     i = -p2+2 ;
 690     if (p1%2 == 0) {
 691       j = (p1+2*p2-4)/2 ;
 692       S6 ;
 693     }
 694     for (p3=1;p3<=-p2;p3++) {
 695       i = -p2+1 ;
 696       if (p1%2 == 0) {
 697         j = (p1+2*p2-2)/2 ;
 698         S5(k = p3) ;
 699       }
 700     }
 701   }
 702   for (p2=ceild(-p1+2*n+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 703     i = -p2+2 ;
 704     if (p1%2 == 0) {
 705       j = (p1+2*p2-4)/2 ;
 706       S6 ;
 707     }
 708   }
 709   if ((p1+3)%4 == 0) {
 710     i = (p1-1)/4 ;
 711     S3 ;
 712   }
 713   for (i=ceild(p1+1,4);i<=n;i++) {
 714     if ((p1+1)%2 == 0) {
 715       j = (p1-2*i-1)/2 ;
 716       S2 ;
 717     }
 718   }
 719 }
 720 for (p1=max(4*n-9,2*n+59);p1<=4*n-8;p1++) {
 721   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 722     if (p1%2 == 0) {
 723       j = (p1+2*p2)/2 ;
 724       S4(i = -p2) ;
 725     }
 726   }
 727   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 728     if (p1%2 == 0) {
 729       j = (p1+2*p2)/2 ;
 730       S4(i = -p2) ;
 731     }
 732     for (p3=1;p3<=-p2;p3++) {
 733       i = -p2+1 ;
 734       if (p1%2 == 0) {
 735         j = (p1+2*p2-2)/2 ;
 736         S5(k = p3) ;
 737       }
 738     }
 739   }
 740   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 741     i = -p2+2 ;
 742     if (p1%2 == 0) {
 743       j = (p1+2*p2-4)/2 ;
 744       S6 ;
 745     }
 746     for (p3=1;p3<=-p2;p3++) {
 747       i = -p2+1 ;
 748       if (p1%2 == 0) {
 749         j = (p1+2*p2-2)/2 ;
 750         S5(k = p3) ;
 751       }
 752     }
 753   }
 754   for (p2=ceild(-p1+2*n+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 755     i = -p2+2 ;
 756     if (p1%2 == 0) {
 757       j = (p1+2*p2-4)/2 ;
 758       S6 ;
 759     }
 760   }
 761   if ((p1+3)%4 == 0) {
 762     i = (p1-1)/4 ;
 763     S3 ;
 764   }
 765   for (i=ceild(p1+1,4);i<=n;i++) {
 766     if ((p1+1)%2 == 0) {
 767       j = (p1-2*i-1)/2 ;
 768       S2 ;
 769     }
 770   }
 771 }
 772 for (p1=max(4*n-7,2*n+59);p1<=4*n-6;p1++) {
 773   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 774     if (p1%2 == 0) {
 775       j = (p1+2*p2)/2 ;
 776       S4(i = -p2) ;
 777     }
 778   }
 779   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 780     if (p1%2 == 0) {
 781       j = (p1+2*p2)/2 ;
 782       S4(i = -p2) ;
 783     }
 784     for (p3=1;p3<=-p2;p3++) {
 785       i = -p2+1 ;
 786       if (p1%2 == 0) {
 787         j = (p1+2*p2-2)/2 ;
 788         S5(k = p3) ;
 789       }
 790     }
 791   }
 792   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+9,4);p2++) {
 793     for (p3=1;p3<=-p2;p3++) {
 794       i = -p2+1 ;
 795       if (p1%2 == 0) {
 796         j = (p1+2*p2-2)/2 ;
 797         S5(k = p3) ;
 798       }
 799     }
 800   }
 801   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n+2,2);p2++) {
 802     i = -p2+2 ;
 803     if (p1%2 == 0) {
 804       j = (p1+2*p2-4)/2 ;
 805       S6 ;
 806     }
 807     for (p3=1;p3<=-p2;p3++) {
 808       i = -p2+1 ;
 809       if (p1%2 == 0) {
 810         j = (p1+2*p2-2)/2 ;
 811         S5(k = p3) ;
 812       }
 813     }
 814   }
 815   for (p2=ceild(-p1+2*n+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 816     i = -p2+2 ;
 817     if (p1%2 == 0) {
 818       j = (p1+2*p2-4)/2 ;
 819       S6 ;
 820     }
 821   }
 822   if ((p1+3)%4 == 0) {
 823     i = (p1-1)/4 ;
 824     S3 ;
 825   }
 826   for (i=ceild(p1+1,4);i<=n;i++) {
 827     if ((p1+1)%2 == 0) {
 828       j = (p1-2*i-1)/2 ;
 829       S2 ;
 830     }
 831   }
 832 }
 833 for (p1=max(4*n-5,2*n+59);p1<=4*n-2;p1++) {
 834   for (p2=ceild(-p1+2,4);p2<=floord(-p1+2*n,2);p2++) {
 835     if (p1%2 == 0) {
 836       j = (p1+2*p2)/2 ;
 837       S4(i = -p2) ;
 838     }
 839   }
 840   for (p2=max(ceild(-p1+2*n+1,2),ceild(-p1+6,4));p2<=floord(-p1+2*n+2,2);p2++) {
 841     for (p3=1;p3<=-p2;p3++) {
 842       i = -p2+1 ;
 843       if (p1%2 == 0) {
 844         j = (p1+2*p2-2)/2 ;
 845         S5(k = p3) ;
 846       }
 847     }
 848   }
 849   for (p2=max(ceild(-p1+10,4),ceild(-p1+2*n+3,2));p2<=floord(-p1+2*n+4,2);p2++) {
 850     i = -p2+2 ;
 851     if (p1%2 == 0) {
 852       j = (p1+2*p2-4)/2 ;
 853       S6 ;
 854     }
 855   }
 856   if ((p1+3)%4 == 0) {
 857     i = (p1-1)/4 ;
 858     S3 ;
 859   }
 860   for (i=ceild(p1+1,4);i<=n;i++) {
 861     if ((p1+1)%2 == 0) {
 862       j = (p1-2*i-1)/2 ;
 863       S2 ;
 864     }
 865   }
 866 }
 867 for (p1=max(4*n,7);p1<=min(4*n+1,2*n+58);p1++) {
 868   if ((p1+3)%4 == 0) {
 869     i = (p1-1)/4 ;
 870     S3 ;
 871   }
 872   if (p1%2 == 0) {
 873     i = (p1-2)/2 ;
 874     S1 ;
 875   }
 876 }
 877 if (n >= 30) {
 878   p1 = 4*n-1 ;
 879   j = n-1 ;
 880   S2(i = n) ;
 881 }
 882 for (p1=max(max(-54*n+4,4*n+2),7);p1<=2*n+58;p1++) {
 883   if (p1%2 == 0) {
 884     i = (p1-2)/2 ;
 885     S1 ;
 886   }
 887 }
 888 for (p1=max(4*n,2*n+59);p1<=4*n+1;p1++) {
 889   if ((p1+3)%4 == 0) {
 890     i = (p1-1)/4 ;
 891     S3 ;
 892   }
 893 }