📄 zlarfx.c
字号:
z__7.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__7.i = v2.r *
c__[i__3].i + v2.i * c__[i__3].r;
z__5.r = z__6.r + z__7.r, z__5.i = z__6.i + z__7.i;
i__4 = j * c_dim1 + 3;
z__8.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__8.i = v3.r *
c__[i__4].i + v3.i * c__[i__4].r;
z__4.r = z__5.r + z__8.r, z__4.i = z__5.i + z__8.i;
i__5 = j * c_dim1 + 4;
z__9.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__9.i = v4.r *
c__[i__5].i + v4.i * c__[i__5].r;
z__3.r = z__4.r + z__9.r, z__3.i = z__4.i + z__9.i;
i__6 = j * c_dim1 + 5;
z__10.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__10.i = v5.r
* c__[i__6].i + v5.i * c__[i__6].r;
z__2.r = z__3.r + z__10.r, z__2.i = z__3.i + z__10.i;
i__7 = j * c_dim1 + 6;
z__11.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__11.i = v6.r
* c__[i__7].i + v6.i * c__[i__7].r;
z__1.r = z__2.r + z__11.r, z__1.i = z__2.i + z__11.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;
/*< C( 6, J ) = C( 6, J ) - SUM*T6 >*/
i__2 = j * c_dim1 + 6;
i__3 = j * c_dim1 + 6;
z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
sum.i * t6.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;
/*< 120 CONTINUE >*/
/* L120: */
}
/*< GO TO 410 >*/
goto L410;
/*< 130 CONTINUE >*/
L130:
/* Special code for 7 x 7 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;
/*< V7 = DCONJG( V( 7 ) ) >*/
d_cnjg(&z__1, &v[7]);
v7.r = z__1.r, v7.i = z__1.i;
/*< T7 = TAU*DCONJG( V7 ) >*/
d_cnjg(&z__2, &v7);
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;
t7.r = z__1.r, t7.i = z__1.i;
/*< DO 140 J = 1, N >*/
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
/*< >*/
i__2 = j * c_dim1 + 1;
z__7.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__7.i = v1.r *
c__[i__2].i + v1.i * c__[i__2].r;
i__3 = j * c_dim1 + 2;
z__8.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__8.i = v2.r *
c__[i__3].i + v2.i * c__[i__3].r;
z__6.r = z__7.r + z__8.r, z__6.i = z__7.i + z__8.i;
i__4 = j * c_dim1 + 3;
z__9.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__9.i = v3.r *
c__[i__4].i + v3.i * c__[i__4].r;
z__5.r = z__6.r + z__9.r, z__5.i = z__6.i + z__9.i;
i__5 = j * c_dim1 + 4;
z__10.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__10.i = v4.r
* c__[i__5].i + v4.i * c__[i__5].r;
z__4.r = z__5.r + z__10.r, z__4.i = z__5.i + z__10.i;
i__6 = j * c_dim1 + 5;
z__11.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__11.i = v5.r
* c__[i__6].i + v5.i * c__[i__6].r;
z__3.r = z__4.r + z__11.r, z__3.i = z__4.i + z__11.i;
i__7 = j * c_dim1 + 6;
z__12.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__12.i = v6.r
* c__[i__7].i + v6.i * c__[i__7].r;
z__2.r = z__3.r + z__12.r, z__2.i = z__3.i + z__12.i;
i__8 = j * c_dim1 + 7;
z__13.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__13.i = v7.r
* c__[i__8].i + v7.i * c__[i__8].r;
z__1.r = z__2.r + z__13.r, z__1.i = z__2.i + z__13.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;
/*< C( 6, J ) = C( 6, J ) - SUM*T6 >*/
i__2 = j * c_dim1 + 6;
i__3 = j * c_dim1 + 6;
z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
sum.i * t6.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( 7, J ) = C( 7, J ) - SUM*T7 >*/
i__2 = j * c_dim1 + 7;
i__3 = j * c_dim1 + 7;
z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
sum.i * t7.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;
/*< 140 CONTINUE >*/
/* L140: */
}
/*< GO TO 410 >*/
goto L410;
/*< 150 CONTINUE >*/
L150:
/* Special code for 8 x 8 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;
/*< V7 = DCONJG( V( 7 ) ) >*/
d_cnjg(&z__1, &v[7]);
v7.r = z__1.r, v7.i = z__1.i;
/*< T7 = TAU*DCONJG( V7 ) >*/
d_cnjg(&z__2, &v7);
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;
t7.r = z__1.r, t7.i = z__1.i;
/*< V8 = DCONJG( V( 8 ) ) >*/
d_cnjg(&z__1, &v[8]);
v8.r = z__1.r, v8.i = z__1.i;
/*< T8 = TAU*DCONJG( V8 ) >*/
d_cnjg(&z__2, &v8);
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;
t8.r = z__1.r, t8.i = z__1.i;
/*< DO 160 J = 1, N >*/
i__1 = *n;
for (j = 1; j <= i__1; ++j) {
/*< >*/
i__2 = j * c_dim1 + 1;
z__8.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__8.i = v1.r *
c__[i__2].i + v1.i * c__[i__2].r;
i__3 = j * c_dim1 + 2;
z__9.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__9.i = v2.r *
c__[i__3].i + v2.i * c__[i__3].r;
z__7.r = z__8.r + z__9.r, z__7.i = z__8.i + z__9.i;
i__4 = j * c_dim1 + 3;
z__10.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__10.i = v3.r
* c__[i__4].i + v3.i * c__[i__4].r;
z__6.r = z__7.r + z__10.r, z__6.i = z__7.i + z__10.i;
i__5 = j * c_dim1 + 4;
z__11.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__11.i = v4.r
* c__[i__5].i + v4.i * c__[i__5].r;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -