CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Functions | Variables
sicif.h File Reference
#include "vdt/vdtMath.h"

Go to the source code of this file.

Functions

float p1evlf (float xx, float *coef, int N)
 
float polevlf (float xx, float *coef, int N)
 
int sicif (float xx, float &si, float &ci)
 

Variables

static float CD []
 
static float CN []
 
static float FD4 []
 
static float FD8 []
 
static float FN4 []
 
static float FN8 []
 
static float GD4 []
 
static float GD8 []
 
static float GN4 []
 
static float GN8 []
 
static float SD []
 
static float SN []
 

Function Documentation

float p1evlf ( float  xx,
float *  coef,
int  N 
)
inline

Definition at line 207 of file sicif.h.

References i, AlCaHLTBitMon_ParallelJobs::p, and x.

Referenced by sicif().

207  {
208 float ans, x;
209 float *p;
210 int i;
211 
212 x = xx;
213 p = coef;
214 ans = x + *p++;
215 i = N-1;
216 
217 do
218  ans = ans * x + *p++;
219 while( --i );
220 
221 return( ans );
222 }
i
int i
Definition: DBlmapReader.cc:9
N
#define N
Definition: blowfish.cc:9
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
x
Definition: DDAxes.h:10
float polevlf ( float  xx,
float *  coef,
int  N 
)
inline

Definition at line 184 of file sicif.h.

References i, N, AlCaHLTBitMon_ParallelJobs::p, and x.

Referenced by sicif().

184  {
185 float ans, x;
186 float *p;
187 int i;
188 
189 x = xx;
190 p = coef;
191 ans = *p++;
192 
193 i = N;
194 do
195  ans = ans * x + *p++;
196 while( --i );
197 
198 return( ans );
199 }
i
int i
Definition: DBlmapReader.cc:9
N
#define N
Definition: blowfish.cc:9
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
x
Definition: DDAxes.h:10
int sicif ( float  xx,
float &  si,
float &  ci 
)
inline

Definition at line 227 of file sicif.h.

References trackerHits::c, CD, CN, f, FD4, FD8, FN4, FN8, g, GD4, GD8, GN4, GN8, p1evlf(), polevlf(), alignCSCRings::s, SD, SN, x, and detailsBasic3DVector::z.

Referenced by VVIObjFDetails::sincosint().

227  {
228  const float MAXNUMF = 1.7014117331926442990585209174225846272e38;
229  const float PIO2F = 1.5707963267948966192;
230  // const float MACHEPF = 5.9604644775390625E-8;
231  const float EUL = 0.57721566490153286061;
232 
233  float x, z, c, s, f, g;
234  int sign;
235 
236  x = xx;
237  if( x < 0.0f )
238  {
239  sign = -1;
240  x = -x;
241  }
242  else
243  sign = 0;
244 
245 
246  if( x == 0.0f )
247  {
248  si = 0.0;
249  ci = -MAXNUMF;
250  return( 0 );
251  }
252 
253 
254  if( x > 1.0e9f )
255  {
256  float su,cu; vdt::fast_sincosf(x,su,cu);
257  si = PIO2F - cu/x;
258  ci = su/x;
259  return( 0 );
260  }
261 
262 
263 
264  if( x > 4.0f )
265  goto asympt;
266 
267  z = x * x;
268  s = x * polevlf( z, SN, 5 ) / polevlf( z, SD, 5 );
269  c = z * polevlf( z, CN, 5 ) / polevlf( z, CD, 5 );
270 
271  if( sign )
272  s = -s;
273  si = s;
274  ci = EUL + vdt::fast_logf(x) + c; /* real part if x < 0 */
275  return(0);
276 
277 
278 
279  /* The auxiliary functions are:
280  *
281  *
282  * *si = *si - PIO2;
283  * c = cos(x);
284  * s = sin(x);
285  *
286  * t = *ci * s - *si * c;
287  * a = *ci * c + *si * s;
288  *
289  * *si = t;
290  * *ci = -a;
291  */
292 
293 
294  asympt:
295  vdt::fast_sincosf(x,s,c);
296  z = 1.0f/(x*x);
297  if( x < 8.0f )
298  {
299  f = polevlf( z, FN4, 6 ) / (x * p1evlf( z, FD4, 7 ));
300  g = z * polevlf( z, GN4, 7 ) / p1evlf( z, GD4, 7 );
301  }
302  else
303  {
304  f = polevlf( z, FN8, 8 ) / (x * p1evlf( z, FD8, 8 ));
305  g = z * polevlf( z, GN8, 8 ) / p1evlf( z, GD8, 9 );
306  }
307  si = PIO2F - f * c - g * s;
308  if( sign )
309  si = -( si );
310  ci = f * s - g * c;
311 
312  return(0);
313 }
FN4
static float FN4[]
Definition: sicif.h:94
SN
static float SN[]
Definition: sicif.h:59
CN
static float CN[]
Definition: sicif.h:76
FN8
static float FN8[]
Definition: sicif.h:115
g
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g
Definition: Activities.doc:4
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15
GD4
static float GD4[]
Definition: sicif.h:148
FD8
static float FD8[]
Definition: sicif.h:126
FD4
static float FD4[]
Definition: sicif.h:103
f
double f[11][100]
Definition: MuScleFitUtils.cc:79
polevlf
float polevlf(float xx, float *coef, int N)
Definition: sicif.h:184
GN8
static float GN8[]
Definition: sicif.h:159
p1evlf
float p1evlf(float xx, float *coef, int N)
Definition: sicif.h:207
GD8
static float GD8[]
Definition: sicif.h:170
SD
static float SD[]
Definition: sicif.h:67
GN4
static float GN4[]
Definition: sicif.h:138
x
Definition: DDAxes.h:10
CD
static float CD[]
Definition: sicif.h:84

Variable Documentation

float CD[]
static
Initial value:
= {
4.07746040061880559506E-12,
3.06780997581887812692E-9,
1.23210355685883423679E-6,
3.17442024775032769882E-4,
5.10028056236446052392E-2,
4.00000000000000000080E0,
}

Definition at line 84 of file sicif.h.

Referenced by sicif().

float CN[]
static
Initial value:
= {
2.02524002389102268789E-11,
-1.35249504915790756375E-8,
3.59325051419993077021E-6,
-4.74007206873407909465E-4,
2.89159652607555242092E-2,
-1.00000000000000000080E0,
}

Definition at line 76 of file sicif.h.

Referenced by EcalLaserAnalyzerYousi::analyze(), RecoTauCleanerImpl< Prod >::produce(), and sicif().

float FD4[]
static
Initial value:
= {
8.16496634205391016773E0,
7.30828822505564552187E0,
1.86792257950184183883E0,
1.78792052963149907262E-1,
7.01710668322789753610E-3,
1.10034357153915731354E-4,
5.48900252756255700982E-7,
}

Definition at line 103 of file sicif.h.

Referenced by sicif().

float FD8[]
static
Initial value:
= {
9.17463611873684053703E-1,
1.78685545332074536321E-1,
1.22253594771971293032E-2,
3.58696481881851580297E-4,
4.92435064317881464393E-6,
3.21956939101046018377E-8,
9.43720590350276732376E-11,
9.70507110881952025725E-14,
}

Definition at line 126 of file sicif.h.

Referenced by sicif().

float FN4[]
static
Initial value:
= {
4.23612862892216586994E0,
5.45937717161812843388E0,
1.62083287701538329132E0,
1.67006611831323023771E-1,
6.81020132472518137426E-3,
1.08936580650328664411E-4,
5.48900223421373614008E-7,
}

Definition at line 94 of file sicif.h.

Referenced by sicif().

float FN8[]
static
Initial value:
= {
4.55880873470465315206E-1,
7.13715274100146711374E-1,
1.60300158222319456320E-1,
1.16064229408124407915E-2,
3.49556442447859055605E-4,
4.86215430826454749482E-6,
3.20092790091004902806E-8,
9.41779576128512936592E-11,
9.70507110881952024631E-14,
}

Definition at line 115 of file sicif.h.

Referenced by sicif().

float GD4[]
static
Initial value:
= {
1.64402202413355338886E0,
6.66296701268987968381E-1,
9.88771761277688796203E-2,
6.22396345441768420760E-3,
1.73221081474177119497E-4,
2.02659182086343991969E-6,
7.82579218933534490868E-9,
}

Definition at line 148 of file sicif.h.

Referenced by sicif().

float GD8[]
static
Initial value:
= {
1.68548898811011640017E0,
4.87852258695304967486E-1,
4.67913194259625806320E-2,
1.90284426674399523638E-3,
3.68475504442561108162E-5,
3.57043223443740838771E-7,
1.72693748966316146736E-9,
3.87830166023954706752E-12,
3.14040098946363335242E-15,
}

Definition at line 170 of file sicif.h.

Referenced by sicif().

float GN4[]
static
Initial value:
= {
8.71001698973114191777E-2,
6.11379109952219284151E-1,
3.97180296392337498885E-1,
7.48527737628469092119E-2,
5.38868681462177273157E-3,
1.61999794598934024525E-4,
1.97963874140963632189E-6,
7.82579040744090311069E-9,
}

Definition at line 138 of file sicif.h.

Referenced by sicif().

float GN8[]
static
Initial value:
= {
6.97359953443276214934E-1,
3.30410979305632063225E-1,
3.84878767649974295920E-2,
1.71718239052347903558E-3,
3.48941165502279436777E-5,
3.47131167084116673800E-7,
1.70404452782044526189E-9,
3.85945925430276600453E-12,
3.14040098946363334640E-15,
}

Definition at line 159 of file sicif.h.

Referenced by sicif().

float SD[]
static
Initial value:
= {
2.03269266195951942049E-12,
1.27997891179943299903E-9,
4.41827842801218905784E-7,
9.96412122043875552487E-5,
1.42085239326149893930E-2,
9.99999999999999996984E-1,
}

Definition at line 67 of file sicif.h.

Referenced by HcalDetDiagLEDMonitor::GetCalib(), HcalDetDiagLaserMonitor::GetCalib(), HcalDetDiagLEDClient::htmlOutput(), and sicif().

float SN[]
static
Initial value:
= {
-8.39167827910303881427E-11,
4.62591714427012837309E-8,
-9.75759303843632795789E-6,
9.76945438170435310816E-4,
-4.13470316229406538752E-2,
1.00000000000000000302E0,
}

Definition at line 59 of file sicif.h.

Referenced by sicif().

Generated for CMSSW Reference Manual by doxygen 1.8.5