friend pvector exc(const pvector& in1); friend pvector excp(const pvector& in1);Return values:
- exc: pvector, each of whose data elements is the result of the application of excVWN() to the components of in1: out.dat[i][j][k]=excVWN(rs), where rs=pow( ((4*pi)/3.)*in1.dat[i][j][k], -1./3.) .
- excp: pvector, each of whose data elements is
for each element of in1, out.dat[i][j][k]==excpVWN(rs)*(-rs/in1.dat[i][j][k]/3.) .
Input:
- in1: input pvector (containing the charge density)
Debugging: These routines will be debugged in Part 3.