ztrsv.c

InsightToolkit-1.4.0(有大量的优化算法程序)

                }
            } else {
                jx = kx + (*n - 1) * *incx;
                for (j = *n - 1; j >= 0; --j) {
                    if (x[jx].r != 0. || x[jx].i != 0.) {
                        if (nounit) {
                            z_div(&x[jx], &x[jx], &a[j + j * *lda]);
                        }
                        temp.r = x[jx].r, temp.i = x[jx].i;
                        ix = jx;
                        for (i = j - 1; i >= 0; --i) {
                            ix -= *incx;
                            i__1 = i + j * *lda;
                            x[ix].r -= temp.r * a[i__1].r - temp.i * a[i__1].i,
                            x[ix].i -= temp.r * a[i__1].i + temp.i * a[i__1].r;
                        }
                    }
                    jx -= *incx;
                }
            }
        } else {
            if (*incx == 1) {
                for (j = 0; j < *n; ++j) {
                    if (x[j].r != 0. || x[j].i != 0.) {
                        if (nounit) {
                            z_div(&x[j], &x[j], &a[j + j * *lda]);
                        }
                        temp.r = x[j].r, temp.i = x[j].i;
                        for (i = j + 1; i < *n; ++i) {
                            i__1 = i + j * *lda;
                            x[i].r -= temp.r * a[i__1].r - temp.i * a[i__1].i,
                            x[i].i -= temp.r * a[i__1].i + temp.i * a[i__1].r;
                        }
                    }
                }
            } else {
                jx = kx;
                for (j = 0; j < *n; ++j) {
                    if (x[jx].r != 0. || x[jx].i != 0.) {
                        if (nounit) {
                            z_div(&x[jx], &x[jx], &a[j + j * *lda]);
                        }
                        temp.r = x[jx].r, temp.i = x[jx].i;
                        ix = jx;
                        for (i = j + 1; i < *n; ++i) {
                            ix += *incx;
                            i__1 = i + j * *lda;
                            x[ix].r -= temp.r * a[i__1].r - temp.i * a[i__1].i,
                            x[ix].i -= temp.r * a[i__1].i + temp.i * a[i__1].r;
                        }
                    }
                    jx += *incx;
                }
            }
        }
    } else {

/*        Form  x := inv( A' )*x  or  x := inv( conjg( A' ) )*x. */

        if (lsame_(uplo, "U")) {
            if (*incx == 1) {
                for (j = 0; j < *n; ++j) {
                    temp.r = x[j].r, temp.i = x[j].i;
                    if (noconj) {
                        for (i = 0; i < j; ++i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[i].r - a[i__1].i * x[i].i,
                            temp.i -= a[i__1].r * x[i].i + a[i__1].i * x[i].r;
                        }
                        if (nounit) {
                            z_div(&temp, &temp, &a[j + j * *lda]);
                        }
                    } else {
                        for (i = 0; i < j; ++i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[i].r + a[i__1].i * x[i].i,
                            temp.i -= a[i__1].r * x[i].i - a[i__1].i * x[i].r;
                        }
                        if (nounit) {
                            d_cnjg(&z__1, &a[j + j * *lda]);
                            z_div(&temp, &temp, &z__1);
                        }
                    }
                    x[j].r = temp.r, x[j].i = temp.i;
                }
            } else {
                jx = kx;
                for (j = 0; j < *n; ++j) {
                    ix = kx;
                    temp.r = x[jx].r, temp.i = x[jx].i;
                    if (noconj) {
                        for (i = 0; i < j; ++i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[ix].r - a[i__1].i * x[ix].i,
                            temp.i -= a[i__1].r * x[ix].i + a[i__1].i * x[ix].r;
                            ix += *incx;
                        }
                        if (nounit) {
                            z_div(&temp, &temp, &a[j + j * *lda]);
                        }
                    } else {
                        for (i = 0; i < j; ++i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[ix].r + a[i__1].i * x[ix].i,
                            temp.i -= a[i__1].r * x[ix].i - a[i__1].i * x[ix].r;
                            ix += *incx;
                        }
                        if (nounit) {
                            d_cnjg(&z__1, &a[j + j * *lda]);
                            z_div(&temp, &temp, &z__1);
                        }
                    }
                    x[jx].r = temp.r, x[jx].i = temp.i;
                    jx += *incx;
                }
            }
        } else {
            if (*incx == 1) {
                for (j = *n - 1; j >= 0; --j) {
                    temp.r = x[j].r, temp.i = x[j].i;
                    if (noconj) {
                        for (i = *n - 1; i > j; --i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[i].r - a[i__1].i * x[i].i,
                            temp.i -= a[i__1].r * x[i].i + a[i__1].i * x[i].r;
                        }
                        if (nounit) {
                            z_div(&temp, &temp, &a[j + j * *lda]);
                        }
                    } else {
                        for (i = *n - 1; i > j; --i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[i].r + a[i__1].i * x[i].i,
                            temp.i -= a[i__1].r * x[i].i - a[i__1].i * x[i].r;
                        }
                        if (nounit) {
                            d_cnjg(&z__1, &a[j + j * *lda]);
                            z_div(&temp, &temp, &z__1);
                        }
                    }
                    x[j].r = temp.r, x[j].i = temp.i;
                }
            } else {
                kx += (*n - 1) * *incx;
                jx = kx;
                for (j = *n - 1; j >= 0; --j) {
                    ix = kx;
                    temp.r = x[jx].r, temp.i = x[jx].i;
                    if (noconj) {
                        for (i = *n - 1; i > j; --i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[ix].r - a[i__1].i * x[ix].i,
                            temp.i -= a[i__1].r * x[ix].i + a[i__1].i * x[ix].r;
                            ix -= *incx;
                        }
                        if (nounit) {
                            z_div(&temp, &temp, &a[j + j * *lda]);
                        }
                    } else {
                        for (i = *n - 1; i > j; --i) {
                            i__1 = i + j * *lda;
                            temp.r -= a[i__1].r * x[ix].r + a[i__1].i * x[ix].i,
                            temp.i -= a[i__1].r * x[ix].i - a[i__1].i * x[ix].r;
                            ix -= *incx;
                        }
                        if (nounit) {
                            d_cnjg(&z__1, &a[j + j * *lda]);
                            z_div(&temp, &temp, &z__1);
                        }
                    }
                    x[jx].r = temp.r, x[jx].i = temp.i;
                    jx -= *incx;
                }
            }
        }
    }
} /* ztrsv_ */