| 343 | |
| 344 | #ifdef KR_headers |
| 345 | VOID pow_ci(p, a, b) /* p = a**b */ |
| 346 | complex *p, *a; integer *b; |
| 347 | #else |
| 348 | extern void pow_zi(doublecomplex*, doublecomplex*, integer*); |
| 349 | void pow_ci(complex *p, complex *a, integer *b) /* p = a**b */ |
| 350 | #endif |
| 351 | { |
| 352 | doublecomplex p1, a1; |
| 353 | |
| 354 | a1.r = a->r; |
| 355 | a1.i = a->i; |
| 356 | |
| 357 | pow_zi(&p1, &a1, b); |
| 358 | |
| 359 | p->r = p1.r; |
| 360 | p->i = p1.i; |
| 361 | } |
| 362 | |
| 363 | /* Unless compiled with -DNO_OVERWRITE, this variant of s_cat allows the |
| 364 | * target of a concatenation to appear on its right-hand side (contrary |
| 365 | * to the Fortran 77 Standard, but in accordance with Fortran 90). |
| 366 | */ |
| 367 | #define NO_OVERWRITE |
| 368 | |
| 369 | |
| 370 | #ifndef NO_OVERWRITE |
| 371 | |
| 372 | #undef abs |
| 373 | #ifdef KR_headers |
| 374 | extern char *F77_aloc(); |
| 375 | extern void free(); |
| 376 | extern void exit_(); |
| 377 | #else |
| 378 | |
| 379 | extern char *F77_aloc(ftnlen, char*); |
| 380 | #endif |
| 381 | |
| 382 | #endif /* NO_OVERWRITE */ |
| 383 | |
| 384 | #ifdef KR_headers |
| 385 | int |
| 386 | s_cat(lp, rpp, rnp, np, ll) char *lp, *rpp[]; ftnlen rnp[], *np, ll; |
| 387 | #else |
| 388 | int |
| 389 | s_cat(char *lp, char *rpp[], ftnlen rnp[], ftnlen *np, ftnlen ll) |
| 390 | #endif |
| 391 | { |
| 392 | ftnlen i, nc; |
| 393 | char *rp; |
| 394 | ftnlen n = *np; |
| 395 | #ifndef NO_OVERWRITE |
| 396 | ftnlen L, m; |
| 397 | char *lp0, *lp1; |
| 398 | |
| 399 | lp0 = 0; |
| 400 | lp1 = lp; |
| 401 | L = ll; |
| 402 | i = 0; |