at.mug.iqm.commons.lau

Class Linear_algebra

java.lang.Object

at.mug.iqm.commons.lau.Linear_algebra

public class Linear_algebra
extends Object

This is from: Hang T. Lau A Numerical Library in Java for Scientists & Engineers Chapman & Hall/Crc 2003

Since:

2010 05

Author:

Helmut Ahammer

Constructor Summary

Constructors

Constructor and Description
Linear_algebra()

Method Summary

Modifier and Type	Method and Description
static void	alljaczer(int n, double alfa, double beta, double[] zer)
static void	alllagzer(int n, double alfa, double[] zer)
static void	allzerortpol(int n, double[] b, double[] c, double[] zer, double[] em)
static void	bakcomhes(double[][] ar, double[][] ai, double[] tr, double[] ti, double[] del, double[][] vr, double[][] vi, int n, int n1, int n2)
static void	bakhrmtri(double[][] a, int n, int n1, int n2, double[][] vecr, double[][] veci, double[] tr, double[] ti)
static void	baklbr(int n, int n1, int n2, double[] d, int[] inter, double[][] vec)
static void	baklbrcom(int n, int n1, int n2, double[] d, int[] inter, double[][] vr, double[][] vi)
static void	bakreahes1(double[][] a, int n, int[] index, double[] v)
static void	bakreahes2(double[][] a, int n, int n1, int n2, int[] index, double[][] vec)
static void	baksymtri1(double[] a, int n, int n1, int n2, double[][] vec)
static void	baksymtri2(double[][] a, int n, int n1, int n2, double[][] vec)
static void	bounds(int n, double[] a, double[] re, double[] im, double rele, double abse, double[] recentre, double[] imcentre, double[] bound)
static void	chldec1(double[] a, int n, double[] aux)
static void	chldec2(double[][] a, int n, double[] aux)
static void	chldecbnd(double[] a, int n, int w, double[] aux)
static void	chldecinv1(double[] a, int n, double[] aux)
static void	chldecinv2(double[][] a, int n, double[] aux)
static void	chldecsol1(double[] a, int n, double[] aux, double[] b)
static void	chldecsol2(double[][] a, int n, double[] aux, double[] b)
static void	chldecsolbnd(double[] a, int n, int w, double[] aux, double[] b)
static double	chldeterm1(double[] a, int n)
static double	chldeterm2(double[][] a, int n)
static double	chldetermbnd(double[] a, int n, int w)
static void	chlinv1(double[] a, int n)
static void	chlinv2(double[][] a, int n)
static void	chlsol1(double[] a, int n, double[] b)
static void	chlsol2(double[][] a, int n, double[] b)
static void	chlsolbnd(double[] a, int n, int w, double[] b)
static int	comeig1(double[][] a, int n, double[] em, double[] re, double[] im, double[][] vec)
static int	comeigval(double[][] a, int n, double[] em, double[] re, double[] im)
static void	comkwd(double pr, double pi, double qr, double qi, double[] gr, double[] gi, double[] kr, double[] ki)
static int	comvalqri(double[][] a, int n, double[] em, double[] re, double[] im)
static void	comveches(double[][] a, int n, double lambda, double mu, double[] em, double[] u, double[] v)
static void	conjgrad(LA_conjgrad_methods method, double[] x, double[] r, int l, int n, int[] iterate, double[] norm2)
static void	dec(double[][] a, int n, double[] aux, int[] p)
static void	decbnd(double[] a, int n, int lw, int rw, double[] aux, double[] m, int[] p)
static void	decinv(double[][] a, int n, double[] aux)
static void	decsol(double[][] a, int n, double[] aux, double[] b)
static void	decsolbnd(double[] a, int n, int lw, int rw, double[] aux, double[] b)
static void	decsolsym2(double[][] a, int n, double[] b, double tol, int[] aux)
static void	decsolsymtri(double[] diag, double[] co, int n, double[] aux, double[] b)
static void	decsoltri(double[] sub, double[] diag, double[] supre, int n, double[] aux, double[] b)
static void	decsoltripiv(double[] sub, double[] diag, double[] supre, int n, double[] aux, double[] b)
static void	decsym2(double[][] a, int n, double tol, int[] aux, int[] p, double[] detaux)
static void	decsymtri(double[] diag, double[] co, int n, double[] aux)
static void	dectri(double[] sub, double[] diag, double[] supre, int n, double[] aux)
static void	dectripiv(double[] sub, double[] diag, double[] supre, int n, double[] aid, double[] aux, boolean[] piv)
static double	determ(double[][] a, int n, int sign)
static double	determbnd(double[] a, int n, int lw, int rw, int sgndet)
static double	determsym2(double[] detaux, int n, int[] aux)
static int	eigcom(double[][] ar, double[][] ai, int n, double[] em, double[] valr, double[] vali, double[][] vr, double[][] vi)
static void	eighrm(double[][] a, int n, int numval, double[] val, double[][] vecr, double[][] veci, double[] em)
static void	eigsym1(double[] a, int n, int numval, double[] val, double[][] vec, double[] em)
static void	eigsym2(double[][] a, int n, int numval, double[] val, double[][] vec, double[] em)
static int	eigvalcom(double[][] ar, double[][] ai, int n, double[] em, double[] valr, double[] vali)
static void	eigvalhrm(double[][] a, int n, int numval, double[] val, double[] em)
static void	eigvalsym1(double[] a, int n, int numval, double[] val, double[] em)
static void	eigvalsym2(double[][] a, int n, int numval, double[] val, double[] em)
static void	eqilbr(double[][] a, int n, double[] em, double[] d, int[] inter)
static void	eqilbrcom(double[][] a1, double[][] a2, int n, double[] em, double[] d, int[] inter)
static void	erbelm(int n, double[] aux, double nrminv)
static void	gsselm(double[][] a, int n, double[] aux, int[] ri, int[] ci)
static void	gsserb(double[][] a, int n, double[] aux, int[] ri, int[] ci)
static void	gssinv(double[][] a, int n, double[] aux)
static void	gssinverb(double[][] a, int n, double[] aux)
static void	gssitisol(double[][] a, int n, double[] aux, double[] b)
static void	gssitisolerb(double[][] a, int n, double[] aux, double[] b)
static void	gssnri(double[][] a, int n, double[] aux, int[] ri, int[] ci)
static void	gsssol(double[][] a, int n, double[] aux, double[] b)
static void	gsssolerb(double[][] a, int n, double[] aux, double[] b)
static void	hestgl2(int n, double[][] a, double[][] b)
static void	hestgl3(int n, double[][] a, double[][] b, double[][] x)
static int	homsol(double[][] a, int m, int n, double[][] v, double[] em)
static void	homsolsvd(double[][] u, double[] val, double[][] v, int m, int n)
static void	hsh2col(int la, int lb, int u, int i, double a1, double a2, double[][] a, double[][] b)
static void	hsh2row2(int l, int ua, int ub, int j, double a1, double a2, double[][] a, double[][] b)
static void	hsh2row3(int l, int ua, int ub, int ux, int j, double a1, double a2, double[][] a, double[][] b, double[][] x)
static void	hsh3col(int la, int lb, int u, int i, double a1, double a2, double a3, double[][] a, double[][] b)
static void	hsh3row2(int l, int u, int j, double a1, double a2, double a3, double[][] a, double[][] b)
static void	hsh3row3(int l, int u, int ux, int j, double a1, double a2, double a3, double[][] a, double[][] b, double[][] x)
static void	hshcomhes(double[][] ar, double[][] ai, int n, double[] em, double[] b, double[] tr, double[] ti, double[] del)
static void	hshdecmul(int n, double[][] a, double[][] b, double dwarf)
static void	hshhrmtri(double[][] a, int n, double[] d, double[] b, double[] bb, double[] em, double[] tr, double[] ti)
static void	hshhrmtrival(double[][] a, int n, double[] d, double[] bb, double[] em)
static void	hshreabid(double[][] a, int m, int n, double[] d, double[] b, double[] em)
static void	inv(double[][] a, int n, int[] p)
static double	inv1(double[][] a, int n, int[] ri, int[] ci, boolean withnorm)
static void	itisol(double[][] a, double[][] lu, int n, double[] aux, int[] ri, int[] ci, double[] b)
static void	itisolerb(double[][] a, double[][] lu, int n, double[] aux, int[] ri, int[] ci, double[] b)
private static void	kcluster(int k, int m, int n, double[] rc, double[] re, double[] im, double[] recentre, double[] imcentre, double[] bound, double[] clust)
static void	lsqdecomp(double[][] a, int n, int m, int n1, double[] aux, double[] aid, int[] ci)
static void	lsqdglinv(double[][] a, int m, double[] aid, int[] ci, double[] diag)
static void	lsqinv(double[][] a, int m, double[] aid, int[] c)
static void	lsqortdec(double[][] a, int n, int m, double[] aux, double[] aid, int[] ci)
static void	lsqortdecsol(double[][] a, int n, int m, double[] aux, double[] diag, double[] b)
static void	lsqrefsol(double[][] a, double[][] qr, int n, int m, int n1, double[] aux, double[] aid, int[] ci, double[] b, double[] ldx, double[] x, double[] res)
static void	lsqsol(double[][] a, int n, int m, double[] aid, int[] ci, double[] b)
static void	lupzerortpol(int n, int m, double[] b, double[] c, double[] zer, double[] em)
private static void	merge(int lo, int ls, int rs, int[] p, double[] a, int[] hp)
static void	mergesort(double[] a, int[] p, int low, int up)
static double	onenrminv(double[][] a, int n)
static void	orthog(int n, int lc, int uc, double[][] x)
static void	pretfmmat(double[][] a, int m, int n, double[] d)
static int	psdinv(double[][] a, int m, int n, double[] em)
static void	psdinvsvd(double[][] u, double[] val, double[][] v, int m, int n, double[] em)
static void	psttfmmat(double[][] a, int n, double[][] v, double[] b)
static int	qricom(double[][] a1, double[][] a2, double[] b, int n, double[] em, double[] val1, double[] val2, double[][] vec1, double[][] vec2)
static int	qrihrm(double[][] a, int n, double[] val, double[][] vr, double[][] vi, double[] em)
static int	qrisngval(double[][] a, int m, int n, double[] val, double[] em)
static int	qrisngvalbid(double[] d, double[] b, int n, double[] em)
static int	qrisngvaldec(double[][] a, int m, int n, double[] val, double[][] v, double[] em)
static int	qrisngvaldecbid(double[] d, double[] b, int m, int n, double[][] u, double[][] v, double[] em)
static int	qrisym(double[][] a, int n, double[] val, double[] em)
static int	qrisymtri(double[][] a, int n, double[] d, double[] b, double[] bb, double[] em)
static int	qrivalhrm(double[][] a, int n, double[] val, double[] em)
static int	qrivalsym1(double[] a, int n, double[] val, double[] em)
static int	qrivalsym2(double[][] a, int n, double[] val, double[] em)
static int	qrivalsymtri(double[] d, double[] bb, int n, double[] em)
static void	qzi(int n, double[][] a, double[][] b, double[][] x, double[] alfr, double[] alfi, double[] beta, int[] iter, double[] em)
static void	qzival(int n, double[][] a, double[][] b, double[] alfr, double[] alfi, double[] beta, int[] iter, double[] em)
static int	reaeig1(double[][] a, int n, double[] em, double[] val, double[][] vec)
static int	reaeig3(double[][] a, int n, double[] em, double[] val, double[][] vec)
static int	reaeigval(double[][] a, int n, double[] em, double[] val)
static int	reaqri(double[][] a, int n, double[] em, double[] val, double[][] vec)
static int	reavalqri(double[][] a, int n, double[] em, double[] val)
static void	reaveches(double[][] a, int n, double lambda, double[] em, double[] v)
static void	rowperm(int[] perm, int low, int upp, int i, double[][] mat)
static void	selzerortpol(int n, int n1, int n2, double[] b, double[] c, double[] zer, double[] em)
static void	sol(double[][] a, int n, int[] p, double[] b)
static void	solbnd(double[] a, int n, int lw, int rw, double[] m, int[] p, double[] b)
static void	solelm(double[][] a, int n, int[] ri, int[] ci, double[] b)
static int	solovr(double[][] a, int m, int n, double[] x, double[] em)
static void	solsvdovr(double[][] u, double[] val, double[][] v, int m, int n, double[] x, double[] em)
static void	solsvdund(double[][] u, double[] val, double[][] v, int m, int n, double[] x, double[] em)
static void	solsym2(double[][] a, int n, double[] b, int[] p, double[] detaux)
static void	solsymtri(double[] diag, double[] co, int n, double[] b)
static void	soltri(double[] sub, double[] diag, double[] supre, int n, double[] b)
static void	soltripiv(double[] sub, double[] diag, double[] supre, int n, double[] aid, boolean[] piv, double[] b)
static int	solund(double[][] a, int m, int n, double[] x, double[] em)
private static double	sturm(double[] d, double[] bb, int n, double x, int k, double machtol, double max, int[] count, double[] lb, double[] ub)
static void	symeigimp(int n, double[][] a, double[][] vec, double[] val, double[] lbound, double[] ubound, double[] aux)
static void	tfmprevec(double[][] a, int n)
static void	tfmreahes(double[][] a, int n, double[] em, int[] index)
static void	tfmsymtri1(double[] a, int n, double[] d, double[] b, double[] bb, double[] em)
static void	tfmsymtri2(double[][] a, int n, double[] d, double[] b, double[] bb, double[] em)
static int	valqricom(double[][] a1, double[][] a2, double[] b, int n, double[] em, double[] val1, double[] val2)
static void	valsymtri(double[] d, double[] bb, int n, int n1, int n2, double[] val, double[] em)
static void	vecperm(int[] perm, int low, int upp, double[] vector)
static void	vecsymtri(double[] d, double[] b, int n, int n1, int n2, double[] val, double[][] vec, double[] em)
static int	zerpol(int n, double[] a, double[] em, double[] re, double[] im, double[] d)
private static boolean	zerpolfunction(int n, double[] d, double[] f, double x, double y, double tol, int[] it, double[] newf)

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

Linear_algebra

public Linear_algebra()

Method Detail

dec

public static void dec(double[][] a,
                       int n,
                       double[] aux,
                       int[] p)

gsselm

public static void gsselm(double[][] a,
                          int n,
                          double[] aux,
                          int[] ri,
                          int[] ci)

onenrminv

public static double onenrminv(double[][] a,
                               int n)

erbelm

public static void erbelm(int n,
                          double[] aux,
                          double nrminv)

gsserb

public static void gsserb(double[][] a,
                          int n,
                          double[] aux,
                          int[] ri,
                          int[] ci)

gssnri

public static void gssnri(double[][] a,
                          int n,
                          double[] aux,
                          int[] ri,
                          int[] ci)

determ

public static double determ(double[][] a,
                            int n,
                            int sign)

sol

public static void sol(double[][] a,
                       int n,
                       int[] p,
                       double[] b)

decsol

public static void decsol(double[][] a,
                          int n,
                          double[] aux,
                          double[] b)

solelm

public static void solelm(double[][] a,
                          int n,
                          int[] ri,
                          int[] ci,
                          double[] b)

gsssol

public static void gsssol(double[][] a,
                          int n,
                          double[] aux,
                          double[] b)

gsssolerb

public static void gsssolerb(double[][] a,
                             int n,
                             double[] aux,
                             double[] b)

inv

public static void inv(double[][] a,
                       int n,
                       int[] p)

decinv

public static void decinv(double[][] a,
                          int n,
                          double[] aux)

inv1

public static double inv1(double[][] a,
                          int n,
                          int[] ri,
                          int[] ci,
                          boolean withnorm)

gssinv

public static void gssinv(double[][] a,
                          int n,
                          double[] aux)

gssinverb

public static void gssinverb(double[][] a,
                             int n,
                             double[] aux)

itisol

public static void itisol(double[][] a,
                          double[][] lu,
                          int n,
                          double[] aux,
                          int[] ri,
                          int[] ci,
                          double[] b)

gssitisol

public static void gssitisol(double[][] a,
                             int n,
                             double[] aux,
                             double[] b)

itisolerb

public static void itisolerb(double[][] a,
                             double[][] lu,
                             int n,
                             double[] aux,
                             int[] ri,
                             int[] ci,
                             double[] b)

gssitisolerb

public static void gssitisolerb(double[][] a,
                                int n,
                                double[] aux,
                                double[] b)

chldec2

public static void chldec2(double[][] a,
                           int n,
                           double[] aux)

chldec1

public static void chldec1(double[] a,
                           int n,
                           double[] aux)

chldeterm2

public static double chldeterm2(double[][] a,
                                int n)

chldeterm1

public static double chldeterm1(double[] a,
                                int n)

chlsol2

public static void chlsol2(double[][] a,
                           int n,
                           double[] b)

chlsol1

public static void chlsol1(double[] a,
                           int n,
                           double[] b)

chldecsol2

public static void chldecsol2(double[][] a,
                              int n,
                              double[] aux,
                              double[] b)

chldecsol1

public static void chldecsol1(double[] a,
                              int n,
                              double[] aux,
                              double[] b)

chlinv2

public static void chlinv2(double[][] a,
                           int n)

chlinv1

public static void chlinv1(double[] a,
                           int n)

chldecinv2

public static void chldecinv2(double[][] a,
                              int n,
                              double[] aux)

chldecinv1

public static void chldecinv1(double[] a,
                              int n,
                              double[] aux)

decsym2

public static void decsym2(double[][] a,
                           int n,
                           double tol,
                           int[] aux,
                           int[] p,
                           double[] detaux)

determsym2

public static double determsym2(double[] detaux,
                                int n,
                                int[] aux)

solsym2

public static void solsym2(double[][] a,
                           int n,
                           double[] b,
                           int[] p,
                           double[] detaux)

decsolsym2

public static void decsolsym2(double[][] a,
                              int n,
                              double[] b,
                              double tol,
                              int[] aux)

lsqortdec

public static void lsqortdec(double[][] a,
                             int n,
                             int m,
                             double[] aux,
                             double[] aid,
                             int[] ci)

lsqdglinv

public static void lsqdglinv(double[][] a,
                             int m,
                             double[] aid,
                             int[] ci,
                             double[] diag)

lsqsol

public static void lsqsol(double[][] a,
                          int n,
                          int m,
                          double[] aid,
                          int[] ci,
                          double[] b)

lsqortdecsol

public static void lsqortdecsol(double[][] a,
                                int n,
                                int m,
                                double[] aux,
                                double[] diag,
                                double[] b)

lsqinv

public static void lsqinv(double[][] a,
                          int m,
                          double[] aid,
                          int[] c)

lsqdecomp

public static void lsqdecomp(double[][] a,
                             int n,
                             int m,
                             int n1,
                             double[] aux,
                             double[] aid,
                             int[] ci)

lsqrefsol

public static void lsqrefsol(double[][] a,
                             double[][] qr,
                             int n,
                             int m,
                             int n1,
                             double[] aux,
                             double[] aid,
                             int[] ci,
                             double[] b,
                             double[] ldx,
                             double[] x,
                             double[] res)

solsvdovr

public static void solsvdovr(double[][] u,
                             double[] val,
                             double[][] v,
                             int m,
                             int n,
                             double[] x,
                             double[] em)

solovr

public static int solovr(double[][] a,
                         int m,
                         int n,
                         double[] x,
                         double[] em)

solsvdund

public static void solsvdund(double[][] u,
                             double[] val,
                             double[][] v,
                             int m,
                             int n,
                             double[] x,
                             double[] em)

solund

public static int solund(double[][] a,
                         int m,
                         int n,
                         double[] x,
                         double[] em)

homsolsvd

public static void homsolsvd(double[][] u,
                             double[] val,
                             double[][] v,
                             int m,
                             int n)

homsol

public static int homsol(double[][] a,
                         int m,
                         int n,
                         double[][] v,
                         double[] em)

psdinvsvd

public static void psdinvsvd(double[][] u,
                             double[] val,
                             double[][] v,
                             int m,
                             int n,
                             double[] em)

psdinv

public static int psdinv(double[][] a,
                         int m,
                         int n,
                         double[] em)

decbnd

public static void decbnd(double[] a,
                          int n,
                          int lw,
                          int rw,
                          double[] aux,
                          double[] m,
                          int[] p)

determbnd

public static double determbnd(double[] a,
                               int n,
                               int lw,
                               int rw,
                               int sgndet)

solbnd

public static void solbnd(double[] a,
                          int n,
                          int lw,
                          int rw,
                          double[] m,
                          int[] p,
                          double[] b)

decsolbnd

public static void decsolbnd(double[] a,
                             int n,
                             int lw,
                             int rw,
                             double[] aux,
                             double[] b)

dectri

public static void dectri(double[] sub,
                          double[] diag,
                          double[] supre,
                          int n,
                          double[] aux)

dectripiv

public static void dectripiv(double[] sub,
                             double[] diag,
                             double[] supre,
                             int n,
                             double[] aid,
                             double[] aux,
                             boolean[] piv)

soltri

public static void soltri(double[] sub,
                          double[] diag,
                          double[] supre,
                          int n,
                          double[] b)

decsoltri

public static void decsoltri(double[] sub,
                             double[] diag,
                             double[] supre,
                             int n,
                             double[] aux,
                             double[] b)

soltripiv

public static void soltripiv(double[] sub,
                             double[] diag,
                             double[] supre,
                             int n,
                             double[] aid,
                             boolean[] piv,
                             double[] b)

decsoltripiv

public static void decsoltripiv(double[] sub,
                                double[] diag,
                                double[] supre,
                                int n,
                                double[] aux,
                                double[] b)

chldecbnd

public static void chldecbnd(double[] a,
                             int n,
                             int w,
                             double[] aux)

chldetermbnd

public static double chldetermbnd(double[] a,
                                  int n,
                                  int w)

chlsolbnd

public static void chlsolbnd(double[] a,
                             int n,
                             int w,
                             double[] b)

chldecsolbnd

public static void chldecsolbnd(double[] a,
                                int n,
                                int w,
                                double[] aux,
                                double[] b)

decsymtri

public static void decsymtri(double[] diag,
                             double[] co,
                             int n,
                             double[] aux)

solsymtri

public static void solsymtri(double[] diag,
                             double[] co,
                             int n,
                             double[] b)

decsolsymtri

public static void decsolsymtri(double[] diag,
                                double[] co,
                                int n,
                                double[] aux,
                                double[] b)

conjgrad

public static void conjgrad(LA_conjgrad_methods method,
                            double[] x,
                            double[] r,
                            int l,
                            int n,
                            int[] iterate,
                            double[] norm2)

eqilbr

public static void eqilbr(double[][] a,
                          int n,
                          double[] em,
                          double[] d,
                          int[] inter)

baklbr

public static void baklbr(int n,
                          int n1,
                          int n2,
                          double[] d,
                          int[] inter,
                          double[][] vec)

eqilbrcom

public static void eqilbrcom(double[][] a1,
                             double[][] a2,
                             int n,
                             double[] em,
                             double[] d,
                             int[] inter)

baklbrcom

public static void baklbrcom(int n,
                             int n1,
                             int n2,
                             double[] d,
                             int[] inter,
                             double[][] vr,
                             double[][] vi)

tfmsymtri2

public static void tfmsymtri2(double[][] a,
                              int n,
                              double[] d,
                              double[] b,
                              double[] bb,
                              double[] em)

baksymtri2

public static void baksymtri2(double[][] a,
                              int n,
                              int n1,
                              int n2,
                              double[][] vec)

tfmprevec

public static void tfmprevec(double[][] a,
                             int n)

tfmsymtri1

public static void tfmsymtri1(double[] a,
                              int n,
                              double[] d,
                              double[] b,
                              double[] bb,
                              double[] em)

baksymtri1

public static void baksymtri1(double[] a,
                              int n,
                              int n1,
                              int n2,
                              double[][] vec)

tfmreahes

public static void tfmreahes(double[][] a,
                             int n,
                             double[] em,
                             int[] index)

bakreahes1

public static void bakreahes1(double[][] a,
                              int n,
                              int[] index,
                              double[] v)

bakreahes2

public static void bakreahes2(double[][] a,
                              int n,
                              int n1,
                              int n2,
                              int[] index,
                              double[][] vec)

hshhrmtri

public static void hshhrmtri(double[][] a,
                             int n,
                             double[] d,
                             double[] b,
                             double[] bb,
                             double[] em,
                             double[] tr,
                             double[] ti)

hshhrmtrival

public static void hshhrmtrival(double[][] a,
                                int n,
                                double[] d,
                                double[] bb,
                                double[] em)

bakhrmtri

public static void bakhrmtri(double[][] a,
                             int n,
                             int n1,
                             int n2,
                             double[][] vecr,
                             double[][] veci,
                             double[] tr,
                             double[] ti)

hshcomhes

public static void hshcomhes(double[][] ar,
                             double[][] ai,
                             int n,
                             double[] em,
                             double[] b,
                             double[] tr,
                             double[] ti,
                             double[] del)

bakcomhes

public static void bakcomhes(double[][] ar,
                             double[][] ai,
                             double[] tr,
                             double[] ti,
                             double[] del,
                             double[][] vr,
                             double[][] vi,
                             int n,
                             int n1,
                             int n2)

hshreabid

public static void hshreabid(double[][] a,
                             int m,
                             int n,
                             double[] d,
                             double[] b,
                             double[] em)

psttfmmat

public static void psttfmmat(double[][] a,
                             int n,
                             double[][] v,
                             double[] b)

pretfmmat

public static void pretfmmat(double[][] a,
                             int m,
                             int n,
                             double[] d)

valsymtri

public static void valsymtri(double[] d,
                             double[] bb,
                             int n,
                             int n1,
                             int n2,
                             double[] val,
                             double[] em)

sturm

private static double sturm(double[] d,
                            double[] bb,
                            int n,
                            double x,
                            int k,
                            double machtol,
                            double max,
                            int[] count,
                            double[] lb,
                            double[] ub)

vecsymtri

public static void vecsymtri(double[] d,
                             double[] b,
                             int n,
                             int n1,
                             int n2,
                             double[] val,
                             double[][] vec,
                             double[] em)

qrivalsymtri

public static int qrivalsymtri(double[] d,
                               double[] bb,
                               int n,
                               double[] em)

qrisymtri

public static int qrisymtri(double[][] a,
                            int n,
                            double[] d,
                            double[] b,
                            double[] bb,
                            double[] em)

eigvalsym2

public static void eigvalsym2(double[][] a,
                              int n,
                              int numval,
                              double[] val,
                              double[] em)

eigsym2

public static void eigsym2(double[][] a,
                           int n,
                           int numval,
                           double[] val,
                           double[][] vec,
                           double[] em)

eigvalsym1

public static void eigvalsym1(double[] a,
                              int n,
                              int numval,
                              double[] val,
                              double[] em)

eigsym1

public static void eigsym1(double[] a,
                           int n,
                           int numval,
                           double[] val,
                           double[][] vec,
                           double[] em)

qrivalsym2

public static int qrivalsym2(double[][] a,
                             int n,
                             double[] val,
                             double[] em)

qrisym

public static int qrisym(double[][] a,
                         int n,
                         double[] val,
                         double[] em)

qrivalsym1

public static int qrivalsym1(double[] a,
                             int n,
                             double[] val,
                             double[] em)

mergesort

public static void mergesort(double[] a,
                             int[] p,
                             int low,
                             int up)

merge

private static void merge(int lo,
                          int ls,
                          int rs,
                          int[] p,
                          double[] a,
                          int[] hp)

vecperm

public static void vecperm(int[] perm,
                           int low,
                           int upp,
                           double[] vector)

rowperm

public static void rowperm(int[] perm,
                           int low,
                           int upp,
                           int i,
                           double[][] mat)

orthog

public static void orthog(int n,
                          int lc,
                          int uc,
                          double[][] x)

symeigimp

public static void symeigimp(int n,
                             double[][] a,
                             double[][] vec,
                             double[] val,
                             double[] lbound,
                             double[] ubound,
                             double[] aux)

reavalqri

public static int reavalqri(double[][] a,
                            int n,
                            double[] em,
                            double[] val)

reaveches

public static void reaveches(double[][] a,
                             int n,
                             double lambda,
                             double[] em,
                             double[] v)

reaqri

public static int reaqri(double[][] a,
                         int n,
                         double[] em,
                         double[] val,
                         double[][] vec)

comvalqri

public static int comvalqri(double[][] a,
                            int n,
                            double[] em,
                            double[] re,
                            double[] im)

comveches

public static void comveches(double[][] a,
                             int n,
                             double lambda,
                             double mu,
                             double[] em,
                             double[] u,
                             double[] v)

reaeigval

public static int reaeigval(double[][] a,
                            int n,
                            double[] em,
                            double[] val)

reaeig1

public static int reaeig1(double[][] a,
                          int n,
                          double[] em,
                          double[] val,
                          double[][] vec)

reaeig3

public static int reaeig3(double[][] a,
                          int n,
                          double[] em,
                          double[] val,
                          double[][] vec)

comeigval

public static int comeigval(double[][] a,
                            int n,
                            double[] em,
                            double[] re,
                            double[] im)

comeig1

public static int comeig1(double[][] a,
                          int n,
                          double[] em,
                          double[] re,
                          double[] im,
                          double[][] vec)

eigvalhrm

public static void eigvalhrm(double[][] a,
                             int n,
                             int numval,
                             double[] val,
                             double[] em)

eighrm

public static void eighrm(double[][] a,
                          int n,
                          int numval,
                          double[] val,
                          double[][] vecr,
                          double[][] veci,
                          double[] em)

qrivalhrm

public static int qrivalhrm(double[][] a,
                            int n,
                            double[] val,
                            double[] em)

qrihrm

public static int qrihrm(double[][] a,
                         int n,
                         double[] val,
                         double[][] vr,
                         double[][] vi,
                         double[] em)

valqricom

public static int valqricom(double[][] a1,
                            double[][] a2,
                            double[] b,
                            int n,
                            double[] em,
                            double[] val1,
                            double[] val2)

qricom

public static int qricom(double[][] a1,
                         double[][] a2,
                         double[] b,
                         int n,
                         double[] em,
                         double[] val1,
                         double[] val2,
                         double[][] vec1,
                         double[][] vec2)

eigvalcom

public static int eigvalcom(double[][] ar,
                            double[][] ai,
                            int n,
                            double[] em,
                            double[] valr,
                            double[] vali)

eigcom

public static int eigcom(double[][] ar,
                         double[][] ai,
                         int n,
                         double[] em,
                         double[] valr,
                         double[] vali,
                         double[][] vr,
                         double[][] vi)

qzival

public static void qzival(int n,
                          double[][] a,
                          double[][] b,
                          double[] alfr,
                          double[] alfi,
                          double[] beta,
                          int[] iter,
                          double[] em)

qzi

public static void qzi(int n,
                       double[][] a,
                       double[][] b,
                       double[][] x,
                       double[] alfr,
                       double[] alfi,
                       double[] beta,
                       int[] iter,
                       double[] em)

hshdecmul

public static void hshdecmul(int n,
                             double[][] a,
                             double[][] b,
                             double dwarf)

hestgl3

public static void hestgl3(int n,
                           double[][] a,
                           double[][] b,
                           double[][] x)

hestgl2

public static void hestgl2(int n,
                           double[][] a,
                           double[][] b)

hsh2col

public static void hsh2col(int la,
                           int lb,
                           int u,
                           int i,
                           double a1,
                           double a2,
                           double[][] a,
                           double[][] b)

hsh3col

public static void hsh3col(int la,
                           int lb,
                           int u,
                           int i,
                           double a1,
                           double a2,
                           double a3,
                           double[][] a,
                           double[][] b)

hsh2row3

public static void hsh2row3(int l,
                            int ua,
                            int ub,
                            int ux,
                            int j,
                            double a1,
                            double a2,
                            double[][] a,
                            double[][] b,
                            double[][] x)

hsh2row2

public static void hsh2row2(int l,
                            int ua,
                            int ub,
                            int j,
                            double a1,
                            double a2,
                            double[][] a,
                            double[][] b)

hsh3row3

public static void hsh3row3(int l,
                            int u,
                            int ux,
                            int j,
                            double a1,
                            double a2,
                            double a3,
                            double[][] a,
                            double[][] b,
                            double[][] x)

hsh3row2

public static void hsh3row2(int l,
                            int u,
                            int j,
                            double a1,
                            double a2,
                            double a3,
                            double[][] a,
                            double[][] b)

qrisngvalbid

public static int qrisngvalbid(double[] d,
                               double[] b,
                               int n,
                               double[] em)

qrisngvaldecbid

public static int qrisngvaldecbid(double[] d,
                                  double[] b,
                                  int m,
                                  int n,
                                  double[][] u,
                                  double[][] v,
                                  double[] em)

qrisngval

public static int qrisngval(double[][] a,
                            int m,
                            int n,
                            double[] val,
                            double[] em)

qrisngvaldec

public static int qrisngvaldec(double[][] a,
                               int m,
                               int n,
                               double[] val,
                               double[][] v,
                               double[] em)

zerpol

public static int zerpol(int n,
                         double[] a,
                         double[] em,
                         double[] re,
                         double[] im,
                         double[] d)

zerpolfunction

private static boolean zerpolfunction(int n,
                                      double[] d,
                                      double[] f,
                                      double x,
                                      double y,
                                      double tol,
                                      int[] it,
                                      double[] newf)

bounds

public static void bounds(int n,
                          double[] a,
                          double[] re,
                          double[] im,
                          double rele,
                          double abse,
                          double[] recentre,
                          double[] imcentre,
                          double[] bound)

kcluster

private static void kcluster(int k,
                             int m,
                             int n,
                             double[] rc,
                             double[] re,
                             double[] im,
                             double[] recentre,
                             double[] imcentre,
                             double[] bound,
                             double[] clust)

allzerortpol

public static void allzerortpol(int n,
                                double[] b,
                                double[] c,
                                double[] zer,
                                double[] em)

lupzerortpol

public static void lupzerortpol(int n,
                                int m,
                                double[] b,
                                double[] c,
                                double[] zer,
                                double[] em)

selzerortpol

public static void selzerortpol(int n,
                                int n1,
                                int n2,
                                double[] b,
                                double[] c,
                                double[] zer,
                                double[] em)

alljaczer

public static void alljaczer(int n,
                             double alfa,
                             double beta,
                             double[] zer)

alllagzer

public static void alllagzer(int n,
                             double alfa,
                             double[] zer)

comkwd

public static void comkwd(double pr,
                          double pi,
                          double qr,
                          double qi,
                          double[] gr,
                          double[] gi,
                          double[] kr,
                          double[] ki)