00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035 #ifndef TEMPLATE_LAPACK_GGBAK_HEADER
00036 #define TEMPLATE_LAPACK_GGBAK_HEADER
00037
00038
00039 template<class Treal>
00040 int template_lapack_ggbak(const char *job, const char *side, const integer *n, const integer *ilo,
00041 const integer *ihi, const Treal *lscale, const Treal *rscale, const integer *m,
00042 Treal *v, const integer *ldv, integer *info)
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
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118 integer v_dim1, v_offset, i__1;
00119
00120 integer i__, k;
00121 logical leftv;
00122 logical rightv;
00123 #define v_ref(a_1,a_2) v[(a_2)*v_dim1 + a_1]
00124
00125 --lscale;
00126 --rscale;
00127 v_dim1 = *ldv;
00128 v_offset = 1 + v_dim1 * 1;
00129 v -= v_offset;
00130
00131
00132 rightv = template_blas_lsame(side, "R");
00133 leftv = template_blas_lsame(side, "L");
00134
00135 *info = 0;
00136 if (! template_blas_lsame(job, "N") && ! template_blas_lsame(job, "P") && ! template_blas_lsame(job, "S")
00137 && ! template_blas_lsame(job, "B")) {
00138 *info = -1;
00139 } else if (! rightv && ! leftv) {
00140 *info = -2;
00141 } else if (*n < 0) {
00142 *info = -3;
00143 } else if (*ilo < 1) {
00144 *info = -4;
00145 } else if (*ihi < *ilo || *ihi > maxMACRO(1,*n)) {
00146 *info = -5;
00147 } else if (*m < 0) {
00148 *info = -6;
00149 } else if (*ldv < maxMACRO(1,*n)) {
00150 *info = -10;
00151 }
00152 if (*info != 0) {
00153 i__1 = -(*info);
00154 template_blas_erbla("GGBAK ", &i__1);
00155 return 0;
00156 }
00157
00158
00159
00160 if (*n == 0) {
00161 return 0;
00162 }
00163 if (*m == 0) {
00164 return 0;
00165 }
00166 if (template_blas_lsame(job, "N")) {
00167 return 0;
00168 }
00169
00170 if (*ilo == *ihi) {
00171 goto L30;
00172 }
00173
00174
00175
00176 if (template_blas_lsame(job, "S") || template_blas_lsame(job, "B")) {
00177
00178
00179
00180 if (rightv) {
00181 i__1 = *ihi;
00182 for (i__ = *ilo; i__ <= i__1; ++i__) {
00183 template_blas_scal(m, &rscale[i__], &v_ref(i__, 1), ldv);
00184
00185 }
00186 }
00187
00188
00189
00190 if (leftv) {
00191 i__1 = *ihi;
00192 for (i__ = *ilo; i__ <= i__1; ++i__) {
00193 template_blas_scal(m, &lscale[i__], &v_ref(i__, 1), ldv);
00194
00195 }
00196 }
00197 }
00198
00199
00200
00201 L30:
00202 if (template_blas_lsame(job, "P") || template_blas_lsame(job, "B")) {
00203
00204
00205
00206 if (rightv) {
00207 if (*ilo == 1) {
00208 goto L50;
00209 }
00210
00211 for (i__ = *ilo - 1; i__ >= 1; --i__) {
00212 k = (integer) rscale[i__];
00213 if (k == i__) {
00214 goto L40;
00215 }
00216 template_blas_swap(m, &v_ref(i__, 1), ldv, &v_ref(k, 1), ldv);
00217 L40:
00218 ;
00219 }
00220
00221 L50:
00222 if (*ihi == *n) {
00223 goto L70;
00224 }
00225 i__1 = *n;
00226 for (i__ = *ihi + 1; i__ <= i__1; ++i__) {
00227 k = (integer) rscale[i__];
00228 if (k == i__) {
00229 goto L60;
00230 }
00231 template_blas_swap(m, &v_ref(i__, 1), ldv, &v_ref(k, 1), ldv);
00232 L60:
00233 ;
00234 }
00235 }
00236
00237
00238
00239 L70:
00240 if (leftv) {
00241 if (*ilo == 1) {
00242 goto L90;
00243 }
00244 for (i__ = *ilo - 1; i__ >= 1; --i__) {
00245 k = (integer) lscale[i__];
00246 if (k == i__) {
00247 goto L80;
00248 }
00249 template_blas_swap(m, &v_ref(i__, 1), ldv, &v_ref(k, 1), ldv);
00250 L80:
00251 ;
00252 }
00253
00254 L90:
00255 if (*ihi == *n) {
00256 goto L110;
00257 }
00258 i__1 = *n;
00259 for (i__ = *ihi + 1; i__ <= i__1; ++i__) {
00260 k = (integer) lscale[i__];
00261 if (k == i__) {
00262 goto L100;
00263 }
00264 template_blas_swap(m, &v_ref(i__, 1), ldv, &v_ref(k, 1), ldv);
00265 L100:
00266 ;
00267 }
00268 }
00269 }
00270
00271 L110:
00272
00273 return 0;
00274
00275
00276
00277 }
00278
00279 #undef v_ref
00280
00281
00282 #endif