📄 zlarfx.c
字号:
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
/*< SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) >*/
i__2 = j * c_dim1 + 1;
z__3.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__3.i = v1.r *
c__[i__2].i + v1.i * c__[i__2].r;
i__3 = j * c_dim1 + 2;
z__4.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__4.i = v2.r *
c__[i__3].i + v2.i * c__[i__3].r;
z__2.r = z__3.r + z__4.r, z__2.i = z__3.i + z__4.i;
i__4 = j * c_dim1 + 3;
z__5.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__5.i = v3.r *
c__[i__4].i + v3.i * c__[i__4].r;
z__1.r = z__2.r + z__5.r, z__1.i = z__2.i + z__5.i;
sum.r = z__1.r, sum.i = z__1.i;
/*< C( 1, J ) = C( 1, J ) - SUM*T1 >*/
i__2 = j * c_dim1 + 1;
i__3 = j * c_dim1 + 1;
z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
sum.i * t1.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 2, J ) = C( 2, J ) - SUM*T2 >*/
i__2 = j * c_dim1 + 2;
i__3 = j * c_dim1 + 2;
z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
sum.i * t2.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 3, J ) = C( 3, J ) - SUM*T3 >*/
i__2 = j * c_dim1 + 3;
i__3 = j * c_dim1 + 3;
z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
sum.i * t3.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< 60 CONTINUE >*/
/* L60: */
}
/*< GO TO 410 >*/
goto L410;
/*< 70 CONTINUE >*/
L70:
/* Special code for 4 x 4 Householder */
/*< V1 = DCONJG( V( 1 ) ) >*/
d_cnjg(&z__1, &v[1]);
v1.r = z__1.r, v1.i = z__1.i;
/*< T1 = TAU*DCONJG( V1 ) >*/
d_cnjg(&z__2, &v1);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t1.r = z__1.r, t1.i = z__1.i;
/*< V2 = DCONJG( V( 2 ) ) >*/
d_cnjg(&z__1, &v[2]);
v2.r = z__1.r, v2.i = z__1.i;
/*< T2 = TAU*DCONJG( V2 ) >*/
d_cnjg(&z__2, &v2);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t2.r = z__1.r, t2.i = z__1.i;
/*< V3 = DCONJG( V( 3 ) ) >*/
d_cnjg(&z__1, &v[3]);
v3.r = z__1.r, v3.i = z__1.i;
/*< T3 = TAU*DCONJG( V3 ) >*/
d_cnjg(&z__2, &v3);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t3.r = z__1.r, t3.i = z__1.i;
/*< V4 = DCONJG( V( 4 ) ) >*/
d_cnjg(&z__1, &v[4]);
v4.r = z__1.r, v4.i = z__1.i;
/*< T4 = TAU*DCONJG( V4 ) >*/
d_cnjg(&z__2, &v4);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t4.r = z__1.r, t4.i = z__1.i;
/*< DO 80 J = 1, N >*/
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
/*< >*/
i__2 = j * c_dim1 + 1;
z__4.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__4.i = v1.r *
c__[i__2].i + v1.i * c__[i__2].r;
i__3 = j * c_dim1 + 2;
z__5.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__5.i = v2.r *
c__[i__3].i + v2.i * c__[i__3].r;
z__3.r = z__4.r + z__5.r, z__3.i = z__4.i + z__5.i;
i__4 = j * c_dim1 + 3;
z__6.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__6.i = v3.r *
c__[i__4].i + v3.i * c__[i__4].r;
z__2.r = z__3.r + z__6.r, z__2.i = z__3.i + z__6.i;
i__5 = j * c_dim1 + 4;
z__7.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__7.i = v4.r *
c__[i__5].i + v4.i * c__[i__5].r;
z__1.r = z__2.r + z__7.r, z__1.i = z__2.i + z__7.i;
sum.r = z__1.r, sum.i = z__1.i;
/*< C( 1, J ) = C( 1, J ) - SUM*T1 >*/
i__2 = j * c_dim1 + 1;
i__3 = j * c_dim1 + 1;
z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
sum.i * t1.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 2, J ) = C( 2, J ) - SUM*T2 >*/
i__2 = j * c_dim1 + 2;
i__3 = j * c_dim1 + 2;
z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
sum.i * t2.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 3, J ) = C( 3, J ) - SUM*T3 >*/
i__2 = j * c_dim1 + 3;
i__3 = j * c_dim1 + 3;
z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
sum.i * t3.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 4, J ) = C( 4, J ) - SUM*T4 >*/
i__2 = j * c_dim1 + 4;
i__3 = j * c_dim1 + 4;
z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
sum.i * t4.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< 80 CONTINUE >*/
/* L80: */
}
/*< GO TO 410 >*/
goto L410;
/*< 90 CONTINUE >*/
L90:
/* Special code for 5 x 5 Householder */
/*< V1 = DCONJG( V( 1 ) ) >*/
d_cnjg(&z__1, &v[1]);
v1.r = z__1.r, v1.i = z__1.i;
/*< T1 = TAU*DCONJG( V1 ) >*/
d_cnjg(&z__2, &v1);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t1.r = z__1.r, t1.i = z__1.i;
/*< V2 = DCONJG( V( 2 ) ) >*/
d_cnjg(&z__1, &v[2]);
v2.r = z__1.r, v2.i = z__1.i;
/*< T2 = TAU*DCONJG( V2 ) >*/
d_cnjg(&z__2, &v2);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t2.r = z__1.r, t2.i = z__1.i;
/*< V3 = DCONJG( V( 3 ) ) >*/
d_cnjg(&z__1, &v[3]);
v3.r = z__1.r, v3.i = z__1.i;
/*< T3 = TAU*DCONJG( V3 ) >*/
d_cnjg(&z__2, &v3);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t3.r = z__1.r, t3.i = z__1.i;
/*< V4 = DCONJG( V( 4 ) ) >*/
d_cnjg(&z__1, &v[4]);
v4.r = z__1.r, v4.i = z__1.i;
/*< T4 = TAU*DCONJG( V4 ) >*/
d_cnjg(&z__2, &v4);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t4.r = z__1.r, t4.i = z__1.i;
/*< V5 = DCONJG( V( 5 ) ) >*/
d_cnjg(&z__1, &v[5]);
v5.r = z__1.r, v5.i = z__1.i;
/*< T5 = TAU*DCONJG( V5 ) >*/
d_cnjg(&z__2, &v5);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t5.r = z__1.r, t5.i = z__1.i;
/*< DO 100 J = 1, N >*/
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
/*< >*/
i__2 = j * c_dim1 + 1;
z__5.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__5.i = v1.r *
c__[i__2].i + v1.i * c__[i__2].r;
i__3 = j * c_dim1 + 2;
z__6.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__6.i = v2.r *
c__[i__3].i + v2.i * c__[i__3].r;
z__4.r = z__5.r + z__6.r, z__4.i = z__5.i + z__6.i;
i__4 = j * c_dim1 + 3;
z__7.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__7.i = v3.r *
c__[i__4].i + v3.i * c__[i__4].r;
z__3.r = z__4.r + z__7.r, z__3.i = z__4.i + z__7.i;
i__5 = j * c_dim1 + 4;
z__8.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__8.i = v4.r *
c__[i__5].i + v4.i * c__[i__5].r;
z__2.r = z__3.r + z__8.r, z__2.i = z__3.i + z__8.i;
i__6 = j * c_dim1 + 5;
z__9.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__9.i = v5.r *
c__[i__6].i + v5.i * c__[i__6].r;
z__1.r = z__2.r + z__9.r, z__1.i = z__2.i + z__9.i;
sum.r = z__1.r, sum.i = z__1.i;
/*< C( 1, J ) = C( 1, J ) - SUM*T1 >*/
i__2 = j * c_dim1 + 1;
i__3 = j * c_dim1 + 1;
z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
sum.i * t1.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 2, J ) = C( 2, J ) - SUM*T2 >*/
i__2 = j * c_dim1 + 2;
i__3 = j * c_dim1 + 2;
z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
sum.i * t2.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 3, J ) = C( 3, J ) - SUM*T3 >*/
i__2 = j * c_dim1 + 3;
i__3 = j * c_dim1 + 3;
z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
sum.i * t3.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 4, J ) = C( 4, J ) - SUM*T4 >*/
i__2 = j * c_dim1 + 4;
i__3 = j * c_dim1 + 4;
z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
sum.i * t4.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< C( 5, J ) = C( 5, J ) - SUM*T5 >*/
i__2 = j * c_dim1 + 5;
i__3 = j * c_dim1 + 5;
z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
sum.i * t5.r;
z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
/*< 100 CONTINUE >*/
/* L100: */
}
/*< GO TO 410 >*/
goto L410;
/*< 110 CONTINUE >*/
L110:
/* Special code for 6 x 6 Householder */
/*< V1 = DCONJG( V( 1 ) ) >*/
d_cnjg(&z__1, &v[1]);
v1.r = z__1.r, v1.i = z__1.i;
/*< T1 = TAU*DCONJG( V1 ) >*/
d_cnjg(&z__2, &v1);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t1.r = z__1.r, t1.i = z__1.i;
/*< V2 = DCONJG( V( 2 ) ) >*/
d_cnjg(&z__1, &v[2]);
v2.r = z__1.r, v2.i = z__1.i;
/*< T2 = TAU*DCONJG( V2 ) >*/
d_cnjg(&z__2, &v2);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t2.r = z__1.r, t2.i = z__1.i;
/*< V3 = DCONJG( V( 3 ) ) >*/
d_cnjg(&z__1, &v[3]);
v3.r = z__1.r, v3.i = z__1.i;
/*< T3 = TAU*DCONJG( V3 ) >*/
d_cnjg(&z__2, &v3);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t3.r = z__1.r, t3.i = z__1.i;
/*< V4 = DCONJG( V( 4 ) ) >*/
d_cnjg(&z__1, &v[4]);
v4.r = z__1.r, v4.i = z__1.i;
/*< T4 = TAU*DCONJG( V4 ) >*/
d_cnjg(&z__2, &v4);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t4.r = z__1.r, t4.i = z__1.i;
/*< V5 = DCONJG( V( 5 ) ) >*/
d_cnjg(&z__1, &v[5]);
v5.r = z__1.r, v5.i = z__1.i;
/*< T5 = TAU*DCONJG( V5 ) >*/
d_cnjg(&z__2, &v5);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t5.r = z__1.r, t5.i = z__1.i;
/*< V6 = DCONJG( V( 6 ) ) >*/
d_cnjg(&z__1, &v[6]);
v6.r = z__1.r, v6.i = z__1.i;
/*< T6 = TAU*DCONJG( V6 ) >*/
d_cnjg(&z__2, &v6);
z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
+ tau->i * z__2.r;
t6.r = z__1.r, t6.i = z__1.i;
/*< DO 120 J = 1, N >*/
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
/*< >*/
i__2 = j * c_dim1 + 1;
z__6.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__6.i = v1.r *
c__[i__2].i + v1.i * c__[i__2].r;
i__3 = j * c_dim1 + 2;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -