[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / f / bakvec.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/04/01 15:02:32  mclareni
* Mathlib gen
*
*
#include "gen/pilot.h"
      SUBROUTINE BAKVEC(NM,N,T,E,M,Z,IERR)
      INTEGER I,J,M,N,NM,IERR
      REAL T(NM,3),E(N),Z(NM,M)
      IERR = 0
      E(1) = 1.0
      IF (N .EQ. 1) GO TO 1001
      DO 100 I = 2, N
         IF (E(I) .NE. 0.0) GO TO 80
         IF (T(I,1) .NE. 0.0 .OR. T(I-1,3) .NE. 0.0) GO TO 1000
         E(I) = 1.0
         GO TO 100
   80    E(I) = E(I-1) * E(I) / T(I-1,3)
  100 CONTINUE
      DO 120 J = 1, M
         DO 120 I = 2, N
         Z(I,J) = Z(I,J) * E(I)
  120 CONTINUE
      GO TO 1001
 1000 IERR = 2 * N + I
 1001 RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/04/01 15:02:32 mclareni
	6	* Mathlib gen
	7	*
	8	*
	9	#include "gen/pilot.h"
	10	SUBROUTINE BAKVEC(NM,N,T,E,M,Z,IERR)
	11	INTEGER I,J,M,N,NM,IERR
	12	REAL T(NM,3),E(N),Z(NM,M)
	13	IERR = 0
	14	E(1) = 1.0
	15	IF (N .EQ. 1) GO TO 1001
	16	DO 100 I = 2, N
	17	IF (E(I) .NE. 0.0) GO TO 80
	18	IF (T(I,1) .NE. 0.0 .OR. T(I-1,3) .NE. 0.0) GO TO 1000
	19	E(I) = 1.0
	20	GO TO 100
	21	80 E(I) = E(I-1) * E(I) / T(I-1,3)
	22	100 CONTINUE
	23	DO 120 J = 1, M
	24	DO 120 I = 2, N
	25	Z(I,J) = Z(I,J) * E(I)
	26	120 CONTINUE
	27	GO TO 1001
	28	1000 IERR = 2 * N + I
	29	1001 RETURN
	30	END