Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros
G4NeutronHPFastLegendre.hh
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27// $Id$
28//
29#ifndef G4NeutronHPFastLegendre_h
30#define G4NeutronHPFastLegendre_h 1
31
32#include "globals.hh"
33
34class G4NeutronHPFastLegendre
35{
36 public:
37
38 G4NeutronHPFastLegendre()
39 {
40 value = new G4double * [31];
41 value[0] = l0;
42 value[1] = l1;
43 value[2] = l2;
44 value[3] = l3;
45 value[4] = l4;
46 value[5] = l5;
47 value[6] = l6;
48 value[7] = l7;
49 value[8] = l8;
50 value[9] = l9;
51 value[10] = l10;
52 value[11] = l11;
53 value[12] = l12;
54 value[13] = l13;
55 value[14] = l14;
56 value[15] = l15;
57 value[16] = l16;
58 value[17] = l17;
59 value[18] = l18;
60 value[19] = l19;
61 value[20] = l20;
62 value[21] = l21;
63 value[22] = l22;
64 value[23] = l23;
65 value[24] = l24;
66 value[25] = l25;
67 value[26] = l26;
68 value[27] = l27;
69 value[28] = l28;
70 value[29] = l29;
71 value[30] = l30;
72 integral = new G4double * [31];
73 integral[0] = i0;
74 integral[1] = i1;
75 integral[2] = i2;
76 integral[3] = i3;
77 integral[4] = i4;
78 integral[5] = i5;
79 integral[6] = i6;
80 integral[7] = i7;
81 integral[8] = i8;
82 integral[9] = i9;
83 integral[10] = i10;
84 integral[11] = i11;
85 integral[12] = i12;
86 integral[13] = i13;
87 integral[14] = i14;
88 integral[15] = i15;
89 integral[16] = i16;
90 integral[17] = i17;
91 integral[18] = i18;
92 integral[19] = i19;
93 integral[20] = i20;
94 integral[21] = i21;
95 integral[22] = i22;
96 integral[23] = i23;
97 integral[24] = i24;
98 integral[25] = i25;
99 integral[26] = i26;
100 integral[27] = i27;
101 integral[28] = i28;
102 integral[29] = i29;
103 integral[30] = i30;
104
105 G4int i;
106 for(i=0;i<31;i++) theNbin[i]=1+200*(i+1);
107 }
108
109 ~G4NeutronHPFastLegendre()
110 {
111 delete [] value;
112 delete [] integral;
113 }
114
115 G4double Integrate(G4int l, G4double costh)
116 {
117 if(l>30) return regularIntegrate(l,costh);
118 G4int bin = GetBin(l, costh);
119 G4double y1, y2;
120// G4cout <<"Testhpw G4NeutronHPFastLegendre::Integrate "<<l<<" "<<bin<<G4endl;
121 y1 = integral[l][bin];
122 y2 = integral[l][bin+1];
123// G4cout <<"Testhpw G4NeutronHPFastLegendre::Integrate exit"<<G4endl;
124 return Interpolate(bin, l, y1, y2, costh);
125 }
126
127 inline G4double Evaluate(G4int l, G4double costh)
128 {
129 if(l>30) return regularEvaluate(l,costh);
130 G4double result;
131 G4int bin = GetBin(l, costh);
132 if(bin != theNbin[l]-1)
133 {
134 G4double y1, y2;
135 y1 = value[l][bin];
136 y2 = value[l][bin+1];
137 result = Interpolate(bin, l, y1, y2, costh);
138 }
139 else
140 {
141 result = value[l][bin];
142 }
143 return result;
144 }
145
146 private:
147
148 G4double regularEvaluate( int l , double x );
149 G4double regularIntegrate( int l , double x );
150
151
152 inline G4int GetBin(G4int l, G4double costh)
153 {
154 G4int bin=0;
155 bin = G4int( (theNbin[l]-1)*(costh+1)/2. );
156 if(bin == theNbin[l]-1) bin--;
157 return bin;
158 }
159
160 inline G4double Interpolate(G4int bin, G4int l, G4double y1, G4double y2, G4double x)
161 {
162 G4double slope = 0, off = 0, x2=0, x1mx2;
163 G4int half = (theNbin[l]-1)/2;
164// x1 = (bin-half)/G4double(half);
165 x2 = (bin+1-half)/G4double(half);
166 x1mx2 = 1./G4double( (theNbin[l]-1)/2 );
167// slope = (y2-y1)/(x2-x1);
168 slope = (y2-y1)/x1mx2;
169 off = y2-x2*slope;
170 return x*slope+off;
171 }
172
173 G4double ** value;
174 G4double ** integral;
175 G4int theNbin[31];
176 static G4double l0[201];
177 static G4double i0[201];
178 static G4double l1[401];
179 static G4double i1[401];
180 static G4double l2[601];
181 static G4double i2[601];
182 static G4double l3[801];
183 static G4double i3[801];
184 static G4double l4[1001];
185 static G4double i4[1001];
186 static G4double l5[1201];
187 static G4double i5[1201];
188 static G4double l6[1401];
189 static G4double i6[1401];
190 static G4double l7[1601];
191 static G4double i7[1601];
192 static G4double l8[1801];
193 static G4double i8[1801];
194 static G4double l9[2001];
195 static G4double i9[2001];
196 static G4double l10[2201];
197 static G4double i10[2201];
198 static G4double l11[2401];
199 static G4double i11[2401];
200 static G4double l12[2601];
201 static G4double i12[2601];
202 static G4double l13[2801];
203 static G4double i13[2801];
204 static G4double l14[3001];
205 static G4double i14[3001];
206 static G4double l15[3201];
207 static G4double i15[3201];
208 static G4double l16[3401];
209 static G4double i16[3401];
210 static G4double l17[3601];
211 static G4double i17[3601];
212 static G4double l18[3801];
213 static G4double i18[3801];
214 static G4double l19[4001];
215 static G4double i19[4001];
216 static G4double l20[4201];
217 static G4double i20[4201];
218 static G4double l21[4401];
219 static G4double i21[4401];
220 static G4double l22[4601];
221 static G4double i22[4601];
222 static G4double l23[4801];
223 static G4double i23[4801];
224 static G4double l24[5001];
225 static G4double i24[5001];
226 static G4double l25[5201];
227 static G4double i25[5201];
228 static G4double l26[5401];
229 static G4double i26[5401];
230 static G4double l27[5601];
231 static G4double i27[5601];
232 static G4double l28[5801];
233 static G4double i28[5801];
234 static G4double l29[6001];
235 static G4double i29[6001];
236 static G4double l30[6201];
237 static G4double i30[6201];
238};
239#endif
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4NeutronHPFastLegendre::Integrate
G4double Integrate(G4int l, G4double costh)
Definition: G4NeutronHPFastLegendre.hh:115
G4NeutronHPFastLegendre::Evaluate
G4double Evaluate(G4int l, G4double costh)
Definition: G4NeutronHPFastLegendre.hh:127