~maddevelopers/mg5amcnlo/simple_unlops : revision 969

tests/input_files/IOTestsComparison/short_ML_SMQCD_default/gg_ttx/CT_interface.f

152

153

END

154

155

SUBROUTINE ML5_0_LOOP_3_4_3( LID, P1, P2, P3, W1, W2, W3, W4, M1

156

$ ,MP_M1, M2,MP_M2, M3,MP_M3, C1,MP_C1, C2,MP_C2, C3,MP_C3, RANK

157

$ , LSYMFACT, LMULTIPLIER, AMPLN, RES, STABLE)

158

159

INTEGER NEXTERNAL

160

PARAMETER (NEXTERNAL=4)

161

INTEGER MAXLCOUPLINGS

162

PARAMETER (MAXLCOUPLINGS=4)

163

INTEGER NLOOPLINE

164

PARAMETER (NLOOPLINE=3)

165

INTEGER NWAVEFUNCS

166

PARAMETER (NWAVEFUNCS=10)

167

INTEGER NCOMB

168

PARAMETER (NCOMB=16)

169

170

C ARGUMENTS

171

172

INTEGER W1, W2, W3, W4

173

COMPLEX*16 M1, M2, M3

174

COMPLEX*32 MP_M1, MP_M2, MP_M3

175

COMPLEX*16 C1, C2, C3

176

COMPLEX*32 MP_C1, MP_C2, MP_C3

177

INTEGER P1, P2, P3

178

COMPLEX*16 RES(3)

179

INTEGER LID, RANK, LSYMFACT, LMULTIPLIER

180

INTEGER AMPLN

181

LOGICAL STABLE

182

183

C LOCAL VARIABLES

184

185

REAL*8 PL(0:3,NLOOPLINE)

186

COMPLEX*16 M2L(NLOOPLINE)

187

INTEGER PAIRING(NLOOPLINE)

188

INTEGER I, J, K, TEMP

189

190

C GLOBAL VARIABLES

191

192

INTEGER WE(NEXTERNAL)

193

INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM

194

COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM

195

196

COMPLEX*16 LC(MAXLCOUPLINGS)

197

COMPLEX*16 ML(NEXTERNAL+2)

198

COMMON/ML5_0_DP_LOOP/LC,ML

199

200

COMPLEX*32 MP_LC(MAXLCOUPLINGS)

201

COMPLEX*32 MP_ML(NEXTERNAL+2)

202

COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML

203

204

COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)

205

INTEGER VALIDH

206

COMMON/ML5_0_WFCTS/W

207

COMMON/ML5_0_VALIDH/VALIDH

208

209

C ----------

210

C BEGIN CODE

211

C ----------

212

213

WE(1)=W1

214

WE(2)=W2

215

WE(3)=W3

216

WE(4)=W4

217

M2L(1)=M3**2

218

M2L(2)=M1**2

219

M2L(3)=M2**2

220

ML(1)=M3

221

ML(2)=M3

222

MP_ML(1)=MP_M3

223

MP_ML(2)=MP_M3

224

ML(3)=M1

225

MP_ML(3)=MP_M1

226

ML(4)=M2

227

MP_ML(4)=MP_M2

228

ML(5)=M3

229

MP_ML(5)=MP_M3

230

PAIRING(1)=P1

231

PAIRING(2)=P2

232

PAIRING(3)=P3

233

LC(1)=C1

234

MP_LC(1)=MP_C1

235

LC(2)=C2

236

MP_LC(2)=MP_C2

237

LC(3)=C3

238

MP_LC(3)=MP_C3

239

AMPLNUM=AMPLN

240

ID=LID

241

SYMFACT=LSYMFACT

242

MULTIPLIER=LMULTIPLIER

243

DO I=0,3

244

TEMP=1

245

DO J=1,NLOOPLINE

246

PL(I,J)=0.D0

247

DO K=TEMP,(TEMP+PAIRING(J)-1)

248

PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))

249

ENDDO

250

TEMP=TEMP+PAIRING(J)

251

ENDDO

252

ENDDO

253

CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)

254

255

END

256

257

155

SUBROUTINE ML5_0_LOOP_2_2( LID, W1, W2, M1,MP_M1, M2,MP_M2, C1

258

156

$ ,MP_C1, C2,MP_C2, RANK, LSYMFACT, LMULTIPLIER, AMPLN, RES,

259

157

$ STABLE)

350

248

351

249

END

352

250

353

SUBROUTINE ML5_0_LOOP_2_3_2( LID, P1, P2, W1, W2, W3, M1,MP_M1,

354

$ M2,MP_M2, C1,MP_C1, C2,MP_C2, RANK, LSYMFACT, LMULTIPLIER,

355

$ AMPLN, RES, STABLE)

356

357

INTEGER NEXTERNAL

358

PARAMETER (NEXTERNAL=4)

359

INTEGER MAXLCOUPLINGS

360

PARAMETER (MAXLCOUPLINGS=4)

361

INTEGER NLOOPLINE

362

PARAMETER (NLOOPLINE=2)

363

INTEGER NWAVEFUNCS

364

PARAMETER (NWAVEFUNCS=10)

365

INTEGER NCOMB

366

PARAMETER (NCOMB=16)

367

368

C ARGUMENTS

369

370

INTEGER W1, W2, W3

371

COMPLEX*16 M1, M2

372

COMPLEX*32 MP_M1, MP_M2

373

COMPLEX*16 C1, C2

374

COMPLEX*32 MP_C1, MP_C2

375

INTEGER P1, P2

376

COMPLEX*16 RES(3)

377

INTEGER LID, RANK, LSYMFACT, LMULTIPLIER

378

INTEGER AMPLN

379

LOGICAL STABLE

380

381

C LOCAL VARIABLES

382

383

REAL*8 PL(0:3,NLOOPLINE)

384

COMPLEX*16 M2L(NLOOPLINE)

385

INTEGER PAIRING(NLOOPLINE)

386

INTEGER I, J, K, TEMP

387

388

C GLOBAL VARIABLES

389

390

INTEGER WE(NEXTERNAL)

391

INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM

392

COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM

393

394

COMPLEX*16 LC(MAXLCOUPLINGS)

395

COMPLEX*16 ML(NEXTERNAL+2)

396

COMMON/ML5_0_DP_LOOP/LC,ML

397

398

COMPLEX*32 MP_LC(MAXLCOUPLINGS)

399

COMPLEX*32 MP_ML(NEXTERNAL+2)

400

COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML

401

402

COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)

403

INTEGER VALIDH

404

COMMON/ML5_0_WFCTS/W

405

COMMON/ML5_0_VALIDH/VALIDH

406

407

C ----------

408

C BEGIN CODE

409

C ----------

410

411

WE(1)=W1

412

WE(2)=W2

413

WE(3)=W3

414

M2L(1)=M2**2

415

M2L(2)=M1**2

416

ML(1)=M2

417

ML(2)=M2

418

MP_ML(1)=MP_M2

419

MP_ML(2)=MP_M2

420

ML(3)=M1

421

MP_ML(3)=MP_M1

422

ML(4)=M2

423

MP_ML(4)=MP_M2

424

PAIRING(1)=P1

425

PAIRING(2)=P2

426

LC(1)=C1

427

MP_LC(1)=MP_C1

428

LC(2)=C2

429

MP_LC(2)=MP_C2

430

AMPLNUM=AMPLN

431

ID=LID

432

SYMFACT=LSYMFACT

433

MULTIPLIER=LMULTIPLIER

434

DO I=0,3

435

TEMP=1

436

DO J=1,NLOOPLINE

437

PL(I,J)=0.D0

438

DO K=TEMP,(TEMP+PAIRING(J)-1)

439

PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))

440

ENDDO

441

TEMP=TEMP+PAIRING(J)

442

ENDDO

443

ENDDO

444

CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)

445

446

END

447

448

251

SUBROUTINE ML5_0_LOOP_3_3( LID, W1, W2, W3, M1,MP_M1, M2,MP_M2,

449

252

$ M3,MP_M3, C1,MP_C1, C2,MP_C2, C3,MP_C3, RANK, LSYMFACT,

450

253

$ LMULTIPLIER, AMPLN, RES, STABLE)

655

458

656

459

END

657

460

461

SUBROUTINE ML5_0_LOOP_2_3_2( LID, P1, P2, W1, W2, W3, M1,MP_M1,

462

$ M2,MP_M2, C1,MP_C1, C2,MP_C2, RANK, LSYMFACT, LMULTIPLIER,

463

$ AMPLN, RES, STABLE)

464

465

INTEGER NEXTERNAL

466

PARAMETER (NEXTERNAL=4)

467

INTEGER MAXLCOUPLINGS

468

PARAMETER (MAXLCOUPLINGS=4)

469

INTEGER NLOOPLINE

470

PARAMETER (NLOOPLINE=2)

471

INTEGER NWAVEFUNCS

472

PARAMETER (NWAVEFUNCS=10)

473

INTEGER NCOMB

474

PARAMETER (NCOMB=16)

475

476

C ARGUMENTS

477

478

INTEGER W1, W2, W3

479

COMPLEX*16 M1, M2

480

COMPLEX*32 MP_M1, MP_M2

481

COMPLEX*16 C1, C2

482

COMPLEX*32 MP_C1, MP_C2

483

INTEGER P1, P2

484

COMPLEX*16 RES(3)

485

INTEGER LID, RANK, LSYMFACT, LMULTIPLIER

486

INTEGER AMPLN

487

LOGICAL STABLE

488

489

C LOCAL VARIABLES

490

491

REAL*8 PL(0:3,NLOOPLINE)

492

COMPLEX*16 M2L(NLOOPLINE)

493

INTEGER PAIRING(NLOOPLINE)

494

INTEGER I, J, K, TEMP

495

496

C GLOBAL VARIABLES

497

498

INTEGER WE(NEXTERNAL)

499

INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM

500

COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM

501

502

COMPLEX*16 LC(MAXLCOUPLINGS)

503

COMPLEX*16 ML(NEXTERNAL+2)

504

COMMON/ML5_0_DP_LOOP/LC,ML

505

506

COMPLEX*32 MP_LC(MAXLCOUPLINGS)

507

COMPLEX*32 MP_ML(NEXTERNAL+2)

508

COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML

509

510

COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)

511

INTEGER VALIDH

512

COMMON/ML5_0_WFCTS/W

513

COMMON/ML5_0_VALIDH/VALIDH

514

515

C ----------

516

C BEGIN CODE

517

C ----------

518

519

WE(1)=W1

520

WE(2)=W2

521

WE(3)=W3

522

M2L(1)=M2**2

523

M2L(2)=M1**2

524

ML(1)=M2

525

ML(2)=M2

526

MP_ML(1)=MP_M2

527

MP_ML(2)=MP_M2

528

ML(3)=M1

529

MP_ML(3)=MP_M1

530

ML(4)=M2

531

MP_ML(4)=MP_M2

532

PAIRING(1)=P1

533

PAIRING(2)=P2

534

LC(1)=C1

535

MP_LC(1)=MP_C1

536

LC(2)=C2

537

MP_LC(2)=MP_C2

538

AMPLNUM=AMPLN

539

ID=LID

540

SYMFACT=LSYMFACT

541

MULTIPLIER=LMULTIPLIER

542

DO I=0,3

543

TEMP=1

544

DO J=1,NLOOPLINE

545

PL(I,J)=0.D0

546

DO K=TEMP,(TEMP+PAIRING(J)-1)

547

PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))

548

ENDDO

549

TEMP=TEMP+PAIRING(J)

550

ENDDO

551

ENDDO

552

CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)

553

554

END

555

556

SUBROUTINE ML5_0_LOOP_3_4_3( LID, P1, P2, P3, W1, W2, W3, W4, M1

557

$ ,MP_M1, M2,MP_M2, M3,MP_M3, C1,MP_C1, C2,MP_C2, C3,MP_C3, RANK

558

$ , LSYMFACT, LMULTIPLIER, AMPLN, RES, STABLE)

559

560

INTEGER NEXTERNAL

561

PARAMETER (NEXTERNAL=4)

562

INTEGER MAXLCOUPLINGS

563

PARAMETER (MAXLCOUPLINGS=4)

564

INTEGER NLOOPLINE

565

PARAMETER (NLOOPLINE=3)

566

INTEGER NWAVEFUNCS

567

PARAMETER (NWAVEFUNCS=10)

568

INTEGER NCOMB

569

PARAMETER (NCOMB=16)

570

571

C ARGUMENTS

572

573

INTEGER W1, W2, W3, W4

574

COMPLEX*16 M1, M2, M3

575

COMPLEX*32 MP_M1, MP_M2, MP_M3

576

COMPLEX*16 C1, C2, C3

577

COMPLEX*32 MP_C1, MP_C2, MP_C3

578

INTEGER P1, P2, P3

579

COMPLEX*16 RES(3)

580

INTEGER LID, RANK, LSYMFACT, LMULTIPLIER

581

INTEGER AMPLN

582

LOGICAL STABLE

583

584

C LOCAL VARIABLES

585

586

REAL*8 PL(0:3,NLOOPLINE)

587

COMPLEX*16 M2L(NLOOPLINE)

588

INTEGER PAIRING(NLOOPLINE)

589

INTEGER I, J, K, TEMP

590

591

C GLOBAL VARIABLES

592

593

INTEGER WE(NEXTERNAL)

594

INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM

595

COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM

596

597

COMPLEX*16 LC(MAXLCOUPLINGS)

598

COMPLEX*16 ML(NEXTERNAL+2)

599

COMMON/ML5_0_DP_LOOP/LC,ML

600

601

COMPLEX*32 MP_LC(MAXLCOUPLINGS)

602

COMPLEX*32 MP_ML(NEXTERNAL+2)

603

COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML

604

605

COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)

606

INTEGER VALIDH

607

COMMON/ML5_0_WFCTS/W

608

COMMON/ML5_0_VALIDH/VALIDH

609

610

C ----------

611

C BEGIN CODE

612

C ----------

613

614

WE(1)=W1

615

WE(2)=W2

616

WE(3)=W3

617

WE(4)=W4

618

M2L(1)=M3**2

619

M2L(2)=M1**2

620

M2L(3)=M2**2

621

ML(1)=M3

622

ML(2)=M3

623

MP_ML(1)=MP_M3

624

MP_ML(2)=MP_M3

625

ML(3)=M1

626

MP_ML(3)=MP_M1

627

ML(4)=M2

628

MP_ML(4)=MP_M2

629

ML(5)=M3

630

MP_ML(5)=MP_M3

631

PAIRING(1)=P1

632

PAIRING(2)=P2

633

PAIRING(3)=P3

634

LC(1)=C1

635

MP_LC(1)=MP_C1

636

LC(2)=C2

637

MP_LC(2)=MP_C2

638

LC(3)=C3

639

MP_LC(3)=MP_C3

640

AMPLNUM=AMPLN

641

ID=LID

642

SYMFACT=LSYMFACT

643

MULTIPLIER=LMULTIPLIER

644

DO I=0,3

645

TEMP=1

646

DO J=1,NLOOPLINE

647

PL(I,J)=0.D0

648

DO K=TEMP,(TEMP+PAIRING(J)-1)

649

PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))

650

ENDDO

651

TEMP=TEMP+PAIRING(J)

652

ENDDO

653

ENDDO

654

CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)

655

656

END

657

Older »