//
|
// Source code recreated from a .class file by IntelliJ IDEA
|
// (powered by FernFlower decompiler)
|
//
|
|
package com.jisuan.Jama;
|
|
//import Jama.util.Maths;
|
import com.jisuan.Jama.util.Maths;
|
|
import java.io.Serializable;
|
|
public class EigenvalueDecomposition implements Serializable {
|
private int n;
|
private boolean issymmetric;
|
private double[] d;
|
private double[] e;
|
private double[][] V;
|
private double[][] H;
|
private double[] ort;
|
private transient double cdivr;
|
private transient double cdivi;
|
private static final long serialVersionUID = 1L;
|
|
private void tred2() {
|
int var1;
|
for(var1 = 0; var1 < this.n; ++var1) {
|
this.d[var1] = this.V[this.n - 1][var1];
|
}
|
|
double var2;
|
double[] var10000;
|
for(var1 = this.n - 1; var1 > 0; --var1) {
|
var2 = 0.0;
|
double var4 = 0.0;
|
|
int var6;
|
for(var6 = 0; var6 < var1; ++var6) {
|
var2 += Math.abs(this.d[var6]);
|
}
|
|
if (var2 == 0.0) {
|
this.e[var1] = this.d[var1 - 1];
|
|
for(var6 = 0; var6 < var1; ++var6) {
|
this.d[var6] = this.V[var1 - 1][var6];
|
this.V[var1][var6] = 0.0;
|
this.V[var6][var1] = 0.0;
|
}
|
} else {
|
for(var6 = 0; var6 < var1; ++var6) {
|
var10000 = this.d;
|
var10000[var6] /= var2;
|
var4 += this.d[var6] * this.d[var6];
|
}
|
|
double var15 = this.d[var1 - 1];
|
double var8 = Math.sqrt(var4);
|
if (var15 > 0.0) {
|
var8 = -var8;
|
}
|
|
this.e[var1] = var2 * var8;
|
var4 -= var15 * var8;
|
this.d[var1 - 1] = var15 - var8;
|
|
int var10;
|
for(var10 = 0; var10 < var1; ++var10) {
|
this.e[var10] = 0.0;
|
}
|
|
var10 = 0;
|
|
label156:
|
while(true) {
|
if (var10 >= var1) {
|
var15 = 0.0;
|
|
for(var10 = 0; var10 < var1; ++var10) {
|
var10000 = this.e;
|
var10000[var10] /= var4;
|
var15 += this.e[var10] * this.d[var10];
|
}
|
|
double var16 = var15 / (var4 + var4);
|
|
int var12;
|
for(var12 = 0; var12 < var1; ++var12) {
|
var10000 = this.e;
|
var10000[var12] -= var16 * this.d[var12];
|
}
|
|
var12 = 0;
|
|
while(true) {
|
if (var12 >= var1) {
|
break label156;
|
}
|
|
var15 = this.d[var12];
|
var8 = this.e[var12];
|
|
for(int var13 = var12; var13 <= var1 - 1; ++var13) {
|
var10000 = this.V[var13];
|
var10000[var12] -= var15 * this.e[var13] + var8 * this.d[var13];
|
}
|
|
this.d[var12] = this.V[var1 - 1][var12];
|
this.V[var1][var12] = 0.0;
|
++var12;
|
}
|
}
|
|
var15 = this.d[var10];
|
this.V[var10][var1] = var15;
|
var8 = this.e[var10] + this.V[var10][var10] * var15;
|
|
for(int var11 = var10 + 1; var11 <= var1 - 1; ++var11) {
|
var8 += this.V[var11][var10] * this.d[var11];
|
var10000 = this.e;
|
var10000[var11] += this.V[var11][var10] * var15;
|
}
|
|
this.e[var10] = var8;
|
++var10;
|
}
|
}
|
|
this.d[var1] = var4;
|
}
|
|
for(var1 = 0; var1 < this.n - 1; ++var1) {
|
this.V[this.n - 1][var1] = this.V[var1][var1];
|
this.V[var1][var1] = 1.0;
|
var2 = this.d[var1 + 1];
|
int var14;
|
if (var2 != 0.0) {
|
for(var14 = 0; var14 <= var1; ++var14) {
|
this.d[var14] = this.V[var14][var1 + 1] / var2;
|
}
|
|
for(var14 = 0; var14 <= var1; ++var14) {
|
double var5 = 0.0;
|
|
int var7;
|
for(var7 = 0; var7 <= var1; ++var7) {
|
var5 += this.V[var7][var1 + 1] * this.V[var7][var14];
|
}
|
|
for(var7 = 0; var7 <= var1; ++var7) {
|
var10000 = this.V[var7];
|
var10000[var14] -= var5 * this.d[var7];
|
}
|
}
|
}
|
|
for(var14 = 0; var14 <= var1; ++var14) {
|
this.V[var14][var1 + 1] = 0.0;
|
}
|
}
|
|
for(var1 = 0; var1 < this.n; ++var1) {
|
this.d[var1] = this.V[this.n - 1][var1];
|
this.V[this.n - 1][var1] = 0.0;
|
}
|
|
this.V[this.n - 1][this.n - 1] = 1.0;
|
this.e[0] = 0.0;
|
}
|
|
private void tql2() {
|
for(int var1 = 1; var1 < this.n; ++var1) {
|
this.e[var1 - 1] = this.e[var1];
|
}
|
|
this.e[this.n - 1] = 0.0;
|
double var34 = 0.0;
|
double var3 = 0.0;
|
double var5 = Math.pow(2.0, -52.0);
|
|
int var7;
|
int var8;
|
for(var7 = 0; var7 < this.n; ++var7) {
|
var3 = Math.max(var3, Math.abs(this.d[var7]) + Math.abs(this.e[var7]));
|
|
for(var8 = var7; var8 < this.n && !(Math.abs(this.e[var8]) <= var5 * var3); ++var8) {
|
}
|
|
if (var8 > var7) {
|
int var9 = 0;
|
|
do {
|
++var9;
|
double var10 = this.d[var7];
|
double var12 = (this.d[var7 + 1] - var10) / (2.0 * this.e[var7]);
|
double var14 = Maths.hypot(var12, 1.0);
|
if (var12 < 0.0) {
|
var14 = -var14;
|
}
|
|
this.d[var7] = this.e[var7] / (var12 + var14);
|
this.d[var7 + 1] = this.e[var7] * (var12 + var14);
|
double var16 = this.d[var7 + 1];
|
double var18 = var10 - this.d[var7];
|
|
for(int var20 = var7 + 2; var20 < this.n; ++var20) {
|
double[] var10000 = this.d;
|
var10000[var20] -= var18;
|
}
|
|
var34 += var18;
|
var12 = this.d[var8];
|
double var36 = 1.0;
|
double var22 = var36;
|
double var24 = var36;
|
double var26 = this.e[var7 + 1];
|
double var28 = 0.0;
|
double var30 = 0.0;
|
|
for(int var32 = var8 - 1; var32 >= var7; --var32) {
|
var24 = var22;
|
var22 = var36;
|
var30 = var28;
|
var10 = var36 * this.e[var32];
|
var18 = var36 * var12;
|
var14 = Maths.hypot(var12, this.e[var32]);
|
this.e[var32 + 1] = var28 * var14;
|
var28 = this.e[var32] / var14;
|
var36 = var12 / var14;
|
var12 = var36 * this.d[var32] - var28 * var10;
|
this.d[var32 + 1] = var18 + var28 * (var36 * var10 + var28 * this.d[var32]);
|
|
for(int var33 = 0; var33 < this.n; ++var33) {
|
var18 = this.V[var33][var32 + 1];
|
this.V[var33][var32 + 1] = var28 * this.V[var33][var32] + var36 * var18;
|
this.V[var33][var32] = var36 * this.V[var33][var32] - var28 * var18;
|
}
|
}
|
|
var12 = -var28 * var30 * var24 * var26 * this.e[var7] / var16;
|
this.e[var7] = var28 * var12;
|
this.d[var7] = var36 * var12;
|
} while(Math.abs(this.e[var7]) > var5 * var3);
|
}
|
|
this.d[var7] += var34;
|
this.e[var7] = 0.0;
|
}
|
|
for(var7 = 0; var7 < this.n - 1; ++var7) {
|
var8 = var7;
|
double var35 = this.d[var7];
|
|
int var11;
|
for(var11 = var7 + 1; var11 < this.n; ++var11) {
|
if (this.d[var11] < var35) {
|
var8 = var11;
|
var35 = this.d[var11];
|
}
|
}
|
|
if (var8 != var7) {
|
this.d[var8] = this.d[var7];
|
this.d[var7] = var35;
|
|
for(var11 = 0; var11 < this.n; ++var11) {
|
var35 = this.V[var11][var7];
|
this.V[var11][var7] = this.V[var11][var8];
|
this.V[var11][var8] = var35;
|
}
|
}
|
}
|
|
}
|
|
private void orthes() {
|
byte var1 = 0;
|
int var2 = this.n - 1;
|
|
int var3;
|
double[] var10000;
|
for(var3 = var1 + 1; var3 <= var2 - 1; ++var3) {
|
double var4 = 0.0;
|
|
for(int var6 = var3; var6 <= var2; ++var6) {
|
var4 += Math.abs(this.H[var6][var3 - 1]);
|
}
|
|
if (var4 != 0.0) {
|
double var15 = 0.0;
|
|
for(int var8 = var2; var8 >= var3; --var8) {
|
this.ort[var8] = this.H[var8][var3 - 1] / var4;
|
var15 += this.ort[var8] * this.ort[var8];
|
}
|
|
double var16 = Math.sqrt(var15);
|
if (this.ort[var3] > 0.0) {
|
var16 = -var16;
|
}
|
|
var15 -= this.ort[var3] * var16;
|
this.ort[var3] -= var16;
|
|
int var10;
|
double var11;
|
int var13;
|
for(var10 = var3; var10 < this.n; ++var10) {
|
var11 = 0.0;
|
|
for(var13 = var2; var13 >= var3; --var13) {
|
var11 += this.ort[var13] * this.H[var13][var10];
|
}
|
|
var11 /= var15;
|
|
for(var13 = var3; var13 <= var2; ++var13) {
|
var10000 = this.H[var13];
|
var10000[var10] -= var11 * this.ort[var13];
|
}
|
}
|
|
for(var10 = 0; var10 <= var2; ++var10) {
|
var11 = 0.0;
|
|
for(var13 = var2; var13 >= var3; --var13) {
|
var11 += this.ort[var13] * this.H[var10][var13];
|
}
|
|
var11 /= var15;
|
|
for(var13 = var3; var13 <= var2; ++var13) {
|
var10000 = this.H[var10];
|
var10000[var13] -= var11 * this.ort[var13];
|
}
|
}
|
|
this.ort[var3] = var4 * this.ort[var3];
|
this.H[var3][var3 - 1] = var4 * var16;
|
}
|
}
|
|
int var14;
|
for(var3 = 0; var3 < this.n; ++var3) {
|
for(var14 = 0; var14 < this.n; ++var14) {
|
this.V[var3][var14] = var3 == var14 ? 1.0 : 0.0;
|
}
|
}
|
|
for(var3 = var2 - 1; var3 >= var1 + 1; --var3) {
|
if (this.H[var3][var3 - 1] != 0.0) {
|
for(var14 = var3 + 1; var14 <= var2; ++var14) {
|
this.ort[var14] = this.H[var14][var3 - 1];
|
}
|
|
for(var14 = var3; var14 <= var2; ++var14) {
|
double var5 = 0.0;
|
|
int var7;
|
for(var7 = var3; var7 <= var2; ++var7) {
|
var5 += this.ort[var7] * this.V[var7][var14];
|
}
|
|
var5 = var5 / this.ort[var3] / this.H[var3][var3 - 1];
|
|
for(var7 = var3; var7 <= var2; ++var7) {
|
var10000 = this.V[var7];
|
var10000[var14] += var5 * this.ort[var7];
|
}
|
}
|
}
|
}
|
|
}
|
|
private void cdiv(double var1, double var3, double var5, double var7) {
|
double var9;
|
double var11;
|
if (Math.abs(var5) > Math.abs(var7)) {
|
var9 = var7 / var5;
|
var11 = var5 + var9 * var7;
|
this.cdivr = (var1 + var9 * var3) / var11;
|
this.cdivi = (var3 - var9 * var1) / var11;
|
} else {
|
var9 = var5 / var7;
|
var11 = var7 + var9 * var5;
|
this.cdivr = (var9 * var1 + var3) / var11;
|
this.cdivi = (var9 * var3 - var1) / var11;
|
}
|
|
}
|
|
private void hqr2() {
|
int var1 = this.n;
|
int var2 = var1 - 1;
|
byte var3 = 0;
|
int var4 = var1 - 1;
|
double var5 = Math.pow(2.0, -52.0);
|
double var7 = 0.0;
|
double var9 = 0.0;
|
double var11 = 0.0;
|
double var13 = 0.0;
|
double var15 = 0.0;
|
double var17 = 0.0;
|
double var27 = 0.0;
|
|
int var29;
|
int var30;
|
for(var29 = 0; var29 < var1; ++var29) {
|
if (var29 < var3 | var29 > var4) {
|
this.d[var29] = this.H[var29][var29];
|
this.e[var29] = 0.0;
|
}
|
|
for(var30 = Math.max(var29 - 1, 0); var30 < var1; ++var30) {
|
var27 += Math.abs(this.H[var29][var30]);
|
}
|
}
|
|
var29 = 0;
|
|
while(true) {
|
double var21;
|
double var23;
|
double var25;
|
int var31;
|
int var32;
|
while(var2 >= var3) {
|
for(var30 = var2; var30 > var3; --var30) {
|
var15 = Math.abs(this.H[var30 - 1][var30 - 1]) + Math.abs(this.H[var30][var30]);
|
if (var15 == 0.0) {
|
var15 = var27;
|
}
|
|
if (Math.abs(this.H[var30][var30 - 1]) < var5 * var15) {
|
break;
|
}
|
}
|
|
if (var30 == var2) {
|
this.H[var2][var2] += var7;
|
this.d[var2] = this.H[var2][var2];
|
this.e[var2] = 0.0;
|
--var2;
|
var29 = 0;
|
} else if (var30 == var2 - 1) {
|
var21 = this.H[var2][var2 - 1] * this.H[var2 - 1][var2];
|
var9 = (this.H[var2 - 1][var2 - 1] - this.H[var2][var2]) / 2.0;
|
var11 = var9 * var9 + var21;
|
var17 = Math.sqrt(Math.abs(var11));
|
this.H[var2][var2] += var7;
|
this.H[var2 - 1][var2 - 1] += var7;
|
var23 = this.H[var2][var2];
|
if (var11 >= 0.0) {
|
if (var9 >= 0.0) {
|
var17 += var9;
|
} else {
|
var17 = var9 - var17;
|
}
|
|
this.d[var2 - 1] = var23 + var17;
|
this.d[var2] = this.d[var2 - 1];
|
if (var17 != 0.0) {
|
this.d[var2] = var23 - var21 / var17;
|
}
|
|
this.e[var2 - 1] = 0.0;
|
this.e[var2] = 0.0;
|
var23 = this.H[var2][var2 - 1];
|
var15 = Math.abs(var23) + Math.abs(var17);
|
var9 = var23 / var15;
|
var11 = var17 / var15;
|
var13 = Math.sqrt(var9 * var9 + var11 * var11);
|
var9 /= var13;
|
var11 /= var13;
|
|
for(var31 = var2 - 1; var31 < var1; ++var31) {
|
var17 = this.H[var2 - 1][var31];
|
this.H[var2 - 1][var31] = var11 * var17 + var9 * this.H[var2][var31];
|
this.H[var2][var31] = var11 * this.H[var2][var31] - var9 * var17;
|
}
|
|
for(var31 = 0; var31 <= var2; ++var31) {
|
var17 = this.H[var31][var2 - 1];
|
this.H[var31][var2 - 1] = var11 * var17 + var9 * this.H[var31][var2];
|
this.H[var31][var2] = var11 * this.H[var31][var2] - var9 * var17;
|
}
|
|
for(var31 = var3; var31 <= var4; ++var31) {
|
var17 = this.V[var31][var2 - 1];
|
this.V[var31][var2 - 1] = var11 * var17 + var9 * this.V[var31][var2];
|
this.V[var31][var2] = var11 * this.V[var31][var2] - var9 * var17;
|
}
|
} else {
|
this.d[var2 - 1] = var23 + var9;
|
this.d[var2] = var23 + var9;
|
this.e[var2 - 1] = var17;
|
this.e[var2] = -var17;
|
}
|
|
var2 -= 2;
|
var29 = 0;
|
} else {
|
var23 = this.H[var2][var2];
|
var25 = 0.0;
|
var21 = 0.0;
|
if (var30 < var2) {
|
var25 = this.H[var2 - 1][var2 - 1];
|
var21 = this.H[var2][var2 - 1] * this.H[var2 - 1][var2];
|
}
|
|
double[] var10000;
|
if (var29 == 10) {
|
var7 += var23;
|
|
for(var31 = var3; var31 <= var2; ++var31) {
|
var10000 = this.H[var31];
|
var10000[var31] -= var23;
|
}
|
|
var15 = Math.abs(this.H[var2][var2 - 1]) + Math.abs(this.H[var2 - 1][var2 - 2]);
|
var23 = var25 = 0.75 * var15;
|
var21 = -0.4375 * var15 * var15;
|
}
|
|
if (var29 == 30) {
|
var15 = (var25 - var23) / 2.0;
|
var15 = var15 * var15 + var21;
|
if (var15 > 0.0) {
|
var15 = Math.sqrt(var15);
|
if (var25 < var23) {
|
var15 = -var15;
|
}
|
|
var15 = var23 - var21 / ((var25 - var23) / 2.0 + var15);
|
|
for(var31 = var3; var31 <= var2; ++var31) {
|
var10000 = this.H[var31];
|
var10000[var31] -= var15;
|
}
|
|
var7 += var15;
|
var21 = 0.964;
|
var25 = 0.964;
|
var23 = 0.964;
|
}
|
}
|
|
++var29;
|
|
for(var31 = var2 - 2; var31 >= var30; --var31) {
|
var17 = this.H[var31][var31];
|
var13 = var23 - var17;
|
var15 = var25 - var17;
|
var9 = (var13 * var15 - var21) / this.H[var31 + 1][var31] + this.H[var31][var31 + 1];
|
var11 = this.H[var31 + 1][var31 + 1] - var17 - var13 - var15;
|
var13 = this.H[var31 + 2][var31 + 1];
|
var15 = Math.abs(var9) + Math.abs(var11) + Math.abs(var13);
|
var9 /= var15;
|
var11 /= var15;
|
var13 /= var15;
|
if (var31 == var30 || Math.abs(this.H[var31][var31 - 1]) * (Math.abs(var11) + Math.abs(var13)) < var5 * Math.abs(var9) * (Math.abs(this.H[var31 - 1][var31 - 1]) + Math.abs(var17) + Math.abs(this.H[var31 + 1][var31 + 1]))) {
|
break;
|
}
|
}
|
|
for(var32 = var31 + 2; var32 <= var2; ++var32) {
|
this.H[var32][var32 - 2] = 0.0;
|
if (var32 > var31 + 2) {
|
this.H[var32][var32 - 3] = 0.0;
|
}
|
}
|
|
for(var32 = var31; var32 <= var2 - 1; ++var32) {
|
boolean var33 = var32 != var2 - 1;
|
if (var32 != var31) {
|
var9 = this.H[var32][var32 - 1];
|
var11 = this.H[var32 + 1][var32 - 1];
|
var13 = var33 ? this.H[var32 + 2][var32 - 1] : 0.0;
|
var23 = Math.abs(var9) + Math.abs(var11) + Math.abs(var13);
|
if (var23 == 0.0) {
|
continue;
|
}
|
|
var9 /= var23;
|
var11 /= var23;
|
var13 /= var23;
|
}
|
|
var15 = Math.sqrt(var9 * var9 + var11 * var11 + var13 * var13);
|
if (var9 < 0.0) {
|
var15 = -var15;
|
}
|
|
if (var15 != 0.0) {
|
if (var32 != var31) {
|
this.H[var32][var32 - 1] = -var15 * var23;
|
} else if (var30 != var31) {
|
this.H[var32][var32 - 1] = -this.H[var32][var32 - 1];
|
}
|
|
var9 += var15;
|
var23 = var9 / var15;
|
var25 = var11 / var15;
|
var17 = var13 / var15;
|
var11 /= var9;
|
var13 /= var9;
|
|
int var34;
|
for(var34 = var32; var34 < var1; ++var34) {
|
var9 = this.H[var32][var34] + var11 * this.H[var32 + 1][var34];
|
if (var33) {
|
var9 += var13 * this.H[var32 + 2][var34];
|
this.H[var32 + 2][var34] -= var9 * var17;
|
}
|
|
this.H[var32][var34] -= var9 * var23;
|
this.H[var32 + 1][var34] -= var9 * var25;
|
}
|
|
for(var34 = 0; var34 <= Math.min(var2, var32 + 3); ++var34) {
|
var9 = var23 * this.H[var34][var32] + var25 * this.H[var34][var32 + 1];
|
if (var33) {
|
var9 += var17 * this.H[var34][var32 + 2];
|
this.H[var34][var32 + 2] -= var9 * var13;
|
}
|
|
this.H[var34][var32] -= var9;
|
this.H[var34][var32 + 1] -= var9 * var11;
|
}
|
|
for(var34 = var3; var34 <= var4; ++var34) {
|
var9 = var23 * this.V[var34][var32] + var25 * this.V[var34][var32 + 1];
|
if (var33) {
|
var9 += var17 * this.V[var34][var32 + 2];
|
this.V[var34][var32 + 2] -= var9 * var13;
|
}
|
|
this.V[var34][var32] -= var9;
|
this.V[var34][var32 + 1] -= var9 * var11;
|
}
|
}
|
}
|
}
|
}
|
|
if (var27 == 0.0) {
|
return;
|
}
|
|
for(var2 = var1 - 1; var2 >= 0; --var2) {
|
var9 = this.d[var2];
|
var11 = this.e[var2];
|
double var19;
|
if (var11 == 0.0) {
|
var30 = var2;
|
this.H[var2][var2] = 1.0;
|
|
for(var31 = var2 - 1; var31 >= 0; --var31) {
|
var21 = this.H[var31][var31] - var9;
|
var13 = 0.0;
|
|
for(var32 = var30; var32 <= var2; ++var32) {
|
var13 += this.H[var31][var32] * this.H[var32][var2];
|
}
|
|
if (this.e[var31] < 0.0) {
|
var17 = var21;
|
var15 = var13;
|
} else {
|
var30 = var31;
|
if (this.e[var31] == 0.0) {
|
if (var21 != 0.0) {
|
this.H[var31][var2] = -var13 / var21;
|
} else {
|
this.H[var31][var2] = -var13 / (var5 * var27);
|
}
|
} else {
|
var23 = this.H[var31][var31 + 1];
|
var25 = this.H[var31 + 1][var31];
|
var11 = (this.d[var31] - var9) * (this.d[var31] - var9) + this.e[var31] * this.e[var31];
|
var19 = (var23 * var15 - var17 * var13) / var11;
|
this.H[var31][var2] = var19;
|
if (Math.abs(var23) > Math.abs(var17)) {
|
this.H[var31 + 1][var2] = (-var13 - var21 * var19) / var23;
|
} else {
|
this.H[var31 + 1][var2] = (-var15 - var25 * var19) / var17;
|
}
|
}
|
|
var19 = Math.abs(this.H[var31][var2]);
|
if (var5 * var19 * var19 > 1.0) {
|
for(var32 = var31; var32 <= var2; ++var32) {
|
this.H[var32][var2] /= var19;
|
}
|
}
|
}
|
}
|
} else if (var11 < 0.0) {
|
var30 = var2 - 1;
|
if (Math.abs(this.H[var2][var2 - 1]) > Math.abs(this.H[var2 - 1][var2])) {
|
this.H[var2 - 1][var2 - 1] = var11 / this.H[var2][var2 - 1];
|
this.H[var2 - 1][var2] = -(this.H[var2][var2] - var9) / this.H[var2][var2 - 1];
|
} else {
|
this.cdiv(0.0, -this.H[var2 - 1][var2], this.H[var2 - 1][var2 - 1] - var9, var11);
|
this.H[var2 - 1][var2 - 1] = this.cdivr;
|
this.H[var2 - 1][var2] = this.cdivi;
|
}
|
|
this.H[var2][var2 - 1] = 0.0;
|
this.H[var2][var2] = 1.0;
|
|
for(var31 = var2 - 2; var31 >= 0; --var31) {
|
double var41 = 0.0;
|
double var42 = 0.0;
|
|
int var40;
|
for(var40 = var30; var40 <= var2; ++var40) {
|
var41 += this.H[var31][var40] * this.H[var40][var2 - 1];
|
var42 += this.H[var31][var40] * this.H[var40][var2];
|
}
|
|
var21 = this.H[var31][var31] - var9;
|
if (this.e[var31] < 0.0) {
|
var17 = var21;
|
var13 = var41;
|
var15 = var42;
|
} else {
|
var30 = var31;
|
if (this.e[var31] == 0.0) {
|
this.cdiv(-var41, -var42, var21, var11);
|
this.H[var31][var2 - 1] = this.cdivr;
|
this.H[var31][var2] = this.cdivi;
|
} else {
|
var23 = this.H[var31][var31 + 1];
|
var25 = this.H[var31 + 1][var31];
|
double var36 = (this.d[var31] - var9) * (this.d[var31] - var9) + this.e[var31] * this.e[var31] - var11 * var11;
|
double var38 = (this.d[var31] - var9) * 2.0 * var11;
|
if (var36 == 0.0 & var38 == 0.0) {
|
var36 = var5 * var27 * (Math.abs(var21) + Math.abs(var11) + Math.abs(var23) + Math.abs(var25) + Math.abs(var17));
|
}
|
|
this.cdiv(var23 * var13 - var17 * var41 + var11 * var42, var23 * var15 - var17 * var42 - var11 * var41, var36, var38);
|
this.H[var31][var2 - 1] = this.cdivr;
|
this.H[var31][var2] = this.cdivi;
|
if (Math.abs(var23) > Math.abs(var17) + Math.abs(var11)) {
|
this.H[var31 + 1][var2 - 1] = (-var41 - var21 * this.H[var31][var2 - 1] + var11 * this.H[var31][var2]) / var23;
|
this.H[var31 + 1][var2] = (-var42 - var21 * this.H[var31][var2] - var11 * this.H[var31][var2 - 1]) / var23;
|
} else {
|
this.cdiv(-var13 - var25 * this.H[var31][var2 - 1], -var15 - var25 * this.H[var31][var2], var17, var11);
|
this.H[var31 + 1][var2 - 1] = this.cdivr;
|
this.H[var31 + 1][var2] = this.cdivi;
|
}
|
}
|
|
var19 = Math.max(Math.abs(this.H[var31][var2 - 1]), Math.abs(this.H[var31][var2]));
|
if (var5 * var19 * var19 > 1.0) {
|
for(var40 = var31; var40 <= var2; ++var40) {
|
this.H[var40][var2 - 1] /= var19;
|
this.H[var40][var2] /= var19;
|
}
|
}
|
}
|
}
|
}
|
}
|
|
for(var30 = 0; var30 < var1; ++var30) {
|
if (var30 < var3 | var30 > var4) {
|
for(var31 = var30; var31 < var1; ++var31) {
|
this.V[var30][var31] = this.H[var30][var31];
|
}
|
}
|
}
|
|
for(var30 = var1 - 1; var30 >= var3; --var30) {
|
for(var31 = var3; var31 <= var4; ++var31) {
|
var17 = 0.0;
|
|
for(var32 = var3; var32 <= Math.min(var30, var4); ++var32) {
|
var17 += this.V[var31][var32] * this.H[var32][var30];
|
}
|
|
this.V[var31][var30] = var17;
|
}
|
}
|
|
return;
|
}
|
}
|
|
public EigenvalueDecomposition(Matrix var1) {
|
double[][] var2 = var1.getArray();
|
this.n = var1.getColumnDimension();
|
this.V = new double[this.n][this.n];
|
this.d = new double[this.n];
|
this.e = new double[this.n];
|
this.issymmetric = true;
|
|
int var3;
|
int var4;
|
for(var3 = 0; var3 < this.n & this.issymmetric; ++var3) {
|
for(var4 = 0; var4 < this.n & this.issymmetric; ++var4) {
|
this.issymmetric = var2[var4][var3] == var2[var3][var4];
|
}
|
}
|
|
if (this.issymmetric) {
|
for(var3 = 0; var3 < this.n; ++var3) {
|
for(var4 = 0; var4 < this.n; ++var4) {
|
this.V[var3][var4] = var2[var3][var4];
|
}
|
}
|
|
this.tred2();
|
this.tql2();
|
} else {
|
this.H = new double[this.n][this.n];
|
this.ort = new double[this.n];
|
|
for(var3 = 0; var3 < this.n; ++var3) {
|
for(var4 = 0; var4 < this.n; ++var4) {
|
this.H[var4][var3] = var2[var4][var3];
|
}
|
}
|
|
this.orthes();
|
this.hqr2();
|
}
|
|
}
|
|
public Matrix getV() {
|
return new Matrix(this.V, this.n, this.n);
|
}
|
|
public double[] getRealEigenvalues() {
|
return this.d;
|
}
|
|
public double[] getImagEigenvalues() {
|
return this.e;
|
}
|
|
public Matrix getD() {
|
Matrix var1 = new Matrix(this.n, this.n);
|
double[][] var2 = var1.getArray();
|
|
for(int var3 = 0; var3 < this.n; ++var3) {
|
for(int var4 = 0; var4 < this.n; ++var4) {
|
var2[var3][var4] = 0.0;
|
}
|
|
var2[var3][var3] = this.d[var3];
|
if (this.e[var3] > 0.0) {
|
var2[var3][var3 + 1] = this.e[var3];
|
} else if (this.e[var3] < 0.0) {
|
var2[var3][var3 - 1] = this.e[var3];
|
}
|
}
|
|
return var1;
|
}
|
}
|
