00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015
00016 int cupgtr_(char *uplo, integer *n, complex *ap, complex *
00017 tau, complex *q, integer *ldq, complex *work, integer *info)
00018 {
00019
00020 integer q_dim1, q_offset, i__1, i__2, i__3, i__4;
00021
00022
00023 integer i__, j, ij;
00024 extern logical lsame_(char *, char *);
00025 integer iinfo;
00026 logical upper;
00027 extern int cung2l_(integer *, integer *, integer *,
00028 complex *, integer *, complex *, complex *, integer *), cung2r_(
00029 integer *, integer *, integer *, complex *, integer *, complex *,
00030 complex *, integer *), xerbla_(char *, integer *);
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102 --ap;
00103 --tau;
00104 q_dim1 = *ldq;
00105 q_offset = 1 + q_dim1;
00106 q -= q_offset;
00107 --work;
00108
00109
00110 *info = 0;
00111 upper = lsame_(uplo, "U");
00112 if (! upper && ! lsame_(uplo, "L")) {
00113 *info = -1;
00114 } else if (*n < 0) {
00115 *info = -2;
00116 } else if (*ldq < max(1,*n)) {
00117 *info = -6;
00118 }
00119 if (*info != 0) {
00120 i__1 = -(*info);
00121 xerbla_("CUPGTR", &i__1);
00122 return 0;
00123 }
00124
00125
00126
00127 if (*n == 0) {
00128 return 0;
00129 }
00130
00131 if (upper) {
00132
00133
00134
00135
00136
00137
00138
00139 ij = 2;
00140 i__1 = *n - 1;
00141 for (j = 1; j <= i__1; ++j) {
00142 i__2 = j - 1;
00143 for (i__ = 1; i__ <= i__2; ++i__) {
00144 i__3 = i__ + j * q_dim1;
00145 i__4 = ij;
00146 q[i__3].r = ap[i__4].r, q[i__3].i = ap[i__4].i;
00147 ++ij;
00148
00149 }
00150 ij += 2;
00151 i__2 = *n + j * q_dim1;
00152 q[i__2].r = 0.f, q[i__2].i = 0.f;
00153
00154 }
00155 i__1 = *n - 1;
00156 for (i__ = 1; i__ <= i__1; ++i__) {
00157 i__2 = i__ + *n * q_dim1;
00158 q[i__2].r = 0.f, q[i__2].i = 0.f;
00159
00160 }
00161 i__1 = *n + *n * q_dim1;
00162 q[i__1].r = 1.f, q[i__1].i = 0.f;
00163
00164
00165
00166 i__1 = *n - 1;
00167 i__2 = *n - 1;
00168 i__3 = *n - 1;
00169 cung2l_(&i__1, &i__2, &i__3, &q[q_offset], ldq, &tau[1], &work[1], &
00170 iinfo);
00171
00172 } else {
00173
00174
00175
00176
00177
00178
00179
00180 i__1 = q_dim1 + 1;
00181 q[i__1].r = 1.f, q[i__1].i = 0.f;
00182 i__1 = *n;
00183 for (i__ = 2; i__ <= i__1; ++i__) {
00184 i__2 = i__ + q_dim1;
00185 q[i__2].r = 0.f, q[i__2].i = 0.f;
00186
00187 }
00188 ij = 3;
00189 i__1 = *n;
00190 for (j = 2; j <= i__1; ++j) {
00191 i__2 = j * q_dim1 + 1;
00192 q[i__2].r = 0.f, q[i__2].i = 0.f;
00193 i__2 = *n;
00194 for (i__ = j + 1; i__ <= i__2; ++i__) {
00195 i__3 = i__ + j * q_dim1;
00196 i__4 = ij;
00197 q[i__3].r = ap[i__4].r, q[i__3].i = ap[i__4].i;
00198 ++ij;
00199
00200 }
00201 ij += 2;
00202
00203 }
00204 if (*n > 1) {
00205
00206
00207
00208 i__1 = *n - 1;
00209 i__2 = *n - 1;
00210 i__3 = *n - 1;
00211 cung2r_(&i__1, &i__2, &i__3, &q[(q_dim1 << 1) + 2], ldq, &tau[1],
00212 &work[1], &iinfo);
00213 }
00214 }
00215 return 0;
00216
00217
00218
00219 }