4 * This source code is part of
8 * GROningen MAchine for Chemical Simulations
11 * Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
12 * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
13 * Copyright (c) 2001-2004, The GROMACS development team,
14 * check out http://www.gromacs.org for more information.
16 * This program is free software; you can redistribute it and/or
17 * modify it under the terms of the GNU General Public License
18 * as published by the Free Software Foundation; either version 2
19 * of the License, or (at your option) any later version.
21 * If you want to redistribute modifications, please consider that
22 * scientific software is very special. Version control is crucial -
23 * bugs must be traceable. We will be happy to consider code for
24 * inclusion in the official distribution, but derived work must not
25 * be called official GROMACS. Details are found in the README & COPYING
26 * files - if they are missing, get the official version at www.gromacs.org.
28 * To help us fund GROMACS development, we humbly ask that you cite
29 * the papers on the package - you can find them in the top README file.
31 * For more info, check our website at http://www.gromacs.org
34 * Gromacs Runs On Most of All Computer Systems
50 typedef t_complex cvec
[DIM
];
52 static t_complex cnul
= { 0.0, 0.0 };
54 static t_complex
rcmul(real r
,t_complex c
)
64 static t_complex
rcexp(real r
)
75 static t_complex
cadd(t_complex a
,t_complex b
)
85 static t_complex
csub(t_complex a
,t_complex b
)
95 static t_complex
cmul(t_complex a
,t_complex b
)
99 c
.re
= a
.re
*b
.re
- a
.im
*b
.im
;
100 c
.im
= a
.re
*b
.im
+ a
.im
*b
.re
;
105 static t_complex
conjugate(t_complex c
)
115 static t_complex
cdiv(t_complex teller
,t_complex noemer
)
117 t_complex res
,anoemer
;
119 anoemer
= cmul(conjugate(noemer
),noemer
);
120 res
= cmul(teller
,conjugate(noemer
));
122 return rcmul(1.0/anoemer
.re
,res
);