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

*
* $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
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/04/01 15:01:46 mclareni
	6	* Mathlib gen
	7	*
	8	*
	9	#include "gen/pilot.h"
	10	SUBROUTINE U100PR (I,K,M)
	11	DIMENSION IPRIME(40),K(40)
	12	DATA IPRIME/2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,
	13	+ 71,73,79,83,89,97,101,103,107,109,113,127,131,137,139,149,151,
	14	+ 157,163,167,173/
	15	II=I
	16	M=1
	17	IF(I.LT.1) GO TO 6
	18	CALL VZERO (K, 40)
	19	IF(I.EQ.1) RETURN
	20	DO 2 KL=1,40
	21	1 IF(MOD(I,IPRIME(KL)).NE.0) GO TO 2
	22	K(KL)=K(KL)+1
	23	I=I/IPRIME(KL)
	24	IF(I.EQ.1) GO TO 3
	25	GO TO 1
	26	2 CONTINUE
	27	IF(I.NE.1) WRITE(6,4) I
	28	4 FORMAT(' PRIME NO. LARGER THAN THE LIST,THE FACTOR LEFT=',I15)
	29	M=I
	30	3 I=II
	31	RETURN
	32	6 WRITE(6,7)
	33	7 FORMAT(' PLEASE OBSERVE THE CONVENTION OF YOUR ASKING NO. BEING',
	34	F'LARGER THAN +0')
	35	RETURN
	36	END