This commit was generated by cvs2svn to compensate for changes in r2,

[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / u / u100pr.F.ori

diff --git a/MINICERN/mathlib/gen/u/u100pr.F.ori b/MINICERN/mathlib/gen/u/u100pr.F.ori
new file mode 100644 (file)
index 0000000..5200cb9
--- /dev/null
+++ b/MINICERN/mathlib/gen/u/u100pr.F.ori
@@ -0,0 +1,36 @@
+*
+* $Id$
+*
+* $Log$
+* Revision 1.1.1.1  1996/04/01 15:01:46  mclareni
+* Mathlib gen
+*
+*
+#include "gen/pilot.h"
+      SUBROUTINE U100PR (I,K,M)
+      DIMENSION IPRIME(40),K(40)
+      DATA IPRIME/2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,
+     + 71,73,79,83,89,97,101,103,107,109,113,127,131,137,139,149,151,
+     + 157,163,167,173/
+      II=I
+      M=1
+      IF(I.LT.1)  GO TO 6
+      CALL VZERO (K, 40)
+      IF(I.EQ.1)  RETURN
+      DO 2 KL=1,40
+    1 IF(MOD(I,IPRIME(KL)).NE.0)   GO TO 2
+      K(KL)=K(KL)+1
+      I=I/IPRIME(KL)
+      IF(I.EQ.1)  GO TO 3
+      GO TO 1
+    2 CONTINUE
+      IF(I.NE.1)    WRITE(6,4)    I
+    4 FORMAT(' PRIME NO. LARGER THAN THE LIST,THE FACTOR LEFT=',I15)
+      M=I
+    3 I=II
+      RETURN
+    6 WRITE(6,7)
+    7 FORMAT(' PLEASE OBSERVE THE CONVENTION OF YOUR ASKING NO. BEING',
+     F'LARGER THAN +0')
+      RETURN
+      END