Changeset 756 in Sophya for trunk/SophyaLib


Ignore:
Timestamp:
Mar 2, 2000, 5:22:08 PM (26 years ago)
Author:
ansari
Message:

Restructuration de Sophya en petits modules - Petites corrections diverses

Reza 2/3/2000

Location:
trunk/SophyaLib/NTools
Files:
2 added
27 deleted
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/SophyaLib/NTools/NTools.o.list

    r750 r756  
    1 arrfitter.o
    21cimage.o
    32cspline.o
    4 cvector.o
    5 datacards.o
    63datatypes.o
    74dates.o
    85datime.o
    96difeq.o
    10 dvlist.o
    117dynccd.o
    128fct1dfit.o
    … …  
    2420integ.o
    2521linfit.o
    26 matrix.o
    2722matxop.o
    2823median.o
    … …  
    3025nbrandom.o
    3126nbtri.o
     27ntoolsinit.o
    3228ntupintf.o
    33 outilsinit.o
     29objfitter.o
    3430perandom.o
    3531poly.o
    … …  
    3733rzimage.o
    3834simplesort.o
    39 strutil.o
    4035tabmath.o
    4136utils.o

  • trunk/SophyaLib/NTools/fftpackc.c

    r753 r756  
    1111
    1212/*  ------ File cfftb.f ------ */
    13 /* Subroutine */ int cfftb_(integer *n, r_4 *c__, r_4 *wsave)
    14 {
    15     extern /* Subroutine */ int cfftb1_(integer *, r_4 *, r_4 *, r_4 *,
     13/* Subroutine */ int cfftb_(int_4 *n, r_4 *c__, r_4 *wsave)
     14{
     15    extern /* Subroutine */ int cfftb1_(int_4 *, r_4 *, r_4 *, r_4 *,
    1616            int_4 *);
    17     static integer iw1, iw2;
     17    static int_4 iw1, iw2;
    1818
    1919    /* Parameter adjustments */
    … …  
    3232
    3333/*  ------ File cfftb1.f ------ */
    34 /* Subroutine */ int cfftb1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa,
     34/* Subroutine */ int cfftb1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa,
    3535        int_4 *ifac)
    3636{
    3737    /* System generated locals */
    38     integer i__1;
     38    int_4 i__1;
    3939
    4040    /* Local variables */
    41     static integer idot, i__;
    42     extern /* Subroutine */ int passb_(integer *, integer *, integer *,
    43             integer *, integer *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *,
     41    static int_4 idot, i__;
     42    extern /* Subroutine */ int passb_(int_4 *, int_4 *, int_4 *,
     43            int_4 *, int_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *,
    4444            r_4 *);
    45     static integer k1, l1, l2, n2;
    46     extern /* Subroutine */ int passb2_(integer *, integer *, r_4 *, r_4 *,
    47             r_4 *), passb3_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    48             r_4 *), passb4_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    49             r_4 *, r_4 *), passb5_(integer *, integer *, r_4 *, r_4 *,
     45    static int_4 k1, l1, l2, n2;
     46    extern /* Subroutine */ int passb2_(int_4 *, int_4 *, r_4 *, r_4 *,
     47            r_4 *), passb3_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     48            r_4 *), passb4_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     49            r_4 *, r_4 *), passb5_(int_4 *, int_4 *, r_4 *, r_4 *,
    5050            r_4 *, r_4 *, r_4 *, r_4 *);
    51     static integer na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
     51    static int_4 na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
    5252
    5353    /* Parameter adjustments */
    … …  
    164164
    165165/*  ------ File cfftf.f ------ */
    166 /* Subroutine */ int cfftf_(integer *n, r_4 *c__, r_4 *wsave)
    167 {
    168     extern /* Subroutine */ int cfftf1_(integer *, r_4 *, r_4 *, r_4 *,
     166/* Subroutine */ int cfftf_(int_4 *n, r_4 *c__, r_4 *wsave)
     167{
     168    extern /* Subroutine */ int cfftf1_(int_4 *, r_4 *, r_4 *, r_4 *,
    169169            int_4 *);
    170     static integer iw1, iw2;
     170    static int_4 iw1, iw2;
    171171
    172172    /* Parameter adjustments */
    … …  
    185185
    186186/*  ------ File cfftf1.f ------ */
    187 /* Subroutine */ int cfftf1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa,
     187/* Subroutine */ int cfftf1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa,
    188188        int_4 *ifac)
    189189{
    190190    /* System generated locals */
    191     integer i__1;
     191    int_4 i__1;
    192192
    193193    /* Local variables */
    194     static integer idot, i__;
    195     extern /* Subroutine */ int passf_(integer *, integer *, integer *,
    196             integer *, integer *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *,
     194    static int_4 idot, i__;
     195    extern /* Subroutine */ int passf_(int_4 *, int_4 *, int_4 *,
     196            int_4 *, int_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *,
    197197            r_4 *);
    198     static integer k1, l1, l2, n2;
    199     extern /* Subroutine */ int passf2_(integer *, integer *, r_4 *, r_4 *,
    200             r_4 *), passf3_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    201             r_4 *), passf4_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    202             r_4 *, r_4 *), passf5_(integer *, integer *, r_4 *, r_4 *,
     198    static int_4 k1, l1, l2, n2;
     199    extern /* Subroutine */ int passf2_(int_4 *, int_4 *, r_4 *, r_4 *,
     200            r_4 *), passf3_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     201            r_4 *), passf4_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     202            r_4 *, r_4 *), passf5_(int_4 *, int_4 *, r_4 *, r_4 *,
    203203            r_4 *, r_4 *, r_4 *, r_4 *);
    204     static integer na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
     204    static int_4 na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
    205205
    206206    /* Parameter adjustments */
    … …  
    317317
    318318/*  ------ File cffti.f ------ */
    319 /* Subroutine */ int cffti_(integer *n, r_4 *wsave)
    320 {
    321     extern /* Subroutine */ int cffti1_(integer *, r_4 *, int_4 *);
    322     static integer iw1, iw2;
     319/* Subroutine */ int cffti_(int_4 *n, r_4 *wsave)
     320{
     321    extern /* Subroutine */ int cffti1_(int_4 *, r_4 *, int_4 *);
     322    static int_4 iw1, iw2;
    323323
    324324    /* Parameter adjustments */
    … …  
    336336
    337337/*  ------ File cffti1.f ------ */
    338 /* Subroutine */ int cffti1_(integer *n, r_4 *wa, int_4 *ifac)
     338/* Subroutine */ int cffti1_(int_4 *n, r_4 *wa, int_4 *ifac)
    339339{
    340340    /* Initialized data */
    341341
    342     static integer ntryh[4] = { 3,4,2,5 };
     342    static int_4 ntryh[4] = { 3,4,2,5 };
    343343
    344344    /* System generated locals */
    345     integer i__1, i__2, i__3;
     345    int_4 i__1, i__2, i__3;
    346346
    347347    /* Builtin functions */
    … …  
    350350    /* Local variables */
    351351    static r_4 argh;
    352     static integer idot, ntry, i__, j;
     352    static int_4 idot, ntry, i__, j;
    353353    static r_4 argld;
    354     static integer i1, k1, l1, l2, ib;
     354    static int_4 i1, k1, l1, l2, ib;
    355355    static r_4 fi;
    356     static integer ld, ii, nf, ip, nl, nq, nr;
     356    static int_4 ld, ii, nf, ip, nl, nq, nr;
    357357    static r_4 arg;
    358     static integer ido, ipm;
     358    static int_4 ido, ipm;
    359359    static r_4 tpi;
    360360
    … …  
    454454
    455455/*  ------ File cosqb.f ------ */
    456 /* Subroutine */ int cosqb_(integer *n, r_4 *x, r_4 *wsave)
     456/* Subroutine */ int cosqb_(int_4 *n, r_4 *x, r_4 *wsave)
    457457{
    458458    /* Initialized data */
    … …  
    461461
    462462    /* System generated locals */
    463     integer i__1;
     463    int_4 i__1;
    464464
    465465    /* Local variables */
    466466    static r_4 x1;
    467     extern /* Subroutine */ int cosqb1_(integer *, r_4 *, r_4 *, r_4 *);
     467    extern /* Subroutine */ int cosqb1_(int_4 *, r_4 *, r_4 *, r_4 *);
    468468
    469469    /* Parameter adjustments */
    … …  
    493493
    494494/*  ------ File cosqb1.f ------ */
    495 /* Subroutine */ int cosqb1_(integer *n, r_4 *x, r_4 *w, r_4 *xh)
     495/* Subroutine */ int cosqb1_(int_4 *n, r_4 *x, r_4 *w, r_4 *xh)
    496496{
    497497    /* System generated locals */
    498     integer i__1;
     498    int_4 i__1;
    499499
    500500    /* Local variables */
    501     static integer modn, i__, k;
    502     extern /* Subroutine */ int rfftb_(integer *, r_4 *, r_4 *);
    503     static integer kc, np2, ns2;
     501    static int_4 modn, i__, k;
     502    extern /* Subroutine */ int rfftb_(int_4 *, r_4 *, r_4 *);
     503    static int_4 kc, np2, ns2;
    504504    static r_4 xim1;
    505505
    … …  
    547547
    548548/*  ------ File cosqf.f ------ */
    549 /* Subroutine */ int cosqf_(integer *n, r_4 *x, r_4 *wsave)
     549/* Subroutine */ int cosqf_(int_4 *n, r_4 *x, r_4 *wsave)
    550550{
    551551    /* Initialized data */
    … …  
    554554
    555555    /* System generated locals */
    556     integer i__1;
     556    int_4 i__1;
    557557
    558558    /* Local variables */
    559559    static r_4 tsqx;
    560     extern /* Subroutine */ int cosqf1_(integer *, r_4 *, r_4 *, r_4 *);
     560    extern /* Subroutine */ int cosqf1_(int_4 *, r_4 *, r_4 *, r_4 *);
    561561
    562562    /* Parameter adjustments */
    … …  
    584584
    585585/*  ------ File cosqf1.f ------ */
    586 /* Subroutine */ int cosqf1_(integer *n, r_4 *x, r_4 *w, r_4 *xh)
     586/* Subroutine */ int cosqf1_(int_4 *n, r_4 *x, r_4 *w, r_4 *xh)
    587587{
    588588    /* System generated locals */
    589     integer i__1;
     589    int_4 i__1;
    590590
    591591    /* Local variables */
    592     static integer modn, i__, k;
    593     extern /* Subroutine */ int rfftf_(integer *, r_4 *, r_4 *);
    594     static integer kc, np2, ns2;
     592    static int_4 modn, i__, k;
     593    extern /* Subroutine */ int rfftf_(int_4 *, r_4 *, r_4 *);
     594    static int_4 kc, np2, ns2;
    595595    static r_4 xim1;
    596596
    … …  
    636636
    637637/*  ------ File cosqi.f ------ */
    638 /* Subroutine */ int cosqi_(integer *n, r_4 *wsave)
     638/* Subroutine */ int cosqi_(int_4 *n, r_4 *wsave)
    639639{
    640640    /* Initialized data */
    … …  
    643643
    644644    /* System generated locals */
    645     integer i__1;
     645    int_4 i__1;
    646646
    647647    /* Builtin functions */
    … …  
    649649
    650650    /* Local variables */
    651     static integer k;
    652     extern /* Subroutine */ int rffti_(integer *, r_4 *);
     651    static int_4 k;
     652    extern /* Subroutine */ int rffti_(int_4 *, r_4 *);
    653653    static r_4 fk, dt;
    654654
    … …  
    670670
    671671/*  ------ File cost.f ------ */
    672 /* Subroutine */ int cost_(integer *n, r_4 *x, r_4 *wsave)
     672/* Subroutine */ int cost_(int_4 *n, r_4 *x, r_4 *wsave)
    673673{
    674674    /* System generated locals */
    675     integer i__1;
     675    int_4 i__1;
    676676
    677677    /* Local variables */
    678     static integer modn, i__, k;
    679     extern /* Subroutine */ int rfftf_(integer *, r_4 *, r_4 *);
     678    static int_4 modn, i__, k;
     679    extern /* Subroutine */ int rfftf_(int_4 *, r_4 *, r_4 *);
    680680    static r_4 c1, t1, t2;
    681     static integer kc;
     681    static int_4 kc;
    682682    static r_4 xi;
    683     static integer nm1, np1;
     683    static int_4 nm1, np1;
    684684    static r_4 x1h;
    685     static integer ns2;
     685    static int_4 ns2;
    686686    static r_4 tx2, x1p3, xim2;
    687687
    … …  
    753753
    754754/*  ------ File costi.f ------ */
    755 /* Subroutine */ int costi_(integer *n, r_4 *wsave)
     755/* Subroutine */ int costi_(int_4 *n, r_4 *wsave)
    756756{
    757757    /* Initialized data */
    … …  
    760760
    761761    /* System generated locals */
    762     integer i__1;
     762    int_4 i__1;
    763763
    764764    /* Builtin functions */
    … …  
    766766
    767767    /* Local variables */
    768     static integer k;
    769     extern /* Subroutine */ int rffti_(integer *, r_4 *);
    770     static integer kc;
     768    static int_4 k;
     769    extern /* Subroutine */ int rffti_(int_4 *, r_4 *);
     770    static int_4 kc;
    771771    static r_4 fk, dt;
    772     static integer nm1, np1, ns2;
     772    static int_4 nm1, np1, ns2;
    773773
    774774    /* Parameter adjustments */
    … …  
    797797
    798798/*  ------ File ezfft1.f ------ */
    799 /* Subroutine */ int ezfft1_(integer *n, r_4 *wa, int_4 *ifac)
     799/* Subroutine */ int ezfft1_(int_4 *n, r_4 *wa, int_4 *ifac)
    800800{
    801801    /* Initialized data */
    802802
    803     static integer ntryh[4] = { 4,2,3,5 };
     803    static int_4 ntryh[4] = { 4,2,3,5 };
    804804    static r_4 tpi = 6.28318530717959f;
    805805
    806806    /* System generated locals */
    807     integer i__1, i__2, i__3;
     807    int_4 i__1, i__2, i__3;
    808808
    809809    /* Builtin functions */
    … …  
    812812    /* Local variables */
    813813    static r_4 argh;
    814     static integer ntry, i__, j, k1, l1, l2, ib, ii, nf, ip, nl, is, nq, nr;
     814    static int_4 ntry, i__, j, k1, l1, l2, ib, ii, nf, ip, nl, is, nq, nr;
    815815    static r_4 ch1, sh1;
    816     static integer ido, ipm;
     816    static int_4 ido, ipm;
    817817    static r_4 dch1, ch1h, arg1, dsh1;
    818     static integer nfm1;
     818    static int_4 nfm1;
    819819
    820820    /* Parameter adjustments */
    … …  
    916916
    917917/*  ------ File ezfftb.f ------ */
    918 /* Subroutine */ int ezfftb_(integer *n, r_4 *r__, r_4 *azero, r_4 *a,
     918/* Subroutine */ int ezfftb_(int_4 *n, r_4 *r__, r_4 *azero, r_4 *a,
    919919        r_4 *b, r_4 *wsave)
    920920{
    921921    /* System generated locals */
    922     integer i__1;
     922    int_4 i__1;
    923923
    924924    /* Local variables */
    925     static integer i__;
    926     extern /* Subroutine */ int rfftb_(integer *, r_4 *, r_4 *);
    927     static integer ns2;
     925    static int_4 i__;
     926    extern /* Subroutine */ int rfftb_(int_4 *, r_4 *, r_4 *);
     927    static int_4 ns2;
    928928
    929929    /* Parameter adjustments */
    … …  
    965965
    966966/*  ------ File ezfftf.f ------ */
    967 /* Subroutine */ int ezfftf_(integer *n, r_4 *r__, r_4 *azero, r_4 *a,
     967/* Subroutine */ int ezfftf_(int_4 *n, r_4 *r__, r_4 *azero, r_4 *a,
    968968        r_4 *b, r_4 *wsave)
    969969{
    970970    /* System generated locals */
    971     integer i__1;
     971    int_4 i__1;
    972972
    973973    /* Local variables */
    974     static integer i__;
    975     extern /* Subroutine */ int rfftf_(integer *, r_4 *, r_4 *);
     974    static int_4 i__;
     975    extern /* Subroutine */ int rfftf_(int_4 *, r_4 *, r_4 *);
    976976    static r_4 cf;
    977     static integer ns2;
     977    static int_4 ns2;
    978978    static r_4 cfm;
    979     static integer ns2m;
     979    static int_4 ns2m;
    980980
    981981
    … …  
    10301030
    10311031/*  ------ File ezffti.f ------ */
    1032 /* Subroutine */ int ezffti_(integer *n, r_4 *wsave)
    1033 {
    1034     extern /* Subroutine */ int ezfft1_(integer *, r_4 *, int_4 *);
     1032/* Subroutine */ int ezffti_(int_4 *n, r_4 *wsave)
     1033{
     1034    extern /* Subroutine */ int ezfft1_(int_4 *, r_4 *, int_4 *);
    10351035
    10361036    /* Parameter adjustments */
    … …  
    10461046
    10471047/*  ------ File passb.f ------ */
    1048 /* Subroutine */ int passb_(integer *nac, integer *ido, integer *ip, integer *
    1049         l1, integer *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2,
     1048/* Subroutine */ int passb_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *
     1049        l1, int_4 *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2,
    10501050        r_4 *wa)
    10511051{
    10521052    /* System generated locals */
    1053     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     1053    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    10541054             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    10551055            i__1, i__2, i__3;
    10561056
    10571057    /* Local variables */
    1058     static integer idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
     1058    static int_4 idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
    10591059            idl, inc, idp;
    10601060    static r_4 wai, war;
    1061     static integer ipp2;
     1061    static int_4 ipp2;
    10621062
    10631063    /* Parameter adjustments */
    … …  
    12871287
    12881288/*  ------ File passb2.f ------ */
    1289 /* Subroutine */ int passb2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     1289/* Subroutine */ int passb2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    12901290        r_4 *wa1)
    12911291{
    12921292    /* System generated locals */
    1293     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     1293    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    12941294
    12951295    /* Local variables */
    1296     static integer i__, k;
     1296    static int_4 i__, k;
    12971297    static r_4 ti2, tr2;
    12981298
    … …  
    13491349
    13501350/*  ------ File passb3.f ------ */
    1351 /* Subroutine */ int passb3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     1351/* Subroutine */ int passb3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    13521352        r_4 *wa1, r_4 *wa2)
    13531353{
    … …  
    13581358
    13591359    /* System generated locals */
    1360     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     1360    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    13611361
    13621362    /* Local variables */
    1363     static integer i__, k;
     1363    static int_4 i__, k;
    13641364    static r_4 ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    13651365
    … …  
    14371437
    14381438/*  ------ File passb4.f ------ */
    1439 /* Subroutine */ int passb4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     1439/* Subroutine */ int passb4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    14401440        r_4 *wa1, r_4 *wa2, r_4 *wa3)
    14411441{
    14421442    /* System generated locals */
    1443     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     1443    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    14441444
    14451445    /* Local variables */
    1446     static integer i__, k;
     1446    static int_4 i__, k;
    14471447    static r_4 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    14481448            tr3, tr4;
    … …  
    15421542
    15431543/*  ------ File passb5.f ------ */
    1544 /* Subroutine */ int passb5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     1544/* Subroutine */ int passb5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    15451545        r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4)
    15461546{
    … …  
    15531553
    15541554    /* System generated locals */
    1555     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     1555    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    15561556
    15571557    /* Local variables */
    1558     static integer i__, k;
     1558    static int_4 i__, k;
    15591559    static r_4 ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
    15601560            ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;
    … …  
    16771677
    16781678/*  ------ File passf.f ------ */
    1679 /* Subroutine */ int passf_(integer *nac, integer *ido, integer *ip, integer *
    1680         l1, integer *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2,
     1679/* Subroutine */ int passf_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *
     1680        l1, int_4 *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2,
    16811681        r_4 *wa)
    16821682{
    16831683    /* System generated locals */
    1684     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     1684    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    16851685             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    16861686            i__1, i__2, i__3;
    16871687
    16881688    /* Local variables */
    1689     static integer idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
     1689    static int_4 idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
    16901690            idl, inc, idp;
    16911691    static r_4 wai, war;
    1692     static integer ipp2;
     1692    static int_4 ipp2;
    16931693
    16941694    /* Parameter adjustments */
    … …  
    19181918
    19191919/*  ------ File passf2.f ------ */
    1920 /* Subroutine */ int passf2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     1920/* Subroutine */ int passf2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    19211921        r_4 *wa1)
    19221922{
    19231923    /* System generated locals */
    1924     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     1924    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    19251925
    19261926    /* Local variables */
    1927     static integer i__, k;
     1927    static int_4 i__, k;
    19281928    static r_4 ti2, tr2;
    19291929
    … …  
    19801980
    19811981/*  ------ File passf3.f ------ */
    1982 /* Subroutine */ int passf3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     1982/* Subroutine */ int passf3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    19831983        r_4 *wa1, r_4 *wa2)
    19841984{
    … …  
    19891989
    19901990    /* System generated locals */
    1991     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     1991    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    19921992
    19931993    /* Local variables */
    1994     static integer i__, k;
     1994    static int_4 i__, k;
    19951995    static r_4 ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    19961996
    … …  
    20682068
    20692069/*  ------ File passf4.f ------ */
    2070 /* Subroutine */ int passf4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     2070/* Subroutine */ int passf4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    20712071        r_4 *wa1, r_4 *wa2, r_4 *wa3)
    20722072{
    20732073    /* System generated locals */
    2074     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     2074    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    20752075
    20762076    /* Local variables */
    2077     static integer i__, k;
     2077    static int_4 i__, k;
    20782078    static r_4 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    20792079            tr3, tr4;
    … …  
    21732173
    21742174/*  ------ File passf5.f ------ */
    2175 /* Subroutine */ int passf5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     2175/* Subroutine */ int passf5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    21762176        r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4)
    21772177{
    … …  
    21842184
    21852185    /* System generated locals */
    2186     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     2186    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    21872187
    21882188    /* Local variables */
    2189     static integer i__, k;
     2189    static int_4 i__, k;
    21902190    static r_4 ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
    21912191            ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;
    … …  
    23082308
    23092309/*  ------ File radb2.f ------ */
    2310 /* Subroutine */ int radb2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     2310/* Subroutine */ int radb2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    23112311        r_4 *wa1)
    23122312{
    23132313    /* System generated locals */
    2314     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     2314    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    23152315
    23162316    /* Local variables */
    2317     static integer i__, k, ic;
     2317    static int_4 i__, k, ic;
    23182318    static r_4 ti2, tr2;
    2319     static integer idp2;
     2319    static int_4 idp2;
    23202320
    23212321    /* Parameter adjustments */
    … …  
    23852385
    23862386/*  ------ File radb3.f ------ */
    2387 /* Subroutine */ int radb3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     2387/* Subroutine */ int radb3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    23882388        r_4 *wa1, r_4 *wa2)
    23892389{
    … …  
    23942394
    23952395    /* System generated locals */
    2396     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     2396    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    23972397
    23982398    /* Local variables */
    2399     static integer i__, k, ic;
     2399    static int_4 i__, k, ic;
    24002400    static r_4 ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    2401     static integer idp2;
     2401    static int_4 idp2;
    24022402
    24032403    /* Parameter adjustments */
    … …  
    24682468
    24692469/*  ------ File radb4.f ------ */
    2470 /* Subroutine */ int radb4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     2470/* Subroutine */ int radb4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    24712471        r_4 *wa1, r_4 *wa2, r_4 *wa3)
    24722472{
    … …  
    24762476
    24772477    /* System generated locals */
    2478     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     2478    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    24792479
    24802480    /* Local variables */
    2481     static integer i__, k, ic;
     2481    static int_4 i__, k, ic;
    24822482    static r_4 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    24832483            tr3, tr4;
    2484     static integer idp2;
     2484    static int_4 idp2;
    24852485
    24862486    /* Parameter adjustments */
    … …  
    25922592
    25932593/*  ------ File radb5.f ------ */
    2594 /* Subroutine */ int radb5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     2594/* Subroutine */ int radb5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    25952595        r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4)
    25962596{
    … …  
    26032603
    26042604    /* System generated locals */
    2605     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     2605    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    26062606
    26072607    /* Local variables */
    2608     static integer i__, k, ic;
     2608    static int_4 i__, k, ic;
    26092609    static r_4 ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
    26102610            ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;
    2611     static integer idp2;
     2611    static int_4 idp2;
    26122612
    26132613    /* Parameter adjustments */
    … …  
    27162716
    27172717/*  ------ File radbg.f ------ */
    2718 /* Subroutine */ int radbg_(integer *ido, integer *ip, integer *l1, integer *
     2718/* Subroutine */ int radbg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *
    27192719        idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa)
    27202720{
    … …  
    27242724
    27252725    /* System generated locals */
    2726     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     2726    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    27272727             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    27282728            i__1, i__2, i__3;
    … …  
    27322732
    27332733    /* Local variables */
    2734     static integer idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
     2734    static int_4 idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
    27352735    static r_4 dc2, ai1, ai2, ar1, ar2, ds2;
    2736     static integer nbd;
     2736    static int_4 nbd;
    27372737    static r_4 dcp, arg, dsp, ar1h, ar2h;
    2738     static integer idp2, ipp2;
     2738    static int_4 idp2, ipp2;
    27392739
    27402740    /* Parameter adjustments */
    … …  
    30573057
    30583058/*  ------ File radf2.f ------ */
    3059 /* Subroutine */ int radf2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     3059/* Subroutine */ int radf2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    30603060        r_4 *wa1)
    30613061{
    30623062    /* System generated locals */
    3063     integer ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
     3063    int_4 ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
    30643064
    30653065    /* Local variables */
    3066     static integer i__, k, ic;
     3066    static int_4 i__, k, ic;
    30673067    static r_4 ti2, tr2;
    3068     static integer idp2;
     3068    static int_4 idp2;
    30693069
    30703070    /* Parameter adjustments */
    … …  
    31353135
    31363136/*  ------ File radf3.f ------ */
    3137 /* Subroutine */ int radf3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     3137/* Subroutine */ int radf3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    31383138        r_4 *wa1, r_4 *wa2)
    31393139{
    … …  
    31443144
    31453145    /* System generated locals */
    3146     integer ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
     3146    int_4 ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
    31473147
    31483148    /* Local variables */
    3149     static integer i__, k, ic;
     3149    static int_4 i__, k, ic;
    31503150    static r_4 ci2, di2, di3, cr2, dr2, dr3, ti2, ti3, tr2, tr3;
    3151     static integer idp2;
     3151    static int_4 idp2;
    31523152
    31533153    /* Parameter adjustments */
    … …  
    32143214
    32153215/*  ------ File radf4.f ------ */
    3216 /* Subroutine */ int radf4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     3216/* Subroutine */ int radf4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    32173217        r_4 *wa1, r_4 *wa2, r_4 *wa3)
    32183218{
    … …  
    32223222
    32233223    /* System generated locals */
    3224     integer cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
     3224    int_4 cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
    32253225
    32263226    /* Local variables */
    3227     static integer i__, k, ic;
     3227    static int_4 i__, k, ic;
    32283228    static r_4 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    32293229            tr3, tr4;
    3230     static integer idp2;
     3230    static int_4 idp2;
    32313231
    32323232    /* Parameter adjustments */
    … …  
    33303330
    33313331/*  ------ File radf5.f ------ */
    3332 /* Subroutine */ int radf5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch,
     3332/* Subroutine */ int radf5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch,
    33333333        r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4)
    33343334{
    … …  
    33413341
    33423342    /* System generated locals */
    3343     integer cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
     3343    int_4 cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
    33443344
    33453345    /* Local variables */
    3346     static integer i__, k, ic;
     3346    static int_4 i__, k, ic;
    33473347    static r_4 ci2, di2, ci4, ci5, di3, di4, di5, ci3, cr2, cr3, dr2, dr3,
    33483348            dr4, dr5, cr5, cr4, ti2, ti3, ti5, ti4, tr2, tr3, tr4, tr5;
    3349     static integer idp2;
     3349    static int_4 idp2;
    33503350
    33513351    /* Parameter adjustments */
    … …  
    34483448
    34493449/*  ------ File radfg.f ------ */
    3450 /* Subroutine */ int radfg_(integer *ido, integer *ip, integer *l1, integer *
     3450/* Subroutine */ int radfg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *
    34513451        idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa)
    34523452{
    … …  
    34563456
    34573457    /* System generated locals */
    3458     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     3458    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    34593459             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    34603460            i__1, i__2, i__3;
    … …  
    34643464
    34653465    /* Local variables */
    3466     static integer idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
     3466    static int_4 idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
    34673467    static r_4 dc2, ai1, ai2, ar1, ar2, ds2;
    3468     static integer nbd;
     3468    static int_4 nbd;
    34693469    static r_4 dcp, arg, dsp, ar1h, ar2h;
    3470     static integer idp2, ipp2;
     3470    static int_4 idp2, ipp2;
    34713471
    34723472    /* Parameter adjustments */
    … …  
    37953795
    37963796/*  ------ File rfftb.f ------ */
    3797 /* Subroutine */ int rfftb_(integer *n, r_4 *r__, r_4 *wsave)
    3798 {
    3799     extern /* Subroutine */ int rfftb1_(integer *, r_4 *, r_4 *, r_4 *,
     3797/* Subroutine */ int rfftb_(int_4 *n, r_4 *r__, r_4 *wsave)
     3798{
     3799    extern /* Subroutine */ int rfftb1_(int_4 *, r_4 *, r_4 *, r_4 *,
    38003800            int_4 *);
    38013801
    … …  
    38133813
    38143814/*  ------ File rfftb1.f ------ */
    3815 /* Subroutine */ int rfftb1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa,
     3815/* Subroutine */ int rfftb1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa,
    38163816        int_4 *ifac)
    38173817{
    38183818    /* System generated locals */
    3819     integer i__1;
     3819    int_4 i__1;
    38203820
    38213821    /* Local variables */
    3822     extern /* Subroutine */ int radb2_(integer *, integer *, r_4 *, r_4 *,
    3823             r_4 *), radb3_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    3824             r_4 *), radb4_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    3825             r_4 *, r_4 *), radb5_(integer *, integer *, r_4 *, r_4 *,
     3822    extern /* Subroutine */ int radb2_(int_4 *, int_4 *, r_4 *, r_4 *,
     3823            r_4 *), radb3_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     3824            r_4 *), radb4_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     3825            r_4 *, r_4 *), radb5_(int_4 *, int_4 *, r_4 *, r_4 *,
    38263826            r_4 *, r_4 *, r_4 *, r_4 *);
    3827     static integer i__;
    3828     extern /* Subroutine */ int radbg_(integer *, integer *, integer *,
    3829             integer *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *);
    3830     static integer k1, l1, l2, na, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
     3827    static int_4 i__;
     3828    extern /* Subroutine */ int radbg_(int_4 *, int_4 *, int_4 *,
     3829            int_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *);
     3830    static int_4 k1, l1, l2, na, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
    38313831
    38323832    /* Parameter adjustments */
    … …  
    39413941
    39423942/*  ------ File rfftf.f ------ */
    3943 /* Subroutine */ int rfftf_(integer *n, r_4 *r__, r_4 *wsave)
    3944 {
    3945     extern /* Subroutine */ int rfftf1_(integer *, r_4 *, r_4 *, r_4 *,
     3943/* Subroutine */ int rfftf_(int_4 *n, r_4 *r__, r_4 *wsave)
     3944{
     3945    extern /* Subroutine */ int rfftf1_(int_4 *, r_4 *, r_4 *, r_4 *,
    39463946            int_4 *);
    39473947
    … …  
    39593959
    39603960/*  ------ File rfftf1.f ------ */
    3961 /* Subroutine */ int rfftf1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa,
     3961/* Subroutine */ int rfftf1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa,
    39623962        int_4 *ifac)
    39633963{
    39643964    /* System generated locals */
    3965     integer i__1;
     3965    int_4 i__1;
    39663966
    39673967    /* Local variables */
    3968     extern /* Subroutine */ int radf2_(integer *, integer *, r_4 *, r_4 *,
    3969             r_4 *), radf3_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    3970             r_4 *), radf4_(integer *, integer *, r_4 *, r_4 *, r_4 *,
    3971             r_4 *, r_4 *), radf5_(integer *, integer *, r_4 *, r_4 *,
     3968    extern /* Subroutine */ int radf2_(int_4 *, int_4 *, r_4 *, r_4 *,
     3969            r_4 *), radf3_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     3970            r_4 *), radf4_(int_4 *, int_4 *, r_4 *, r_4 *, r_4 *,
     3971            r_4 *, r_4 *), radf5_(int_4 *, int_4 *, r_4 *, r_4 *,
    39723972            r_4 *, r_4 *, r_4 *, r_4 *);
    3973     static integer i__;
    3974     extern /* Subroutine */ int radfg_(integer *, integer *, integer *,
    3975             integer *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *);
    3976     static integer k1, l1, l2, na, kh, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
     3973    static int_4 i__;
     3974    extern /* Subroutine */ int radfg_(int_4 *, int_4 *, int_4 *,
     3975            int_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *, r_4 *);
     3976    static int_4 k1, l1, l2, na, kh, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
    39773977
    39783978    /* Parameter adjustments */
    … …  
    40824082
    40834083/*  ------ File rffti.f ------ */
    4084 /* Subroutine */ int rffti_(integer *n, r_4 *wsave)
    4085 {
    4086     extern /* Subroutine */ int rffti1_(integer *, r_4 *, int_4 *);
     4084/* Subroutine */ int rffti_(int_4 *n, r_4 *wsave)
     4085{
     4086    extern /* Subroutine */ int rffti1_(int_4 *, r_4 *, int_4 *);
    40874087
    40884088    /* Parameter adjustments */
    … …  
    40984098
    40994099/*  ------ File rffti1.f ------ */
    4100 /* Subroutine */ int rffti1_(integer *n, r_4 *wa, int_4 *ifac)
     4100/* Subroutine */ int rffti1_(int_4 *n, r_4 *wa, int_4 *ifac)
    41014101{
    41024102    /* Initialized data */
    41034103
    4104     static integer ntryh[4] = { 4,2,3,5 };
     4104    static int_4 ntryh[4] = { 4,2,3,5 };
    41054105
    41064106    /* System generated locals */
    4107     integer i__1, i__2, i__3;
     4107    int_4 i__1, i__2, i__3;
    41084108
    41094109    /* Builtin functions */
    … …  
    41124112    /* Local variables */
    41134113    static r_4 argh;
    4114     static integer ntry, i__, j;
     4114    static int_4 ntry, i__, j;
    41154115    static r_4 argld;
    4116     static integer k1, l1, l2, ib;
     4116    static int_4 k1, l1, l2, ib;
    41174117    static r_4 fi;
    4118     static integer ld, ii, nf, ip, nl, is, nq, nr;
     4118    static int_4 ld, ii, nf, ip, nl, is, nq, nr;
    41194119    static r_4 arg;
    4120     static integer ido, ipm;
     4120    static int_4 ido, ipm;
    41214121    static r_4 tpi;
    4122     static integer nfm1;
     4122    static int_4 nfm1;
    41234123
    41244124    /* Parameter adjustments */
    … …  
    42134213
    42144214/*  ------ File sinqb.f ------ */
    4215 /* Subroutine */ int sinqb_(integer *n, r_4 *x, r_4 *wsave)
     4215/* Subroutine */ int sinqb_(int_4 *n, r_4 *x, r_4 *wsave)
    42164216{
    42174217    /* System generated locals */
    4218     integer i__1;
     4218    int_4 i__1;
    42194219
    42204220    /* Local variables */
    4221     static integer k;
    4222     extern /* Subroutine */ int cosqb_(integer *, r_4 *, r_4 *);
     4221    static int_4 k;
     4222    extern /* Subroutine */ int cosqb_(int_4 *, r_4 *, r_4 *);
    42234223    static r_4 xhold;
    4224     static integer kc, ns2;
     4224    static int_4 kc, ns2;
    42254225
    42264226    /* Parameter adjustments */
    … …  
    42544254
    42554255/*  ------ File sinqf.f ------ */
    4256 /* Subroutine */ int sinqf_(integer *n, r_4 *x, r_4 *wsave)
     4256/* Subroutine */ int sinqf_(int_4 *n, r_4 *x, r_4 *wsave)
    42574257{
    42584258    /* System generated locals */
    4259     integer i__1;
     4259    int_4 i__1;
    42604260
    42614261    /* Local variables */
    4262     static integer k;
    4263     extern /* Subroutine */ int cosqf_(integer *, r_4 *, r_4 *);
     4262    static int_4 k;
     4263    extern /* Subroutine */ int cosqf_(int_4 *, r_4 *, r_4 *);
    42644264    static r_4 xhold;
    4265     static integer kc, ns2;
     4265    static int_4 kc, ns2;
    42664266
    42674267    /* Parameter adjustments */
    … …  
    42924292
    42934293/*  ------ File sinqi.f ------ */
    4294 /* Subroutine */ int sinqi_(integer *n, r_4 *wsave)
    4295 {
    4296     extern /* Subroutine */ int cosqi_(integer *, r_4 *);
     4294/* Subroutine */ int sinqi_(int_4 *n, r_4 *wsave)
     4295{
     4296    extern /* Subroutine */ int cosqi_(int_4 *, r_4 *);
    42974297
    42984298    /* Parameter adjustments */
    … …  
    43054305
    43064306/*  ------ File sint.f ------ */
    4307 /* Subroutine */ int sint_(integer *n, r_4 *x, r_4 *wsave)
    4308 {
    4309     extern /* Subroutine */ int sint1_(integer *, r_4 *, r_4 *, r_4 *,
     4307/* Subroutine */ int sint_(int_4 *n, r_4 *x, r_4 *wsave)
     4308{
     4309    extern /* Subroutine */ int sint1_(int_4 *, r_4 *, r_4 *, r_4 *,
    43104310            r_4 *, int_4 *);
    4311     static integer np1, iw1, iw2, iw3;
     4311    static int_4 np1, iw1, iw2, iw3;
    43124312
    43134313    /* Parameter adjustments */
    … …  
    43254325
    43264326/*  ------ File sint1.f ------ */
    4327 /* Subroutine */ int sint1_(integer *n, r_4 *war, r_4 *was, r_4 *xh, r_4 *
     4327/* Subroutine */ int sint1_(int_4 *n, r_4 *war, r_4 *was, r_4 *xh, r_4 *
    43284328        x, int_4 *ifac)
    43294329{
    … …  
    43334333
    43344334    /* System generated locals */
    4335     integer i__1;
     4335    int_4 i__1;
    43364336
    43374337    /* Local variables */
    4338     static integer modn, i__, k;
     4338    static int_4 modn, i__, k;
    43394339    static r_4 xhold, t1, t2;
    4340     extern /* Subroutine */ int rfftf1_(integer *, r_4 *, r_4 *, r_4 *,
     4340    extern /* Subroutine */ int rfftf1_(int_4 *, r_4 *, r_4 *, r_4 *,
    43414341            int_4 *);
    4342     static integer kc, np1, ns2;
     4342    static int_4 kc, np1, ns2;
    43434343
    43444344    /* Parameter adjustments */
    … …  
    44114411
    44124412/*  ------ File sinti.f ------ */
    4413 /* Subroutine */ int sinti_(integer *n, r_4 *wsave)
     4413/* Subroutine */ int sinti_(int_4 *n, r_4 *wsave)
    44144414{
    44154415    /* Initialized data */
    … …  
    44184418
    44194419    /* System generated locals */
    4420     integer i__1;
     4420    int_4 i__1;
    44214421
    44224422    /* Builtin functions */
    … …  
    44244424
    44254425    /* Local variables */
    4426     static integer k;
    4427     extern /* Subroutine */ int rffti_(integer *, r_4 *);
     4426    static int_4 k;
     4427    extern /* Subroutine */ int rffti_(int_4 *, r_4 *);
    44284428    static r_4 dt;
    4429     static integer np1, ns2;
     4429    static int_4 np1, ns2;
    44304430
    44314431    /* Parameter adjustments */
    … …  
    44524452/* --------------------------------------------------------------- */
    44534453
    4454 /* Subroutine */ int cdfftb_(integer *n, r_8 *c__, r_8 *wsave)
    4455 {
    4456     extern /* Subroutine */ int cdfftb1_(integer *, r_8 *, r_8 *, r_8 *,
     4454/* Subroutine */ int cdfftb_(int_4 *n, r_8 *c__, r_8 *wsave)
     4455{
     4456    extern /* Subroutine */ int cdfftb1_(int_4 *, r_8 *, r_8 *, r_8 *,
    44574457            int_8 *);
    4458     static integer iw1, iw2;
     4458    static int_4 iw1, iw2;
    44594459
    44604460    /* Parameter adjustments */
    … …  
    44734473
    44744474/*  ------ File cdfftb1.f ------ */
    4475 /* Subroutine */ int cdfftb1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa,
     4475/* Subroutine */ int cdfftb1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa,
    44764476        int_8 *ifac)
    44774477{
    44784478    /* System generated locals */
    4479     integer i__1;
     4479    int_4 i__1;
    44804480
    44814481    /* Local variables */
    4482     static integer idot, i__;
    4483     extern /* Subroutine */ int dpassb_(integer *, integer *, integer *,
    4484             integer *, integer *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *,
     4482    static int_4 idot, i__;
     4483    extern /* Subroutine */ int dpassb_(int_4 *, int_4 *, int_4 *,
     4484            int_4 *, int_4 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *,
    44854485            r_8 *);
    4486     static integer k1, l1, l2, n2;
    4487     extern /* Subroutine */ int dpassb2_(integer *, integer *, r_8 *, r_8 *,
    4488             r_8 *), dpassb3_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    4489             r_8 *), dpassb4_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    4490             r_8 *, r_8 *), dpassb5_(integer *, integer *, r_8 *, r_8 *,
     4486    static int_4 k1, l1, l2, n2;
     4487    extern /* Subroutine */ int dpassb2_(int_4 *, int_4 *, r_8 *, r_8 *,
     4488            r_8 *), dpassb3_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     4489            r_8 *), dpassb4_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     4490            r_8 *, r_8 *), dpassb5_(int_4 *, int_4 *, r_8 *, r_8 *,
    44914491            r_8 *, r_8 *, r_8 *, r_8 *);
    4492     static integer na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
     4492    static int_4 na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
    44934493
    44944494    /* Parameter adjustments */
    … …  
    46054605
    46064606/*  ------ File cdfftf.f ------ */
    4607 /* Subroutine */ int cdfftf_(integer *n, r_8 *c__, r_8 *wsave)
    4608 {
    4609     extern /* Subroutine */ int cdfftf1_(integer *, r_8 *, r_8 *, r_8 *,
     4607/* Subroutine */ int cdfftf_(int_4 *n, r_8 *c__, r_8 *wsave)
     4608{
     4609    extern /* Subroutine */ int cdfftf1_(int_4 *, r_8 *, r_8 *, r_8 *,
    46104610            int_8 *);
    4611     static integer iw1, iw2;
     4611    static int_4 iw1, iw2;
    46124612
    46134613    /* Parameter adjustments */
    … …  
    46264626
    46274627/*  ------ File cdfftf1.f ------ */
    4628 /* Subroutine */ int cdfftf1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa,
     4628/* Subroutine */ int cdfftf1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa,
    46294629        int_8 *ifac)
    46304630{
    46314631    /* System generated locals */
    4632     integer i__1;
     4632    int_4 i__1;
    46334633
    46344634    /* Local variables */
    4635     static integer idot, i__;
    4636     extern /* Subroutine */ int dpassf_(integer *, integer *, integer *,
    4637             integer *, integer *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *,
     4635    static int_4 idot, i__;
     4636    extern /* Subroutine */ int dpassf_(int_4 *, int_4 *, int_4 *,
     4637            int_4 *, int_4 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *,
    46384638            r_8 *);
    4639     static integer k1, l1, l2, n2;
    4640     extern /* Subroutine */ int dpassf2_(integer *, integer *, r_8 *, r_8 *,
    4641             r_8 *), dpassf3_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    4642             r_8 *), dpassf4_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    4643             r_8 *, r_8 *), dpassf5_(integer *, integer *, r_8 *, r_8 *,
     4639    static int_4 k1, l1, l2, n2;
     4640    extern /* Subroutine */ int dpassf2_(int_4 *, int_4 *, r_8 *, r_8 *,
     4641            r_8 *), dpassf3_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     4642            r_8 *), dpassf4_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     4643            r_8 *, r_8 *), dpassf5_(int_4 *, int_4 *, r_8 *, r_8 *,
    46444644            r_8 *, r_8 *, r_8 *, r_8 *);
    4645     static integer na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
     4645    static int_4 na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
    46464646
    46474647    /* Parameter adjustments */
    … …  
    47584758
    47594759/*  ------ File cdffti.f ------ */
    4760 /* Subroutine */ int cdffti_(integer *n, r_8 *wsave)
    4761 {
    4762     extern /* Subroutine */ int cdffti1_(integer *, r_8 *, int_8 *);
    4763     static integer iw1, iw2;
     4760/* Subroutine */ int cdffti_(int_4 *n, r_8 *wsave)
     4761{
     4762    extern /* Subroutine */ int cdffti1_(int_4 *, r_8 *, int_8 *);
     4763    static int_4 iw1, iw2;
    47644764
    47654765    /* Parameter adjustments */
    … …  
    47774777
    47784778/*  ------ File cdffti1.f ------ */
    4779 /* Subroutine */ int cdffti1_(integer *n, r_8 *wa, int_8 *ifac)
     4779/* Subroutine */ int cdffti1_(int_4 *n, r_8 *wa, int_8 *ifac)
    47804780{
    47814781    /* Initialized data */
    47824782
    4783     static integer ntryh[4] = { 3,4,2,5 };
     4783    static int_4 ntryh[4] = { 3,4,2,5 };
    47844784
    47854785    /* System generated locals */
    4786     integer i__1, i__2, i__3;
     4786    int_4 i__1, i__2, i__3;
    47874787
    47884788    /* Builtin functions */
    … …  
    47914791    /* Local variables */
    47924792    static r_8 argh;
    4793     static integer idot, ntry, i__, j;
     4793    static int_4 idot, ntry, i__, j;
    47944794    static r_8 argld;
    4795     static integer i1, k1, l1, l2, ib;
     4795    static int_4 i1, k1, l1, l2, ib;
    47964796    static r_8 fi;
    4797     static integer ld, ii, nf, ip, nl, nq, nr;
     4797    static int_4 ld, ii, nf, ip, nl, nq, nr;
    47984798    static r_8 arg;
    4799     static integer ido, ipm;
     4799    static int_4 ido, ipm;
    48004800    static r_8 tpi;
    48014801
    … …  
    48954895
    48964896/*  ------ File dcosqb.f ------ */
    4897 /* Subroutine */ int dcosqb_(integer *n, r_8 *x, r_8 *wsave)
     4897/* Subroutine */ int dcosqb_(int_4 *n, r_8 *x, r_8 *wsave)
    48984898{
    48994899    /* Initialized data */
    … …  
    49024902
    49034903    /* System generated locals */
    4904     integer i__1;
     4904    int_4 i__1;
    49054905
    49064906    /* Local variables */
    49074907    static r_8 x1;
    4908     extern /* Subroutine */ int dcosqb1_(integer *, r_8 *, r_8 *, r_8 *);
     4908    extern /* Subroutine */ int dcosqb1_(int_4 *, r_8 *, r_8 *, r_8 *);
    49094909
    49104910    /* Parameter adjustments */
    … …  
    49344934
    49354935/*  ------ File dcosqb1.f ------ */
    4936 /* Subroutine */ int dcosqb1_(integer *n, r_8 *x, r_8 *w, r_8 *xh)
     4936/* Subroutine */ int dcosqb1_(int_4 *n, r_8 *x, r_8 *w, r_8 *xh)
    49374937{
    49384938    /* System generated locals */
    4939     integer i__1;
     4939    int_4 i__1;
    49404940
    49414941    /* Local variables */
    4942     static integer modn, i__, k;
    4943     extern /* Subroutine */ int dfftb_(integer *, r_8 *, r_8 *);
    4944     static integer kc, np2, ns2;
     4942    static int_4 modn, i__, k;
     4943    extern /* Subroutine */ int dfftb_(int_4 *, r_8 *, r_8 *);
     4944    static int_4 kc, np2, ns2;
    49454945    static r_8 xim1;
    49464946
    … …  
    49884988
    49894989/*  ------ File dcosqf.f ------ */
    4990 /* Subroutine */ int dcosqf_(integer *n, r_8 *x, r_8 *wsave)
     4990/* Subroutine */ int dcosqf_(int_4 *n, r_8 *x, r_8 *wsave)
    49914991{
    49924992    /* Initialized data */
    … …  
    49954995
    49964996    /* System generated locals */
    4997     integer i__1;
     4997    int_4 i__1;
    49984998
    49994999    /* Local variables */
    50005000    static r_8 tsqx;
    5001     extern /* Subroutine */ int dcosqf1_(integer *, r_8 *, r_8 *, r_8 *);
     5001    extern /* Subroutine */ int dcosqf1_(int_4 *, r_8 *, r_8 *, r_8 *);
    50025002
    50035003    /* Parameter adjustments */
    … …  
    50255025
    50265026/*  ------ File dcosqf1.f ------ */
    5027 /* Subroutine */ int dcosqf1_(integer *n, r_8 *x, r_8 *w, r_8 *xh)
     5027/* Subroutine */ int dcosqf1_(int_4 *n, r_8 *x, r_8 *w, r_8 *xh)
    50285028{
    50295029    /* System generated locals */
    5030     integer i__1;
     5030    int_4 i__1;
    50315031
    50325032    /* Local variables */
    5033     static integer modn, i__, k;
    5034     extern /* Subroutine */ int dfftf_(integer *, r_8 *, r_8 *);
    5035     static integer kc, np2, ns2;
     5033    static int_4 modn, i__, k;
     5034    extern /* Subroutine */ int dfftf_(int_4 *, r_8 *, r_8 *);
     5035    static int_4 kc, np2, ns2;
    50365036    static r_8 xim1;
    50375037
    … …  
    50775077
    50785078/*  ------ File dcosqi.f ------ */
    5079 /* Subroutine */ int dcosqi_(integer *n, r_8 *wsave)
     5079/* Subroutine */ int dcosqi_(int_4 *n, r_8 *wsave)
    50805080{
    50815081    /* Initialized data */
    … …  
    50845084
    50855085    /* System generated locals */
    5086     integer i__1;
     5086    int_4 i__1;
    50875087
    50885088    /* Builtin functions */
    … …  
    50905090
    50915091    /* Local variables */
    5092     static integer k;
    5093     extern /* Subroutine */ int dffti_(integer *, r_8 *);
     5092    static int_4 k;
     5093    extern /* Subroutine */ int dffti_(int_4 *, r_8 *);
    50945094    static r_8 fk, dt;
    50955095
    … …  
    51115111
    51125112/*  ------ File dcost.f ------ */
    5113 /* Subroutine */ int dcost_(integer *n, r_8 *x, r_8 *wsave)
     5113/* Subroutine */ int dcost_(int_4 *n, r_8 *x, r_8 *wsave)
    51145114{
    51155115    /* System generated locals */
    5116     integer i__1;
     5116    int_4 i__1;
    51175117
    51185118    /* Local variables */
    5119     static integer modn, i__, k;
    5120     extern /* Subroutine */ int dfftf_(integer *, r_8 *, r_8 *);
     5119    static int_4 modn, i__, k;
     5120    extern /* Subroutine */ int dfftf_(int_4 *, r_8 *, r_8 *);
    51215121    static r_8 c1, t1, t2;
    5122     static integer kc;
     5122    static int_4 kc;
    51235123    static r_8 xi;
    5124     static integer nm1, np1;
     5124    static int_4 nm1, np1;
    51255125    static r_8 x1h;
    5126     static integer ns2;
     5126    static int_4 ns2;
    51275127    static r_8 tx2, x1p3, xim2;
    51285128
    … …  
    51945194
    51955195/*  ------ File dcosti.f ------ */
    5196 /* Subroutine */ int dcosti_(integer *n, r_8 *wsave)
     5196/* Subroutine */ int dcosti_(int_4 *n, r_8 *wsave)
    51975197{
    51985198    /* Initialized data */
    … …  
    52015201
    52025202    /* System generated locals */
    5203     integer i__1;
     5203    int_4 i__1;
    52045204
    52055205    /* Builtin functions */
    … …  
    52075207
    52085208    /* Local variables */
    5209     static integer k;
    5210     extern /* Subroutine */ int dffti_(integer *, r_8 *);
    5211     static integer kc;
     5209    static int_4 k;
     5210    extern /* Subroutine */ int dffti_(int_4 *, r_8 *);
     5211    static int_4 kc;
    52125212    static r_8 fk, dt;
    5213     static integer nm1, np1, ns2;
     5213    static int_4 nm1, np1, ns2;
    52145214
    52155215    /* Parameter adjustments */
    … …  
    52385238
    52395239/*  ------ File dezfft1.f ------ */
    5240 /* Subroutine */ int dezfft1_(integer *n, r_8 *wa, int_8 *ifac)
     5240/* Subroutine */ int dezfft1_(int_4 *n, r_8 *wa, int_8 *ifac)
    52415241{
    52425242    /* Initialized data */
    52435243
    5244     static integer ntryh[4] = { 4,2,3,5 };
     5244    static int_4 ntryh[4] = { 4,2,3,5 };
    52455245    static r_8 tpi = 6.28318530717959f;
    52465246
    52475247    /* System generated locals */
    5248     integer i__1, i__2, i__3;
     5248    int_4 i__1, i__2, i__3;
    52495249
    52505250    /* Builtin functions */
    … …  
    52535253    /* Local variables */
    52545254    static r_8 argh;
    5255     static integer ntry, i__, j, k1, l1, l2, ib, ii, nf, ip, nl, is, nq, nr;
     5255    static int_4 ntry, i__, j, k1, l1, l2, ib, ii, nf, ip, nl, is, nq, nr;
    52565256    static r_8 ch1, sh1;
    5257     static integer ido, ipm;
     5257    static int_4 ido, ipm;
    52585258    static r_8 dch1, ch1h, arg1, dsh1;
    5259     static integer nfm1;
     5259    static int_4 nfm1;
    52605260
    52615261    /* Parameter adjustments */
    … …  
    53575357
    53585358/*  ------ File dezfftb.f ------ */
    5359 /* Subroutine */ int dezfftb_(integer *n, r_8 *r__, r_8 *azero, r_8 *a,
     5359/* Subroutine */ int dezfftb_(int_4 *n, r_8 *r__, r_8 *azero, r_8 *a,
    53605360        r_8 *b, r_8 *wsave)
    53615361{
    53625362    /* System generated locals */
    5363     integer i__1;
     5363    int_4 i__1;
    53645364
    53655365    /* Local variables */
    5366     static integer i__;
    5367     extern /* Subroutine */ int dfftb_(integer *, r_8 *, r_8 *);
    5368     static integer ns2;
     5366    static int_4 i__;
     5367    extern /* Subroutine */ int dfftb_(int_4 *, r_8 *, r_8 *);
     5368    static int_4 ns2;
    53695369
    53705370    /* Parameter adjustments */
    … …  
    54065406
    54075407/*  ------ File dezfftf.f ------ */
    5408 /* Subroutine */ int dezfftf_(integer *n, r_8 *r__, r_8 *azero, r_8 *a,
     5408/* Subroutine */ int dezfftf_(int_4 *n, r_8 *r__, r_8 *azero, r_8 *a,
    54095409        r_8 *b, r_8 *wsave)
    54105410{
    54115411    /* System generated locals */
    5412     integer i__1;
     5412    int_4 i__1;
    54135413
    54145414    /* Local variables */
    5415     static integer i__;
    5416     extern /* Subroutine */ int dfftf_(integer *, r_8 *, r_8 *);
     5415    static int_4 i__;
     5416    extern /* Subroutine */ int dfftf_(int_4 *, r_8 *, r_8 *);
    54175417    static r_8 cf;
    5418     static integer ns2;
     5418    static int_4 ns2;
    54195419    static r_8 cfm;
    5420     static integer ns2m;
     5420    static int_4 ns2m;
    54215421
    54225422
    … …  
    54715471
    54725472/*  ------ File dezffti.f ------ */
    5473 /* Subroutine */ int dezffti_(integer *n, r_8 *wsave)
    5474 {
    5475     extern /* Subroutine */ int dezfft1_(integer *, r_8 *, int_8 *);
     5473/* Subroutine */ int dezffti_(int_4 *n, r_8 *wsave)
     5474{
     5475    extern /* Subroutine */ int dezfft1_(int_4 *, r_8 *, int_8 *);
    54765476
    54775477    /* Parameter adjustments */
    … …  
    54875487
    54885488/*  ------ File dpassb.f ------ */
    5489 /* Subroutine */ int dpassb_(integer *nac, integer *ido, integer *ip, integer *
    5490         l1, integer *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2,
     5489/* Subroutine */ int dpassb_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *
     5490        l1, int_4 *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2,
    54915491        r_8 *wa)
    54925492{
    54935493    /* System generated locals */
    5494     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     5494    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    54955495             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    54965496            i__1, i__2, i__3;
    54975497
    54985498    /* Local variables */
    5499     static integer idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
     5499    static int_4 idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
    55005500            idl, inc, idp;
    55015501    static r_8 wai, war;
    5502     static integer ipp2;
     5502    static int_4 ipp2;
    55035503
    55045504    /* Parameter adjustments */
    … …  
    57285728
    57295729/*  ------ File dpassb2.f ------ */
    5730 /* Subroutine */ int dpassb2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     5730/* Subroutine */ int dpassb2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    57315731        r_8 *wa1)
    57325732{
    57335733    /* System generated locals */
    5734     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     5734    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    57355735
    57365736    /* Local variables */
    5737     static integer i__, k;
     5737    static int_4 i__, k;
    57385738    static r_8 ti2, tr2;
    57395739
    … …  
    57905790
    57915791/*  ------ File dpassb3.f ------ */
    5792 /* Subroutine */ int dpassb3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     5792/* Subroutine */ int dpassb3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    57935793        r_8 *wa1, r_8 *wa2)
    57945794{
    … …  
    57995799
    58005800    /* System generated locals */
    5801     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     5801    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    58025802
    58035803    /* Local variables */
    5804     static integer i__, k;
     5804    static int_4 i__, k;
    58055805    static r_8 ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    58065806
    … …  
    58785878
    58795879/*  ------ File dpassb4.f ------ */
    5880 /* Subroutine */ int dpassb4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     5880/* Subroutine */ int dpassb4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    58815881        r_8 *wa1, r_8 *wa2, r_8 *wa3)
    58825882{
    58835883    /* System generated locals */
    5884     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     5884    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    58855885
    58865886    /* Local variables */
    5887     static integer i__, k;
     5887    static int_4 i__, k;
    58885888    static r_8 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    58895889            tr3, tr4;
    … …  
    59835983
    59845984/*  ------ File dpassb5.f ------ */
    5985 /* Subroutine */ int dpassb5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     5985/* Subroutine */ int dpassb5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    59865986        r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4)
    59875987{
    … …  
    59945994
    59955995    /* System generated locals */
    5996     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     5996    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    59975997
    59985998    /* Local variables */
    5999     static integer i__, k;
     5999    static int_4 i__, k;
    60006000    static r_8 ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
    60016001            ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;
    … …  
    61186118
    61196119/*  ------ File dpassf.f ------ */
    6120 /* Subroutine */ int dpassf_(integer *nac, integer *ido, integer *ip, integer *
    6121         l1, integer *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2,
     6120/* Subroutine */ int dpassf_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *
     6121        l1, int_4 *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2,
    61226122        r_8 *wa)
    61236123{
    61246124    /* System generated locals */
    6125     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     6125    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    61266126             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    61276127            i__1, i__2, i__3;
    61286128
    61296129    /* Local variables */
    6130     static integer idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
     6130    static int_4 idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
    61316131            idl, inc, idp;
    61326132    static r_8 wai, war;
    6133     static integer ipp2;
     6133    static int_4 ipp2;
    61346134
    61356135    /* Parameter adjustments */
    … …  
    63596359
    63606360/*  ------ File dpassf2.f ------ */
    6361 /* Subroutine */ int dpassf2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     6361/* Subroutine */ int dpassf2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    63626362        r_8 *wa1)
    63636363{
    63646364    /* System generated locals */
    6365     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     6365    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    63666366
    63676367    /* Local variables */
    6368     static integer i__, k;
     6368    static int_4 i__, k;
    63696369    static r_8 ti2, tr2;
    63706370
    … …  
    64216421
    64226422/*  ------ File dpassf3.f ------ */
    6423 /* Subroutine */ int dpassf3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     6423/* Subroutine */ int dpassf3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    64246424        r_8 *wa1, r_8 *wa2)
    64256425{
    … …  
    64306430
    64316431    /* System generated locals */
    6432     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     6432    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    64336433
    64346434    /* Local variables */
    6435     static integer i__, k;
     6435    static int_4 i__, k;
    64366436    static r_8 ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    64376437
    … …  
    65096509
    65106510/*  ------ File dpassf4.f ------ */
    6511 /* Subroutine */ int dpassf4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     6511/* Subroutine */ int dpassf4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    65126512        r_8 *wa1, r_8 *wa2, r_8 *wa3)
    65136513{
    65146514    /* System generated locals */
    6515     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     6515    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    65166516
    65176517    /* Local variables */
    6518     static integer i__, k;
     6518    static int_4 i__, k;
    65196519    static r_8 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    65206520            tr3, tr4;
    … …  
    66146614
    66156615/*  ------ File dpassf5.f ------ */
    6616 /* Subroutine */ int dpassf5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     6616/* Subroutine */ int dpassf5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    66176617        r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4)
    66186618{
    … …  
    66256625
    66266626    /* System generated locals */
    6627     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     6627    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    66286628
    66296629    /* Local variables */
    6630     static integer i__, k;
     6630    static int_4 i__, k;
    66316631    static r_8 ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
    66326632            ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;
    … …  
    67496749
    67506750/*  ------ File dadb2.f ------ */
    6751 /* Subroutine */ int dadb2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     6751/* Subroutine */ int dadb2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    67526752        r_8 *wa1)
    67536753{
    67546754    /* System generated locals */
    6755     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     6755    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    67566756
    67576757    /* Local variables */
    6758     static integer i__, k, ic;
     6758    static int_4 i__, k, ic;
    67596759    static r_8 ti2, tr2;
    6760     static integer idp2;
     6760    static int_4 idp2;
    67616761
    67626762    /* Parameter adjustments */
    … …  
    68266826
    68276827/*  ------ File dadb3.f ------ */
    6828 /* Subroutine */ int dadb3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     6828/* Subroutine */ int dadb3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    68296829        r_8 *wa1, r_8 *wa2)
    68306830{
    … …  
    68356835
    68366836    /* System generated locals */
    6837     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     6837    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    68386838
    68396839    /* Local variables */
    6840     static integer i__, k, ic;
     6840    static int_4 i__, k, ic;
    68416841    static r_8 ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    6842     static integer idp2;
     6842    static int_4 idp2;
    68436843
    68446844    /* Parameter adjustments */
    … …  
    69096909
    69106910/*  ------ File dadb4.f ------ */
    6911 /* Subroutine */ int dadb4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     6911/* Subroutine */ int dadb4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    69126912        r_8 *wa1, r_8 *wa2, r_8 *wa3)
    69136913{
    … …  
    69176917
    69186918    /* System generated locals */
    6919     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     6919    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    69206920
    69216921    /* Local variables */
    6922     static integer i__, k, ic;
     6922    static int_4 i__, k, ic;
    69236923    static r_8 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    69246924            tr3, tr4;
    6925     static integer idp2;
     6925    static int_4 idp2;
    69266926
    69276927    /* Parameter adjustments */
    … …  
    70337033
    70347034/*  ------ File dadb5.f ------ */
    7035 /* Subroutine */ int dadb5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     7035/* Subroutine */ int dadb5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    70367036        r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4)
    70377037{
    … …  
    70447044
    70457045    /* System generated locals */
    7046     integer cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
     7046    int_4 cc_dim1, cc_offset, ch_dim1, ch_dim2, ch_offset, i__1, i__2;
    70477047
    70487048    /* Local variables */
    7049     static integer i__, k, ic;
     7049    static int_4 i__, k, ic;
    70507050    static r_8 ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
    70517051            ti2, ti3, ti4, ti5, dr3, dr4, dr5, dr2, tr2, tr3, tr4, tr5;
    7052     static integer idp2;
     7052    static int_4 idp2;
    70537053
    70547054    /* Parameter adjustments */
    … …  
    71577157
    71587158/*  ------ File dadbg.f ------ */
    7159 /* Subroutine */ int dadbg_(integer *ido, integer *ip, integer *l1, integer *
     7159/* Subroutine */ int dadbg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *
    71607160        idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa)
    71617161{
    … …  
    71657165
    71667166    /* System generated locals */
    7167     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     7167    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    71687168             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    71697169            i__1, i__2, i__3;
    … …  
    71737173
    71747174    /* Local variables */
    7175     static integer idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
     7175    static int_4 idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
    71767176    static r_8 dc2, ai1, ai2, ar1, ar2, ds2;
    7177     static integer nbd;
     7177    static int_4 nbd;
    71787178    static r_8 dcp, arg, dsp, ar1h, ar2h;
    7179     static integer idp2, ipp2;
     7179    static int_4 idp2, ipp2;
    71807180
    71817181    /* Parameter adjustments */
    … …  
    74987498
    74997499/*  ------ File dadf2.f ------ */
    7500 /* Subroutine */ int dadf2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     7500/* Subroutine */ int dadf2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    75017501        r_8 *wa1)
    75027502{
    75037503    /* System generated locals */
    7504     integer ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
     7504    int_4 ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
    75057505
    75067506    /* Local variables */
    7507     static integer i__, k, ic;
     7507    static int_4 i__, k, ic;
    75087508    static r_8 ti2, tr2;
    7509     static integer idp2;
     7509    static int_4 idp2;
    75107510
    75117511    /* Parameter adjustments */
    … …  
    75767576
    75777577/*  ------ File dadf3.f ------ */
    7578 /* Subroutine */ int dadf3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     7578/* Subroutine */ int dadf3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    75797579        r_8 *wa1, r_8 *wa2)
    75807580{
    … …  
    75857585
    75867586    /* System generated locals */
    7587     integer ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
     7587    int_4 ch_dim1, ch_offset, cc_dim1, cc_dim2, cc_offset, i__1, i__2;
    75887588
    75897589    /* Local variables */
    7590     static integer i__, k, ic;
     7590    static int_4 i__, k, ic;
    75917591    static r_8 ci2, di2, di3, cr2, dr2, dr3, ti2, ti3, tr2, tr3;
    7592     static integer idp2;
     7592    static int_4 idp2;
    75937593
    75947594    /* Parameter adjustments */
    … …  
    76557655
    76567656/*  ------ File dadf4.f ------ */
    7657 /* Subroutine */ int dadf4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     7657/* Subroutine */ int dadf4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    76587658        r_8 *wa1, r_8 *wa2, r_8 *wa3)
    76597659{
    … …  
    76637663
    76647664    /* System generated locals */
    7665     integer cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
     7665    int_4 cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
    76667666
    76677667    /* Local variables */
    7668     static integer i__, k, ic;
     7668    static int_4 i__, k, ic;
    76697669    static r_8 ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    76707670            tr3, tr4;
    7671     static integer idp2;
     7671    static int_4 idp2;
    76727672
    76737673    /* Parameter adjustments */
    … …  
    77717771
    77727772/*  ------ File dadf5.f ------ */
    7773 /* Subroutine */ int dadf5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch,
     7773/* Subroutine */ int dadf5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch,
    77747774        r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4)
    77757775{
    … …  
    77827782
    77837783    /* System generated locals */
    7784     integer cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
     7784    int_4 cc_dim1, cc_dim2, cc_offset, ch_dim1, ch_offset, i__1, i__2;
    77857785
    77867786    /* Local variables */
    7787     static integer i__, k, ic;
     7787    static int_4 i__, k, ic;
    77887788    static r_8 ci2, di2, ci4, ci5, di3, di4, di5, ci3, cr2, cr3, dr2, dr3,
    77897789            dr4, dr5, cr5, cr4, ti2, ti3, ti5, ti4, tr2, tr3, tr4, tr5;
    7790     static integer idp2;
     7790    static int_4 idp2;
    77917791
    77927792    /* Parameter adjustments */
    … …  
    78897889
    78907890/*  ------ File dadfg.f ------ */
    7891 /* Subroutine */ int dadfg_(integer *ido, integer *ip, integer *l1, integer *
     7891/* Subroutine */ int dadfg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *
    78927892        idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa)
    78937893{
    … …  
    78977897
    78987898    /* System generated locals */
    7899     integer ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
     7899    int_4 ch_dim1, ch_dim2, ch_offset, cc_dim1, cc_dim2, cc_offset, c1_dim1,
    79007900             c1_dim2, c1_offset, c2_dim1, c2_offset, ch2_dim1, ch2_offset,
    79017901            i__1, i__2, i__3;
    … …  
    79057905
    79067906    /* Local variables */
    7907     static integer idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
     7907    static int_4 idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
    79087908    static r_8 dc2, ai1, ai2, ar1, ar2, ds2;
    7909     static integer nbd;
     7909    static int_4 nbd;
    79107910    static r_8 dcp, arg, dsp, ar1h, ar2h;
    7911     static integer idp2, ipp2;
     7911    static int_4 idp2, ipp2;
    79127912
    79137913    /* Parameter adjustments */
    … …  
    82368236
    82378237/*  ------ File dfftb.f ------ */
    8238 /* Subroutine */ int dfftb_(integer *n, r_8 *r__, r_8 *wsave)
    8239 {
    8240     extern /* Subroutine */ int dfftb1_(integer *, r_8 *, r_8 *, r_8 *,
     8238/* Subroutine */ int dfftb_(int_4 *n, r_8 *r__, r_8 *wsave)
     8239{
     8240    extern /* Subroutine */ int dfftb1_(int_4 *, r_8 *, r_8 *, r_8 *,
    82418241            int_8 *);
    82428242
    … …  
    82548254
    82558255/*  ------ File dfftb1.f ------ */
    8256 /* Subroutine */ int dfftb1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa,
     8256/* Subroutine */ int dfftb1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa,
    82578257        int_8 *ifac)
    82588258{
    82598259    /* System generated locals */
    8260     integer i__1;
     8260    int_4 i__1;
    82618261
    82628262    /* Local variables */
    8263     extern /* Subroutine */ int dadb2_(integer *, integer *, r_8 *, r_8 *,
    8264             r_8 *), dadb3_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    8265             r_8 *), dadb4_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    8266             r_8 *, r_8 *), dadb5_(integer *, integer *, r_8 *, r_8 *,
     8263    extern /* Subroutine */ int dadb2_(int_4 *, int_4 *, r_8 *, r_8 *,
     8264            r_8 *), dadb3_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     8265            r_8 *), dadb4_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     8266            r_8 *, r_8 *), dadb5_(int_4 *, int_4 *, r_8 *, r_8 *,
    82678267            r_8 *, r_8 *, r_8 *, r_8 *);
    8268     static integer i__;
    8269     extern /* Subroutine */ int dadbg_(integer *, integer *, integer *,
    8270             integer *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *);
    8271     static integer k1, l1, l2, na, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
     8268    static int_4 i__;
     8269    extern /* Subroutine */ int dadbg_(int_4 *, int_4 *, int_4 *,
     8270            int_4 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *);
     8271    static int_4 k1, l1, l2, na, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
    82728272
    82738273    /* Parameter adjustments */
    … …  
    83828382
    83838383/*  ------ File dfftf.f ------ */
    8384 /* Subroutine */ int dfftf_(integer *n, r_8 *r__, r_8 *wsave)
    8385 {
    8386     extern /* Subroutine */ int dfftf1_(integer *, r_8 *, r_8 *, r_8 *,
     8384/* Subroutine */ int dfftf_(int_4 *n, r_8 *r__, r_8 *wsave)
     8385{
     8386    extern /* Subroutine */ int dfftf1_(int_4 *, r_8 *, r_8 *, r_8 *,
    83878387            int_8 *);
    83888388
    … …  
    84008400
    84018401/*  ------ File dfftf1.f ------ */
    8402 /* Subroutine */ int dfftf1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa,
     8402/* Subroutine */ int dfftf1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa,
    84038403        int_8 *ifac)
    84048404{
    84058405    /* System generated locals */
    8406     integer i__1;
     8406    int_4 i__1;
    84078407
    84088408    /* Local variables */
    8409     extern /* Subroutine */ int dadf2_(integer *, integer *, r_8 *, r_8 *,
    8410             r_8 *), dadf3_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    8411             r_8 *), dadf4_(integer *, integer *, r_8 *, r_8 *, r_8 *,
    8412             r_8 *, r_8 *), dadf5_(integer *, integer *, r_8 *, r_8 *,
     8409    extern /* Subroutine */ int dadf2_(int_4 *, int_4 *, r_8 *, r_8 *,
     8410            r_8 *), dadf3_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     8411            r_8 *), dadf4_(int_4 *, int_4 *, r_8 *, r_8 *, r_8 *,
     8412            r_8 *, r_8 *), dadf5_(int_4 *, int_4 *, r_8 *, r_8 *,
    84138413            r_8 *, r_8 *, r_8 *, r_8 *);
    8414     static integer i__;
    8415     extern /* Subroutine */ int dadfg_(integer *, integer *, integer *,
    8416             integer *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *);
    8417     static integer k1, l1, l2, na, kh, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
     8414    static int_4 i__;
     8415    extern /* Subroutine */ int dadfg_(int_4 *, int_4 *, int_4 *,
     8416            int_4 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *, r_8 *);
     8417    static int_4 k1, l1, l2, na, kh, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
    84188418
    84198419    /* Parameter adjustments */
    … …  
    85238523
    85248524/*  ------ File dffti.f ------ */
    8525 /* Subroutine */ int dffti_(integer *n, r_8 *wsave)
    8526 {
    8527     extern /* Subroutine */ int dffti1_(integer *, r_8 *, int_8 *);
     8525/* Subroutine */ int dffti_(int_4 *n, r_8 *wsave)
     8526{
     8527    extern /* Subroutine */ int dffti1_(int_4 *, r_8 *, int_8 *);
    85288528
    85298529    /* Parameter adjustments */
    … …  
    85398539
    85408540/*  ------ File dffti1.f ------ */
    8541 /* Subroutine */ int dffti1_(integer *n, r_8 *wa, int_8 *ifac)
     8541/* Subroutine */ int dffti1_(int_4 *n, r_8 *wa, int_8 *ifac)
    85428542{
    85438543    /* Initialized data */
    85448544
    8545     static integer ntryh[4] = { 4,2,3,5 };
     8545    static int_4 ntryh[4] = { 4,2,3,5 };
    85468546
    85478547    /* System generated locals */
    8548     integer i__1, i__2, i__3;
     8548    int_4 i__1, i__2, i__3;
    85498549
    85508550    /* Builtin functions */
    … …  
    85538553    /* Local variables */
    85548554    static r_8 argh;
    8555     static integer ntry, i__, j;
     8555    static int_4 ntry, i__, j;
    85568556    static r_8 argld;
    8557     static integer k1, l1, l2, ib;
     8557    static int_4 k1, l1, l2, ib;
    85588558    static r_8 fi;
    8559     static integer ld, ii, nf, ip, nl, is, nq, nr;
     8559    static int_4 ld, ii, nf, ip, nl, is, nq, nr;
    85608560    static r_8 arg;
    8561     static integer ido, ipm;
     8561    static int_4 ido, ipm;
    85628562    static r_8 tpi;
    8563     static integer nfm1;
     8563    static int_4 nfm1;
    85648564
    85658565    /* Parameter adjustments */
    … …  
    86548654
    86558655/*  ------ File dsinqb.f ------ */
    8656 /* Subroutine */ int dsinqb_(integer *n, r_8 *x, r_8 *wsave)
     8656/* Subroutine */ int dsinqb_(int_4 *n, r_8 *x, r_8 *wsave)
    86578657{
    86588658    /* System generated locals */
    8659     integer i__1;
     8659    int_4 i__1;
    86608660
    86618661    /* Local variables */
    8662     static integer k;
    8663     extern /* Subroutine */ int dcosqb_(integer *, r_8 *, r_8 *);
     8662    static int_4 k;
     8663    extern /* Subroutine */ int dcosqb_(int_4 *, r_8 *, r_8 *);
    86648664    static r_8 xhold;
    8665     static integer kc, ns2;
     8665    static int_4 kc, ns2;
    86668666
    86678667    /* Parameter adjustments */
    … …  
    86958695
    86968696/*  ------ File dsinqf.f ------ */
    8697 /* Subroutine */ int dsinqf_(integer *n, r_8 *x, r_8 *wsave)
     8697/* Subroutine */ int dsinqf_(int_4 *n, r_8 *x, r_8 *wsave)
    86988698{
    86998699    /* System generated locals */
    8700     integer i__1;
     8700    int_4 i__1;
    87018701
    87028702    /* Local variables */
    8703     static integer k;
    8704     extern /* Subroutine */ int dcosqf_(integer *, r_8 *, r_8 *);
     8703    static int_4 k;
     8704    extern /* Subroutine */ int dcosqf_(int_4 *, r_8 *, r_8 *);
    87058705    static r_8 xhold;
    8706     static integer kc, ns2;
     8706    static int_4 kc, ns2;
    87078707
    87088708    /* Parameter adjustments */
    … …  
    87338733
    87348734/*  ------ File dsinqi.f ------ */
    8735 /* Subroutine */ int dsinqi_(integer *n, r_8 *wsave)
    8736 {
    8737     extern /* Subroutine */ int dcosqi_(integer *, r_8 *);
     8735/* Subroutine */ int dsinqi_(int_4 *n, r_8 *wsave)
     8736{
     8737    extern /* Subroutine */ int dcosqi_(int_4 *, r_8 *);
    87388738
    87398739    /* Parameter adjustments */
    … …  
    87468746
    87478747/*  ------ File dsint.f ------ */
    8748 /* Subroutine */ int dsint_(integer *n, r_8 *x, r_8 *wsave)
    8749 {
    8750     extern /* Subroutine */ int dsint1_(integer *, r_8 *, r_8 *, r_8 *,
     8748/* Subroutine */ int dsint_(int_4 *n, r_8 *x, r_8 *wsave)
     8749{
     8750    extern /* Subroutine */ int dsint1_(int_4 *, r_8 *, r_8 *, r_8 *,
    87518751            r_8 *, int_8 *);
    8752     static integer np1, iw1, iw2, iw3;
     8752    static int_4 np1, iw1, iw2, iw3;
    87538753
    87548754    /* Parameter adjustments */
    … …  
    87668766
    87678767/*  ------ File dsint1.f ------ */
    8768 /* Subroutine */ int dsint1_(integer *n, r_8 *war, r_8 *was, r_8 *xh, r_8 *
     8768/* Subroutine */ int dsint1_(int_4 *n, r_8 *war, r_8 *was, r_8 *xh, r_8 *
    87698769        x, int_8 *ifac)
    87708770{
    … …  
    87748774
    87758775    /* System generated locals */
    8776     integer i__1;
     8776    int_4 i__1;
    87778777
    87788778    /* Local variables */
    8779     static integer modn, i__, k;
     8779    static int_4 modn, i__, k;
    87808780    static r_8 xhold, t1, t2;
    8781     extern /* Subroutine */ int dfftf1_(integer *, r_8 *, r_8 *, r_8 *,
     8781    extern /* Subroutine */ int dfftf1_(int_4 *, r_8 *, r_8 *, r_8 *,
    87828782            int_8 *);
    8783     static integer kc, np1, ns2;
     8783    static int_4 kc, np1, ns2;
    87848784
    87858785    /* Parameter adjustments */
    … …  
    88528852
    88538853/*  ------ File dsinti.f ------ */
    8854 /* Subroutine */ int dsinti_(integer *n, r_8 *wsave)
     8854/* Subroutine */ int dsinti_(int_4 *n, r_8 *wsave)
    88558855{
    88568856    /* Initialized data */
    … …  
    88598859
    88608860    /* System generated locals */
    8861     integer i__1;
     8861    int_4 i__1;
    88628862
    88638863    /* Builtin functions */
    … …  
    88658865
    88668866    /* Local variables */
    8867     static integer k;
    8868     extern /* Subroutine */ int dffti_(integer *, r_8 *);
     8867    static int_4 k;
     8868    extern /* Subroutine */ int dffti_(int_4 *, r_8 *);
    88698869    static r_8 dt;
    8870     static integer np1, ns2;
     8870    static int_4 np1, ns2;
    88718871
    88728872    /* Parameter adjustments */

  • trunk/SophyaLib/NTools/fftpackc.h

    r717 r756  
    88#endif
    99
    10 typedef int integer;
    11 typedef unsigned uinteger;
    12 typedef char *address;
    13 typedef short int shortint;
    14   // typedef float real;
    15   // typedef double doublereal;
    16   //typedef struct { float r, i; } complex;
    17   //typedef struct { double r, i; } doublecomplex;
     10  /* typedef int integer;
     11  typedef unsigned uinteger;
     12  typedef char *address;
     13  typedef short int shortint;
     14  typedef float real;
     15  typedef double doublereal;
     16  typedef struct { float r, i; } complex;
     17  typedef struct { double r, i; } doublecomplex;  */
    1818
    1919/* typedef int_4 logical;
    … …  
    2424/* ---------------------- */
    2525
    26 extern int cfftb_(integer *n, r_4 *c__, r_4 *wsave);
    27 extern int cfftb1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
    28 extern int cfftf_(integer *n, r_4 *c__, r_4 *wsave);
    29 extern int cfftf1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
    30 extern int cffti_(integer *n, r_4 *wsave);
    31 extern int cffti1_(integer *n, r_4 *wa, int_4 *ifac);
    32 extern int cosqb_(integer *n, r_4 *x, r_4 *wsave);
    33 extern int cosqb1_(integer *n, r_4 *x, r_4 *w, r_4 *xh);
    34 extern int cosqf_(integer *n, r_4 *x, r_4 *wsave);
    35 extern int cosqf1_(integer *n, r_4 *x, r_4 *w, r_4 *xh);
    36 extern int cosqi_(integer *n, r_4 *wsave);
    37 extern int cost_(integer *n, r_4 *x, r_4 *wsave);
    38 extern int costi_(integer *n, r_4 *wsave);
    39 extern int ezfft1_(integer *n, r_4 *wa, int_4 *ifac);
    40 extern int ezfftb_(integer *n, r_4 *r__, r_4 *azero, r_4 *a, r_4 *b, r_4 *wsave);
    41 extern int ezfftf_(integer *n, r_4 *r__, r_4 *azero, r_4 *a, r_4 *b, r_4 *wsave);
    42 extern int ezffti_(integer *n, r_4 *wsave);
    43 extern int passb_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
    44 extern int passb2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
    45 extern int passb3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
    46 extern int passb4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
    47 extern int passb5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
    48 extern int passf_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
    49 extern int passf2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
    50 extern int passf3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
    51 extern int passf4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
    52 extern int passf5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
    53 extern int radb2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
    54 extern int radb3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
    55 extern int radb4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
    56 extern int radb5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
    57 extern int radbg_(integer *ido, integer *ip, integer *l1, integer *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
    58 extern int radf2_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
    59 extern int radf3_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
    60 extern int radf4_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
    61 extern int radf5_(integer *ido, integer *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
    62 extern int radfg_(integer *ido, integer *ip, integer *l1, integer *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
    63 extern int rfftb_(integer *n, r_4 *r__, r_4 *wsave);
    64 extern int rfftb1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
    65 extern int rfftf_(integer *n, r_4 *r__, r_4 *wsave);
    66 extern int rfftf1_(integer *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
    67 extern int rffti_(integer *n, r_4 *wsave);
    68 extern int rffti1_(integer *n, r_4 *wa, int_4 *ifac);
    69 extern int sinqb_(integer *n, r_4 *x, r_4 *wsave);
    70 extern int sinqf_(integer *n, r_4 *x, r_4 *wsave);
    71 extern int sinqi_(integer *n, r_4 *wsave);
    72 extern int sint_(integer *n, r_4 *x, r_4 *wsave);
    73 extern int sint1_(integer *n, r_4 *war, r_4 *was, r_4 *xh, r_4 *x, int_4 *ifac);
    74 extern int sinti_(integer *n, r_4 *wsave);
     26extern int cfftb_(int_4 *n, r_4 *c__, r_4 *wsave);
     27extern int cfftb1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
     28extern int cfftf_(int_4 *n, r_4 *c__, r_4 *wsave);
     29extern int cfftf1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
     30extern int cffti_(int_4 *n, r_4 *wsave);
     31extern int cffti1_(int_4 *n, r_4 *wa, int_4 *ifac);
     32extern int cosqb_(int_4 *n, r_4 *x, r_4 *wsave);
     33extern int cosqb1_(int_4 *n, r_4 *x, r_4 *w, r_4 *xh);
     34extern int cosqf_(int_4 *n, r_4 *x, r_4 *wsave);
     35extern int cosqf1_(int_4 *n, r_4 *x, r_4 *w, r_4 *xh);
     36extern int cosqi_(int_4 *n, r_4 *wsave);
     37extern int cost_(int_4 *n, r_4 *x, r_4 *wsave);
     38extern int costi_(int_4 *n, r_4 *wsave);
     39extern int ezfft1_(int_4 *n, r_4 *wa, int_4 *ifac);
     40extern int ezfftb_(int_4 *n, r_4 *r__, r_4 *azero, r_4 *a, r_4 *b, r_4 *wsave);
     41extern int ezfftf_(int_4 *n, r_4 *r__, r_4 *azero, r_4 *a, r_4 *b, r_4 *wsave);
     42extern int ezffti_(int_4 *n, r_4 *wsave);
     43extern int passb_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
     44extern int passb2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
     45extern int passb3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
     46extern int passb4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
     47extern int passb5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
     48extern int passf_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
     49extern int passf2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
     50extern int passf3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
     51extern int passf4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
     52extern int passf5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
     53extern int radb2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
     54extern int radb3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
     55extern int radb4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
     56extern int radb5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
     57extern int radbg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
     58extern int radf2_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1);
     59extern int radf3_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2);
     60extern int radf4_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3);
     61extern int radf5_(int_4 *ido, int_4 *l1, r_4 *cc, r_4 *ch, r_4 *wa1, r_4 *wa2, r_4 *wa3, r_4 *wa4);
     62extern int radfg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_4 *cc, r_4 *c1, r_4 *c2, r_4 *ch, r_4 *ch2, r_4 *wa);
     63extern int rfftb_(int_4 *n, r_4 *r__, r_4 *wsave);
     64extern int rfftb1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
     65extern int rfftf_(int_4 *n, r_4 *r__, r_4 *wsave);
     66extern int rfftf1_(int_4 *n, r_4 *c__, r_4 *ch, r_4 *wa, int_4 *ifac);
     67extern int rffti_(int_4 *n, r_4 *wsave);
     68extern int rffti1_(int_4 *n, r_4 *wa, int_4 *ifac);
     69extern int sinqb_(int_4 *n, r_4 *x, r_4 *wsave);
     70extern int sinqf_(int_4 *n, r_4 *x, r_4 *wsave);
     71extern int sinqi_(int_4 *n, r_4 *wsave);
     72extern int sint_(int_4 *n, r_4 *x, r_4 *wsave);
     73extern int sint1_(int_4 *n, r_4 *war, r_4 *was, r_4 *xh, r_4 *x, int_4 *ifac);
     74extern int sinti_(int_4 *n, r_4 *wsave);
    7575
    7676/* ---------------------- */
    7777
    78 extern int cdfftb_(integer *n, r_8 *c__, r_8 *wsave);
    79 extern int cdfftb1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
    80 extern int cdfftf_(integer *n, r_8 *c__, r_8 *wsave);
    81 extern int cdfftf1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
    82 extern int cdffti_(integer *n, r_8 *wsave);
    83 extern int cdffti1_(integer *n, r_8 *wa, int_8 *ifac);
    84 extern int dcosqb_(integer *n, r_8 *x, r_8 *wsave);
    85 extern int dcosqb1_(integer *n, r_8 *x, r_8 *w, r_8 *xh);
    86 extern int dcosqf_(integer *n, r_8 *x, r_8 *wsave);
    87 extern int dcosqf1_(integer *n, r_8 *x, r_8 *w, r_8 *xh);
    88 extern int dcosqi_(integer *n, r_8 *wsave);
    89 extern int dcost_(integer *n, r_8 *x, r_8 *wsave);
    90 extern int dcosti_(integer *n, r_8 *wsave);
    91 extern int dezfft1_(integer *n, r_8 *wa, int_8 *ifac);
    92 extern int dezfftb_(integer *n, r_8 *r__, r_8 *azero, r_8 *a, r_8 *b, r_8 *wsave);
    93 extern int dezfftf_(integer *n, r_8 *r__, r_8 *azero, r_8 *a, r_8 *b, r_8 *wsave);
    94 extern int dezffti_(integer *n, r_8 *wsave);
    95 extern int dpassb_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
    96 extern int dpassb2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
    97 extern int dpassb3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
    98 extern int dpassb4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
    99 extern int dpassb5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
    100 extern int dpassf_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
    101 extern int dpassf2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
    102 extern int dpassf3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
    103 extern int dpassf4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
    104 extern int dpassf5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
    105 extern int dadb2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
    106 extern int dadb3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
    107 extern int dadb4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
    108 extern int dadb5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
    109 extern int dadbg_(integer *ido, integer *ip, integer *l1, integer *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
    110 extern int dadf2_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
    111 extern int dadf3_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
    112 extern int dadf4_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
    113 extern int dadf5_(integer *ido, integer *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
    114 extern int dadfg_(integer *ido, integer *ip, integer *l1, integer *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
    115 extern int dfftb_(integer *n, r_8 *r__, r_8 *wsave);
    116 extern int dfftb1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
    117 extern int dfftf_(integer *n, r_8 *r__, r_8 *wsave);
    118 extern int dfftf1_(integer *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
    119 extern int dffti_(integer *n, r_8 *wsave);
    120 extern int dffti1_(integer *n, r_8 *wa, int_8 *ifac);
    121 extern int dsinqb_(integer *n, r_8 *x, r_8 *wsave);
    122 extern int dsinqf_(integer *n, r_8 *x, r_8 *wsave);
    123 extern int dsinqi_(integer *n, r_8 *wsave);
    124 extern int dsint_(integer *n, r_8 *x, r_8 *wsave);
    125 extern int dsint1_(integer *n, r_8 *war, r_8 *was, r_8 *xh, r_8 *x, int_8 *ifac);
    126 extern int dsinti_(integer *n, r_8 *wsave);
     78extern int cdfftb_(int_4 *n, r_8 *c__, r_8 *wsave);
     79extern int cdfftb1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
     80extern int cdfftf_(int_4 *n, r_8 *c__, r_8 *wsave);
     81extern int cdfftf1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
     82extern int cdffti_(int_4 *n, r_8 *wsave);
     83extern int cdffti1_(int_4 *n, r_8 *wa, int_8 *ifac);
     84extern int dcosqb_(int_4 *n, r_8 *x, r_8 *wsave);
     85extern int dcosqb1_(int_4 *n, r_8 *x, r_8 *w, r_8 *xh);
     86extern int dcosqf_(int_4 *n, r_8 *x, r_8 *wsave);
     87extern int dcosqf1_(int_4 *n, r_8 *x, r_8 *w, r_8 *xh);
     88extern int dcosqi_(int_4 *n, r_8 *wsave);
     89extern int dcost_(int_4 *n, r_8 *x, r_8 *wsave);
     90extern int dcosti_(int_4 *n, r_8 *wsave);
     91extern int dezfft1_(int_4 *n, r_8 *wa, int_8 *ifac);
     92extern int dezfftb_(int_4 *n, r_8 *r__, r_8 *azero, r_8 *a, r_8 *b, r_8 *wsave);
     93extern int dezfftf_(int_4 *n, r_8 *r__, r_8 *azero, r_8 *a, r_8 *b, r_8 *wsave);
     94extern int dezffti_(int_4 *n, r_8 *wsave);
     95extern int dpassb_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
     96extern int dpassb2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
     97extern int dpassb3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
     98extern int dpassb4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
     99extern int dpassb5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
     100extern int dpassf_(int_4 *nac, int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
     101extern int dpassf2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
     102extern int dpassf3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
     103extern int dpassf4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
     104extern int dpassf5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
     105extern int dadb2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
     106extern int dadb3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
     107extern int dadb4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
     108extern int dadb5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
     109extern int dadbg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
     110extern int dadf2_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1);
     111extern int dadf3_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2);
     112extern int dadf4_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3);
     113extern int dadf5_(int_4 *ido, int_4 *l1, r_8 *cc, r_8 *ch, r_8 *wa1, r_8 *wa2, r_8 *wa3, r_8 *wa4);
     114extern int dadfg_(int_4 *ido, int_4 *ip, int_4 *l1, int_4 *idl1, r_8 *cc, r_8 *c1, r_8 *c2, r_8 *ch, r_8 *ch2, r_8 *wa);
     115extern int dfftb_(int_4 *n, r_8 *r__, r_8 *wsave);
     116extern int dfftb1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
     117extern int dfftf_(int_4 *n, r_8 *r__, r_8 *wsave);
     118extern int dfftf1_(int_4 *n, r_8 *c__, r_8 *ch, r_8 *wa, int_8 *ifac);
     119extern int dffti_(int_4 *n, r_8 *wsave);
     120extern int dffti1_(int_4 *n, r_8 *wa, int_8 *ifac);
     121extern int dsinqb_(int_4 *n, r_8 *x, r_8 *wsave);
     122extern int dsinqf_(int_4 *n, r_8 *x, r_8 *wsave);
     123extern int dsinqi_(int_4 *n, r_8 *wsave);
     124extern int dsint_(int_4 *n, r_8 *x, r_8 *wsave);
     125extern int dsint1_(int_4 *n, r_8 *war, r_8 *was, r_8 *xh, r_8 *x, int_8 *ifac);
     126extern int dsinti_(int_4 *n, r_8 *wsave);
    127127
    128128/* Rerunning f2c -P may change prototypes or declarations. */

  • trunk/SophyaLib/NTools/nbmath.h

    r244 r756  
    55
    66#include "machdefs.h"
     7#include "smathconst.h"
    78
    89#ifdef __cplusplus
    … …  
    1011#endif
    1112
    12 /* EA conflit entre les diverses copies de nbmath.h, nbconst.h... A nettoyer... */
    13 #undef  Pi   
    14 #undef  Pis2
    15 #undef  DeuxPi
    16 #undef  SPi   
    17 #undef  S2Pi   
    18 #undef  Rac2   
    19 #undef  Log2   
    20 #undef  LnPi   
    21 #undef  LgPi   
    22 #undef  Ln10   
    23 #undef  DftoDm
    24 #undef  Hln2pi
    25 #undef  JourSec
    26 #undef  AnSec 
    27 
    28 #undef  GRAND 
    29 #undef  GRAND2
    30 #undef  IGRAND
    31 
    32 #define  Pi     (double) ( 3.14159265358979310e+00)  /* c'est Pi */
    33 #define  Pis2   (double) (1.57079632679489655)       /* c'est Pi/2 */
    34 #define  DeuxPi (double) ( 6.28318530717958620e+00)  /* c'est 2*Pi */
    35 #define  SPi    (double) ( 1.77245385090551590e+00)  /* c'est sqrt(Pi) */
    36 #define  S2Pi   (double) ( 2.50662827463100020e+00)  /* c'est sqrt(2*Pi) */
    37 #define  Rac2   (double) ( 1.41421356237309510e+00)  /* c'est sqrt(2) */
    38 #define  Log2   (double) ( 3.01029995663981200e-01)  /* c'est log10(2) */
    39 #define  LnPi   (double) ( 1.14472988584940020e+00)  /* c'est ln(Pi) */
    40 #define  LgPi   (double) ( 4.97149872694133850e-01)  /* c'est log10(Pi) */
    41 #define  Ln10   (double) ( 2.30258509299404590e+00)  /* c'est ln(10) */
    42 #define  DftoDm (double) ( 1.08573620475812940e+00)  /* c'est 2.5/ln(10) */
    43 #define  Hln2pi (double) ( 9.18938533204672670e-01)  /* c'est Ln(2*pi)/2 */
    44 
    45 #define  JourSec (float) (86400.000000) /* nombre de secondes dans 24H */
    46 #define  AnSec  (int) (31557600) /* nombre de secondes dans 365.25 jours */
    47 
    48 #define  GRAND  (float)  (1.e+35)
    49 #define  GRAND2 (double) (1.e+35)
    50 #define  IGRAND (int_4) (2147483647)
    5113
    5214/* pour le fit 1-dim */

  • trunk/SophyaLib/NTools/rzimage.h

    r552 r756  
    103103  virtual void  WriteSelf(POutPersist&) const;
    104104  virtual void  ReadSelf(PInPersist&);
     105  virtual AnyDataObj* DataObj() { return(this); }
     106  virtual void       SetDataObj(AnyDataObj &)
     107  { throw  ForbiddenError("RzImage::SetDataObj() not allowed ! "); }
    105108
    106109//  Fonctions d'acces aux pixels avec conversion
Note: See TracChangeset for help on using the changeset viewer.