1(*
2 Copyright (C) 2017 M.A.L. Marques
3
4 This Source Code Form is subject to the terms of the Mozilla Public
5 License, v. 2.0. If a copy of the MPL was not distributed with this
6 file, You can obtain one at http://mozilla.org/MPL/2.0/.
7*)
8
9(* type: gga_exc *)
10
11lag_a1 :=   0.041106:
12lag_a2 :=   2.626712:
13lag_a3 :=   0.092070:
14lag_a4 :=   0.657946:
15
16lag_f0 := s-> lag_a1 * s^lag_a2/(1 + lag_a3 * s^lag_a2)^lag_a4:
17lag_f  := x-> lag_f0(X2S*x):
18
19f := (rs, zeta, xt, xs0, xs1) -> gga_exchange(lag_f, rs, zeta, xs0, xs1):
20