21 lines
370 B
C
21 lines
370 B
C
|
#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;
|
||
|
|
}
|