[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / minfccod.inc

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/04/01 15:02:20  mclareni
* Mathlib gen
*
*
*
* minfccod.inc
*
      PARAMETER (W5 = 2.23606 79774 99790D0)
      PARAMETER (HV = (3-W5)/2, HW = (W5-1)/2, R1 = 1, HF = R1/2)

      N=-1
      IF(A .NE. B) N=NINT(2.08*LOG(ABS((A-B)/EPS)))
      C=A
      D=B
      IF(A .GT. B) THEN
       C=B
       D=A
      ENDIF
      LLT=.TRUE.
      LGE=.TRUE.

    1 H=D-C
      IF(N .LT. 0) THEN
       X=HF*(C+D)
       Y=F(X)
       LLM=ABS(X-A) .GT. DELTA .AND. ABS(X-B) .GT. DELTA
       RETURN
      ENDIF
      IF(LLT) THEN
       V=C+HV*H
       FV=F(V)
      ENDIF
      IF(LGE) THEN
       W=C+HW*H
       FW=F(W)
      ENDIF
      IF(FV .LT. FW) THEN
       LLT=.TRUE.
       LGE=.FALSE.
       D=W
       W=V
       FW=FV
      ELSE
       LLT=.FALSE.
       LGE=.TRUE.
       C=V
       V=W
       FV=FW
      ENDIF
      N=N-1
      GO TO 1
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/04/01 15:02:20 mclareni
	6	* Mathlib gen
	7	*
	8	*
	9	*
	10	* minfccod.inc
	11	*
	12	PARAMETER (W5 = 2.23606 79774 99790D0)
	13	PARAMETER (HV = (3-W5)/2, HW = (W5-1)/2, R1 = 1, HF = R1/2)
	14
	15	N=-1
	16	IF(A .NE. B) N=NINT(2.08*LOG(ABS((A-B)/EPS)))
	17	C=A
	18	D=B
	19	IF(A .GT. B) THEN
	20	C=B
	21	D=A
	22	ENDIF
	23	LLT=.TRUE.
	24	LGE=.TRUE.
	25
	26	1 H=D-C
	27	IF(N .LT. 0) THEN
	28	X=HF*(C+D)
	29	Y=F(X)
	30	LLM=ABS(X-A) .GT. DELTA .AND. ABS(X-B) .GT. DELTA
	31	RETURN
	32	ENDIF
	33	IF(LLT) THEN
	34	V=C+HV*H
	35	FV=F(V)
	36	ENDIF
	37	IF(LGE) THEN
	38	W=C+HW*H
	39	FW=F(W)
	40	ENDIF
	41	IF(FV .LT. FW) THEN
	42	LLT=.TRUE.
	43	LGE=.FALSE.
	44	D=W
	45	W=V
	46	FW=FV
	47	ELSE
	48	LLT=.FALSE.
	49	LGE=.TRUE.
	50	C=V
	51	V=W
	52	FV=FW
	53	ENDIF
	54	N=N-1
	55	GO TO 1
	56	END