475 |
|
for (i = sd.dim[0].ne; i--; ) |
476 |
|
rvec[i] += c0*sd.arr.d[i] |
477 |
|
+ c1*sd.arr.d[i+stride]; |
478 |
< |
return(0.); |
479 |
< |
} |
480 |
< |
if (dp->type == SPECTY) { |
478 |
> |
} else if (dp->type == SPECTY) { |
479 |
|
double f; |
480 |
|
sd.arr.s = dp->arr.s + i*stride; |
481 |
|
f = ldexp(1.0, (int)sd.arr.s[sd.dim[0].ne] |
487 |
|
- (COLXS+8)); |
488 |
|
for (i = sd.dim[0].ne; i--; ) |
489 |
|
rvec[i] += c1*f*(sd.arr.s[i] + 0.5); |
490 |
< |
return(0.); |
490 |
> |
} else { |
491 |
> |
sd.arr.c = dp->arr.c + i*stride; |
492 |
> |
for (i = sd.dim[0].ne; i--; ) |
493 |
> |
rvec[i] += c0*colrval(sd.arr.c[i],sd.type) |
494 |
> |
+ c1*colrval(sd.arr.c[i+stride],sd.type); |
495 |
|
} |
496 |
< |
sd.arr.c = dp->arr.c + i*stride; |
495 |
< |
for (i = sd.dim[0].ne; i--; ) |
496 |
< |
rvec[i] += c0*colrval(sd.arr.c[i],sd.type) |
497 |
< |
+ c1*colrval(sd.arr.c[i+stride],sd.type); |
498 |
< |
return(0.); |
496 |
> |
return(0.); /* return value ignored */ |
497 |
|
} |
498 |
|
/* get dependent variable */ |
499 |
|
if (dp->nd > 1) { |