← Back to branch summary

~maddevelopers/mg5amcnlo/find_granny_inverted_border_enhanced

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/polynomial.f

  • Committer: Rikkert Frederix
  • Date: 2016-03-09 13:16:26 UTC
  • mfrom: (253.48.95 2.3.4)
  • Revision ID: frederix@physik.uzh.ch-20160309131626-1j7l4vnn3kw5mjhs
merge with latest 2.3.4

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/polynomial.f
15
15
      PARAMETER (IMAG1=(ZERO,ONE))
16
16
      COMPLEX*16 CMPLX_ZERO
17
17
      PARAMETER (CMPLX_ZERO=(ZERO,ZERO))
18
 
      INTEGER MAXLWFSIZE
19
 
      PARAMETER (MAXLWFSIZE=4)
20
 
      INTEGER LOOPMAXCOEFS
21
 
      PARAMETER (LOOPMAXCOEFS=35)
 
18
      INCLUDE 'loop_max_coefs.inc'
 
19
      INCLUDE 'coef_specs.inc'
22
20
      INTEGER    NCOLORROWS
23
21
      PARAMETER (NCOLORROWS=129)
24
22
      INTEGER    NLOOPGROUPS
98
96
      END
99
97
 
100
98
 
 
99
 
101
100
C     THE SUBROUTINE TO CREATE THE COEFFICIENTS FROM LAST LOOP WF AND 
102
101
C     MULTIPLY BY THE BORN
103
102
 
115
114
      PARAMETER (IMAG1=(ZERO,ONE))
116
115
      COMPLEX*32 CMPLX_ZERO
117
116
      PARAMETER (CMPLX_ZERO=(ZERO,ZERO))
118
 
      INTEGER MAXLWFSIZE
119
 
      PARAMETER (MAXLWFSIZE=4)
120
 
      INTEGER LOOPMAXCOEFS
121
 
      PARAMETER (LOOPMAXCOEFS=35)
 
117
      INCLUDE 'loop_max_coefs.inc'
 
118
      INCLUDE 'coef_specs.inc'
122
119
      INTEGER    NCOLORROWS
123
120
      PARAMETER (NCOLORROWS=129)
124
121
      INTEGER    NLOOPGROUPS
198
195
      END
199
196
 
200
197
 
 
198
 
201
199
      SUBROUTINE ML5_0_EVAL_POLY(C,R,Q,OUT)
202
200
      INCLUDE 'coef_specs.inc'
203
 
      COMPLEX*16 C(0:LOOP_MAXCOEFS-1)
 
201
      INCLUDE 'loop_max_coefs.inc'
 
202
      COMPLEX*16 C(0:LOOPMAXCOEFS-1)
204
203
      INTEGER R
205
204
      COMPLEX*16 Q(0:3)
206
205
      COMPLEX*16 OUT
228
227
 
229
228
      SUBROUTINE MP_ML5_0_EVAL_POLY(C,R,Q,OUT)
230
229
      INCLUDE 'coef_specs.inc'
231
 
      COMPLEX*32 C(0:LOOP_MAXCOEFS-1)
 
230
      INCLUDE 'loop_max_coefs.inc'
 
231
      COMPLEX*32 C(0:LOOPMAXCOEFS-1)
232
232
      INTEGER R
233
233
      COMPLEX*32 Q(0:3)
234
234
      COMPLEX*32 OUT
256
256
 
257
257
      SUBROUTINE ML5_0_ADD_COEFS(A,RA,B,RB)
258
258
      INCLUDE 'coef_specs.inc'
 
259
      INCLUDE 'loop_max_coefs.inc'
259
260
      INTEGER I
260
 
      COMPLEX*16 A(0:LOOP_MAXCOEFS-1),B(0:LOOP_MAXCOEFS-1)
 
261
      COMPLEX*16 A(0:LOOPMAXCOEFS-1),B(0:LOOPMAXCOEFS-1)
261
262
      INTEGER RA,RB
262
263
 
263
264
      INTEGER NCOEF_R(0:3)
270
271
 
271
272
      SUBROUTINE MP_ML5_0_ADD_COEFS(A,RA,B,RB)
272
273
      INCLUDE 'coef_specs.inc'
 
274
      INCLUDE 'loop_max_coefs.inc'
273
275
      INTEGER I
274
 
      COMPLEX*32 A(0:LOOP_MAXCOEFS-1),B(0:LOOP_MAXCOEFS-1)
 
276
      COMPLEX*32 A(0:LOOPMAXCOEFS-1),B(0:LOOPMAXCOEFS-1)
275
277
      INTEGER RA,RB
276
278
 
277
279
      INTEGER NCOEF_R(0:3)
284
286
 
285
287
      SUBROUTINE ML5_0_MERGE_WL(WL,R,LCUT_SIZE,CONST,OUT)
286
288
      INCLUDE 'coef_specs.inc'
 
289
      INCLUDE 'loop_max_coefs.inc'
287
290
      INTEGER I,J
288
 
      COMPLEX*16 WL(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
291
      COMPLEX*16 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
289
292
      INTEGER R,LCUT_SIZE
290
293
      COMPLEX*16 CONST
291
 
      COMPLEX*16 OUT(0:LOOP_MAXCOEFS-1)
 
294
      COMPLEX*16 OUT(0:LOOPMAXCOEFS-1)
292
295
 
293
296
      INTEGER NCOEF_R(0:3)
294
297
      DATA NCOEF_R/1,5,15,35/
302
305
 
303
306
      SUBROUTINE MP_ML5_0_MERGE_WL(WL,R,LCUT_SIZE,CONST,OUT)
304
307
      INCLUDE 'coef_specs.inc'
 
308
      INCLUDE 'loop_max_coefs.inc'
305
309
      INTEGER I,J
306
 
      COMPLEX*32 WL(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
310
      COMPLEX*32 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
307
311
      INTEGER R,LCUT_SIZE
308
312
      COMPLEX*32 CONST
309
 
      COMPLEX*32 OUT(0:LOOP_MAXCOEFS-1)
 
313
      COMPLEX*32 OUT(0:LOOPMAXCOEFS-1)
310
314
 
311
315
      INTEGER NCOEF_R(0:3)
312
316
      DATA NCOEF_R/1,5,15,35/
321
325
      SUBROUTINE ML5_0_UPDATE_WL_0_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
322
326
     $ ,OUT)
323
327
      INCLUDE 'coef_specs.inc'
 
328
      INCLUDE 'loop_max_coefs.inc'
324
329
      INTEGER I,J,K
325
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
330
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
326
331
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
327
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
332
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
328
333
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
329
334
 
330
335
      DO I=1,LCUT_SIZE
346
351
      SUBROUTINE MP_ML5_0_UPDATE_WL_0_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
347
352
     $ ,OUT)
348
353
      INCLUDE 'coef_specs.inc'
 
354
      INCLUDE 'loop_max_coefs.inc'
349
355
      INTEGER I,J,K
350
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
356
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
351
357
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
352
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
358
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
353
359
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
354
360
 
355
361
      DO I=1,LCUT_SIZE
371
377
      SUBROUTINE ML5_0_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
372
378
     $ ,OUT)
373
379
      INCLUDE 'coef_specs.inc'
 
380
      INCLUDE 'loop_max_coefs.inc'
374
381
      INTEGER I,J,K
375
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
382
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
376
383
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
377
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
384
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
378
385
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
379
386
 
380
387
      DO I=1,LCUT_SIZE
392
399
      SUBROUTINE MP_ML5_0_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
393
400
     $ ,OUT)
394
401
      INCLUDE 'coef_specs.inc'
 
402
      INCLUDE 'loop_max_coefs.inc'
395
403
      INTEGER I,J,K
396
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
404
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
397
405
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
398
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
406
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
399
407
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
400
408
 
401
409
      DO I=1,LCUT_SIZE
413
421
      SUBROUTINE ML5_0_UPDATE_WL_3_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
414
422
     $ ,OUT)
415
423
      INCLUDE 'coef_specs.inc'
 
424
      INCLUDE 'loop_max_coefs.inc'
416
425
      INTEGER I,J,K
417
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
426
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
418
427
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
419
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
428
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
420
429
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
421
430
 
422
431
      DO I=1,LCUT_SIZE
468
477
      SUBROUTINE MP_ML5_0_UPDATE_WL_3_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
469
478
     $ ,OUT)
470
479
      INCLUDE 'coef_specs.inc'
 
480
      INCLUDE 'loop_max_coefs.inc'
471
481
      INTEGER I,J,K
472
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
482
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
473
483
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
474
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
484
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
475
485
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
476
486
 
477
487
      DO I=1,LCUT_SIZE
523
533
      SUBROUTINE ML5_0_UPDATE_WL_2_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
524
534
     $ ,OUT)
525
535
      INCLUDE 'coef_specs.inc'
 
536
      INCLUDE 'loop_max_coefs.inc'
526
537
      INTEGER I,J,K
527
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
538
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
528
539
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
529
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
540
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
530
541
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
531
542
 
532
543
      DO I=1,LCUT_SIZE
598
609
      SUBROUTINE MP_ML5_0_UPDATE_WL_2_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
599
610
     $ ,OUT)
600
611
      INCLUDE 'coef_specs.inc'
 
612
      INCLUDE 'loop_max_coefs.inc'
601
613
      INTEGER I,J,K
602
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
614
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
603
615
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
604
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
616
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
605
617
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
606
618
 
607
619
      DO I=1,LCUT_SIZE
673
685
      SUBROUTINE ML5_0_UPDATE_WL_2_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
674
686
     $ ,OUT)
675
687
      INCLUDE 'coef_specs.inc'
 
688
      INCLUDE 'loop_max_coefs.inc'
676
689
      INTEGER I,J,K
677
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
690
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
678
691
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
679
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
692
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
680
693
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
681
694
 
682
695
      DO I=1,LCUT_SIZE
708
721
      SUBROUTINE MP_ML5_0_UPDATE_WL_2_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
709
722
     $ ,OUT)
710
723
      INCLUDE 'coef_specs.inc'
 
724
      INCLUDE 'loop_max_coefs.inc'
711
725
      INTEGER I,J,K
712
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
726
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
713
727
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
714
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
728
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
715
729
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
716
730
 
717
731
      DO I=1,LCUT_SIZE
743
757
      SUBROUTINE ML5_0_UPDATE_WL_1_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
744
758
     $ ,OUT)
745
759
      INCLUDE 'coef_specs.inc'
 
760
      INCLUDE 'loop_max_coefs.inc'
746
761
      INTEGER I,J,K
747
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
762
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
748
763
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
749
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
764
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
750
765
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
751
766
 
752
767
      DO I=1,LCUT_SIZE
768
783
      SUBROUTINE MP_ML5_0_UPDATE_WL_1_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
769
784
     $ ,OUT)
770
785
      INCLUDE 'coef_specs.inc'
 
786
      INCLUDE 'loop_max_coefs.inc'
771
787
      INTEGER I,J,K
772
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
788
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
773
789
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
774
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
790
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
775
791
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
776
792
 
777
793
      DO I=1,LCUT_SIZE
793
809
      SUBROUTINE ML5_0_UPDATE_WL_1_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
794
810
     $ ,OUT)
795
811
      INCLUDE 'coef_specs.inc'
 
812
      INCLUDE 'loop_max_coefs.inc'
796
813
      INTEGER I,J,K
797
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
814
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
798
815
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
799
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
816
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
800
817
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
801
818
 
802
819
      DO I=1,LCUT_SIZE
828
845
      SUBROUTINE MP_ML5_0_UPDATE_WL_1_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
829
846
     $ ,OUT)
830
847
      INCLUDE 'coef_specs.inc'
 
848
      INCLUDE 'loop_max_coefs.inc'
831
849
      INTEGER I,J,K
832
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
850
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
833
851
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
834
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOP_MAXCOEFS-1,MAXLWFSIZE)
 
852
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
835
853
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
836
854
 
837
855
      DO I=1,LCUT_SIZE