37 #ifndef TEMPLATE_BLAS_TPSV_HEADER
38 #define TEMPLATE_BLAS_TPSV_HEADER
43 const Treal *ap, Treal *x,
const integer *incx)
135 }
else if (*incx == 0) {
150 kx = 1 - (*n - 1) * *incx;
151 }
else if (*incx != 1) {
159 kk = *n * (*n + 1) / 2;
161 for (j = *n; j >= 1; --j) {
168 for (i__ = j - 1; i__ >= 1; --i__) {
169 x[i__] -= temp * ap[k];
178 jx = kx + (*n - 1) * *incx;
179 for (j = *n; j >= 1; --j) {
187 for (k = kk - 1; k >= i__1; --k) {
189 x[ix] -= temp * ap[k];
202 for (j = 1; j <= i__1; ++j) {
210 for (i__ = j + 1; i__ <= i__2; ++i__) {
211 x[i__] -= temp * ap[k];
222 for (j = 1; j <= i__1; ++j) {
230 for (k = kk + 1; k <= i__2; ++k) {
232 x[ix] -= temp * ap[k];
248 for (j = 1; j <= i__1; ++j) {
252 for (i__ = 1; i__ <= i__2; ++i__) {
253 temp -= ap[k] * x[i__];
258 temp /= ap[kk + j - 1];
267 for (j = 1; j <= i__1; ++j) {
271 for (k = kk; k <= i__2; ++k) {
272 temp -= ap[k] * x[ix];
277 temp /= ap[kk + j - 1];
286 kk = *n * (*n + 1) / 2;
288 for (j = *n; j >= 1; --j) {
292 for (i__ = *n; i__ >= i__1; --i__) {
293 temp -= ap[k] * x[i__];
298 temp /= ap[kk - *n + j];
305 kx += (*n - 1) * *incx;
307 for (j = *n; j >= 1; --j) {
310 i__1 = kk - (*n - (j + 1));
311 for (k = kk; k >= i__1; --k) {
312 temp -= ap[k] * x[ix];
317 temp /= ap[kk - *n + j];