public static void dlaein(
boolean paramBoolean1,
boolean paramBoolean2,
int paramInt1,
double[] paramArrayOfDouble1,
int paramInt2,
int paramInt3,
double paramDouble1,
double paramDouble2,
double[] paramArrayOfDouble2,
int paramInt4,
double[] paramArrayOfDouble3,
int paramInt5,
double[] paramArrayOfDouble4,
int paramInt6,
int paramInt7,
double[] paramArrayOfDouble5,
int paramInt8,
double paramDouble3,
double paramDouble4,
double paramDouble5,
intW paramintW) {
String str1 = new String(" ");
String str2 = new String(" ");
int i = 0;
int j = 0;
int k = 0;
int m = 0;
intW localintW = new intW(0);
int n = 0;
int i1 = 0;
double d1 = 0.0D;
double d2 = 0.0D;
double d3 = 0.0D;
double d4 = 0.0D;
double d5 = 0.0D;
double d6 = 0.0D;
double d7 = 0.0D;
double d8 = 0.0D;
double d9 = 0.0D;
doubleW localdoubleW = new doubleW(0.0D);
double d10 = 0.0D;
double d11 = 0.0D;
double d12 = 0.0D;
double d13 = 0.0D;
double d14 = 0.0D;
double d15 = 0.0D;
double d16 = 0.0D;
double d17 = 0.0D;
double d18 = 0.0D;
double d19 = 0.0D;
paramintW.val = 0;
d9 = Math.sqrt(paramInt1);
d5 = 0.1D / d9;
d7 = Math.max(1.0D, paramDouble3 * d9) * paramDouble4;
i1 = 1;
int i3;
for (int i2 = paramInt1 - 1 + 1; i2 > 0; i2--) {
i = 1;
for (i3 = i1 - 1 - 1 + 1; i3 > 0; i3--) {
paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] =
paramArrayOfDouble1[(i - 1 + (i1 - 1) * paramInt3 + paramInt2)];
i += 1;
}
paramArrayOfDouble1[(i1 - 1 + (i1 - 1) * paramInt3 + paramInt2)] -= paramDouble1;
i1 += 1;
}
if ((paramDouble2 != 0.0D ? 0 : 1) != 0) {
if (paramBoolean2) {
i = 1;
for (i2 = paramInt1 - 1 + 1; i2 > 0; i2--) {
paramArrayOfDouble2[(i - 1 + paramInt4)] = paramDouble3;
i += 1;
}
} else {
d13 = Dnrm2.dnrm2(paramInt1, paramArrayOfDouble2, paramInt4, 1);
Dscal.dscal(
paramInt1, paramDouble3 * d9 / Math.max(d13, d7), paramArrayOfDouble2, paramInt4, 1);
}
if (paramBoolean1) {
i = 1;
for (i2 = paramInt1 - 1 - 1 + 1; i2 > 0; i2--) {
d3 = paramArrayOfDouble1[(i + 1 - 1 + (i - 1) * paramInt3 + paramInt2)];
if ((Math.abs(paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)])
>= Math.abs(d3)
? 0
: 1)
!= 0) {
d16 = paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] / d3;
paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] = d3;
i1 = i + 1;
for (i3 = paramInt1 - (i + 1) + 1; i3 > 0; i3--) {
d10 = paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)];
paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] - d16 * d10);
paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] = d10;
i1 += 1;
}
} else {
if ((paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] != 0.0D ? 0 : 1)
!= 0) {
paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] = paramDouble3;
}
d16 = d3 / paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)];
if ((d16 == 0.0D ? 0 : 1) != 0) {
i1 = i + 1;
for (i3 = paramInt1 - (i + 1) + 1; i3 > 0; i3--) {
paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] -=
d16 * paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)];
i1 += 1;
}
}
}
i += 1;
}
if ((paramArrayOfDouble4[(paramInt1 - 1 + (paramInt1 - 1) * paramInt7 + paramInt6)] != 0.0D
? 0
: 1)
!= 0) {
paramArrayOfDouble4[(paramInt1 - 1 + (paramInt1 - 1) * paramInt7 + paramInt6)] =
paramDouble3;
}
str2 = "N";
} else {
i1 = paramInt1;
for (i2 = (2 - paramInt1 + -1) / -1; i2 > 0; i2--) {
d4 = paramArrayOfDouble1[(i1 - 1 + (i1 - 1 - 1) * paramInt3 + paramInt2)];
if ((Math.abs(paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)])
>= Math.abs(d4)
? 0
: 1)
!= 0) {
d16 = paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] / d4;
paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] = d4;
i = 1;
for (i3 = i1 - 1 - 1 + 1; i3 > 0; i3--) {
d10 = paramArrayOfDouble4[(i - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)];
paramArrayOfDouble4[(i - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] - d16 * d10);
paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] = d10;
i += 1;
}
} else {
if ((paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] != 0.0D ? 0 : 1)
!= 0) {
paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] = paramDouble3;
}
d16 = d4 / paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)];
if ((d16 == 0.0D ? 0 : 1) != 0) {
i = 1;
for (i3 = i1 - 1 - 1 + 1; i3 > 0; i3--) {
paramArrayOfDouble4[(i - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] -=
d16 * paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)];
i += 1;
}
}
}
i1 += -1;
}
if ((paramArrayOfDouble4[(1 - 1 + (1 - 1) * paramInt7 + paramInt6)] != 0.0D ? 0 : 1) != 0) {
paramArrayOfDouble4[(1 - 1 + (1 - 1) * paramInt7 + paramInt6)] = paramDouble3;
}
str2 = "T";
}
str1 = "N";
n = 1;
for (i2 = paramInt1 - 1 + 1; i2 > 0; i2--) {
Dlatrs.dlatrs(
"Upper",
str2,
"Nonunit",
str1,
paramInt1,
paramArrayOfDouble4,
paramInt6,
paramInt7,
paramArrayOfDouble2,
paramInt4,
localdoubleW,
paramArrayOfDouble5,
paramInt8,
localintW);
str1 = "Y";
d13 = Dasum.dasum(paramInt1, paramArrayOfDouble2, paramInt4, 1);
if ((d13 < d5 * localdoubleW.val ? 0 : 1) != 0) {
break;
}
d10 = paramDouble3 / (d9 + 1.0D);
paramArrayOfDouble2[(1 - 1 + paramInt4)] = paramDouble3;
i = 2;
for (i3 = paramInt1 - 2 + 1; i3 > 0; i3--) {
paramArrayOfDouble2[(i - 1 + paramInt4)] = d10;
i += 1;
}
paramArrayOfDouble2[(paramInt1 - n + 1 - 1 + paramInt4)] -= paramDouble3 * d9;
n += 1;
}
paramintW.val = 1;
i = Idamax.idamax(paramInt1, paramArrayOfDouble2, paramInt4, 1);
Dscal.dscal(
paramInt1,
1.0D / Math.abs(paramArrayOfDouble2[(i - 1 + paramInt4)]),
paramArrayOfDouble2,
paramInt4,
1);
} else {
if (paramBoolean2) {
i = 1;
for (i2 = paramInt1 - 1 + 1; i2 > 0; i2--) {
paramArrayOfDouble2[(i - 1 + paramInt4)] = paramDouble3;
paramArrayOfDouble3[(i - 1 + paramInt5)] = 0.0D;
i += 1;
}
} else {
d6 =
Dlapy2.dlapy2(
Dnrm2.dnrm2(paramInt1, paramArrayOfDouble2, paramInt4, 1),
Dnrm2.dnrm2(paramInt1, paramArrayOfDouble3, paramInt5, 1));
d8 = paramDouble3 * d9 / Math.max(d6, d7);
Dscal.dscal(paramInt1, d8, paramArrayOfDouble2, paramInt4, 1);
Dscal.dscal(paramInt1, d8, paramArrayOfDouble3, paramInt5, 1);
}
if (paramBoolean1) {
paramArrayOfDouble4[(2 - 1 + (1 - 1) * paramInt7 + paramInt6)] = (-paramDouble2);
i = 2;
for (i2 = paramInt1 - 2 + 1; i2 > 0; i2--) {
paramArrayOfDouble4[(i + 1 - 1 + (1 - 1) * paramInt7 + paramInt6)] = 0.0D;
i += 1;
}
i = 1;
for (i2 = paramInt1 - 1 - 1 + 1; i2 > 0; i2--) {
d1 =
Dlapy2.dlapy2(
paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)],
paramArrayOfDouble4[(i + 1 - 1 + (i - 1) * paramInt7 + paramInt6)]);
d3 = paramArrayOfDouble1[(i + 1 - 1 + (i - 1) * paramInt3 + paramInt2)];
if ((d1 >= Math.abs(d3) ? 0 : 1) != 0) {
d18 = paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] / d3;
d17 = paramArrayOfDouble4[(i + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] / d3;
paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] = d3;
paramArrayOfDouble4[(i + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] = 0.0D;
i1 = i + 1;
for (i3 = paramInt1 - (i + 1) + 1; i3 > 0; i3--) {
d10 = paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)];
paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] - d18 * d10);
paramArrayOfDouble4[(i1 + 1 - 1 + (i + 1 - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] - d17 * d10);
paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] = d10;
paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] = 0.0D;
i1 += 1;
}
paramArrayOfDouble4[(i + 2 - 1 + (i - 1) * paramInt7 + paramInt6)] = (-paramDouble2);
paramArrayOfDouble4[(i + 1 - 1 + (i + 1 - 1) * paramInt7 + paramInt6)] -=
d17 * paramDouble2;
paramArrayOfDouble4[(i + 2 - 1 + (i + 1 - 1) * paramInt7 + paramInt6)] +=
d18 * paramDouble2;
} else {
if ((d1 != 0.0D ? 0 : 1) != 0) {
paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] = paramDouble3;
paramArrayOfDouble4[(i + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] = 0.0D;
d1 = paramDouble3;
}
d3 = d3 / d1 / d1;
d18 = paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)] * d3;
d17 = -(paramArrayOfDouble4[(i + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] * d3);
i1 = i + 1;
for (i3 = paramInt1 - (i + 1) + 1; i3 > 0; i3--) {
paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)]
- d18 * paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)]
+ d17 * paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)]);
paramArrayOfDouble4[(i1 + 1 - 1 + (i + 1 - 1) * paramInt7 + paramInt6)] =
(-(d18 * paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)])
- d17 * paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)]);
i1 += 1;
}
paramArrayOfDouble4[(i + 2 - 1 + (i + 1 - 1) * paramInt7 + paramInt6)] -= paramDouble2;
}
paramArrayOfDouble5[(i - 1 + paramInt8)] =
(Dasum.dasum(
paramInt1 - i,
paramArrayOfDouble4,
i - 1 + (i + 1 - 1) * paramInt7 + paramInt6,
paramInt7)
+ Dasum.dasum(
paramInt1 - i,
paramArrayOfDouble4,
i + 2 - 1 + (i - 1) * paramInt7 + paramInt6,
1));
i += 1;
}
if ((paramArrayOfDouble4[(paramInt1 - 1 + (paramInt1 - 1) * paramInt7 + paramInt6)] != 0.0D
? 0
: 1)
!= 0) {}
if (((paramArrayOfDouble4[(paramInt1 + 1 - 1 + (paramInt1 - 1) * paramInt7 + paramInt6)]
!= 0.0D
? 0
: 1)
!= 0
? 1
: 0)
!= 0) {
paramArrayOfDouble4[(paramInt1 - 1 + (paramInt1 - 1) * paramInt7 + paramInt6)] =
paramDouble3;
}
paramArrayOfDouble5[(paramInt1 - 1 + paramInt8)] = 0.0D;
j = paramInt1;
k = 1;
m = -1;
} else {
paramArrayOfDouble4[(paramInt1 + 1 - 1 + (paramInt1 - 1) * paramInt7 + paramInt6)] =
paramDouble2;
i1 = 1;
for (i2 = paramInt1 - 1 - 1 + 1; i2 > 0; i2--) {
paramArrayOfDouble4[(paramInt1 + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] = 0.0D;
i1 += 1;
}
i1 = paramInt1;
for (i2 = (2 - paramInt1 + -1) / -1; i2 > 0; i2--) {
d4 = paramArrayOfDouble1[(i1 - 1 + (i1 - 1 - 1) * paramInt3 + paramInt2)];
d2 =
Dlapy2.dlapy2(
paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)],
paramArrayOfDouble4[(i1 + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)]);
if ((d2 >= Math.abs(d4) ? 0 : 1) != 0) {
d18 = paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] / d4;
d17 = paramArrayOfDouble4[(i1 + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] / d4;
paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] = d4;
paramArrayOfDouble4[(i1 + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] = 0.0D;
i = 1;
for (i3 = i1 - 1 - 1 + 1; i3 > 0; i3--) {
d10 = paramArrayOfDouble4[(i - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)];
paramArrayOfDouble4[(i - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] - d18 * d10);
paramArrayOfDouble4[(i1 - 1 + (i - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] - d17 * d10);
paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)] = d10;
paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] = 0.0D;
i += 1;
}
paramArrayOfDouble4[(i1 + 1 - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] = paramDouble2;
paramArrayOfDouble4[(i1 - 1 - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] +=
d17 * paramDouble2;
paramArrayOfDouble4[(i1 - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] -=
d18 * paramDouble2;
} else {
if ((d2 != 0.0D ? 0 : 1) != 0) {
paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] = paramDouble3;
paramArrayOfDouble4[(i1 + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] = 0.0D;
d2 = paramDouble3;
}
d4 = d4 / d2 / d2;
d18 = paramArrayOfDouble4[(i1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] * d4;
d17 = -(paramArrayOfDouble4[(i1 + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)] * d4);
i = 1;
for (i3 = i1 - 1 - 1 + 1; i3 > 0; i3--) {
paramArrayOfDouble4[(i - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] =
(paramArrayOfDouble4[(i - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)]
- d18 * paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)]
+ d17 * paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)]);
paramArrayOfDouble4[(i1 - 1 + (i - 1) * paramInt7 + paramInt6)] =
(-(d18 * paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)])
- d17 * paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)]);
i += 1;
}
paramArrayOfDouble4[(i1 - 1 + (i1 - 1 - 1) * paramInt7 + paramInt6)] += paramDouble2;
}
paramArrayOfDouble5[(i1 - 1 + paramInt8)] =
(Dasum.dasum(i1 - 1, paramArrayOfDouble4, 1 - 1 + (i1 - 1) * paramInt7 + paramInt6, 1)
+ Dasum.dasum(
i1 - 1,
paramArrayOfDouble4,
i1 + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7));
i1 += -1;
}
if ((paramArrayOfDouble4[(1 - 1 + (1 - 1) * paramInt7 + paramInt6)] != 0.0D ? 0 : 1)
!= 0) {}
if (((paramArrayOfDouble4[(2 - 1 + (1 - 1) * paramInt7 + paramInt6)] != 0.0D ? 0 : 1) != 0
? 1
: 0)
!= 0) {
paramArrayOfDouble4[(1 - 1 + (1 - 1) * paramInt7 + paramInt6)] = paramDouble3;
}
paramArrayOfDouble5[(1 - 1 + paramInt8)] = 0.0D;
j = 1;
k = paramInt1;
m = 1;
}
n = 1;
for (i2 = paramInt1 - 1 + 1; i2 > 0; i2--) {
localdoubleW.val = 1.0D;
d12 = 1.0D;
d11 = paramDouble5;
i = j;
for (i3 = (k - j + m) / m; i3 > 0; i3--) {
if ((paramArrayOfDouble5[(i - 1 + paramInt8)] <= d11 ? 0 : 1) != 0) {
d8 = 1.0D / d12;
Dscal.dscal(paramInt1, d8, paramArrayOfDouble2, paramInt4, 1);
Dscal.dscal(paramInt1, d8, paramArrayOfDouble3, paramInt5, 1);
localdoubleW.val *= d8;
d12 = 1.0D;
d11 = paramDouble5;
}
d18 = paramArrayOfDouble2[(i - 1 + paramInt4)];
d17 = paramArrayOfDouble3[(i - 1 + paramInt5)];
int i4;
if (paramBoolean1) {
i1 = i + 1;
for (i4 = paramInt1 - (i + 1) + 1; i4 > 0; i4--) {
d18 =
d18
- paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble2[(i1 - 1 + paramInt4)]
+ paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble3[(i1 - 1 + paramInt5)];
d17 =
d17
- paramArrayOfDouble4[(i - 1 + (i1 - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble3[(i1 - 1 + paramInt5)]
- paramArrayOfDouble4[(i1 + 1 - 1 + (i - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble2[(i1 - 1 + paramInt4)];
i1 += 1;
}
} else {
i1 = 1;
for (i4 = i - 1 - 1 + 1; i4 > 0; i4--) {
d18 =
d18
- paramArrayOfDouble4[(i1 - 1 + (i - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble2[(i1 - 1 + paramInt4)]
+ paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble3[(i1 - 1 + paramInt5)];
d17 =
d17
- paramArrayOfDouble4[(i1 - 1 + (i - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble3[(i1 - 1 + paramInt5)]
- paramArrayOfDouble4[(i + 1 - 1 + (i1 - 1) * paramInt7 + paramInt6)]
* paramArrayOfDouble2[(i1 - 1 + paramInt4)];
i1 += 1;
}
}
d14 =
Math.abs(paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)])
+ Math.abs(paramArrayOfDouble4[(i + 1 - 1 + (i - 1) * paramInt7 + paramInt6)]);
if ((d14 <= paramDouble4 ? 0 : 1) != 0) {
if ((d14 >= 1.0D ? 0 : 1) != 0) {
d15 = Math.abs(d18) + Math.abs(d17);
if ((d15 <= d14 * paramDouble5 ? 0 : 1) != 0) {
d8 = 1.0D / d15;
Dscal.dscal(paramInt1, d8, paramArrayOfDouble2, paramInt4, 1);
Dscal.dscal(paramInt1, d8, paramArrayOfDouble3, paramInt5, 1);
d18 = paramArrayOfDouble2[(i - 1 + paramInt4)];
d17 = paramArrayOfDouble3[(i - 1 + paramInt5)];
localdoubleW.val *= d8;
d12 *= d8;
}
}
dladiv_adapter(
d18,
d17,
paramArrayOfDouble4[(i - 1 + (i - 1) * paramInt7 + paramInt6)],
paramArrayOfDouble4[(i + 1 - 1 + (i - 1) * paramInt7 + paramInt6)],
paramArrayOfDouble2,
i - 1 + paramInt4,
paramArrayOfDouble3,
i - 1 + paramInt5);
d12 =
Math.max(
Math.abs(paramArrayOfDouble2[(i - 1 + paramInt4)])
+ Math.abs(paramArrayOfDouble3[(i - 1 + paramInt5)]),
d12);
d11 = paramDouble5 / d12;
} else {
i1 = 1;
for (i4 = paramInt1 - 1 + 1; i4 > 0; i4--) {
paramArrayOfDouble2[(i1 - 1 + paramInt4)] = 0.0D;
paramArrayOfDouble3[(i1 - 1 + paramInt5)] = 0.0D;
i1 += 1;
}
paramArrayOfDouble2[(i - 1 + paramInt4)] = 1.0D;
paramArrayOfDouble3[(i - 1 + paramInt5)] = 1.0D;
localdoubleW.val = 0.0D;
d12 = 1.0D;
d11 = paramDouble5;
}
i += m;
}
d13 =
Dasum.dasum(paramInt1, paramArrayOfDouble2, paramInt4, 1)
+ Dasum.dasum(paramInt1, paramArrayOfDouble3, paramInt5, 1);
if ((d13 < d5 * localdoubleW.val ? 0 : 1) != 0) {
break;
}
d19 = paramDouble3 / (d9 + 1.0D);
paramArrayOfDouble2[(1 - 1 + paramInt4)] = paramDouble3;
paramArrayOfDouble3[(1 - 1 + paramInt5)] = 0.0D;
i = 2;
for (i3 = paramInt1 - 2 + 1; i3 > 0; i3--) {
paramArrayOfDouble2[(i - 1 + paramInt4)] = d19;
paramArrayOfDouble3[(i - 1 + paramInt5)] = 0.0D;
i += 1;
}
paramArrayOfDouble2[(paramInt1 - n + 1 - 1 + paramInt4)] -= paramDouble3 * d9;
n += 1;
}
paramintW.val = 1;
d13 = 0.0D;
i = 1;
for (i2 = paramInt1 - 1 + 1; i2 > 0; i2--) {
d13 =
Math.max(
d13,
Math.abs(paramArrayOfDouble2[(i - 1 + paramInt4)])
+ Math.abs(paramArrayOfDouble3[(i - 1 + paramInt5)]));
i += 1;
}
Dscal.dscal(paramInt1, 1.0D / d13, paramArrayOfDouble2, paramInt4, 1);
Dscal.dscal(paramInt1, 1.0D / d13, paramArrayOfDouble3, paramInt5, 1);
}
}
public static void dpbstf(
String paramString,
int paramInt1,
int paramInt2,
double[] paramArrayOfDouble,
int paramInt3,
int paramInt4,
intW paramintW) {
boolean bool = false;
int i = 0;
int j = 0;
int k = 0;
int m = 0;
double d = 0.0D;
paramintW.val = 0;
bool = Lsame.lsame(paramString, "U");
if ((((bool ^ true)) && ((Lsame.lsame(paramString, "L") ^ true)) ? 1 : 0) != 0) {
paramintW.val = -1;
} else if ((paramInt1 >= 0 ? 0 : 1) != 0) {
paramintW.val = -2;
} else if ((paramInt2 >= 0 ? 0 : 1) != 0) {
paramintW.val = -3;
} else if ((paramInt4 >= paramInt2 + 1 ? 0 : 1) != 0) {
paramintW.val = -5;
}
if ((paramintW.val == 0 ? 0 : 1) != 0) {
Xerbla.xerbla("DPBSTF", -paramintW.val);
return;
}
if ((paramInt1 != 0 ? 0 : 1) != 0) {
return;
}
j = Math.max(1, paramInt4 - 1);
m = (paramInt1 + paramInt2) / 2;
int n;
if (bool) {
i = paramInt1;
for (n = (m + 1 - paramInt1 + -1) / -1; n > 0; n--) {
d = paramArrayOfDouble[(paramInt2 + 1 - 1 + (i - 1) * paramInt4 + paramInt3)];
if ((d > 0.0D ? 0 : 1) != 0) {
break;
}
d = Math.sqrt(d);
paramArrayOfDouble[(paramInt2 + 1 - 1 + (i - 1) * paramInt4 + paramInt3)] = d;
k = Math.min(i - 1, paramInt2);
Dscal.dscal(
k,
1.0D / d,
paramArrayOfDouble,
paramInt2 + 1 - k - 1 + (i - 1) * paramInt4 + paramInt3,
1);
Dsyr.dsyr(
"Upper",
k,
-1.0D,
paramArrayOfDouble,
paramInt2 + 1 - k - 1 + (i - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble,
paramInt2 + 1 - 1 + (i - k - 1) * paramInt4 + paramInt3,
j);
i += -1;
}
i = 1;
for (n = m - 1 + 1; n > 0; n--) {
d = paramArrayOfDouble[(paramInt2 + 1 - 1 + (i - 1) * paramInt4 + paramInt3)];
if ((d > 0.0D ? 0 : 1) != 0) {
break;
}
d = Math.sqrt(d);
paramArrayOfDouble[(paramInt2 + 1 - 1 + (i - 1) * paramInt4 + paramInt3)] = d;
k = Math.min(paramInt2, m - i);
if ((k <= 0 ? 0 : 1) != 0) {
Dscal.dscal(
k,
1.0D / d,
paramArrayOfDouble,
paramInt2 - 1 + (i + 1 - 1) * paramInt4 + paramInt3,
j);
Dsyr.dsyr(
"Upper",
k,
-1.0D,
paramArrayOfDouble,
paramInt2 - 1 + (i + 1 - 1) * paramInt4 + paramInt3,
j,
paramArrayOfDouble,
paramInt2 + 1 - 1 + (i + 1 - 1) * paramInt4 + paramInt3,
j);
}
i += 1;
}
} else {
i = paramInt1;
for (n = (m + 1 - paramInt1 + -1) / -1; n > 0; n--) {
d = paramArrayOfDouble[(1 - 1 + (i - 1) * paramInt4 + paramInt3)];
if ((d > 0.0D ? 0 : 1) != 0) {
break;
}
d = Math.sqrt(d);
paramArrayOfDouble[(1 - 1 + (i - 1) * paramInt4 + paramInt3)] = d;
k = Math.min(i - 1, paramInt2);
Dscal.dscal(
k, 1.0D / d, paramArrayOfDouble, k + 1 - 1 + (i - k - 1) * paramInt4 + paramInt3, j);
Dsyr.dsyr(
"Lower",
k,
-1.0D,
paramArrayOfDouble,
k + 1 - 1 + (i - k - 1) * paramInt4 + paramInt3,
j,
paramArrayOfDouble,
1 - 1 + (i - k - 1) * paramInt4 + paramInt3,
j);
i += -1;
}
i = 1;
for (n = m - 1 + 1; n > 0; n--) {
d = paramArrayOfDouble[(1 - 1 + (i - 1) * paramInt4 + paramInt3)];
if ((d > 0.0D ? 0 : 1) != 0) {
break;
}
d = Math.sqrt(d);
paramArrayOfDouble[(1 - 1 + (i - 1) * paramInt4 + paramInt3)] = d;
k = Math.min(paramInt2, m - i);
if ((k <= 0 ? 0 : 1) != 0) {
Dscal.dscal(k, 1.0D / d, paramArrayOfDouble, 2 - 1 + (i - 1) * paramInt4 + paramInt3, 1);
Dsyr.dsyr(
"Lower",
k,
-1.0D,
paramArrayOfDouble,
2 - 1 + (i - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble,
1 - 1 + (i + 1 - 1) * paramInt4 + paramInt3,
j);
}
i += 1;
}
}
return;
paramintW.val = i;
}
public static void dstev(
String paramString,
int paramInt1,
double[] paramArrayOfDouble1,
int paramInt2,
double[] paramArrayOfDouble2,
int paramInt3,
double[] paramArrayOfDouble3,
int paramInt4,
int paramInt5,
double[] paramArrayOfDouble4,
int paramInt6,
intW paramintW) {
boolean bool = false;
int i = 0;
int j = 0;
double d1 = 0.0D;
double d2 = 0.0D;
double d3 = 0.0D;
double d4 = 0.0D;
double d5 = 0.0D;
double d6 = 0.0D;
double d7 = 0.0D;
double d8 = 0.0D;
bool = Lsame.lsame(paramString, "V");
paramintW.val = 0;
if ((((!bool) && (!Lsame.lsame(paramString, "N")) ? 0 : 1) ^ 0x1) != 0) {
paramintW.val = -1;
} else if ((paramInt1 >= 0 ? 0 : 1) != 0) {
paramintW.val = -2;
} else {
if ((paramInt5 >= 1 ? 0 : 1) == 0) {
if (!bool) {}
}
if ((((paramInt5 >= paramInt1 ? 0 : 1) != 0 ? 1 : 0) == 0 ? 0 : 1) != 0) {
paramintW.val = -6;
}
}
if ((paramintW.val == 0 ? 0 : 1) != 0) {
Xerbla.xerbla("DSTEV ", -paramintW.val);
return;
}
if ((paramInt1 != 0 ? 0 : 1) != 0) {
return;
}
if ((paramInt1 != 1 ? 0 : 1) != 0) {
if (bool) {
paramArrayOfDouble3[(1 - 1 + (1 - 1) * paramInt5 + paramInt4)] = 1.0D;
}
return;
}
d5 = Dlamch.dlamch("Safe minimum");
d2 = Dlamch.dlamch("Precision");
d7 = d5 / d2;
d1 = 1.0D / d7;
d4 = Math.sqrt(d7);
d3 = Math.sqrt(d1);
j = 0;
d8 =
Dlanst.dlanst(
"M", paramInt1, paramArrayOfDouble1, paramInt2, paramArrayOfDouble2, paramInt3);
if ((d8 <= 0.0D ? 0 : 1) != 0) {}
if (((d8 >= d4 ? 0 : 1) != 0 ? 1 : 0) != 0) {
j = 1;
d6 = d4 / d8;
} else if ((d8 <= d3 ? 0 : 1) != 0) {
j = 1;
d6 = d3 / d8;
}
if ((j != 1 ? 0 : 1) != 0) {
Dscal.dscal(paramInt1, d6, paramArrayOfDouble1, paramInt2, 1);
Dscal.dscal(paramInt1 - 1, d6, paramArrayOfDouble2, 1 - 1 + paramInt3, 1);
}
if ((bool ^ true)) {
Dsterf.dsterf(
paramInt1, paramArrayOfDouble1, paramInt2, paramArrayOfDouble2, paramInt3, paramintW);
} else {
Dsteqr.dsteqr(
"I",
paramInt1,
paramArrayOfDouble1,
paramInt2,
paramArrayOfDouble2,
paramInt3,
paramArrayOfDouble3,
paramInt4,
paramInt5,
paramArrayOfDouble4,
paramInt6,
paramintW);
}
if ((j != 1 ? 0 : 1) != 0) {
if ((paramintW.val != 0 ? 0 : 1) != 0) {
i = paramInt1;
} else {
i = paramintW.val - 1;
}
Dscal.dscal(i, 1.0D / d6, paramArrayOfDouble1, paramInt2, 1);
}
}
public static void dsytrs(
String paramString,
int paramInt1,
int paramInt2,
double[] paramArrayOfDouble1,
int paramInt3,
int paramInt4,
int[] paramArrayOfInt,
int paramInt5,
double[] paramArrayOfDouble2,
int paramInt6,
int paramInt7,
intW paramintW) {
boolean bool = false;
int i = 0;
int j = 0;
int k = 0;
double d1 = 0.0D;
double d2 = 0.0D;
double d3 = 0.0D;
double d4 = 0.0D;
double d5 = 0.0D;
double d6 = 0.0D;
paramintW.val = 0;
bool = Lsame.lsame(paramString, "U");
if ((((bool ^ true)) && ((Lsame.lsame(paramString, "L") ^ true)) ? 1 : 0) != 0) {
paramintW.val = -1;
} else if ((paramInt1 >= 0 ? 0 : 1) != 0) {
paramintW.val = -2;
} else if ((paramInt2 >= 0 ? 0 : 1) != 0) {
paramintW.val = -3;
} else if ((paramInt4 >= Math.max(1, paramInt1) ? 0 : 1) != 0) {
paramintW.val = -5;
} else if ((paramInt7 >= Math.max(1, paramInt1) ? 0 : 1) != 0) {
paramintW.val = -8;
}
if ((paramintW.val == 0 ? 0 : 1) != 0) {
Xerbla.xerbla("DSYTRS", -paramintW.val);
return;
}
if ((paramInt1 != 0 ? 0 : 1) == 0) {}
if (((paramInt2 != 0 ? 0 : 1) == 0 ? 0 : 1) != 0) {
return;
}
int m;
if (bool) {
j = paramInt1;
while ((j >= 1 ? 0 : 1) == 0) {
if ((paramArrayOfInt[(j - 1 + paramInt5)] <= 0 ? 0 : 1) != 0) {
k = paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
Dger.dger(
j - 1,
paramInt2,
-1.0D,
paramArrayOfDouble1,
1 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
Dscal.dscal(
paramInt2,
1.0D / paramArrayOfDouble1[(j - 1 + (j - 1) * paramInt4 + paramInt3)],
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
j -= 1;
} else {
k = -paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j - 1 ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j - 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
Dger.dger(
j - 2,
paramInt2,
-1.0D,
paramArrayOfDouble1,
1 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
Dger.dger(
j - 2,
paramInt2,
-1.0D,
paramArrayOfDouble1,
1 - 1 + (j - 1 - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble2,
j - 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
d3 = paramArrayOfDouble1[(j - 1 - 1 + (j - 1) * paramInt4 + paramInt3)];
d2 = paramArrayOfDouble1[(j - 1 - 1 + (j - 1 - 1) * paramInt4 + paramInt3)] / d3;
d1 = paramArrayOfDouble1[(j - 1 + (j - 1) * paramInt4 + paramInt3)] / d3;
d6 = d2 * d1 - 1.0D;
i = 1;
for (m = paramInt2 - 1 + 1; m > 0; m--) {
d5 = paramArrayOfDouble2[(j - 1 - 1 + (i - 1) * paramInt7 + paramInt6)] / d3;
d4 = paramArrayOfDouble2[(j - 1 + (i - 1) * paramInt7 + paramInt6)] / d3;
paramArrayOfDouble2[(j - 1 - 1 + (i - 1) * paramInt7 + paramInt6)] =
((d1 * d5 - d4) / d6);
paramArrayOfDouble2[(j - 1 + (i - 1) * paramInt7 + paramInt6)] = ((d2 * d4 - d5) / d6);
i += 1;
}
j -= 2;
}
}
j = 1;
while ((j <= paramInt1 ? 0 : 1) == 0) {
if ((paramArrayOfInt[(j - 1 + paramInt5)] <= 0 ? 0 : 1) != 0) {
Dgemv.dgemv(
"Transpose",
j - 1,
paramInt2,
-1.0D,
paramArrayOfDouble2,
paramInt6,
paramInt7,
paramArrayOfDouble1,
1 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
1.0D,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
k = paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
j += 1;
} else {
Dgemv.dgemv(
"Transpose",
j - 1,
paramInt2,
-1.0D,
paramArrayOfDouble2,
paramInt6,
paramInt7,
paramArrayOfDouble1,
1 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
1.0D,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
Dgemv.dgemv(
"Transpose",
j - 1,
paramInt2,
-1.0D,
paramArrayOfDouble2,
paramInt6,
paramInt7,
paramArrayOfDouble1,
1 - 1 + (j + 1 - 1) * paramInt4 + paramInt3,
1,
1.0D,
paramArrayOfDouble2,
j + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
k = -paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
j += 2;
}
}
} else {
j = 1;
while ((j <= paramInt1 ? 0 : 1) == 0) {
if ((paramArrayOfInt[(j - 1 + paramInt5)] <= 0 ? 0 : 1) != 0) {
k = paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
if ((j >= paramInt1 ? 0 : 1) != 0) {
Dger.dger(
paramInt1 - j,
paramInt2,
-1.0D,
paramArrayOfDouble1,
j + 1 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
j + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
Dscal.dscal(
paramInt2,
1.0D / paramArrayOfDouble1[(j - 1 + (j - 1) * paramInt4 + paramInt3)],
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
j += 1;
} else {
k = -paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j + 1 ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
if ((j >= paramInt1 - 1 ? 0 : 1) != 0) {
Dger.dger(
paramInt1 - j - 1,
paramInt2,
-1.0D,
paramArrayOfDouble1,
j + 2 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
j + 2 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
Dger.dger(
paramInt1 - j - 1,
paramInt2,
-1.0D,
paramArrayOfDouble1,
j + 2 - 1 + (j + 1 - 1) * paramInt4 + paramInt3,
1,
paramArrayOfDouble2,
j + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
j + 2 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
d3 = paramArrayOfDouble1[(j + 1 - 1 + (j - 1) * paramInt4 + paramInt3)];
d2 = paramArrayOfDouble1[(j - 1 + (j - 1) * paramInt4 + paramInt3)] / d3;
d1 = paramArrayOfDouble1[(j + 1 - 1 + (j + 1 - 1) * paramInt4 + paramInt3)] / d3;
d6 = d2 * d1 - 1.0D;
i = 1;
for (m = paramInt2 - 1 + 1; m > 0; m--) {
d5 = paramArrayOfDouble2[(j - 1 + (i - 1) * paramInt7 + paramInt6)] / d3;
d4 = paramArrayOfDouble2[(j + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] / d3;
paramArrayOfDouble2[(j - 1 + (i - 1) * paramInt7 + paramInt6)] = ((d1 * d5 - d4) / d6);
paramArrayOfDouble2[(j + 1 - 1 + (i - 1) * paramInt7 + paramInt6)] =
((d2 * d4 - d5) / d6);
i += 1;
}
j += 2;
}
}
j = paramInt1;
while ((j >= 1 ? 0 : 1) == 0) {
if ((paramArrayOfInt[(j - 1 + paramInt5)] <= 0 ? 0 : 1) != 0) {
if ((j >= paramInt1 ? 0 : 1) != 0) {
Dgemv.dgemv(
"Transpose",
paramInt1 - j,
paramInt2,
-1.0D,
paramArrayOfDouble2,
j + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble1,
j + 1 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
1.0D,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
k = paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
j -= 1;
} else {
if ((j >= paramInt1 ? 0 : 1) != 0) {
Dgemv.dgemv(
"Transpose",
paramInt1 - j,
paramInt2,
-1.0D,
paramArrayOfDouble2,
j + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble1,
j + 1 - 1 + (j - 1) * paramInt4 + paramInt3,
1,
1.0D,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
Dgemv.dgemv(
"Transpose",
paramInt1 - j,
paramInt2,
-1.0D,
paramArrayOfDouble2,
j + 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble1,
j + 1 - 1 + (j - 1 - 1) * paramInt4 + paramInt3,
1,
1.0D,
paramArrayOfDouble2,
j - 1 - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
k = -paramArrayOfInt[(j - 1 + paramInt5)];
if ((k == j ? 0 : 1) != 0) {
Dswap.dswap(
paramInt2,
paramArrayOfDouble2,
j - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7,
paramArrayOfDouble2,
k - 1 + (1 - 1) * paramInt7 + paramInt6,
paramInt7);
}
j -= 2;
}
}
}
}