[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / v / rnmnml.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/04/01 15:02:56  mclareni
* Mathlib gen
*
*
#include "gen/pilot.h"
      SUBROUTINE RNMNML(N,NSUM,PCUM,NVEC,IERR)
 
      PARAMETER (NRV = 100)
      DIMENSION PCUM(*),NVEC(*),RV(NRV)
 
      IERR=2
      IF(PCUM(N) .NE. 1) GO TO 9
      IERR=1
      DO 1 I = 2,N
      IF(PCUM(I) .LT. PCUM(I-1)) GO TO 9
    1 CONTINUE
      IERR=0
      DO 2 K = 1,N
    2 NVEC(K)=0
      J0=NSUM/NRV+1
      I0=NRV
      DO 3 J = 1,J0
      IF(J .EQ. J0) I0=MOD(NSUM,NRV)
      CALL RANLUX(RV,I0)
      DO 4 I = 1,I0
      DO 5 K = 1,N
      IF(PCUM(K) .GE. RV(I)) THEN
       NVEC(K)=NVEC(K)+1
       GO TO 4
      ENDIF
    5 CONTINUE
    4 CONTINUE
    3 CONTINUE
    9 RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/04/01 15:02:56 mclareni
	6	* Mathlib gen
	7	*
	8	*
	9	#include "gen/pilot.h"
	10	SUBROUTINE RNMNML(N,NSUM,PCUM,NVEC,IERR)
	11
	12	PARAMETER (NRV = 100)
	13	DIMENSION PCUM(),NVEC(),RV(NRV)
	14
	15	IERR=2
	16	IF(PCUM(N) .NE. 1) GO TO 9
	17	IERR=1
	18	DO 1 I = 2,N
	19	IF(PCUM(I) .LT. PCUM(I-1)) GO TO 9
	20	1 CONTINUE
	21	IERR=0
	22	DO 2 K = 1,N
	23	2 NVEC(K)=0
	24	J0=NSUM/NRV+1
	25	I0=NRV
	26	DO 3 J = 1,J0
	27	IF(J .EQ. J0) I0=MOD(NSUM,NRV)
	28	CALL RANLUX(RV,I0)
	29	DO 4 I = 1,I0
	30	DO 5 K = 1,N
	31	IF(PCUM(K) .GE. RV(I)) THEN
	32	NVEC(K)=NVEC(K)+1
	33	GO TO 4
	34	ENDIF
	35	5 CONTINUE
	36	4 CONTINUE
	37	3 CONTINUE
	38	9 RETURN
	39	END