1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
|
/* Library libcerf:
* Compute complex error functions, based on a new implementation of
* Faddeeva's w_of_z. Also provide Dawson and Voigt functions.
*
* File runvoigt.c:
* Interactive evaluation of Voigt's function.
*
* Copyright:
* (C) 2013 Forschungszentrum Jülich GmbH
*
* Licence:
* Public domain.
*
* Author:
* Joachim Wuttke, Forschungszentrum Jülich, 2013
*
* Website:
* http://apps.jcns.fz-juelich.de/libcerf
*/
#include <stdio.h>
#include <stdlib.h>
#include "cerf.h"
#include "defs.h"
#if _WIN32
#define IMPORT __declspec(dllimport)
#else
#define IMPORT
#endif
IMPORT extern int faddeeva_algorithm;
IMPORT extern int faddeeva_nofterms;
int main( int argc, char **argv )
{
double x, y;
if( argc!=3 ){
fprintf( stderr, "usage:\n" );
fprintf( stderr, " run_w_of_z <x> <y>\n" );
exit(-1);
}
x = atof( argv[1] );
y = atof( argv[2] );
_cerf_cmplx w = w_of_z( C(x,y) );
double v[2][2];
v[0][0] = creal(w);
v[0][1] = cimag(w);
printf( "%25.19g %25.19g %3i %3i\n", v[0][0], v[0][1],
faddeeva_algorithm, faddeeva_nofterms );
return 0;
}
|
