MadWeightAnalyzer: main_code.f

File main_code.f, 4.0 KB (added by anonymous, 7 years ago)
Line 
1C+-----------------------------------------------------------------------+
2C|                         MAIN CODE FOR MADWEIGHT                       |
3C|                                                                       |
4C|     Author: Pierre Artoisenet (UCL-CP3)                               |
5C|             Olivier Mattelaer (UCL-CP3)                               |
6C+-----------------------------------------------------------------------+
7C|     This file is generated automaticly by MADWEIGHT-ANALYZER          |
8C+-----------------------------------------------------------------------+     
9      subroutine main_code(x,n_var,num_sol)
10C+-----------------------------------------------------------------------+
11C|     Central Routine for the change of variable choice                 |
12C|       - num_sol: number of the solution to charge                     |
13C|       - x      : random number from vegas                             |
14C+-----------------------------------------------------------------------+
15
16        integer num_sol
17        double precision x(20)
18        integer n_var
19        double precision jac_loc
20
21C+-----------------------------------------------------------------------+
22C|     Scedullar Part                                                    |
23C+-----------------------------------------------------------------------+   
24       if (num_sol.eq.1) then
25
26C+-----------------------------------------------------------------------+
27C|                                                                       |
28C|        ** Enlarged Contraint Sector global information **             |
29C|                                                                       |
30C|    Class: E                                                           |
31C|    particle in ECS : 4(missing)    6(missing)    3(visible)           |
32C|                      5(visible)                                       |
33C|    blob linked are generated by :                                     |
34C|                                                                       |
35C|                                                                       |
36C+-----------------------------------------------------------------------+
37
38
39
40
41C+-----------------------------------------------------------------------+
42C|                                                                       |
43C|     ENLARGED CONTRAINT SECTOR      CLASS E                            |
44C|                                                                       |
45C+-----------------------------------------------------------------------+
46
47
48        call class_e(x,n_var,jac_loc,4,6,3,5,-1,-2,-3)
49        call block_stat(1,'e-4')
50        if (jac_loc.lt.0d0) return
51
52       elseif (num_sol.eq.2) then
53
54C+-----------------------------------------------------------------------+
55C|                                                                       |
56C|        ** Enlarged Contraint Sector global information **             |
57C|                                                                       |
58C|    Class: E                                                           |
59C|    particle in ECS : 4(missing)    6(missing)    3(visible)           |
60C|                      5(visible)                                       |
61C|    blob linked are generated by :                                     |
62C|                                                                       |
63C|                                                                       |
64C+-----------------------------------------------------------------------+
65
66
67
68
69C+-----------------------------------------------------------------------+
70C|                                                                       |
71C|     ENLARGED CONTRAINT SECTOR      CLASS E                            |
72C|                                                                       |
73C+-----------------------------------------------------------------------+
74
75
76        call class_e(x,n_var,jac_loc,4,6,3,5,-1,-2,-3)
77        call block_stat(1,'e-4')
78        if (jac_loc.lt.0d0) return
79
80        endif