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

« back to all changes in this revision

Viewing changes to src/lud_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/lud_hack.c

#ifdef HAVE_STDLIB_H

#include <stdlib.h>

#endif

#include <math.h>

#define NRANSI

#include "nrutil.h"

/* #include "yagi.h" */

#define TINY 1.0e-20;

void ludcmp(double **a, int n, int *indx, double *d)

{

int i,imax,j,k;

double big,dum,sum,temp;

double *vv;

vv=dvector(1,n);

*d=1.0;

for (i=1;i<=n;i++) {

big=0.0;

for (j=1;j<=n;j++)

if ((temp=fabs(a[i][j])) > big) big=temp;

if (big == 0.0) nrerror("Singular matrix in routine ludcmp");

vv[i]=1.0/big;

}

for (j=1;j<=n;j++) {

for (i=1;i<j;i++) {

sum=a[i][j];

for (k=1;k<i;k++) sum -= a[i][k]*a[k][j];

a[i][j]=sum;

}

big=0.0;

for (i=j;i<=n;i++) {

sum=a[i][j];

for (k=1;k<j;k++)

sum -= a[i][k]*a[k][j];

a[i][j]=sum;

if ( (dum=vv[i]*fabs(sum)) >= big) {

big=dum;

imax=i;

}

if (j != imax) {

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

dum=a[imax][k];

a[imax][k]=a[j][k];

a[j][k]=dum;

}

*d = -(*d);

vv[imax]=vv[j];

}

indx[j]=imax;

if (a[j][j] == 0.0) a[j][j]=TINY;

if (j != n) {

dum=1.0/(a[j][j]);

for (i=j+1;i<=n;i++) a[i][j] *= dum;

}

free_dvector(vv,1,n);

}

#undef TINY

#undef NRANSI

/* (C) Copr. 1986-92 Numerical Recipes Software 5#,. */

Older »