test/vivien2.c

   1 /* Generated from vivien2.cloog by CLooG v0.14.0 gmp bits in 0.51s. */
   2 /* CLooG asked for 1188 KBytes. */
   3 for (p1=-54*n+4;p1<=4;p1++) {
   4   if (p1%2 == 0) {
   5     i = (p1-2)/2 ;
   6     S1 ;
   7   }
   8 }
   9 S3(i = 1) ;
  10 S4(i = 1,j = 2) ;
  11 S1(i = 2) ;
  12 S6(i = 1,j = 2) ;
  13 for (p1=7;p1<=9;p1++) {
  14   for (p2=ceild(-p1+2,4);p2<=-1;p2++) {
  15     if (p1%2 == 0) {
  16       j = (p1+2*p2)/2 ;
  17       S4(i = -p2) ;
  18     }
  19   }
  20   if ((p1+3)%4 == 0) {
  21     i = (p1-1)/4 ;
  22     S3 ;
  23   }
  24   for (i=ceild(p1+1,4);i<=floord(p1-3,2);i++) {
  25     if ((p1+1)%2 == 0) {
  26       j = (p1-2*i-1)/2 ;
  27       S2 ;
  28     }
  29   }
  30   if (p1%2 == 0) {
  31     i = (p1-2)/2 ;
  32     S1 ;
  33   }
  34   if (p1%2 == 0) {
  35     j = (p1-2)/2 ;
  36     S6(i = 1) ;
  37   }
  38 }
  39 for (p1=10;p1<=12;p1++) {
  40   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
  41     if (p1%2 == 0) {
  42       j = (p1+2*p2)/2 ;
  43       S4(i = -p2) ;
  44     }
  45   }
  46   for (p2=ceild(-p1+6,4);p2<=-1;p2++) {
  47     if (p1%2 == 0) {
  48       j = (p1+2*p2)/2 ;
  49       S4(i = -p2) ;
  50     }
  51     for (p3=1;p3<=-p2;p3++) {
  52       i = -p2+1 ;
  53       if (p1%2 == 0) {
  54         j = (p1+2*p2-2)/2 ;
  55         S5(k = p3) ;
  56       }
  57     }
  58   }
  59   if (p1%2 == 0) {
  60     j = (p1-4)/2 ;
  61     S6(i = 2) ;
  62   }
  63   for (i=max(ceild(p1+1,4),3);i<=floord(p1-3,2);i++) {
  64     if ((p1+1)%2 == 0) {
  65       j = (p1-2*i-1)/2 ;
  66       S2 ;
  67     }
  68   }
  69   if (p1%2 == 0) {
  70     i = (p1-2)/2 ;
  71     S1 ;
  72   }
  73   if (p1%2 == 0) {
  74     j = (p1-2)/2 ;
  75     S6(i = 1) ;
  76   }
  77 }
  78 S3(i = 3) ;
  79 for (i=4;i<=5;i++) {
  80   j = -i+6 ;
  81   S2 ;
  82 }
  83 for (p1=14;p1<=2*n+2;p1++) {
  84   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
  85     if (p1%2 == 0) {
  86       j = (p1+2*p2)/2 ;
  87       S4(i = -p2) ;
  88     }
  89   }
  90   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
  91     if (p1%2 == 0) {
  92       j = (p1+2*p2)/2 ;
  93       S4(i = -p2) ;
  94     }
  95     for (p3=1;p3<=-p2;p3++) {
  96       i = -p2+1 ;
  97       if (p1%2 == 0) {
  98         j = (p1+2*p2-2)/2 ;
  99         S5(k = p3) ;
 100       }
 101     }
 102   }
 103   for (p2=ceild(-p1+10,4);p2<=-1;p2++) {
 104     if (p1%2 == 0) {
 105       j = (p1+2*p2)/2 ;
 106       S4(i = -p2) ;
 107     }
 108     i = -p2+2 ;
 109     if (p1%2 == 0) {
 110       j = (p1+2*p2-4)/2 ;
 111       S6 ;
 112     }
 113     for (p3=1;p3<=-p2;p3++) {
 114       i = -p2+1 ;
 115       if (p1%2 == 0) {
 116         j = (p1+2*p2-2)/2 ;
 117         S5(k = p3) ;
 118       }
 119     }
 120   }
 121   if (p1%2 == 0) {
 122     j = (p1-4)/2 ;
 123     S6(i = 2) ;
 124   }
 125   if ((p1+3)%4 == 0) {
 126     i = (p1-1)/4 ;
 127     S3 ;
 128   }
 129   for (i=ceild(p1+1,4);i<=floord(p1-3,2);i++) {
 130     if ((p1+1)%2 == 0) {
 131       j = (p1-2*i-1)/2 ;
 132       S2 ;
 133     }
 134   }
 135   if (p1%2 == 0) {
 136     i = (p1-2)/2 ;
 137     S1 ;
 138   }
 139   if (p1%2 == 0) {
 140     j = (p1-2)/2 ;
 141     S6(i = 1) ;
 142   }
 143 }
 144 for (p1=2*n+3;p1<=2*n+4;p1++) {
 145   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 146     if (p1%2 == 0) {
 147       j = (p1+2*p2)/2 ;
 148       S4(i = -p2) ;
 149     }
 150   }
 151   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
 152     if (p1%2 == 0) {
 153       j = (p1+2*p2)/2 ;
 154       S4(i = -p2) ;
 155     }
 156     for (p3=1;p3<=-p2;p3++) {
 157       i = -p2+1 ;
 158       if (p1%2 == 0) {
 159         j = (p1+2*p2-2)/2 ;
 160         S5(k = p3) ;
 161       }
 162     }
 163   }
 164   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n,2);p2++) {
 165     if (p1%2 == 0) {
 166       j = (p1+2*p2)/2 ;
 167       S4(i = -p2) ;
 168     }
 169     i = -p2+2 ;
 170     if (p1%2 == 0) {
 171       j = (p1+2*p2-4)/2 ;
 172       S6 ;
 173     }
 174     for (p3=1;p3<=-p2;p3++) {
 175       i = -p2+1 ;
 176       if (p1%2 == 0) {
 177         j = (p1+2*p2-2)/2 ;
 178         S5(k = p3) ;
 179       }
 180     }
 181   }
 182   for (p2=ceild(-p1+2*n+1,2);p2<=-1;p2++) {
 183     i = -p2+2 ;
 184     if (p1%2 == 0) {
 185       j = (p1+2*p2-4)/2 ;
 186       S6 ;
 187     }
 188     for (p3=1;p3<=-p2;p3++) {
 189       i = -p2+1 ;
 190       if (p1%2 == 0) {
 191         j = (p1+2*p2-2)/2 ;
 192         S5(k = p3) ;
 193       }
 194     }
 195   }
 196   if (p1%2 == 0) {
 197     j = (p1-4)/2 ;
 198     S6(i = 2) ;
 199   }
 200   if ((p1+3)%4 == 0) {
 201     i = (p1-1)/4 ;
 202     S3 ;
 203   }
 204   for (i=ceild(p1+1,4);i<=n;i++) {
 205     if ((p1+1)%2 == 0) {
 206       j = (p1-2*i-1)/2 ;
 207       S2 ;
 208     }
 209   }
 210   if (p1%2 == 0) {
 211     i = (p1-2)/2 ;
 212     S1 ;
 213   }
 214 }
 215 for (p1=2*n+5;p1<=min(4*n-10,2*n+58);p1++) {
 216   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 217     if (p1%2 == 0) {
 218       j = (p1+2*p2)/2 ;
 219       S4(i = -p2) ;
 220     }
 221   }
 222   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
 223     if (p1%2 == 0) {
 224       j = (p1+2*p2)/2 ;
 225       S4(i = -p2) ;
 226     }
 227     for (p3=1;p3<=-p2;p3++) {
 228       i = -p2+1 ;
 229       if (p1%2 == 0) {
 230         j = (p1+2*p2-2)/2 ;
 231         S5(k = p3) ;
 232       }
 233     }
 234   }
 235   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n,2);p2++) {
 236     if (p1%2 == 0) {
 237       j = (p1+2*p2)/2 ;
 238       S4(i = -p2) ;
 239     }
 240     i = -p2+2 ;
 241     if (p1%2 == 0) {
 242       j = (p1+2*p2-4)/2 ;
 243       S6 ;
 244     }
 245     for (p3=1;p3<=-p2;p3++) {
 246       i = -p2+1 ;
 247       if (p1%2 == 0) {
 248         j = (p1+2*p2-2)/2 ;
 249         S5(k = p3) ;
 250       }
 251     }
 252   }
 253   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 254     i = -p2+2 ;
 255     if (p1%2 == 0) {
 256       j = (p1+2*p2-4)/2 ;
 257       S6 ;
 258     }
 259     for (p3=1;p3<=-p2;p3++) {
 260       i = -p2+1 ;
 261       if (p1%2 == 0) {
 262         j = (p1+2*p2-2)/2 ;
 263         S5(k = p3) ;
 264       }
 265     }
 266   }
 267   for (p2=ceild(-p1+2*n+3,2);p2<=min(floord(-p1+2*n+4,2),-1);p2++) {
 268     i = -p2+2 ;
 269     if (p1%2 == 0) {
 270       j = (p1+2*p2-4)/2 ;
 271       S6 ;
 272     }
 273   }
 274   if ((p1+3)%4 == 0) {
 275     i = (p1-1)/4 ;
 276     S3 ;
 277   }
 278   for (i=ceild(p1+1,4);i<=n;i++) {
 279     if ((p1+1)%2 == 0) {
 280       j = (p1-2*i-1)/2 ;
 281       S2 ;
 282     }
 283   }
 284   if (p1%2 == 0) {
 285     i = (p1-2)/2 ;
 286     S1 ;
 287   }
 288 }
 289 for (p1=4*n-9;p1<=min(4*n-8,2*n+58);p1++) {
 290   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 291     if (p1%2 == 0) {
 292       j = (p1+2*p2)/2 ;
 293       S4(i = -p2) ;
 294     }
 295   }
 296   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 297     if (p1%2 == 0) {
 298       j = (p1+2*p2)/2 ;
 299       S4(i = -p2) ;
 300     }
 301     for (p3=1;p3<=-p2;p3++) {
 302       i = -p2+1 ;
 303       if (p1%2 == 0) {
 304         j = (p1+2*p2-2)/2 ;
 305         S5(k = p3) ;
 306       }
 307     }
 308   }
 309   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 310     i = -p2+2 ;
 311     if (p1%2 == 0) {
 312       j = (p1+2*p2-4)/2 ;
 313       S6 ;
 314     }
 315     for (p3=1;p3<=-p2;p3++) {
 316       i = -p2+1 ;
 317       if (p1%2 == 0) {
 318         j = (p1+2*p2-2)/2 ;
 319         S5(k = p3) ;
 320       }
 321     }
 322   }
 323   for (p2=ceild(-p1+2*n+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 324     i = -p2+2 ;
 325     if (p1%2 == 0) {
 326       j = (p1+2*p2-4)/2 ;
 327       S6 ;
 328     }
 329   }
 330   if ((p1+3)%4 == 0) {
 331     i = (p1-1)/4 ;
 332     S3 ;
 333   }
 334   for (i=ceild(p1+1,4);i<=n;i++) {
 335     if ((p1+1)%2 == 0) {
 336       j = (p1-2*i-1)/2 ;
 337       S2 ;
 338     }
 339   }
 340   if (p1%2 == 0) {
 341     i = (p1-2)/2 ;
 342     S1 ;
 343   }
 344 }
 345 for (p1=4*n-7;p1<=min(4*n-6,2*n+58);p1++) {
 346   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 347     if (p1%2 == 0) {
 348       j = (p1+2*p2)/2 ;
 349       S4(i = -p2) ;
 350     }
 351   }
 352   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 353     if (p1%2 == 0) {
 354       j = (p1+2*p2)/2 ;
 355       S4(i = -p2) ;
 356     }
 357     for (p3=1;p3<=-p2;p3++) {
 358       i = -p2+1 ;
 359       if (p1%2 == 0) {
 360         j = (p1+2*p2-2)/2 ;
 361         S5(k = p3) ;
 362       }
 363     }
 364   }
 365   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+9,4);p2++) {
 366     for (p3=1;p3<=-p2;p3++) {
 367       i = -p2+1 ;
 368       if (p1%2 == 0) {
 369         j = (p1+2*p2-2)/2 ;
 370         S5(k = p3) ;
 371       }
 372     }
 373   }
 374   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n+2,2);p2++) {
 375     i = -p2+2 ;
 376     if (p1%2 == 0) {
 377       j = (p1+2*p2-4)/2 ;
 378       S6 ;
 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+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 389     i = -p2+2 ;
 390     if (p1%2 == 0) {
 391       j = (p1+2*p2-4)/2 ;
 392       S6 ;
 393     }
 394   }
 395   if ((p1+3)%4 == 0) {
 396     i = (p1-1)/4 ;
 397     S3 ;
 398   }
 399   for (i=ceild(p1+1,4);i<=n;i++) {
 400     if ((p1+1)%2 == 0) {
 401       j = (p1-2*i-1)/2 ;
 402       S2 ;
 403     }
 404   }
 405   if (p1%2 == 0) {
 406     i = (p1-2)/2 ;
 407     S1 ;
 408   }
 409 }
 410 for (p1=4*n-5;p1<=2*n+58;p1++) {
 411   for (p2=ceild(-p1+2,4);p2<=floord(-p1+2*n,2);p2++) {
 412     if (p1%2 == 0) {
 413       j = (p1+2*p2)/2 ;
 414       S4(i = -p2) ;
 415     }
 416   }
 417   for (p2=max(ceild(-p1+2*n+1,2),ceild(-p1+6,4));p2<=floord(-p1+2*n+2,2);p2++) {
 418     for (p3=1;p3<=-p2;p3++) {
 419       i = -p2+1 ;
 420       if (p1%2 == 0) {
 421         j = (p1+2*p2-2)/2 ;
 422         S5(k = p3) ;
 423       }
 424     }
 425   }
 426   for (p2=max(ceild(-p1+2*n+3,2),ceild(-p1+10,4));p2<=floord(-p1+2*n+4,2);p2++) {
 427     i = -p2+2 ;
 428     if (p1%2 == 0) {
 429       j = (p1+2*p2-4)/2 ;
 430       S6 ;
 431     }
 432   }
 433   if ((p1+3)%4 == 0) {
 434     i = (p1-1)/4 ;
 435     S3 ;
 436   }
 437   for (i=ceild(p1+1,4);i<=n;i++) {
 438     if ((p1+1)%2 == 0) {
 439       j = (p1-2*i-1)/2 ;
 440       S2 ;
 441     }
 442   }
 443   if (p1%2 == 0) {
 444     i = (p1-2)/2 ;
 445     S1 ;
 446   }
 447 }
 448 for (p1=2*n+59;p1<=4*n-10;p1++) {
 449   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 450     if (p1%2 == 0) {
 451       j = (p1+2*p2)/2 ;
 452       S4(i = -p2) ;
 453     }
 454   }
 455   for (p2=ceild(-p1+6,4);p2<=floord(-p1+9,4);p2++) {
 456     if (p1%2 == 0) {
 457       j = (p1+2*p2)/2 ;
 458       S4(i = -p2) ;
 459     }
 460     for (p3=1;p3<=-p2;p3++) {
 461       i = -p2+1 ;
 462       if (p1%2 == 0) {
 463         j = (p1+2*p2-2)/2 ;
 464         S5(k = p3) ;
 465       }
 466     }
 467   }
 468   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n,2);p2++) {
 469     if (p1%2 == 0) {
 470       j = (p1+2*p2)/2 ;
 471       S4(i = -p2) ;
 472     }
 473     i = -p2+2 ;
 474     if (p1%2 == 0) {
 475       j = (p1+2*p2-4)/2 ;
 476       S6 ;
 477     }
 478     for (p3=1;p3<=-p2;p3++) {
 479       i = -p2+1 ;
 480       if (p1%2 == 0) {
 481         j = (p1+2*p2-2)/2 ;
 482         S5(k = p3) ;
 483       }
 484     }
 485   }
 486   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 487     i = -p2+2 ;
 488     if (p1%2 == 0) {
 489       j = (p1+2*p2-4)/2 ;
 490       S6 ;
 491     }
 492     for (p3=1;p3<=-p2;p3++) {
 493       i = -p2+1 ;
 494       if (p1%2 == 0) {
 495         j = (p1+2*p2-2)/2 ;
 496         S5(k = p3) ;
 497       }
 498     }
 499   }
 500   for (p2=ceild(-p1+2*n+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 501     i = -p2+2 ;
 502     if (p1%2 == 0) {
 503       j = (p1+2*p2-4)/2 ;
 504       S6 ;
 505     }
 506   }
 507   if ((p1+3)%4 == 0) {
 508     i = (p1-1)/4 ;
 509     S3 ;
 510   }
 511   for (i=ceild(p1+1,4);i<=n;i++) {
 512     if ((p1+1)%2 == 0) {
 513       j = (p1-2*i-1)/2 ;
 514       S2 ;
 515     }
 516   }
 517 }
 518 for (p1=max(4*n-9,2*n+59);p1<=4*n-8;p1++) {
 519   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 520     if (p1%2 == 0) {
 521       j = (p1+2*p2)/2 ;
 522       S4(i = -p2) ;
 523     }
 524   }
 525   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 526     if (p1%2 == 0) {
 527       j = (p1+2*p2)/2 ;
 528       S4(i = -p2) ;
 529     }
 530     for (p3=1;p3<=-p2;p3++) {
 531       i = -p2+1 ;
 532       if (p1%2 == 0) {
 533         j = (p1+2*p2-2)/2 ;
 534         S5(k = p3) ;
 535       }
 536     }
 537   }
 538   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+2*n+2,2);p2++) {
 539     i = -p2+2 ;
 540     if (p1%2 == 0) {
 541       j = (p1+2*p2-4)/2 ;
 542       S6 ;
 543     }
 544     for (p3=1;p3<=-p2;p3++) {
 545       i = -p2+1 ;
 546       if (p1%2 == 0) {
 547         j = (p1+2*p2-2)/2 ;
 548         S5(k = p3) ;
 549       }
 550     }
 551   }
 552   for (p2=ceild(-p1+2*n+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 553     i = -p2+2 ;
 554     if (p1%2 == 0) {
 555       j = (p1+2*p2-4)/2 ;
 556       S6 ;
 557     }
 558   }
 559   if ((p1+3)%4 == 0) {
 560     i = (p1-1)/4 ;
 561     S3 ;
 562   }
 563   for (i=ceild(p1+1,4);i<=n;i++) {
 564     if ((p1+1)%2 == 0) {
 565       j = (p1-2*i-1)/2 ;
 566       S2 ;
 567     }
 568   }
 569 }
 570 for (p1=max(4*n-7,2*n+59);p1<=4*n-6;p1++) {
 571   for (p2=ceild(-p1+2,4);p2<=floord(-p1+5,4);p2++) {
 572     if (p1%2 == 0) {
 573       j = (p1+2*p2)/2 ;
 574       S4(i = -p2) ;
 575     }
 576   }
 577   for (p2=ceild(-p1+6,4);p2<=floord(-p1+2*n,2);p2++) {
 578     if (p1%2 == 0) {
 579       j = (p1+2*p2)/2 ;
 580       S4(i = -p2) ;
 581     }
 582     for (p3=1;p3<=-p2;p3++) {
 583       i = -p2+1 ;
 584       if (p1%2 == 0) {
 585         j = (p1+2*p2-2)/2 ;
 586         S5(k = p3) ;
 587       }
 588     }
 589   }
 590   for (p2=ceild(-p1+2*n+1,2);p2<=floord(-p1+9,4);p2++) {
 591     for (p3=1;p3<=-p2;p3++) {
 592       i = -p2+1 ;
 593       if (p1%2 == 0) {
 594         j = (p1+2*p2-2)/2 ;
 595         S5(k = p3) ;
 596       }
 597     }
 598   }
 599   for (p2=ceild(-p1+10,4);p2<=floord(-p1+2*n+2,2);p2++) {
 600     i = -p2+2 ;
 601     if (p1%2 == 0) {
 602       j = (p1+2*p2-4)/2 ;
 603       S6 ;
 604     }
 605     for (p3=1;p3<=-p2;p3++) {
 606       i = -p2+1 ;
 607       if (p1%2 == 0) {
 608         j = (p1+2*p2-2)/2 ;
 609         S5(k = p3) ;
 610       }
 611     }
 612   }
 613   for (p2=ceild(-p1+2*n+3,2);p2<=floord(-p1+2*n+4,2);p2++) {
 614     i = -p2+2 ;
 615     if (p1%2 == 0) {
 616       j = (p1+2*p2-4)/2 ;
 617       S6 ;
 618     }
 619   }
 620   if ((p1+3)%4 == 0) {
 621     i = (p1-1)/4 ;
 622     S3 ;
 623   }
 624   for (i=ceild(p1+1,4);i<=n;i++) {
 625     if ((p1+1)%2 == 0) {
 626       j = (p1-2*i-1)/2 ;
 627       S2 ;
 628     }
 629   }
 630 }
 631 for (p1=max(4*n-5,2*n+59);p1<=4*n-2;p1++) {
 632   for (p2=ceild(-p1+2,4);p2<=floord(-p1+2*n,2);p2++) {
 633     if (p1%2 == 0) {
 634       j = (p1+2*p2)/2 ;
 635       S4(i = -p2) ;
 636     }
 637   }
 638   for (p2=max(ceild(-p1+2*n+1,2),ceild(-p1+6,4));p2<=floord(-p1+2*n+2,2);p2++) {
 639     for (p3=1;p3<=-p2;p3++) {
 640       i = -p2+1 ;
 641       if (p1%2 == 0) {
 642         j = (p1+2*p2-2)/2 ;
 643         S5(k = p3) ;
 644       }
 645     }
 646   }
 647   for (p2=max(ceild(-p1+10,4),ceild(-p1+2*n+3,2));p2<=floord(-p1+2*n+4,2);p2++) {
 648     i = -p2+2 ;
 649     if (p1%2 == 0) {
 650       j = (p1+2*p2-4)/2 ;
 651       S6 ;
 652     }
 653   }
 654   if ((p1+3)%4 == 0) {
 655     i = (p1-1)/4 ;
 656     S3 ;
 657   }
 658   for (i=ceild(p1+1,4);i<=n;i++) {
 659     if ((p1+1)%2 == 0) {
 660       j = (p1-2*i-1)/2 ;
 661       S2 ;
 662     }
 663   }
 664 }
 665 p1 = 4*n-1 ;
 666 j = n-1 ;
 667 S2(i = n) ;
 668 for (p1=4*n;p1<=4*n+1;p1++) {
 669   if ((p1+3)%4 == 0) {
 670     i = (p1-1)/4 ;
 671     S3 ;
 672   }
 673 }