]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1996/04/01 15:02:22 mclareni | |
6 | * Mathlib gen | |
7 | * | |
8 | * | |
9 | #include "gen/pilot.h" | |
10 | #if defined(CERNLIB_DOUBLE) | |
11 | FUNCTION DFRDH2(F,G,X,WS,IDIM,N) | |
12 | C | |
13 | #include "gen/imp64.inc" | |
14 | C | |
15 | CHARACTER*(*) NAME | |
16 | PARAMETER(NAME='DFRDH2') | |
17 | #endif | |
18 | #if !defined(CERNLIB_DOUBLE) | |
19 | FUNCTION RFRDH2(F,G,X,WS,IDIM,N) | |
20 | C | |
21 | CHARACTER*(*) NAME | |
22 | PARAMETER(NAME='RFRDH2') | |
23 | #endif | |
24 | C | |
25 | DIMENSION WS(IDIM,*) | |
26 | ||
27 | SUM=F(X) | |
28 | DO 2 J = 1,N | |
29 | 2 SUM=SUM+WS(J,N+4)*G(X,WS(J,N+3))*WS(J,N+1) | |
30 | GO TO 9 | |
31 | ||
32 | #if defined(CERNLIB_DOUBLE) | |
33 | ENTRY DFRDH3(H,WS,IDIM,N) | |
34 | #endif | |
35 | #if !defined(CERNLIB_DOUBLE) | |
36 | ENTRY RFRDH3(H,WS,IDIM,N) | |
37 | #endif | |
38 | ||
39 | SUM=0 | |
40 | DO 3 J = 1,N | |
41 | 3 SUM=SUM+WS(J,N+4)*H(WS(J,N+3))*WS(J,N+1) | |
42 | #if defined(CERNLIB_DOUBLE) | |
43 | 9 DFRDH2=SUM | |
44 | #endif | |
45 | #if !defined(CERNLIB_DOUBLE) | |
46 | 9 RFRDH2=SUM | |
47 | #endif | |
48 | RETURN | |
49 | END |