Changeset 710 in Sophya for trunk/SophyaLib


Ignore:
Timestamp:
Jan 21, 2000, 6:00:45 PM (26 years ago)
Author:
ansari
Message:

Introduction FFTServerInterface et FFTPackServer - Reza 21/01/2000

Location:
trunk/SophyaLib/NTools
Files:
4 added
4 edited

Legend:

Unmodified
Added
Removed

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

    r608 r710  
    1212fct2dfit.o
    1313fftpackc.o
     14fftpserver.o
    1415fftserver.o
     16fftservintf.o
    1517generaldata.o
    1618generalfit.o

  • trunk/SophyaLib/NTools/cspline.cc

    r520 r710  
    44#include <stdio.h>
    55#include <math.h>
     6#include <iostream.h>
     7
    68#ifdef __MWERKS__
    79   #include <iostream>

  • trunk/SophyaLib/NTools/fftpackc.c

    r663 r710  
    1111
    1212/*  ------ File cfftb.f ------ */
    13 /* Subroutine */ int cfftb_(integer *n, real *c__, real *wsave)
    14 {
    15     extern /* Subroutine */ int cfftb1_(integer *, real *, real *, real *,
     13/* Subroutine */ int cfftb_(integer *n, float *c__, float *wsave)
     14{
     15    extern /* Subroutine */ int cfftb1_(integer *, float *, float *, float *,
    1616            integer *);
    1717    static integer iw1, iw2;
    … …  
    3232
    3333/*  ------ File cfftb1.f ------ */
    34 /* Subroutine */ int cfftb1_(integer *n, real *c__, real *ch, real *wa,
     34/* Subroutine */ int cfftb1_(integer *n, float *c__, float *ch, float *wa,
    3535        integer *ifac)
    3636{
    … …  
    4141    static integer idot, i__;
    4242    extern /* Subroutine */ int passb_(integer *, integer *, integer *,
    43             integer *, integer *, real *, real *, real *, real *, real *,
    44             real *);
     43            integer *, integer *, float *, float *, float *, float *, float *,
     44            float *);
    4545    static integer k1, l1, l2, n2;
    46     extern /* Subroutine */ int passb2_(integer *, integer *, real *, real *,
    47             real *), passb3_(integer *, integer *, real *, real *, real *,
    48             real *), passb4_(integer *, integer *, real *, real *, real *,
    49             real *, real *), passb5_(integer *, integer *, real *, real *,
    50             real *, real *, real *, real *);
     46    extern /* Subroutine */ int passb2_(integer *, integer *, float *, float *,
     47            float *), passb3_(integer *, integer *, float *, float *, float *,
     48            float *), passb4_(integer *, integer *, float *, float *, float *,
     49            float *, float *), passb5_(integer *, integer *, float *, float *,
     50            float *, float *, float *, float *);
    5151    static integer na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
    5252
    … …  
    164164
    165165/*  ------ File cfftf.f ------ */
    166 /* Subroutine */ int cfftf_(integer *n, real *c__, real *wsave)
    167 {
    168     extern /* Subroutine */ int cfftf1_(integer *, real *, real *, real *,
     166/* Subroutine */ int cfftf_(integer *n, float *c__, float *wsave)
     167{
     168    extern /* Subroutine */ int cfftf1_(integer *, float *, float *, float *,
    169169            integer *);
    170170    static integer iw1, iw2;
    … …  
    185185
    186186/*  ------ File cfftf1.f ------ */
    187 /* Subroutine */ int cfftf1_(integer *n, real *c__, real *ch, real *wa,
     187/* Subroutine */ int cfftf1_(integer *n, float *c__, float *ch, float *wa,
    188188        integer *ifac)
    189189{
    … …  
    194194    static integer idot, i__;
    195195    extern /* Subroutine */ int passf_(integer *, integer *, integer *,
    196             integer *, integer *, real *, real *, real *, real *, real *,
    197             real *);
     196            integer *, integer *, float *, float *, float *, float *, float *,
     197            float *);
    198198    static integer k1, l1, l2, n2;
    199     extern /* Subroutine */ int passf2_(integer *, integer *, real *, real *,
    200             real *), passf3_(integer *, integer *, real *, real *, real *,
    201             real *), passf4_(integer *, integer *, real *, real *, real *,
    202             real *, real *), passf5_(integer *, integer *, real *, real *,
    203             real *, real *, real *, real *);
     199    extern /* Subroutine */ int passf2_(integer *, integer *, float *, float *,
     200            float *), passf3_(integer *, integer *, float *, float *, float *,
     201            float *), passf4_(integer *, integer *, float *, float *, float *,
     202            float *, float *), passf5_(integer *, integer *, float *, float *,
     203            float *, float *, float *, float *);
    204204    static integer na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
    205205
    … …  
    317317
    318318/*  ------ File cffti.f ------ */
    319 /* Subroutine */ int cffti_(integer *n, real *wsave)
    320 {
    321     extern /* Subroutine */ int cffti1_(integer *, real *, integer *);
     319/* Subroutine */ int cffti_(integer *n, float *wsave)
     320{
     321    extern /* Subroutine */ int cffti1_(integer *, float *, integer *);
    322322    static integer iw1, iw2;
    323323
    … …  
    336336
    337337/*  ------ File cffti1.f ------ */
    338 /* Subroutine */ int cffti1_(integer *n, real *wa, integer *ifac)
     338/* Subroutine */ int cffti1_(integer *n, float *wa, integer *ifac)
    339339{
    340340    /* Initialized data */
    … …  
    346346
    347347    /* Builtin functions */
    348 /*    double cos(doublereal), sin(doublereal);   commente, remplace par math.h - Reza  29/11/99 */
     348/*    double cos(doubletruc), sin(doubletruc);   commente, remplace par math.h - Reza  29/11/99 */
    349349
    350350    /* Local variables */
    351     static real argh;
     351    static float argh;
    352352    static integer idot, ntry, i__, j;
    353     static real argld;
     353    static float argld;
    354354    static integer i1, k1, l1, l2, ib;
    355     static real fi;
     355    static float fi;
    356356    static integer ld, ii, nf, ip, nl, nq, nr;
    357     static real arg;
     357    static float arg;
    358358    static integer ido, ipm;
    359     static real tpi;
     359    static float tpi;
    360360
    361361    /* Parameter adjustments */
    … …  
    411411    ifac[2] = nf;
    412412    tpi = 6.28318530717959f;
    413     argh = tpi / (real) (*n);
     413    argh = tpi / (float) (*n);
    414414    i__ = 2;
    415415    l1 = 1;
    … …  
    429429            ld += l1;
    430430            fi = 0.f;
    431             argld = (real) ld * argh;
     431            argld = (float) ld * argh;
    432432            i__3 = idot;
    433433            for (ii = 4; ii <= i__3; ii += 2) {
    … …  
    454454
    455455/*  ------ File cosqb.f ------ */
    456 /* Subroutine */ int cosqb_(integer *n, real *x, real *wsave)
     456/* Subroutine */ int cosqb_(integer *n, float *x, float *wsave)
    457457{
    458458    /* Initialized data */
    459459
    460     static real tsqrt2 = 2.82842712474619f;
     460    static float tsqrt2 = 2.82842712474619f;
    461461
    462462    /* System generated locals */
    … …  
    464464
    465465    /* Local variables */
    466     static real x1;
    467     extern /* Subroutine */ int cosqb1_(integer *, real *, real *, real *);
     466    static float x1;
     467    extern /* Subroutine */ int cosqb1_(integer *, float *, float *, float *);
    468468
    469469    /* Parameter adjustments */
    … …  
    493493
    494494/*  ------ File cosqb1.f ------ */
    495 /* Subroutine */ int cosqb1_(integer *n, real *x, real *w, real *xh)
     495/* Subroutine */ int cosqb1_(integer *n, float *x, float *w, float *xh)
    496496{
    497497    /* System generated locals */
    … …  
    500500    /* Local variables */
    501501    static integer modn, i__, k;
    502     extern /* Subroutine */ int rfftb_(integer *, real *, real *);
     502    extern /* Subroutine */ int rfftb_(integer *, float *, float *);
    503503    static integer kc, np2, ns2;
    504     static real xim1;
     504    static float xim1;
    505505
    506506    /* Parameter adjustments */
    … …  
    547547
    548548/*  ------ File cosqf.f ------ */
    549 /* Subroutine */ int cosqf_(integer *n, real *x, real *wsave)
     549/* Subroutine */ int cosqf_(integer *n, float *x, float *wsave)
    550550{
    551551    /* Initialized data */
    552552
    553     static real sqrt2 = 1.4142135623731f;
     553    static float sqrt2 = 1.4142135623731f;
    554554
    555555    /* System generated locals */
    … …  
    557557
    558558    /* Local variables */
    559     static real tsqx;
    560     extern /* Subroutine */ int cosqf1_(integer *, real *, real *, real *);
     559    static float tsqx;
     560    extern /* Subroutine */ int cosqf1_(integer *, float *, float *, float *);
    561561
    562562    /* Parameter adjustments */
    … …  
    584584
    585585/*  ------ File cosqf1.f ------ */
    586 /* Subroutine */ int cosqf1_(integer *n, real *x, real *w, real *xh)
     586/* Subroutine */ int cosqf1_(integer *n, float *x, float *w, float *xh)
    587587{
    588588    /* System generated locals */
    … …  
    591591    /* Local variables */
    592592    static integer modn, i__, k;
    593     extern /* Subroutine */ int rfftf_(integer *, real *, real *);
     593    extern /* Subroutine */ int rfftf_(integer *, float *, float *);
    594594    static integer kc, np2, ns2;
    595     static real xim1;
     595    static float xim1;
    596596
    597597    /* Parameter adjustments */
    … …  
    636636
    637637/*  ------ File cosqi.f ------ */
    638 /* Subroutine */ int cosqi_(integer *n, real *wsave)
     638/* Subroutine */ int cosqi_(integer *n, float *wsave)
    639639{
    640640    /* Initialized data */
    641641
    642     static real pih = 1.57079632679491f;
     642    static float pih = 1.57079632679491f;
    643643
    644644    /* System generated locals */
    … …  
    646646
    647647    /* Builtin functions */
    648 /*    double cos(doublereal);    commente - Remplace par math.h  Reza 29/11/99 */
     648/*    double cos(doubletruc);    commente - Remplace par math.h  Reza 29/11/99 */
    649649
    650650    /* Local variables */
    651651    static integer k;
    652     extern /* Subroutine */ int rffti_(integer *, real *);
    653     static real fk, dt;
     652    extern /* Subroutine */ int rffti_(integer *, float *);
     653    static float fk, dt;
    654654
    655655    /* Parameter adjustments */
    … …  
    657657
    658658    /* Function Body */
    659     dt = pih / (real) (*n);
     659    dt = pih / (float) (*n);
    660660    fk = 0.f;
    661661    i__1 = *n;
    … …  
    670670
    671671/*  ------ File cost.f ------ */
    672 /* Subroutine */ int cost_(integer *n, real *x, real *wsave)
     672/* Subroutine */ int cost_(integer *n, float *x, float *wsave)
    673673{
    674674    /* System generated locals */
    … …  
    677677    /* Local variables */
    678678    static integer modn, i__, k;
    679     extern /* Subroutine */ int rfftf_(integer *, real *, real *);
    680     static real c1, t1, t2;
     679    extern /* Subroutine */ int rfftf_(integer *, float *, float *);
     680    static float c1, t1, t2;
    681681    static integer kc;
    682     static real xi;
     682    static float xi;
    683683    static integer nm1, np1;
    684     static real x1h;
     684    static float x1h;
    685685    static integer ns2;
    686     static real tx2, x1p3, xim2;
     686    static float tx2, x1p3, xim2;
    687687
    688688    /* Parameter adjustments */
    … …  
    753753
    754754/*  ------ File costi.f ------ */
    755 /* Subroutine */ int costi_(integer *n, real *wsave)
     755/* Subroutine */ int costi_(integer *n, float *wsave)
    756756{
    757757    /* Initialized data */
    758758
    759     static real pi = 3.14159265358979f;
     759    static float pi = 3.14159265358979f;
    760760
    761761    /* System generated locals */
    … …  
    763763
    764764    /* Builtin functions */
    765 /*    double sin(doublereal), cos(doublereal);  commente - Remplace par math.h  Reza 29/11/99 */
     765/*    double sin(doubletruc), cos(doubletruc);  commente - Remplace par math.h  Reza 29/11/99 */
    766766
    767767    /* Local variables */
    768768    static integer k;
    769     extern /* Subroutine */ int rffti_(integer *, real *);
     769    extern /* Subroutine */ int rffti_(integer *, float *);
    770770    static integer kc;
    771     static real fk, dt;
     771    static float fk, dt;
    772772    static integer nm1, np1, ns2;
    773773
    … …  
    782782    np1 = *n + 1;
    783783    ns2 = *n / 2;
    784     dt = pi / (real) nm1;
     784    dt = pi / (float) nm1;
    785785    fk = 0.f;
    786786    i__1 = ns2;
    … …  
    797797
    798798/*  ------ File ezfft1.f ------ */
    799 /* Subroutine */ int ezfft1_(integer *n, real *wa, integer *ifac)
     799/* Subroutine */ int ezfft1_(integer *n, float *wa, integer *ifac)
    800800{
    801801    /* Initialized data */
    802802
    803803    static integer ntryh[4] = { 4,2,3,5 };
    804     static real tpi = 6.28318530717959f;
     804    static float tpi = 6.28318530717959f;
    805805
    806806    /* System generated locals */
    … …  
    808808
    809809    /* Builtin functions */
    810 /*    double cos(doublereal), sin(doublereal);  commente - Remplace par math.h  Reza 29/11/99 */
     810/*    double cos(doubletruc), sin(doubletruc);  commente - Remplace par math.h  Reza 29/11/99 */
    811811
    812812    /* Local variables */
    813     static real argh;
     813    static float argh;
    814814    static integer ntry, i__, j, k1, l1, l2, ib, ii, nf, ip, nl, is, nq, nr;
    815     static real ch1, sh1;
     815    static float ch1, sh1;
    816816    static integer ido, ipm;
    817     static real dch1, ch1h, arg1, dsh1;
     817    static float dch1, ch1h, arg1, dsh1;
    818818    static integer nfm1;
    819819
    … …  
    869869    ifac[1] = *n;
    870870    ifac[2] = nf;
    871     argh = tpi / (real) (*n);
     871    argh = tpi / (float) (*n);
    872872    is = 0;
    873873    nfm1 = nf - 1;
    … …  
    882882        ido = *n / l2;
    883883        ipm = ip - 1;
    884         arg1 = (real) l1 * argh;
     884        arg1 = (float) l1 * argh;
    885885        ch1 = 1.f;
    886886        sh1 = 0.f;
    … …  
    916916
    917917/*  ------ File ezfftb.f ------ */
    918 /* Subroutine */ int ezfftb_(integer *n, real *r__, real *azero, real *a,
    919         real *b, real *wsave)
     918/* Subroutine */ int ezfftb_(integer *n, float *r__, float *azero, float *a,
     919        float *b, float *wsave)
    920920{
    921921    /* System generated locals */
    … …  
    924924    /* Local variables */
    925925    static integer i__;
    926     extern /* Subroutine */ int rfftb_(integer *, real *, real *);
     926    extern /* Subroutine */ int rfftb_(integer *, float *, float *);
    927927    static integer ns2;
    928928
    … …  
    965965
    966966/*  ------ File ezfftf.f ------ */
    967 /* Subroutine */ int ezfftf_(integer *n, real *r__, real *azero, real *a,
    968         real *b, real *wsave)
     967/* Subroutine */ int ezfftf_(integer *n, float *r__, float *azero, float *a,
     968        float *b, float *wsave)
    969969{
    970970    /* System generated locals */
    … …  
    973973    /* Local variables */
    974974    static integer i__;
    975     extern /* Subroutine */ int rfftf_(integer *, real *, real *);
    976     static real cf;
     975    extern /* Subroutine */ int rfftf_(integer *, float *, float *);
     976    static float cf;
    977977    static integer ns2;
    978     static real cfm;
     978    static float cfm;
    979979    static integer ns2m;
    980980
    … …  
    10101010    }
    10111011    rfftf_(n, &wsave[1], &wsave[*n + 1]);
    1012     cf = 2.f / (real) (*n);
     1012    cf = 2.f / (float) (*n);
    10131013    cfm = -cf;
    10141014    *azero = cf * .5f * wsave[1];
    … …  
    10301030
    10311031/*  ------ File ezffti.f ------ */
    1032 /* Subroutine */ int ezffti_(integer *n, real *wsave)
    1033 {
    1034     extern /* Subroutine */ int ezfft1_(integer *, real *, integer *);
     1032/* Subroutine */ int ezffti_(integer *n, float *wsave)
     1033{
     1034    extern /* Subroutine */ int ezfft1_(integer *, float *, integer *);
    10351035
    10361036    /* Parameter adjustments */
    … …  
    10471047/*  ------ File passb.f ------ */
    10481048/* Subroutine */ int passb_(integer *nac, integer *ido, integer *ip, integer *
    1049         l1, integer *idl1, real *cc, real *c1, real *c2, real *ch, real *ch2,
    1050         real *wa)
     1049        l1, integer *idl1, float *cc, float *c1, float *c2, float *ch, float *ch2,
     1050        float *wa)
    10511051{
    10521052    /* System generated locals */
    … …  
    10581058    static integer idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
    10591059            idl, inc, idp;
    1060     static real wai, war;
     1060    static float wai, war;
    10611061    static integer ipp2;
    10621062
    … …  
    12871287
    12881288/*  ------ File passb2.f ------ */
    1289 /* Subroutine */ int passb2_(integer *ido, integer *l1, real *cc, real *ch,
    1290         real *wa1)
     1289/* Subroutine */ int passb2_(integer *ido, integer *l1, float *cc, float *ch,
     1290        float *wa1)
    12911291{
    12921292    /* System generated locals */
    … …  
    12951295    /* Local variables */
    12961296    static integer i__, k;
    1297     static real ti2, tr2;
     1297    static float ti2, tr2;
    12981298
    12991299    /* Parameter adjustments */
    … …  
    13491349
    13501350/*  ------ File passb3.f ------ */
    1351 /* Subroutine */ int passb3_(integer *ido, integer *l1, real *cc, real *ch,
    1352         real *wa1, real *wa2)
     1351/* Subroutine */ int passb3_(integer *ido, integer *l1, float *cc, float *ch,
     1352        float *wa1, float *wa2)
    13531353{
    13541354    /* Initialized data */
    13551355
    1356     static real taur = -.5f;
    1357     static real taui = .866025403784439f;
     1356    static float taur = -.5f;
     1357    static float taui = .866025403784439f;
    13581358
    13591359    /* System generated locals */
    … …  
    13621362    /* Local variables */
    13631363    static integer i__, k;
    1364     static real ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
     1364    static float ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    13651365
    13661366    /* Parameter adjustments */
    … …  
    14371437
    14381438/*  ------ File passb4.f ------ */
    1439 /* Subroutine */ int passb4_(integer *ido, integer *l1, real *cc, real *ch,
    1440         real *wa1, real *wa2, real *wa3)
     1439/* Subroutine */ int passb4_(integer *ido, integer *l1, float *cc, float *ch,
     1440        float *wa1, float *wa2, float *wa3)
    14411441{
    14421442    /* System generated locals */
    … …  
    14451445    /* Local variables */
    14461446    static integer i__, k;
    1447     static real ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
     1447    static float ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    14481448            tr3, tr4;
    14491449
    … …  
    15421542
    15431543/*  ------ File passb5.f ------ */
    1544 /* Subroutine */ int passb5_(integer *ido, integer *l1, real *cc, real *ch,
    1545         real *wa1, real *wa2, real *wa3, real *wa4)
     1544/* Subroutine */ int passb5_(integer *ido, integer *l1, float *cc, float *ch,
     1545        float *wa1, float *wa2, float *wa3, float *wa4)
    15461546{
    15471547    /* Initialized data */
    15481548
    1549     static real tr11 = .309016994374947f;
    1550     static real ti11 = .951056516295154f;
    1551     static real tr12 = -.809016994374947f;
    1552     static real ti12 = .587785252292473f;
     1549    static float tr11 = .309016994374947f;
     1550    static float ti11 = .951056516295154f;
     1551    static float tr12 = -.809016994374947f;
     1552    static float ti12 = .587785252292473f;
    15531553
    15541554    /* System generated locals */
    … …  
    15571557    /* Local variables */
    15581558    static integer i__, k;
    1559     static real ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
     1559    static float 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;
    15611561
    … …  
    16781678/*  ------ File passf.f ------ */
    16791679/* Subroutine */ int passf_(integer *nac, integer *ido, integer *ip, integer *
    1680         l1, integer *idl1, real *cc, real *c1, real *c2, real *ch, real *ch2,
    1681         real *wa)
     1680        l1, integer *idl1, float *cc, float *c1, float *c2, float *ch, float *ch2,
     1681        float *wa)
    16821682{
    16831683    /* System generated locals */
    … …  
    16891689    static integer idij, idlj, idot, ipph, i__, j, k, l, jc, lc, ik, nt, idj,
    16901690            idl, inc, idp;
    1691     static real wai, war;
     1691    static float wai, war;
    16921692    static integer ipp2;
    16931693
    … …  
    19181918
    19191919/*  ------ File passf2.f ------ */
    1920 /* Subroutine */ int passf2_(integer *ido, integer *l1, real *cc, real *ch,
    1921         real *wa1)
     1920/* Subroutine */ int passf2_(integer *ido, integer *l1, float *cc, float *ch,
     1921        float *wa1)
    19221922{
    19231923    /* System generated locals */
    … …  
    19261926    /* Local variables */
    19271927    static integer i__, k;
    1928     static real ti2, tr2;
     1928    static float ti2, tr2;
    19291929
    19301930    /* Parameter adjustments */
    … …  
    19801980
    19811981/*  ------ File passf3.f ------ */
    1982 /* Subroutine */ int passf3_(integer *ido, integer *l1, real *cc, real *ch,
    1983         real *wa1, real *wa2)
     1982/* Subroutine */ int passf3_(integer *ido, integer *l1, float *cc, float *ch,
     1983        float *wa1, float *wa2)
    19841984{
    19851985    /* Initialized data */
    19861986
    1987     static real taur = -.5f;
    1988     static real taui = -.866025403784439f;
     1987    static float taur = -.5f;
     1988    static float taui = -.866025403784439f;
    19891989
    19901990    /* System generated locals */
    … …  
    19931993    /* Local variables */
    19941994    static integer i__, k;
    1995     static real ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
     1995    static float ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    19961996
    19971997    /* Parameter adjustments */
    … …  
    20682068
    20692069/*  ------ File passf4.f ------ */
    2070 /* Subroutine */ int passf4_(integer *ido, integer *l1, real *cc, real *ch,
    2071         real *wa1, real *wa2, real *wa3)
     2070/* Subroutine */ int passf4_(integer *ido, integer *l1, float *cc, float *ch,
     2071        float *wa1, float *wa2, float *wa3)
    20722072{
    20732073    /* System generated locals */
    … …  
    20762076    /* Local variables */
    20772077    static integer i__, k;
    2078     static real ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
     2078    static float ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    20792079            tr3, tr4;
    20802080
    … …  
    21732173
    21742174/*  ------ File passf5.f ------ */
    2175 /* Subroutine */ int passf5_(integer *ido, integer *l1, real *cc, real *ch,
    2176         real *wa1, real *wa2, real *wa3, real *wa4)
     2175/* Subroutine */ int passf5_(integer *ido, integer *l1, float *cc, float *ch,
     2176        float *wa1, float *wa2, float *wa3, float *wa4)
    21772177{
    21782178    /* Initialized data */
    21792179
    2180     static real tr11 = .309016994374947f;
    2181     static real ti11 = -.951056516295154f;
    2182     static real tr12 = -.809016994374947f;
    2183     static real ti12 = -.587785252292473f;
     2180    static float tr11 = .309016994374947f;
     2181    static float ti11 = -.951056516295154f;
     2182    static float tr12 = -.809016994374947f;
     2183    static float ti12 = -.587785252292473f;
    21842184
    21852185    /* System generated locals */
    … …  
    21882188    /* Local variables */
    21892189    static integer i__, k;
    2190     static real ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
     2190    static float 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;
    21922192
    … …  
    23082308
    23092309/*  ------ File radb2.f ------ */
    2310 /* Subroutine */ int radb2_(integer *ido, integer *l1, real *cc, real *ch,
    2311         real *wa1)
     2310/* Subroutine */ int radb2_(integer *ido, integer *l1, float *cc, float *ch,
     2311        float *wa1)
    23122312{
    23132313    /* System generated locals */
    … …  
    23162316    /* Local variables */
    23172317    static integer i__, k, ic;
    2318     static real ti2, tr2;
     2318    static float ti2, tr2;
    23192319    static integer idp2;
    23202320
    … …  
    23852385
    23862386/*  ------ File radb3.f ------ */
    2387 /* Subroutine */ int radb3_(integer *ido, integer *l1, real *cc, real *ch,
    2388         real *wa1, real *wa2)
     2387/* Subroutine */ int radb3_(integer *ido, integer *l1, float *cc, float *ch,
     2388        float *wa1, float *wa2)
    23892389{
    23902390    /* Initialized data */
    23912391
    2392     static real taur = -.5f;
    2393     static real taui = .866025403784439f;
     2392    static float taur = -.5f;
     2393    static float taui = .866025403784439f;
    23942394
    23952395    /* System generated locals */
    … …  
    23982398    /* Local variables */
    23992399    static integer i__, k, ic;
    2400     static real ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
     2400    static float ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
    24012401    static integer idp2;
    24022402
    … …  
    24682468
    24692469/*  ------ File radb4.f ------ */
    2470 /* Subroutine */ int radb4_(integer *ido, integer *l1, real *cc, real *ch,
    2471         real *wa1, real *wa2, real *wa3)
     2470/* Subroutine */ int radb4_(integer *ido, integer *l1, float *cc, float *ch,
     2471        float *wa1, float *wa2, float *wa3)
    24722472{
    24732473    /* Initialized data */
    24742474
    2475     static real sqrt2 = 1.414213562373095f;
     2475    static float sqrt2 = 1.414213562373095f;
    24762476
    24772477    /* System generated locals */
    … …  
    24802480    /* Local variables */
    24812481    static integer i__, k, ic;
    2482     static real ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
     2482    static float ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    24832483            tr3, tr4;
    24842484    static integer idp2;
    … …  
    25922592
    25932593/*  ------ File radb5.f ------ */
    2594 /* Subroutine */ int radb5_(integer *ido, integer *l1, real *cc, real *ch,
    2595         real *wa1, real *wa2, real *wa3, real *wa4)
     2594/* Subroutine */ int radb5_(integer *ido, integer *l1, float *cc, float *ch,
     2595        float *wa1, float *wa2, float *wa3, float *wa4)
    25962596{
    25972597    /* Initialized data */
    25982598
    2599     static real tr11 = .309016994374947f;
    2600     static real ti11 = .951056516295154f;
    2601     static real tr12 = -.809016994374947f;
    2602     static real ti12 = .587785252292473f;
     2599    static float tr11 = .309016994374947f;
     2600    static float ti11 = .951056516295154f;
     2601    static float tr12 = -.809016994374947f;
     2602    static float ti12 = .587785252292473f;
    26032603
    26042604    /* System generated locals */
    … …  
    26072607    /* Local variables */
    26082608    static integer i__, k, ic;
    2609     static real ci2, ci3, ci4, ci5, di3, di4, di5, di2, cr2, cr3, cr5, cr4,
     2609    static float 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;
    26112611    static integer idp2;
    … …  
    27172717/*  ------ File radbg.f ------ */
    27182718/* Subroutine */ int radbg_(integer *ido, integer *ip, integer *l1, integer *
    2719         idl1, real *cc, real *c1, real *c2, real *ch, real *ch2, real *wa)
     2719        idl1, float *cc, float *c1, float *c2, float *ch, float *ch2, float *wa)
    27202720{
    27212721    /* Initialized data */
    27222722
    2723     static real tpi = 6.28318530717959f;
     2723    static float tpi = 6.28318530717959f;
    27242724
    27252725    /* System generated locals */
    … …  
    27292729
    27302730    /* Builtin functions */
    2731 /*    double cos(doublereal), sin(doublereal);  commente - Remplace par math.h  Reza 29/11/99 */
     2731/*    double cos(doubletruc), sin(doubletruc);  commente - Remplace par math.h  Reza 29/11/99 */
    27322732
    27332733    /* Local variables */
    27342734    static integer idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
    2735     static real dc2, ai1, ai2, ar1, ar2, ds2;
     2735    static float dc2, ai1, ai2, ar1, ar2, ds2;
    27362736    static integer nbd;
    2737     static real dcp, arg, dsp, ar1h, ar2h;
     2737    static float dcp, arg, dsp, ar1h, ar2h;
    27382738    static integer idp2, ipp2;
    27392739
    … …  
    27602760
    27612761    /* Function Body */
    2762     arg = tpi / (real) (*ip);
     2762    arg = tpi / (float) (*ip);
    27632763    dcp = cos(arg);
    27642764    dsp = sin(arg);
    … …  
    30573057
    30583058/*  ------ File radf2.f ------ */
    3059 /* Subroutine */ int radf2_(integer *ido, integer *l1, real *cc, real *ch,
    3060         real *wa1)
     3059/* Subroutine */ int radf2_(integer *ido, integer *l1, float *cc, float *ch,
     3060        float *wa1)
    30613061{
    30623062    /* System generated locals */
    … …  
    30653065    /* Local variables */
    30663066    static integer i__, k, ic;
    3067     static real ti2, tr2;
     3067    static float ti2, tr2;
    30683068    static integer idp2;
    30693069
    … …  
    31353135
    31363136/*  ------ File radf3.f ------ */
    3137 /* Subroutine */ int radf3_(integer *ido, integer *l1, real *cc, real *ch,
    3138         real *wa1, real *wa2)
     3137/* Subroutine */ int radf3_(integer *ido, integer *l1, float *cc, float *ch,
     3138        float *wa1, float *wa2)
    31393139{
    31403140    /* Initialized data */
    31413141
    3142     static real taur = -.5f;
    3143     static real taui = .866025403784439f;
     3142    static float taur = -.5f;
     3143    static float taui = .866025403784439f;
    31443144
    31453145    /* System generated locals */
    … …  
    31483148    /* Local variables */
    31493149    static integer i__, k, ic;
    3150     static real ci2, di2, di3, cr2, dr2, dr3, ti2, ti3, tr2, tr3;
     3150    static float ci2, di2, di3, cr2, dr2, dr3, ti2, ti3, tr2, tr3;
    31513151    static integer idp2;
    31523152
    … …  
    32143214
    32153215/*  ------ File radf4.f ------ */
    3216 /* Subroutine */ int radf4_(integer *ido, integer *l1, real *cc, real *ch,
    3217         real *wa1, real *wa2, real *wa3)
     3216/* Subroutine */ int radf4_(integer *ido, integer *l1, float *cc, float *ch,
     3217        float *wa1, float *wa2, float *wa3)
    32183218{
    32193219    /* Initialized data */
    32203220
    3221     static real hsqt2 = .7071067811865475f;
     3221    static float hsqt2 = .7071067811865475f;
    32223222
    32233223    /* System generated locals */
    … …  
    32263226    /* Local variables */
    32273227    static integer i__, k, ic;
    3228     static real ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
     3228    static float ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2,
    32293229            tr3, tr4;
    32303230    static integer idp2;
    … …  
    33303330
    33313331/*  ------ File radf5.f ------ */
    3332 /* Subroutine */ int radf5_(integer *ido, integer *l1, real *cc, real *ch,
    3333         real *wa1, real *wa2, real *wa3, real *wa4)
     3332/* Subroutine */ int radf5_(integer *ido, integer *l1, float *cc, float *ch,
     3333        float *wa1, float *wa2, float *wa3, float *wa4)
    33343334{
    33353335    /* Initialized data */
    33363336
    3337     static real tr11 = .309016994374947f;
    3338     static real ti11 = .951056516295154f;
    3339     static real tr12 = -.809016994374947f;
    3340     static real ti12 = .587785252292473f;
     3337    static float tr11 = .309016994374947f;
     3338    static float ti11 = .951056516295154f;
     3339    static float tr12 = -.809016994374947f;
     3340    static float ti12 = .587785252292473f;
    33413341
    33423342    /* System generated locals */
    … …  
    33453345    /* Local variables */
    33463346    static integer i__, k, ic;
    3347     static real ci2, di2, ci4, ci5, di3, di4, di5, ci3, cr2, cr3, dr2, dr3,
     3347    static float 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;
    33493349    static integer idp2;
    … …  
    34493449/*  ------ File radfg.f ------ */
    34503450/* Subroutine */ int radfg_(integer *ido, integer *ip, integer *l1, integer *
    3451         idl1, real *cc, real *c1, real *c2, real *ch, real *ch2, real *wa)
     3451        idl1, float *cc, float *c1, float *c2, float *ch, float *ch2, float *wa)
    34523452{
    34533453    /* Initialized data */
    34543454
    3455     static real tpi = 6.28318530717959f;
     3455    static float tpi = 6.28318530717959f;
    34563456
    34573457    /* System generated locals */
    … …  
    34613461
    34623462    /* Builtin functions */
    3463 /*    double cos(doublereal), sin(doublereal);    */
     3463/*    double cos(doubletruc), sin(doubletruc);    */
    34643464
    34653465    /* Local variables */
    34663466    static integer idij, ipph, i__, j, k, l, j2, ic, jc, lc, ik, is;
    3467     static real dc2, ai1, ai2, ar1, ar2, ds2;
     3467    static float dc2, ai1, ai2, ar1, ar2, ds2;
    34683468    static integer nbd;
    3469     static real dcp, arg, dsp, ar1h, ar2h;
     3469    static float dcp, arg, dsp, ar1h, ar2h;
    34703470    static integer idp2, ipp2;
    34713471
    … …  
    34923492
    34933493    /* Function Body */
    3494     arg = tpi / (real) (*ip);
     3494    arg = tpi / (float) (*ip);
    34953495    dcp = cos(arg);
    34963496    dsp = sin(arg);
    … …  
    37953795
    37963796/*  ------ File rfftb.f ------ */
    3797 /* Subroutine */ int rfftb_(integer *n, real *r__, real *wsave)
    3798 {
    3799     extern /* Subroutine */ int rfftb1_(integer *, real *, real *, real *,
     3797/* Subroutine */ int rfftb_(integer *n, float *r__, float *wsave)
     3798{
     3799    extern /* Subroutine */ int rfftb1_(integer *, float *, float *, float *,
    38003800            integer *);
    38013801
    … …  
    38133813
    38143814/*  ------ File rfftb1.f ------ */
    3815 /* Subroutine */ int rfftb1_(integer *n, real *c__, real *ch, real *wa,
     3815/* Subroutine */ int rfftb1_(integer *n, float *c__, float *ch, float *wa,
    38163816        integer *ifac)
    38173817{
    … …  
    38203820
    38213821    /* Local variables */
    3822     extern /* Subroutine */ int radb2_(integer *, integer *, real *, real *,
    3823             real *), radb3_(integer *, integer *, real *, real *, real *,
    3824             real *), radb4_(integer *, integer *, real *, real *, real *,
    3825             real *, real *), radb5_(integer *, integer *, real *, real *,
    3826             real *, real *, real *, real *);
     3822    extern /* Subroutine */ int radb2_(integer *, integer *, float *, float *,
     3823            float *), radb3_(integer *, integer *, float *, float *, float *,
     3824            float *), radb4_(integer *, integer *, float *, float *, float *,
     3825            float *, float *), radb5_(integer *, integer *, float *, float *,
     3826            float *, float *, float *, float *);
    38273827    static integer i__;
    38283828    extern /* Subroutine */ int radbg_(integer *, integer *, integer *,
    3829             integer *, real *, real *, real *, real *, real *, real *);
     3829            integer *, float *, float *, float *, float *, float *, float *);
    38303830    static integer k1, l1, l2, na, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
    38313831
    … …  
    39413941
    39423942/*  ------ File rfftf.f ------ */
    3943 /* Subroutine */ int rfftf_(integer *n, real *r__, real *wsave)
    3944 {
    3945     extern /* Subroutine */ int rfftf1_(integer *, real *, real *, real *,
     3943/* Subroutine */ int rfftf_(integer *n, float *r__, float *wsave)
     3944{
     3945    extern /* Subroutine */ int rfftf1_(integer *, float *, float *, float *,
    39463946            integer *);
    39473947
    … …  
    39593959
    39603960/*  ------ File rfftf1.f ------ */
    3961 /* Subroutine */ int rfftf1_(integer *n, real *c__, real *ch, real *wa,
     3961/* Subroutine */ int rfftf1_(integer *n, float *c__, float *ch, float *wa,
    39623962        integer *ifac)
    39633963{
    … …  
    39663966
    39673967    /* Local variables */
    3968     extern /* Subroutine */ int radf2_(integer *, integer *, real *, real *,
    3969             real *), radf3_(integer *, integer *, real *, real *, real *,
    3970             real *), radf4_(integer *, integer *, real *, real *, real *,
    3971             real *, real *), radf5_(integer *, integer *, real *, real *,
    3972             real *, real *, real *, real *);
     3968    extern /* Subroutine */ int radf2_(integer *, integer *, float *, float *,
     3969            float *), radf3_(integer *, integer *, float *, float *, float *,
     3970            float *), radf4_(integer *, integer *, float *, float *, float *,
     3971            float *, float *), radf5_(integer *, integer *, float *, float *,
     3972            float *, float *, float *, float *);
    39733973    static integer i__;
    39743974    extern /* Subroutine */ int radfg_(integer *, integer *, integer *,
    3975             integer *, real *, real *, real *, real *, real *, real *);
     3975            integer *, float *, float *, float *, float *, float *, float *);
    39763976    static integer k1, l1, l2, na, kh, nf, ip, iw, ix2, ix3, ix4, ido, idl1;
    39773977
    … …  
    40824082
    40834083/*  ------ File rffti.f ------ */
    4084 /* Subroutine */ int rffti_(integer *n, real *wsave)
    4085 {
    4086     extern /* Subroutine */ int rffti1_(integer *, real *, integer *);
     4084/* Subroutine */ int rffti_(integer *n, float *wsave)
     4085{
     4086    extern /* Subroutine */ int rffti1_(integer *, float *, integer *);
    40874087
    40884088    /* Parameter adjustments */
    … …  
    40984098
    40994099/*  ------ File rffti1.f ------ */
    4100 /* Subroutine */ int rffti1_(integer *n, real *wa, integer *ifac)
     4100/* Subroutine */ int rffti1_(integer *n, float *wa, integer *ifac)
    41014101{
    41024102    /* Initialized data */
    … …  
    41084108
    41094109    /* Builtin functions */
    4110 /*   double cos(doublereal), sin(doublereal);  Remplace par math.h , Reza 29/11/99 */
     4110/*   double cos(doubletruc), sin(doubletruc);  Remplace par math.h , Reza 29/11/99 */
    41114111
    41124112    /* Local variables */
    4113     static real argh;
     4113    static float argh;
    41144114    static integer ntry, i__, j;
    4115     static real argld;
     4115    static float argld;
    41164116    static integer k1, l1, l2, ib;
    4117     static real fi;
     4117    static float fi;
    41184118    static integer ld, ii, nf, ip, nl, is, nq, nr;
    4119     static real arg;
     4119    static float arg;
    41204120    static integer ido, ipm;
    4121     static real tpi;
     4121    static float tpi;
    41224122    static integer nfm1;
    41234123
    … …  
    41744174    ifac[2] = nf;
    41754175    tpi = 6.28318530717959f;
    4176     argh = tpi / (real) (*n);
     4176    argh = tpi / (float) (*n);
    41774177    is = 0;
    41784178    nfm1 = nf - 1;
    … …  
    41924192            ld += l1;
    41934193            i__ = is;
    4194             argld = (real) ld * argh;
     4194            argld = (float) ld * argh;
    41954195            fi = 0.f;
    41964196            i__3 = ido;
    … …  
    42134213
    42144214/*  ------ File sinqb.f ------ */
    4215 /* Subroutine */ int sinqb_(integer *n, real *x, real *wsave)
     4215/* Subroutine */ int sinqb_(integer *n, float *x, float *wsave)
    42164216{
    42174217    /* System generated locals */
    … …  
    42204220    /* Local variables */
    42214221    static integer k;
    4222     extern /* Subroutine */ int cosqb_(integer *, real *, real *);
    4223     static real xhold;
     4222    extern /* Subroutine */ int cosqb_(integer *, float *, float *);
     4223    static float xhold;
    42244224    static integer kc, ns2;
    42254225
    … …  
    42544254
    42554255/*  ------ File sinqf.f ------ */
    4256 /* Subroutine */ int sinqf_(integer *n, real *x, real *wsave)
     4256/* Subroutine */ int sinqf_(integer *n, float *x, float *wsave)
    42574257{
    42584258    /* System generated locals */
    … …  
    42614261    /* Local variables */
    42624262    static integer k;
    4263     extern /* Subroutine */ int cosqf_(integer *, real *, real *);
    4264     static real xhold;
     4263    extern /* Subroutine */ int cosqf_(integer *, float *, float *);
     4264    static float xhold;
    42654265    static integer kc, ns2;
    42664266
    … …  
    42924292
    42934293/*  ------ File sinqi.f ------ */
    4294 /* Subroutine */ int sinqi_(integer *n, real *wsave)
    4295 {
    4296     extern /* Subroutine */ int cosqi_(integer *, real *);
     4294/* Subroutine */ int sinqi_(integer *n, float *wsave)
     4295{
     4296    extern /* Subroutine */ int cosqi_(integer *, float *);
    42974297
    42984298    /* Parameter adjustments */
    … …  
    43054305
    43064306/*  ------ File sint.f ------ */
    4307 /* Subroutine */ int sint_(integer *n, real *x, real *wsave)
    4308 {
    4309     extern /* Subroutine */ int sint1_(integer *, real *, real *, real *,
    4310             real *, integer *);
     4307/* Subroutine */ int sint_(integer *n, float *x, float *wsave)
     4308{
     4309    extern /* Subroutine */ int sint1_(integer *, float *, float *, float *,
     4310            float *, integer *);
    43114311    static integer np1, iw1, iw2, iw3;
    43124312
    … …  
    43254325
    43264326/*  ------ File sint1.f ------ */
    4327 /* Subroutine */ int sint1_(integer *n, real *war, real *was, real *xh, real *
     4327/* Subroutine */ int sint1_(integer *n, float *war, float *was, float *xh, float *
    43284328        x, integer *ifac)
    43294329{
    43304330    /* Initialized data */
    43314331
    4332     static real sqrt3 = 1.73205080756888f;
     4332    static float sqrt3 = 1.73205080756888f;
    43334333
    43344334    /* System generated locals */
    … …  
    43374337    /* Local variables */
    43384338    static integer modn, i__, k;
    4339     static real xhold, t1, t2;
    4340     extern /* Subroutine */ int rfftf1_(integer *, real *, real *, real *,
     4339    static float xhold, t1, t2;
     4340    extern /* Subroutine */ int rfftf1_(integer *, float *, float *, float *,
    43414341            integer *);
    43424342    static integer kc, np1, ns2;
    … …  
    44114411
    44124412/*  ------ File sinti.f ------ */
    4413 /* Subroutine */ int sinti_(integer *n, real *wsave)
     4413/* Subroutine */ int sinti_(integer *n, float *wsave)
    44144414{
    44154415    /* Initialized data */
    44164416
    4417     static real pi = 3.14159265358979f;
     4417    static float pi = 3.14159265358979f;
    44184418
    44194419    /* System generated locals */
    … …  
    44214421
    44224422    /* Builtin functions */
    4423 /*    double sin(doublereal);  remplace par math.h , Reza 29/11/99 */
     4423/*    double sin(doubletruc);  remplace par math.h , Reza 29/11/99 */
    44244424
    44254425    /* Local variables */
    44264426    static integer k;
    4427     extern /* Subroutine */ int rffti_(integer *, real *);
    4428     static real dt;
     4427    extern /* Subroutine */ int rffti_(integer *, float *);
     4428    static float dt;
    44294429    static integer np1, ns2;
    44304430
    … …  
    44384438    ns2 = *n / 2;
    44394439    np1 = *n + 1;
    4440     dt = pi / (real) np1;
     4440    dt = pi / (float) np1;
    44414441    i__1 = ns2;
    44424442    for (k = 1; k <= i__1; ++k) {
    … …  
    47854785
    47864786    /* Builtin functions */
    4787 /*    double cos(doublereal), sin(doublereal);  remplace par math.h  Reza 29/11/99 */
     4787/*    double cos(doubletruc), sin(doubletruc);  remplace par math.h  Reza 29/11/99 */
    47884788
    47894789    /* Local variables */
    … …  
    50855085
    50865086    /* Builtin functions */
    5087 /*    double cos(doublereal); remplace par math.h  Reza 29/11/99 */
     5087/*    double cos(doubletruc); remplace par math.h  Reza 29/11/99 */
    50885088
    50895089    /* Local variables */
    … …  
    52025202
    52035203    /* Builtin functions */
    5204 /*    double sin(doublereal), cos(doublereal);  remplace par math.h  Reza 29/11/99 */
     5204/*    double sin(doubletruc), cos(doubletruc);  remplace par math.h  Reza 29/11/99 */
    52055205
    52065206    /* Local variables */
    … …  
    52475247
    52485248    /* Builtin functions */
    5249 /*    double cos(doublereal), sin(doublereal);  remplace par math.h  Reza 29/11/99  */
     5249/*    double cos(doubletruc), sin(doubletruc);  remplace par math.h  Reza 29/11/99  */
    52505250
    52515251    /* Local variables */
    … …  
    71687168
    71697169    /* Builtin functions */
    7170 /*    double cos(doublereal), sin(doublereal);  remplace par math.h  Reza 29/11/99 */
     7170/*    double cos(doubletruc), sin(doubletruc);  remplace par math.h  Reza 29/11/99 */
    71717171
    71727172    /* Local variables */
    … …  
    79007900
    79017901    /* Builtin functions */
    7902 /*    double cos(doublereal), sin(doublereal);   */
     7902/*    double cos(doubletruc), sin(doubletruc);   */
    79037903
    79047904    /* Local variables */
    … …  
    85478547
    85488548    /* Builtin functions */
    8549 /*    double cos(doublereal), sin(doublereal); remplace par math.h  Reza 29/11/99 */
     8549/*    double cos(doubletruc), sin(doubletruc); remplace par math.h  Reza 29/11/99 */
    85508550
    85518551    /* Local variables */
    … …  
    88608860
    88618861    /* Builtin functions */
    8862 /*    double sin(doublereal);  remplace par math.h  Reza 29/11/99 */
     8862/*    double sin(doubletruc);  remplace par math.h  Reza 29/11/99 */
    88638863
    88648864    /* Local variables */

  • trunk/SophyaLib/NTools/fftpackc.h

    r459 r710  
    44#include "machdefs.h"
    55/* #include "f2c.h" -- Remplace par ce qui suit */
     6#ifdef __cplusplus
     7extern "C" {
     8#endif
    69
    710typedef int integer;
    … …  
    912typedef char *address;
    1013typedef short int shortint;
    11 typedef float real;
    12 typedef double doublereal;
    13 typedef struct { float r, i; } complex;
    14 typedef struct { double r, i; } doublecomplex;
     14  // typedef float real;
     15  // typedef double doublereal;
     16  //typedef struct { float r, i; } complex;
     17  //typedef struct { double r, i; } doublecomplex;
    1518
    1619/* typedef int_4 logical;
    … …  
    2124/* ---------------------- */
    2225
    23 extern int cfftb_(integer *n, real *c__, real *wsave);
    24 extern int cfftb1_(integer *n, real *c__, real *ch, real *wa, integer *ifac);
    25 extern int cfftf_(integer *n, real *c__, real *wsave);
    26 extern int cfftf1_(integer *n, real *c__, real *ch, real *wa, integer *ifac);
    27 extern int cffti_(integer *n, real *wsave);
    28 extern int cffti1_(integer *n, real *wa, integer *ifac);
    29 extern int cosqb_(integer *n, real *x, real *wsave);
    30 extern int cosqb1_(integer *n, real *x, real *w, real *xh);
    31 extern int cosqf_(integer *n, real *x, real *wsave);
    32 extern int cosqf1_(integer *n, real *x, real *w, real *xh);
    33 extern int cosqi_(integer *n, real *wsave);
    34 extern int cost_(integer *n, real *x, real *wsave);
    35 extern int costi_(integer *n, real *wsave);
    36 extern int ezfft1_(integer *n, real *wa, integer *ifac);
    37 extern int ezfftb_(integer *n, real *r__, real *azero, real *a, real *b, real *wsave);
    38 extern int ezfftf_(integer *n, real *r__, real *azero, real *a, real *b, real *wsave);
    39 extern int ezffti_(integer *n, real *wsave);
    40 extern int passb_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, real *cc, real *c1, real *c2, real *ch, real *ch2, real *wa);
    41 extern int passb2_(integer *ido, integer *l1, real *cc, real *ch, real *wa1);
    42 extern int passb3_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2);
    43 extern int passb4_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3);
    44 extern int passb5_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3, real *wa4);
    45 extern int passf_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, real *cc, real *c1, real *c2, real *ch, real *ch2, real *wa);
    46 extern int passf2_(integer *ido, integer *l1, real *cc, real *ch, real *wa1);
    47 extern int passf3_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2);
    48 extern int passf4_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3);
    49 extern int passf5_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3, real *wa4);
    50 extern int radb2_(integer *ido, integer *l1, real *cc, real *ch, real *wa1);
    51 extern int radb3_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2);
    52 extern int radb4_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3);
    53 extern int radb5_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3, real *wa4);
    54 extern int radbg_(integer *ido, integer *ip, integer *l1, integer *idl1, real *cc, real *c1, real *c2, real *ch, real *ch2, real *wa);
    55 extern int radf2_(integer *ido, integer *l1, real *cc, real *ch, real *wa1);
    56 extern int radf3_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2);
    57 extern int radf4_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3);
    58 extern int radf5_(integer *ido, integer *l1, real *cc, real *ch, real *wa1, real *wa2, real *wa3, real *wa4);
    59 extern int radfg_(integer *ido, integer *ip, integer *l1, integer *idl1, real *cc, real *c1, real *c2, real *ch, real *ch2, real *wa);
    60 extern int rfftb_(integer *n, real *r__, real *wsave);
    61 extern int rfftb1_(integer *n, real *c__, real *ch, real *wa, integer *ifac);
    62 extern int rfftf_(integer *n, real *r__, real *wsave);
    63 extern int rfftf1_(integer *n, real *c__, real *ch, real *wa, integer *ifac);
    64 extern int rffti_(integer *n, real *wsave);
    65 extern int rffti1_(integer *n, real *wa, integer *ifac);
    66 extern int sinqb_(integer *n, real *x, real *wsave);
    67 extern int sinqf_(integer *n, real *x, real *wsave);
    68 extern int sinqi_(integer *n, real *wsave);
    69 extern int sint_(integer *n, real *x, real *wsave);
    70 extern int sint1_(integer *n, real *war, real *was, real *xh, real *x, integer *ifac);
    71 extern int sinti_(integer *n, real *wsave);
     26extern int cfftb_(integer *n, float *c__, float *wsave);
     27extern int cfftb1_(integer *n, float *c__, float *ch, float *wa, integer *ifac);
     28extern int cfftf_(integer *n, float *c__, float *wsave);
     29extern int cfftf1_(integer *n, float *c__, float *ch, float *wa, integer *ifac);
     30extern int cffti_(integer *n, float *wsave);
     31extern int cffti1_(integer *n, float *wa, integer *ifac);
     32extern int cosqb_(integer *n, float *x, float *wsave);
     33extern int cosqb1_(integer *n, float *x, float *w, float *xh);
     34extern int cosqf_(integer *n, float *x, float *wsave);
     35extern int cosqf1_(integer *n, float *x, float *w, float *xh);
     36extern int cosqi_(integer *n, float *wsave);
     37extern int cost_(integer *n, float *x, float *wsave);
     38extern int costi_(integer *n, float *wsave);
     39extern int ezfft1_(integer *n, float *wa, integer *ifac);
     40extern int ezfftb_(integer *n, float *r__, float *azero, float *a, float *b, float *wsave);
     41extern int ezfftf_(integer *n, float *r__, float *azero, float *a, float *b, float *wsave);
     42extern int ezffti_(integer *n, float *wsave);
     43extern int passb_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, float *cc, float *c1, float *c2, float *ch, float *ch2, float *wa);
     44extern int passb2_(integer *ido, integer *l1, float *cc, float *ch, float *wa1);
     45extern int passb3_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2);
     46extern int passb4_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3);
     47extern int passb5_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3, float *wa4);
     48extern int passf_(integer *nac, integer *ido, integer *ip, integer *l1, integer *idl1, float *cc, float *c1, float *c2, float *ch, float *ch2, float *wa);
     49extern int passf2_(integer *ido, integer *l1, float *cc, float *ch, float *wa1);
     50extern int passf3_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2);
     51extern int passf4_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3);
     52extern int passf5_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3, float *wa4);
     53extern int radb2_(integer *ido, integer *l1, float *cc, float *ch, float *wa1);
     54extern int radb3_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2);
     55extern int radb4_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3);
     56extern int radb5_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3, float *wa4);
     57extern int radbg_(integer *ido, integer *ip, integer *l1, integer *idl1, float *cc, float *c1, float *c2, float *ch, float *ch2, float *wa);
     58extern int radf2_(integer *ido, integer *l1, float *cc, float *ch, float *wa1);
     59extern int radf3_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2);
     60extern int radf4_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3);
     61extern int radf5_(integer *ido, integer *l1, float *cc, float *ch, float *wa1, float *wa2, float *wa3, float *wa4);
     62extern int radfg_(integer *ido, integer *ip, integer *l1, integer *idl1, float *cc, float *c1, float *c2, float *ch, float *ch2, float *wa);
     63extern int rfftb_(integer *n, float *r__, float *wsave);
     64extern int rfftb1_(integer *n, float *c__, float *ch, float *wa, integer *ifac);
     65extern int rfftf_(integer *n, float *r__, float *wsave);
     66extern int rfftf1_(integer *n, float *c__, float *ch, float *wa, integer *ifac);
     67extern int rffti_(integer *n, float *wsave);
     68extern int rffti1_(integer *n, float *wa, integer *ifac);
     69extern int sinqb_(integer *n, float *x, float *wsave);
     70extern int sinqf_(integer *n, float *x, float *wsave);
     71extern int sinqi_(integer *n, float *wsave);
     72extern int sint_(integer *n, float *x, float *wsave);
     73extern int sint1_(integer *n, float *war, float *was, float *xh, float *x, integer *ifac);
     74extern int sinti_(integer *n, float *wsave);
    7275
    7376/* ---------------------- */
    … …  
    124127
    125128/* Rerunning f2c -P may change prototypes or declarations. */
     129#ifdef __cplusplus
     130}
     131#endif
    126132
    127133#endif /*  FFTPACK_H_SEEN */
Note: See TracChangeset for help on using the changeset viewer.