23 integer i__,
j, k, kk, ix, iy, jx, jy, kx, ky,
info;
145 }
else if (*incx == 0) {
147 }
else if (*incy == 0) {
157 if (*n == 0 || (*alpha == 0. && *beta == 1.)) {
166 kx = 1 - (*n - 1) * *incx;
171 ky = 1 - (*n - 1) * *incy;
183 for (i__ = 1; i__ <= i__1; ++i__) {
189 for (i__ = 1; i__ <= i__1; ++i__) {
190 y[i__] = *beta * y[i__];
198 for (i__ = 1; i__ <= i__1; ++i__) {
205 for (i__ = 1; i__ <= i__1; ++i__) {
206 y[iy] = *beta * y[iy];
221 if (*incx == 1 && *incy == 1) {
223 for (j = 1; j <= i__1; ++
j) {
224 temp1 = *alpha * x[
j];
228 for (i__ = 1; i__ <= i__2; ++i__) {
229 y[i__] += temp1 * ap[k];
230 temp2 += ap[k] * x[i__];
234 y[
j] = y[
j] + temp1 * ap[kk + j - 1] + *alpha * temp2;
242 for (j = 1; j <= i__1; ++
j) {
243 temp1 = *alpha * x[jx];
248 for (k = kk; k <= i__2; ++k) {
249 y[iy] += temp1 * ap[k];
250 temp2 += ap[k] * x[ix];
255 y[jy] = y[jy] + temp1 * ap[kk + j - 1] + *alpha * temp2;
266 if (*incx == 1 && *incy == 1) {
268 for (j = 1; j <= i__1; ++
j) {
269 temp1 = *alpha * x[
j];
271 y[
j] += temp1 * ap[kk];
274 for (i__ = j + 1; i__ <= i__2; ++i__) {
275 y[i__] += temp1 * ap[k];
276 temp2 += ap[k] * x[i__];
280 y[
j] += *alpha * temp2;
288 for (j = 1; j <= i__1; ++
j) {
289 temp1 = *alpha * x[jx];
291 y[jy] += temp1 * ap[kk];
295 for (k = kk + 1; k <= i__2; ++k) {
298 y[iy] += temp1 * ap[k];
299 temp2 += ap[k] * x[ix];
302 y[jy] += *alpha * temp2;
int dspmv_(char *uplo, integer *n, doublereal *alpha, doublereal *ap, doublereal *x, integer *incx, doublereal *beta, doublereal *y, integer *incy, ftnlen uplo_len)
RealScalar RealScalar int * incx
EIGEN_WEAK_LINKING int xerbla_(const char *msg, int *info, int)
logical lsame_(char *ca, char *cb, ftnlen ca_len, ftnlen cb_len)
set noclip points set clip one set noclip two set bar set border lt lw set xdata set ydata set zdata set x2data set y2data set boxwidth set dummy x
int RealScalar int RealScalar int * incy