~ubuntu-branches/ubuntu/maverick/yagiuda/maverick

« back to all changes in this revision

Viewing changes to src/cis_hack.c

Committer: Bazaar Package Importer
Author(s): Joop Stakenborg
Date: 2005-08-22 20:20:13 UTC
Revision ID: james.westby@ubuntu.com-20050822202013-mhhxp4xirdxrdfx1

Tags: upstream-1.19

Import upstream version 1.19

files added:

AUTHORS

BUGS

COPYING

ChangeLog

INSTALL

Makefile.am

Makefile.in

NEWS

README

TODO

aclocal.m4

config.guess

config.sub

configure

configure.ac

depcomp

doc/Makefile.am

doc/Makefile.in

doc/dipole.1.html

doc/first.1.html

doc/index.html

doc/input.1.html

doc/optimise.1.html

doc/output.1.html

doc/pattern.jpg

doc/yagi.1.html

doc/yagi.jpg

install-sh

man/Makefile.am

man/Makefile.in

man/man1

man/man1/Makefile.am

man/man1/Makefile.in

man/man1/dipole.1

man/man1/first.1

man/man1/input.1

man/man1/optimise.1

man/man1/output.1

man/man1/yagi.1

man/man5

man/man5/Makefile.am

man/man5/Makefile.in

man/man5/first.5

man/man5/input.5

man/man5/optimise.5

man/man5/output.5

man/man5/yagi.5

missing

mkinstalldirs

src/3-element-yagi

src/Makefile.am

src/Makefile.in

src/auto.c

src/better.c

src/change_max_percentage_changes.c

src/check_flags.c

src/ci.c

src/cin.c

src/cis_hack.c

src/com_hack.c

src/com_hack.h

src/copy_matrix.c

src/dipole.c

src/display_antenna_currents.c

src/dobetter.c

src/dynamic.c

src/end_stop.c

src/error_message.c

src/fill_v_vector.c

src/fill_z_matrix.c

src/first.c

src/fitness.c

src/gain.c

src/gaussian.c

src/genetic.c

src/genetic_algorithm_lib.c

src/get_command_line_options.c

src/get_number_of_elements.c

src/getfiles.c

src/getopt.c

src/globals.h

src/input.c

src/install-sh

src/linear.c

src/lub_hack.c

src/lud_hack.c

src/max_gain.c

src/max_side.c

src/missing

src/mkinstalldirs

src/mpr_hack.c

src/mutual.c

src/new_length.c

src/nr_hack.h

src/nrutil.c

src/nrutil.h

src/optimise.c

src/optimising_for.c

src/output.c

src/perform.c

src/power_input.c

src/print_relavent_performance_data.c

src/print_z_matrix.c

src/random.c

src/randomise.c

src/rantest.c

src/read_header.c

src/read_yagi_data.c

src/reflection_coefficient.c

src/self.c

src/self2.c

src/selftest.c

src/sens.c

src/set_mean_structure.c

src/set_performance_structures.c

src/show_all_optimise_parameters.c

src/solve.c

src/string.c

src/subtract_structures.c

src/test2.c

src/test_stop.c

src/usage_first.c

src/usage_input.c

src/usage_optimise.c

src/usage_output.c

src/usage_yagi.c

src/version.c

src/vswr.c

src/write_coordinates.c

src/write_gain.c

src/write_gnuplot_header.c

src/write_header_to_disk.c

src/write_input_data.c

src/yagi.c

src/yagi.h

src/z_input.c

src/zbr_hack.c

Show diffs side-by-side

added added

removed removed

src/cis_hack.c

#ifdef HAVE_STDLIB_H

#include <stdlib.h>

#endif

#include <math.h>

#include "yagi.h"

#define EPS 6.0e-8

#define EULER 0.57721566

#define MAXIT 100

#define PIBY2 1.5707963

#define FPMIN 1.0e-30

#define TMIN 2.0

#define TRUE 1

#define ONE Complex(1.0,0.0)

void cisi(double x, double *ci, double *si)

{

int i,k,odd;

double a,err,fact,sign,sum,sumc,sums,t,term;

fcomplex h,b,c,d,del;

t=fabs(x);

if (t == 0.0) {

*si=0.0;

*ci = -1.0/FPMIN;

return;

}

if (t > TMIN) {

b=Complex(1.0,t);

c=Complex(1.0/FPMIN,0.0);

d=h=Cdiv(ONE,b);

for (i=2;i<=MAXIT;i++) {

a = -(i-1)*(i-1);

b=Cadd(b,Complex(2.0,0.0));

d=Cdiv(ONE,Cadd(RCmul(a,d),b));

c=Cadd(b,Cdiv(Complex(a,0.0),c));

del=Cmul(c,d);

h=Cmul(h,del);

if (fabs(del.r-1.0)+fabs(del.i) < EPS) break;

}

if (i > MAXIT) nrerror("cf failed in cisi_hack.c");

h=Cmul(Complex(cos(t),-sin(t)),h);

*ci = -h.r;

*si=PIBY2+h.i;

} else {

if (t < sqrt(FPMIN)) {

sumc=0.0;

sums=t;

} else {

sum=sums=sumc=0.0;

sign=fact=1.0;

odd=TRUE;

for (k=1;k<=MAXIT;k++) {

fact *= t/k;

term=fact/k;

sum += sign*term;

err=term/fabs(sum);

if (odd) {

sign = -sign;

sums=sum;

sum=sumc;

} else {

sumc=sum;

sum=sums;

}

if (err < EPS) break;

odd=!odd;

}

if (k > MAXIT) nrerror("maxits exceeded in cisi");

}

*si=sums;

*ci=sumc+log(t)+EULER;

}

if (x < 0.0) *si = -(*si);

}

#undef EPS

#undef EULER

#undef MAXIT

#undef PIBY2

#undef FPMIN

#undef TMIN

#undef TRUE

#undef ONE

Older »