MINICERN/packlib/kernlib/kernnum/f011fort/rfact.F

   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          RFACT(N,A,IDIM,IR,IFAIL,DET,JFAIL)
  11           INTEGER             IR(*),    IPAIRF
  12           REAL                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_NUMRDBLE)
  17           REAL                S11,      S12,      DOTF
  18           DOTF(X,Y,S11)  =  X*Y + S11
  19 #endif
  20 #if defined(CERNLIB_NUMRDBLE)
  21           DOUBLE PRECISION    S11,      S12,      DOTF
  22           DOTF(X,Y,S11)  =  DBLE(X)*DBLE(Y) + S11
  23 #endif
  24           IPAIRF(J,K)  =  J*2**12 + K
  25           PIVOTF(X)    =  ABS(X)
  26           SIZEF(X)     =  ABS(X)
  27 #if defined(CERNLIB_NUME2465)
  28           DATA      G1, G2              /  1.E-1232, 1.E1232  /
  29 #endif
  30 #if defined(CERNLIB_NUME293)
  31           DATA      G1, G2              /  1.E-146, 1.E146  /
  32 #endif
  33 #if defined(CERNLIB_NUME75)
  34           DATA      G1, G2              /  1.E-37,  1.E37  /
  35 #endif
  36 #if defined(CERNLIB_NUME38)
  37           DATA      G1, G2              /  1.E-19,  1.E19  /
  38 #endif
  39 #if defined(CERNLIB_NUME999)
  40           DATA      ?????  G1, G2 NOT DEFINED  ?????
  41 #endif
  42           DATA      HNAME               /  ' RFACT'  /
  43           DATA      ZERO, ONE           /  0., 1.  /
  44           DATA      NORMAL, IMPOSS      /  0, -1  /
  45           DATA      JRANGE, JOVER, JUNDER  /  0, +1, -1  /
  46 #include "fact.inc"
  47           RETURN
  48           END