1c     Correlation energy and potential as determined via the local
2c     density functional of Vosko, Wilk and Nusair (Can. J. of Phys. 58,
3c     1200 (1980)), functional I, the RPA parameterization.
4*
5* $Id$
6*
7#define FNL_VWN_1_RPA
8#include "xc_vwn_a.fh"
9c
10c     Compute the function deltaEc(rs,zeta) function and its derivatives
11c     wrt rs and zeta for the spin-unrestricted case - the rest has the
12c     same form for all VWN functionals and is handled in the header files.
13c
14            dec_rsz = dec_rs1*fz
15            d1dec_rsz(1) = d1dec_rs1*fz
16            d1dec_rsz(2) = dec_rs1*d1fz
17#ifdef SECOND_DERIV
18            d2dec_rsz(1) = d2dec_rs1*fz
19            d2dec_rsz(2) = d1dec_rs1*d1fz
20            d2dec_rsz(3) = dec_rs1*d2fz
21#endif
22#ifdef THIRD_DERIV
23            d3dec_rsz(1) = d3dec_rs1*fz
24            d3dec_rsz(2) = d2dec_rs1*d1fz
25            d3dec_rsz(3) = d1dec_rs1*d2fz
26            d3dec_rsz(4) = dec_rs1*d3fz
27#endif
28c
29#include "xc_vwn_b.fh"
30c
31#ifndef SECOND_DERIV
32#define SECOND_DERIV
33c
34c     Compile source again for the 2nd derivative case
35c
36#include "xc_vwn_1_rpa.F"
37#endif
38#ifndef THIRD_DERIV
39#define THIRD_DERIV
40c
41c     Compile source again for the 3rd derivative case
42c
43#include "xc_vwn_1_rpa.F"
44#endif
45