* * $Id: rsinv.F,v 1.1.1.1 1996/02/15 17:49:05 mclareni Exp $ * * $Log: rsinv.F,v $ * Revision 1.1.1.1 1996/02/15 17:49:05 mclareni * Kernlib * * #include "kernnum/pilot.h" SUBROUTINE RSINV(N,A,IDIM,IFAIL) REAL A(IDIM,*), ZERO, ONE, X, Y REAL PIVOTF CHARACTER*6 HNAME #if !defined(CERNLIB_NUMRDBLE) REAL S1, S31, S32, S33, DOTF DOTF(X,Y,S1) = X*Y + S1 #endif #if defined(CERNLIB_NUMRDBLE) DOUBLE PRECISION S1, S31, S32, S33, DOTF DOTF(X,Y,S1) = DBLE(X)*DBLE(Y) + S1 #endif PIVOTF(X) = X DATA HNAME / 'RSINV ' / DATA ZERO, ONE / 0., 1. / IF(IDIM .LT. N .OR. N .LE. 0) GOTO 900 #include "sfact.inc" #include "sfinv.inc" RETURN 900 CALL TMPRNT(HNAME,N,IDIM,0) RETURN END