]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1996/02/15 17:49:03 mclareni | |
6 | * Kernlib | |
7 | * | |
8 | * | |
9 | #include "kernnum/pilot.h" | |
10 | SUBROUTINE CFACT(N,A,IDIM,IR,IFAIL,DET,JFAIL) | |
11 | INTEGER IR(*), IPAIRF | |
12 | COMPLEX A(IDIM,*),DET, ZERO, ONE,X,Y,TF | |
13 | REAL G1, G2 | |
14 | REAL PIVOTF, P, Q, SIZEF, T | |
15 | CHARACTER*6 HNAME | |
16 | #if !defined(CERNLIB_NUMCDBLE) | |
17 | COMPLEX S11, S12, DOTF | |
18 | DOTF(X,Y,S11) = X*Y + S11 | |
19 | #endif | |
20 | #if defined(CERNLIB_NUMCDBLE) | |
21 | COMPLEX*16 S11, S12, DC, DOTF | |
22 | DC(X) = DCMPLX(DBLE(REAL(X)),DBLE(AIMAG(X))) | |
23 | DOTF(X,Y,S11) = DC(X)*DC(Y) + S11 | |
24 | #endif | |
25 | IPAIRF(J,K) = J*2**12 + K | |
26 | PIVOTF(X) = MAX(ABS(REAL(X)),ABS(AIMAG(X))) | |
27 | SIZEF(X) = MAX(ABS(REAL(X)),ABS(AIMAG(X))) | |
28 | #if defined(CERNLIB_NUME2465) | |
29 | DATA G1, G2 / 1.E-1232, 1.E1232 / | |
30 | #endif | |
31 | #if defined(CERNLIB_NUME293) | |
32 | DATA G1, G2 / 1.E-146, 1.E146 / | |
33 | #endif | |
34 | #if defined(CERNLIB_NUME75) | |
35 | DATA G1, G2 / 1.E-37, 1.E37 / | |
36 | #endif | |
37 | #if defined(CERNLIB_NUME38) | |
38 | DATA G1, G2 / 1.E-19, 1.E19 / | |
39 | #endif | |
40 | #if defined(CERNLIB_NUME999) | |
41 | DATA ????? G1, G2 NOT DEFINED ????? | |
42 | #endif | |
43 | DATA HNAME / ' CFACT' / | |
44 | DATA ZERO, ONE / (0., 0.), (1., 0.) / | |
45 | DATA NORMAL, IMPOSS / 0, -1 / | |
46 | DATA JRANGE, JOVER, JUNDER / 0, +1, -1 / | |
47 | #include "fact.inc" | |
48 | RETURN | |
49 | END |