blob: 48afca63d6d65acaebb2b383458d462b18a6eb83 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
#include "f2c.h"
#ifdef KR_headers
double log(), f__cabs(), atan2();
VOID z_log(resx, z) doublecomplex *resx, *z;
#else
#undef abs
#include <math.h>
extern double f__cabs(double, double);
void z_log(doublecomplex *resx, doublecomplex *z)
#endif
{
doublecomplex res;
res.i = atan2(z->i, z->r);
res.r = log( f__cabs( z->r, z->i ) );
resx->r = res.r;
resx->i = res.i;
}
|
