6 subroutine k_tf_sp(n,rhoup,rhodn,dtdru,dtdrd)
9 integer,
intent(in) :: n
10 real(8),
intent(in) :: rhoup(n),rhodn(n)
11 real(8),
intent(out) :: dtdru(n),dtdrd(n)
15 call k_tf1(2.d0*rhoup(i),dtdru(i))
18 call k_tf1(2.d0*rhodn(i),dtdrd(i))
subroutine k_tf_sp(n, rhoup, rhodn, dtdru, dtdrd)
elemental subroutine k_tf1(rho, dtdr)