ruby-dcl 1.6.2
Sign up to get free protection for your applications and to get access to all the features.
- data/ChangeLog +321 -0
- data/GenWrapper/Makefile +14 -0
- data/GenWrapper/cproto2init.rb +41 -0
- data/GenWrapper/dcl_narrayed_funcs.rb +181 -0
- data/GenWrapper/dcl_rb_footing +29 -0
- data/GenWrapper/dcl_rb_heading +129 -0
- data/GenWrapper/def.rb +20 -0
- data/GenWrapper/elim_ary_size.rb +398 -0
- data/GenWrapper/etc/Mk_proto +375 -0
- data/GenWrapper/etc/dcl_gen +13 -0
- data/GenWrapper/etc/p_header +63 -0
- data/GenWrapper/etc/p_init +410 -0
- data/GenWrapper/pparse.rb +137 -0
- data/GenWrapper/proto/Makefile +35 -0
- data/GenWrapper/proto/grph1_00 +17 -0
- data/GenWrapper/proto/grph1_csgi.fp +27 -0
- data/GenWrapper/proto/grph1_scpack.fp +442 -0
- data/GenWrapper/proto/grph1_sgpack.fp +1406 -0
- data/GenWrapper/proto/grph1_slpack.fp +316 -0
- data/GenWrapper/proto/grph1_stpack.fp +411 -0
- data/GenWrapper/proto/grph1_swpack.fp +746 -0
- data/GenWrapper/proto/grph2_00 +16 -0
- data/GenWrapper/proto/grph2_grpack.fp +139 -0
- data/GenWrapper/proto/grph2_ucpack.fp +522 -0
- data/GenWrapper/proto/grph2_udpack.fp +692 -0
- data/GenWrapper/proto/grph2_uepack.fp +601 -0
- data/GenWrapper/proto/grph2_ugpack.fp +436 -0
- data/GenWrapper/proto/grph2_uhpack.fp +246 -0
- data/GenWrapper/proto/grph2_ulpack.fp +516 -0
- data/GenWrapper/proto/grph2_umpack.fp +676 -0
- data/GenWrapper/proto/grph2_uspack.fp +989 -0
- data/GenWrapper/proto/grph2_uupack.fp +663 -0
- data/GenWrapper/proto/grph2_uvpack.fp +246 -0
- data/GenWrapper/proto/grph2_uwpack.fp +203 -0
- data/GenWrapper/proto/grph2_uxpack.fp +295 -0
- data/GenWrapper/proto/grph2_uypack.fp +295 -0
- data/GenWrapper/proto/grph2_uzpack.fp +574 -0
- data/GenWrapper/proto/math1_00 +12 -0
- data/GenWrapper/proto/math1_blklib.fp +58 -0
- data/GenWrapper/proto/math1_chrlib.fp +83 -0
- data/GenWrapper/proto/math1_fnclib.fp +79 -0
- data/GenWrapper/proto/math1_gnmlib.fp +118 -0
- data/GenWrapper/proto/math1_gt2dlib.fp +144 -0
- data/GenWrapper/proto/math1_ifalib.fp +123 -0
- data/GenWrapper/proto/math1_indxlib.fp +222 -0
- data/GenWrapper/proto/math1_intlib.fp +46 -0
- data/GenWrapper/proto/math1_lrllib.fp +276 -0
- data/GenWrapper/proto/math1_maplib.fp +24 -0
- data/GenWrapper/proto/math1_oslib.fp +28 -0
- data/GenWrapper/proto/math1_rfalib.fp +420 -0
- data/GenWrapper/proto/math1_rfblib.fp +51 -0
- data/GenWrapper/proto/math1_sublib.fp +60 -0
- data/GenWrapper/proto/math1_syslib.fp +767 -0
- data/GenWrapper/proto/math1_vialib.fp +339 -0
- data/GenWrapper/proto/math1_viblib.fp +264 -0
- data/GenWrapper/proto/math1_vralib.fp +339 -0
- data/GenWrapper/proto/math1_vrblib.fp +264 -0
- data/GenWrapper/proto/math1_xfclib.fp +68 -0
- data/GenWrapper/proto/math2_00 +12 -0
- data/GenWrapper/proto/math2_fftlib.fp +803 -0
- data/GenWrapper/proto/math2_intrlib.fp +32 -0
- data/GenWrapper/proto/math2_odelib.fp_notused +568 -0
- data/GenWrapper/proto/math2_rnmlib.fp +54 -0
- data/GenWrapper/proto/math2_shtlib.fp +1292 -0
- data/GenWrapper/proto/math2_shtlib.fp_old +1294 -0
- data/GenWrapper/proto/math2_vstlib.fp +84 -0
- data/GenWrapper/proto/misc1_00 +17 -0
- data/GenWrapper/proto/misc1_chnlib.fp +47 -0
- data/GenWrapper/proto/misc1_datelib.fp +352 -0
- data/GenWrapper/proto/misc1_fmtlib.fp +19 -0
- data/GenWrapper/proto/misc1_misclib.fp +26 -0
- data/GenWrapper/proto/misc1_randlib.fp +36 -0
- data/GenWrapper/proto/misc1_timelib.fp +151 -0
- data/GenWrapper/proto2c.rb +13 -0
- data/GenWrapper/prototype.rb +209 -0
- data/GenWrapper/util.rb +9 -0
- data/GenWrapper/variable.rb +671 -0
- data/README +29 -0
- data/Rakefile +44 -0
- data/ToDo +2 -0
- data/dcl_cary2obj.c +354 -0
- data/dcl_narrayed_funcs.c +518 -0
- data/dcl_obj2cary.c +512 -0
- data/dcl_rubydcloriginal.c +129 -0
- data/demo/gokuraku/hop/hop.rb +26 -0
- data/demo/gokuraku/jump/jump1.rb +64 -0
- data/demo/gokuraku/jump/jump2.rb +50 -0
- data/demo/gokuraku/layout/lay1.rb +46 -0
- data/demo/gokuraku/layout/lay2.rb +44 -0
- data/demo/gokuraku/step/step0.rb +34 -0
- data/demo/gokuraku/step/step1.rb +84 -0
- data/demo/gokuraku/step/step2.rb +62 -0
- data/demo/gokuraku/u1d/u1d1.rb +37 -0
- data/demo/gokuraku/u1d/u1d2.rb +50 -0
- data/demo/gokuraku/u2d/u2d1.rb +46 -0
- data/demo/gokuraku/u2d/u2d2.rb +49 -0
- data/demo/gokuraku/u2d/u2d3.rb +45 -0
- data/demo/gokuraku/u2d/u2d4.rb +61 -0
- data/demo/grph1/scpack/scpkt2.rb +40 -0
- data/demo/grph1/scpack/scpkt3.rb +102 -0
- data/demo/grph1/scpack/scpkt6.rb +84 -0
- data/demo/grph1/scpack/scpkt7.rb +59 -0
- data/demo/grph1/scpack/scpkt8.rb +60 -0
- data/demo/grph1/scpack/t810630.dat +296 -0
- data/demo/grph1/scpack/t811231.dat +296 -0
- data/demo/grph1/sgpack/sgfont.rb +52 -0
- data/demo/grph1/sgpack/sgfonz.rb +53 -0
- data/demo/grph1/sgpack/sgksx1.rb +44 -0
- data/demo/grph1/sgpack/sgksx2.rb +90 -0
- data/demo/grph1/sgpack/sgksx3.rb +75 -0
- data/demo/grph1/sgpack/sglidx.rb +40 -0
- data/demo/grph1/sgpack/sgltyp.rb +48 -0
- data/demo/grph1/sgpack/sgpk01.rb +39 -0
- data/demo/grph1/sgpack/sgpk02.rb +79 -0
- data/demo/grph1/sgpack/sgpk03.rb +107 -0
- data/demo/grph1/sgpack/sgpk04.rb +86 -0
- data/demo/grph1/sgpack/sgpk05.rb +66 -0
- data/demo/grph1/sgpack/sgpk06.rb +91 -0
- data/demo/grph1/sgpack/sgpk07.rb +66 -0
- data/demo/grph1/sgpack/sgpk08.rb +67 -0
- data/demo/grph1/sgpack/sgpk09.rb +61 -0
- data/demo/grph1/sgpack/sgpk10.rb +58 -0
- data/demo/grph1/sgpack/sgtclr.rb +59 -0
- data/demo/grph1/sgpack/sgtone.rb +70 -0
- data/demo/grph1/sgpack/sgtonz.rb +66 -0
- data/demo/grph1/slpack/slpk01.rb +37 -0
- data/demo/grph1/slpack/slpk02.rb +29 -0
- data/demo/grph1/slpack/slpk03.rb +36 -0
- data/demo/grph1/slpack/slpk04.rb +36 -0
- data/demo/grph2/g2pack/g2pk01.rb +78 -0
- data/demo/grph2/g2pack/g2pk02.rb +75 -0
- data/demo/grph2/grpack/grpk01.rb +65 -0
- data/demo/grph2/ucpack/ucpk01.rb +95 -0
- data/demo/grph2/ucpack/ucpk02.rb +95 -0
- data/demo/grph2/udegpk/u2df01.rb +50 -0
- data/demo/grph2/udegpk/u2df02.rb +76 -0
- data/demo/grph2/udegpk/u2df03.rb +80 -0
- data/demo/grph2/udegpk/u2df04.rb +52 -0
- data/demo/grph2/udegpk/u2df05.rb +70 -0
- data/demo/grph2/udegpk/u2df06.rb +59 -0
- data/demo/grph2/udegpk/u2df07.rb +68 -0
- data/demo/grph2/udegpk/u2df08.rb +91 -0
- data/demo/grph2/udegpk/u2df09b.rb +88 -0
- data/demo/grph2/udegpk/u2df09c.rb +89 -0
- data/demo/grph2/udegpk/u2df09e.rb +88 -0
- data/demo/grph2/udegpk/u2df09f.rb +89 -0
- data/demo/grph2/ulpack/ulpk01.rb +111 -0
- data/demo/grph2/ulpack/ulpk01n.rb +111 -0
- data/demo/grph2/ulpack/ulpk02.rb +111 -0
- data/demo/grph2/ulpack/ulpk02n.rb +111 -0
- data/demo/grph2/umpack/t811231.dat +296 -0
- data/demo/grph2/umpack/test01.rb +69 -0
- data/demo/grph2/umpack/test02.rb +64 -0
- data/demo/grph2/umpack/test03.rb +84 -0
- data/demo/grph2/umpack/test04.rb +54 -0
- data/demo/grph2/umpack/test05.rb +85 -0
- data/demo/grph2/umpack/test06.rb +44 -0
- data/demo/grph2/umpack/test07.rb +43 -0
- data/demo/grph2/umpack/test08.rb +37 -0
- data/demo/grph2/umpack/test09.rb +86 -0
- data/demo/grph2/umpack/test10.rb +52 -0
- data/demo/grph2/umpack/umpk01.rb +65 -0
- data/demo/grph2/umpack/umpk02.rb +54 -0
- data/demo/grph2/umpack/umpk03.rb +46 -0
- data/demo/grph2/umpack/umpk04.rb +65 -0
- data/demo/grph2/umpack/umpk05.rb +84 -0
- data/demo/grph2/uspack/uspk01.rb +39 -0
- data/demo/grph2/uspack/uspk02.rb +47 -0
- data/demo/grph2/uspack/uspk03.rb +38 -0
- data/demo/grph2/uspack/uspk04.rb +64 -0
- data/demo/grph2/uspack/uspk05.rb +68 -0
- data/demo/grph2/uspack/uspk06.rb +43 -0
- data/demo/grph2/uspack/uspk07.rb +62 -0
- data/demo/grph2/uspack/uspk08.rb +46 -0
- data/demo/grph2/uspack/uspk09.rb +89 -0
- data/demo/grph2/uspack/uspk10.rb +48 -0
- data/demo/grph2/uspack/uspk11.rb +71 -0
- data/demo/grph2/uspack/uspk12.rb +51 -0
- data/demo/grph2/uupack/uupk01.rb +50 -0
- data/demo/grph2/uupack/uupk02.rb +74 -0
- data/demo/grph2/uupack/uupk03.rb +56 -0
- data/demo/grph2/uupack/uupk04.rb +84 -0
- data/demo/grph2/uupack/uupk05.rb +74 -0
- data/demo/grph2/uupack/uupk06.rb +77 -0
- data/demo/grph2/uupack/uupk07.rb +88 -0
- data/demo/grph2/uxyzpk/uxyz01.rb +38 -0
- data/demo/grph2/uxyzpk/uxyz02.rb +36 -0
- data/demo/grph2/uxyzpk/uxyz03.rb +37 -0
- data/demo/grph2/uxyzpk/uxyz04.rb +54 -0
- data/demo/grph2/uxyzpk/uxyz05.rb +41 -0
- data/demo/grph2/uxyzpk/uxyz06.rb +42 -0
- data/demo/grph2/uxyzpk/uxyz07.rb +38 -0
- data/demo/grph2/uxyzpk/uxyz08.rb +60 -0
- data/demo/grph2/uxyzpk/uxyz09.rb +52 -0
- data/demo/grph2/uxyzpk/uxyz10.rb +57 -0
- data/demo/grph2/ximage/MEMO +23 -0
- data/demo/grph2/ximage/tomsclm.dat +1344 -0
- data/demo/grph2/ximage/ximg01.rb +145 -0
- data/demo/math1/gt2dlib/gt2d01.rb +33 -0
- data/demo/math2/fftlib/fftl01.rb +28 -0
- data/demo/math2/fftlib/fftl02.rb +224 -0
- data/demo/math2/intrlib/intr01.rb +22 -0
- data/demo/math2/rnmlib/rnml01.rb +21 -0
- data/demo/math2/vstlib/vstl01.rb +29 -0
- data/demo/math2/vstlib/vstl02.rb +38 -0
- data/demo/rakuraku/color/color1.rb +61 -0
- data/demo/rakuraku/color/color2.rb +57 -0
- data/demo/rakuraku/color/color3.rb +81 -0
- data/demo/rakuraku/kihon/kihon1.rb +39 -0
- data/demo/rakuraku/kihon/kihon2.rb +56 -0
- data/demo/rakuraku/kihon/kihon3.rb +57 -0
- data/demo/rakuraku/kihon/kihon4.rb +70 -0
- data/demo/rakuraku/kihon/kihon5.rb +56 -0
- data/demo/rakuraku/kihon/kihon6.rb +40 -0
- data/demo/rakuraku/kihon/kihon7.rb +98 -0
- data/demo/rakuraku/kihon/kihon8.rb +93 -0
- data/demo/rakuraku/kihon/kihon9.rb +77 -0
- data/demo/rakuraku/kihon/kihona.rb +123 -0
- data/demo/rakuraku/kihon/kihonb.rb +73 -0
- data/demo/rakuraku/kihon/kihonc.rb +110 -0
- data/demo/rakuraku/layout/lay1.rb +29 -0
- data/demo/rakuraku/layout/lay2.rb +28 -0
- data/demo/rakuraku/layout/lay3.rb +29 -0
- data/demo/rakuraku/map3d/map3d1.rb +68 -0
- data/demo/rakuraku/map3d/map3d2.rb +67 -0
- data/demo/rakuraku/map3d/map3d3.rb +66 -0
- data/demo/rakuraku/map3d/map3d4.rb +107 -0
- data/demo/rakuraku/map3d/map3d5.rb +118 -0
- data/demo/rakuraku/map3d/map3d6.rb +96 -0
- data/demo/rakuraku/map3d/map3d7.rb +98 -0
- data/demo/rakuraku/miss/miss1.rb +76 -0
- data/demo/rakuraku/miss/miss2.rb +60 -0
- data/demo/rakuraku/quick/quick1.rb +35 -0
- data/demo/rakuraku/quick/quick2.rb +33 -0
- data/demo/rakuraku/quick/quick3.rb +46 -0
- data/demo/rakuraku/quick/quick4.rb +47 -0
- data/demo/rakuraku/quick/quick5.rb +46 -0
- data/demo/rakuraku/u2d/u2d1.rb +50 -0
- data/demo/rakuraku/u2d/u2d2.rb +51 -0
- data/demo/rakuraku/u2d/u2d3.rb +70 -0
- data/demo/rakuraku/u2d/u2d4.rb +59 -0
- data/demo/rakuraku/u2d/u2d5.rb +61 -0
- data/demo/rakuraku/u2d/u2d6.rb +57 -0
- data/demo/rakuraku/u2d/u2d7.rb +103 -0
- data/demo/rakuraku/uspack/uspac1.rb +35 -0
- data/demo/rakuraku/uspack/uspac2.rb +66 -0
- data/demo/rakuraku/uspack/uspac3.rb +37 -0
- data/demo/rakuraku/uspack/uspac4.rb +49 -0
- data/demo/rakuraku/uxyz/uxyz1.rb +37 -0
- data/demo/rakuraku/uxyz/uxyz2.rb +35 -0
- data/demo/rakuraku/uxyz/uxyz3.rb +38 -0
- data/demo/rakuraku/uxyz/uxyz4.rb +36 -0
- data/demo/rakuraku/uxyz/uxyz5.rb +57 -0
- data/demo/rakuraku/uxyz/uxyz6.rb +43 -0
- data/demo/rakuraku/uxyz/uxyz7.rb +39 -0
- data/demo/rakuraku/uxyz/uxyz8.rb +56 -0
- data/demo/rubydcloriginal/uemrkz1.rb +37 -0
- data/dummy.c +9 -0
- data/extconf.rb +133 -0
- data/grph1_csgi.c +119 -0
- data/grph1_scpack.c +1388 -0
- data/grph1_sgpack.c +3876 -0
- data/grph1_slpack.c +432 -0
- data/grph1_stpack.c +1044 -0
- data/grph1_swpack.c +1922 -0
- data/grph1_zgpack.c.org +141 -0
- data/grph2_grpack.c +368 -0
- data/grph2_ucpack.c +1232 -0
- data/grph2_udpack.c +1404 -0
- data/grph2_uepack.c +1443 -0
- data/grph2_ugpack.c +1073 -0
- data/grph2_uhpack.c +943 -0
- data/grph2_ulpack.c +1304 -0
- data/grph2_umpack.c +1540 -0
- data/grph2_uspack.c +2481 -0
- data/grph2_uupack.c +1689 -0
- data/grph2_uvpack.c +943 -0
- data/grph2_uwpack.c +552 -0
- data/grph2_uxpack.c +555 -0
- data/grph2_uypack.c +555 -0
- data/grph2_uzpack.c +1455 -0
- data/init.c.default +187 -0
- data/init.c.gtk +189 -0
- data/lib/dcl.rb +5084 -0
- data/math1_blklib.c +227 -0
- data/math1_chrlib.c +239 -0
- data/math1_fnclib.c +254 -0
- data/math1_gnmlib.c +356 -0
- data/math1_gt2dlib.c +514 -0
- data/math1_ifalib.c +427 -0
- data/math1_indxlib.c +709 -0
- data/math1_intlib.c +167 -0
- data/math1_lrllib.c +817 -0
- data/math1_maplib.c +172 -0
- data/math1_oslib.c +111 -0
- data/math1_rfalib.c +1417 -0
- data/math1_rfblib.c +232 -0
- data/math1_sublib.c +221 -0
- data/math1_syslib.c +2025 -0
- data/math1_vialib.c +1156 -0
- data/math1_viblib.c +1027 -0
- data/math1_vralib.c +1156 -0
- data/math1_vrblib.c +1027 -0
- data/math1_xfclib.c +142 -0
- data/math2_fftlib.c +833 -0
- data/math2_intrlib.c +138 -0
- data/math2_rnmlib.c +226 -0
- data/math2_shtlib.c +1568 -0
- data/math2_vstlib.c +366 -0
- data/misc1_chnlib.c +179 -0
- data/misc1_datelib.c +953 -0
- data/misc1_fmtlib.c +99 -0
- data/misc1_misclib.c +118 -0
- data/misc1_randlib.c +142 -0
- data/misc1_timelib.c +380 -0
- data/obsolete/README +2 -0
- data/obsolete/dcl_ext.html +112 -0
- data/obsolete/dcl_ext.rb +313 -0
- metadata +433 -0
@@ -0,0 +1,68 @@
|
|
1
|
+
/*
|
2
|
+
* math1_xfclib.fp
|
3
|
+
*
|
4
|
+
* $Id: math1_xfclib.fp,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
|
5
|
+
*/
|
6
|
+
|
7
|
+
/*****************************************************
|
8
|
+
* math1/xfclib/jfromc.f
|
9
|
+
****************************************************/
|
10
|
+
/*
|
11
|
+
* INTEGER FUNCTION JFROMC(CX)
|
12
|
+
* CHARACTER CX*(*)
|
13
|
+
*/
|
14
|
+
/*
|
15
|
+
INTEGER FUNCTION JFROMC(
|
16
|
+
CHARACTER*(*) CX
|
17
|
+
)
|
18
|
+
*/
|
19
|
+
|
20
|
+
/*****************************************************
|
21
|
+
* math1/xfclib/ifromc.f
|
22
|
+
****************************************************/
|
23
|
+
/* NB: (IFROMC) CX size unknown. */
|
24
|
+
/*
|
25
|
+
* INTEGER FUNCTION IFROMC(CX)
|
26
|
+
* CHARACTER CX*(*)
|
27
|
+
*/
|
28
|
+
INTEGER FUNCTION IFROMC(
|
29
|
+
CHARACTER*(*) CX
|
30
|
+
)
|
31
|
+
|
32
|
+
/*****************************************************
|
33
|
+
* math1/xfclib/lfromc.f
|
34
|
+
****************************************************/
|
35
|
+
/* NB: (LFROMC) CX size unknown. */
|
36
|
+
/*
|
37
|
+
* LOGICAL FUNCTION LFROMC(CX)
|
38
|
+
* CHARACTER CX*(*)
|
39
|
+
*/
|
40
|
+
LOGICAL FUNCTION LFROMC(
|
41
|
+
CHARACTER*(*) CX
|
42
|
+
)
|
43
|
+
|
44
|
+
/*****************************************************
|
45
|
+
* math1/xfclib/ffromc.f
|
46
|
+
****************************************************/
|
47
|
+
/*
|
48
|
+
* REAL FUNCTION FFROMC(CX)
|
49
|
+
* CHARACTER CX*(*)
|
50
|
+
*/
|
51
|
+
/*
|
52
|
+
REAL FUNCTION FFROMC(
|
53
|
+
CHARACTER*(*) CX
|
54
|
+
)
|
55
|
+
*/
|
56
|
+
|
57
|
+
/*****************************************************
|
58
|
+
* math1/xfclib/rfromc.f
|
59
|
+
****************************************************/
|
60
|
+
/* NB: (RFROMC) CX size unknown. */
|
61
|
+
/*
|
62
|
+
* REAL FUNCTION RFROMC(CX)
|
63
|
+
* CHARACTER CX*(*)
|
64
|
+
*/
|
65
|
+
REAL FUNCTION RFROMC(
|
66
|
+
CHARACTER*(*) CX
|
67
|
+
)
|
68
|
+
|
@@ -0,0 +1,803 @@
|
|
1
|
+
/*
|
2
|
+
* math2_fftlib.fp
|
3
|
+
*
|
4
|
+
* $Id: math2_fftlib.fp,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
|
5
|
+
*/
|
6
|
+
|
7
|
+
/*****************************************************
|
8
|
+
* math2/fftlib/cfftb.f
|
9
|
+
****************************************************/
|
10
|
+
/*
|
11
|
+
* SUBROUTINE CFFTB (N,C,WSAVE)
|
12
|
+
* DIMENSION C(*) ,WSAVE(*)
|
13
|
+
*/
|
14
|
+
SUBROUTINE CFFTB(
|
15
|
+
INTEGER N
|
16
|
+
REAL (2*N) io C
|
17
|
+
REAL (4*N+15) WSAVE
|
18
|
+
)
|
19
|
+
|
20
|
+
/*****************************************************
|
21
|
+
* math2/fftlib/cfftb1.f
|
22
|
+
****************************************************/
|
23
|
+
/*
|
24
|
+
* SUBROUTINE CFFTB1 (N,C,CH,WA,IFAC)
|
25
|
+
* DIMENSION CH(*) ,C(*) ,WA(*) ,IFAC(*)
|
26
|
+
*/
|
27
|
+
/*
|
28
|
+
SUBROUTINE CFFTB1(
|
29
|
+
INTEGER N
|
30
|
+
REAL (*) io C
|
31
|
+
REAL (*) t CH
|
32
|
+
REAL (*) t WA
|
33
|
+
INTEGER (*) t IFAC
|
34
|
+
)
|
35
|
+
*/
|
36
|
+
|
37
|
+
/*****************************************************
|
38
|
+
* math2/fftlib/cfftf.f
|
39
|
+
****************************************************/
|
40
|
+
/*
|
41
|
+
* SUBROUTINE CFFTF (N,C,WSAVE)
|
42
|
+
* DIMENSION C(*) ,WSAVE(*)
|
43
|
+
*/
|
44
|
+
SUBROUTINE CFFTF(
|
45
|
+
INTEGER N
|
46
|
+
REAL (2*N) io C
|
47
|
+
REAL (4*N+15) WSAVE
|
48
|
+
)
|
49
|
+
|
50
|
+
/*****************************************************
|
51
|
+
* math2/fftlib/cfftf1.f
|
52
|
+
****************************************************/
|
53
|
+
/*
|
54
|
+
* SUBROUTINE CFFTF1 (N,C,CH,WA,IFAC)
|
55
|
+
* DIMENSION CH(*) ,C(*) ,WA(*) ,IFAC(*)
|
56
|
+
*/
|
57
|
+
/*
|
58
|
+
SUBROUTINE CFFTF1(
|
59
|
+
INTEGER N
|
60
|
+
REAL (*) io C
|
61
|
+
REAL (*) t CH
|
62
|
+
REAL (*) t WA
|
63
|
+
INTEGER (*) t IFAC
|
64
|
+
)
|
65
|
+
*/
|
66
|
+
|
67
|
+
/*****************************************************
|
68
|
+
* math2/fftlib/cffti.f
|
69
|
+
****************************************************/
|
70
|
+
/*
|
71
|
+
* SUBROUTINE CFFTI (N,WSAVE)
|
72
|
+
* DIMENSION WSAVE(*)
|
73
|
+
*/
|
74
|
+
SUBROUTINE CFFTI(
|
75
|
+
INTEGER N
|
76
|
+
REAL (4*N+15) o WSAVE
|
77
|
+
)
|
78
|
+
|
79
|
+
/*****************************************************
|
80
|
+
* math2/fftlib/cffti1.f
|
81
|
+
****************************************************/
|
82
|
+
/*
|
83
|
+
* SUBROUTINE CFFTI1 (N,WA,IFAC)
|
84
|
+
* DIMENSION WA(*) ,IFAC(*) ,NTRYH(4)
|
85
|
+
*/
|
86
|
+
/*
|
87
|
+
SUBROUTINE CFFTI1(
|
88
|
+
INTEGER N
|
89
|
+
REAL (*) o WA
|
90
|
+
INTEGER (*) o IFAC
|
91
|
+
)
|
92
|
+
*/
|
93
|
+
|
94
|
+
/*****************************************************
|
95
|
+
* math2/fftlib/cosqb.f
|
96
|
+
****************************************************/
|
97
|
+
/*
|
98
|
+
* SUBROUTINE COSQB (N,X,WSAVE)
|
99
|
+
* DIMENSION X(*) ,WSAVE(*)
|
100
|
+
*/
|
101
|
+
SUBROUTINE COSQB(
|
102
|
+
INTEGER N
|
103
|
+
REAL (N) io X
|
104
|
+
REAL (3*N+15) WSAVE
|
105
|
+
)
|
106
|
+
|
107
|
+
/*****************************************************
|
108
|
+
* math2/fftlib/cosqb1.f
|
109
|
+
****************************************************/
|
110
|
+
/*
|
111
|
+
* SUBROUTINE COSQB1 (N,X,W,XH)
|
112
|
+
* DIMENSION X(*) ,W(*) ,XH(*)
|
113
|
+
*/
|
114
|
+
/*
|
115
|
+
SUBROUTINE COSQB1(
|
116
|
+
INTEGER N
|
117
|
+
REAL io X
|
118
|
+
REAL (*) t W
|
119
|
+
REAL (*) t XH
|
120
|
+
)
|
121
|
+
*/
|
122
|
+
|
123
|
+
/*****************************************************
|
124
|
+
* math2/fftlib/cosqf.f
|
125
|
+
****************************************************/
|
126
|
+
/*
|
127
|
+
* SUBROUTINE COSQF (N,X,WSAVE)
|
128
|
+
* DIMENSION X(*) ,WSAVE(*)
|
129
|
+
*/
|
130
|
+
SUBROUTINE COSQF(
|
131
|
+
INTEGER N
|
132
|
+
REAL (N) io X
|
133
|
+
REAL (3*N+15) WSAVE
|
134
|
+
)
|
135
|
+
|
136
|
+
/*****************************************************
|
137
|
+
* math2/fftlib/cosqf1.f
|
138
|
+
****************************************************/
|
139
|
+
/*
|
140
|
+
* SUBROUTINE COSQF1 (N,X,W,XH)
|
141
|
+
* DIMENSION X(*) ,W(*) ,XH(*)
|
142
|
+
*/
|
143
|
+
/*
|
144
|
+
SUBROUTINE COSQF1(
|
145
|
+
INTEGER N
|
146
|
+
REAL io X
|
147
|
+
REAL (*) t W
|
148
|
+
REAL (*) t XH
|
149
|
+
)
|
150
|
+
*/
|
151
|
+
|
152
|
+
/*****************************************************
|
153
|
+
* math2/fftlib/cosqi.f
|
154
|
+
****************************************************/
|
155
|
+
/*
|
156
|
+
* SUBROUTINE COSQI (N,WSAVE)
|
157
|
+
* DIMENSION WSAVE(*)
|
158
|
+
*/
|
159
|
+
SUBROUTINE COSQI(
|
160
|
+
INTEGER N
|
161
|
+
REAL (3*N+15) o WSAVE
|
162
|
+
)
|
163
|
+
|
164
|
+
/*****************************************************
|
165
|
+
* math2/fftlib/cost.f
|
166
|
+
****************************************************/
|
167
|
+
/*
|
168
|
+
* SUBROUTINE COST (N,X,WSAVE)
|
169
|
+
* DIMENSION X(*) ,WSAVE(*)
|
170
|
+
*/
|
171
|
+
SUBROUTINE COST(
|
172
|
+
INTEGER N
|
173
|
+
REAL (N) io X
|
174
|
+
REAL (3*N+15) WSAVE
|
175
|
+
)
|
176
|
+
|
177
|
+
/*****************************************************
|
178
|
+
* math2/fftlib/costi.f
|
179
|
+
****************************************************/
|
180
|
+
/*
|
181
|
+
* SUBROUTINE COSTI (N,WSAVE)
|
182
|
+
* DIMENSION WSAVE(*)
|
183
|
+
*/
|
184
|
+
SUBROUTINE COSTI(
|
185
|
+
INTEGER N
|
186
|
+
REAL (3*N+15) o WSAVE
|
187
|
+
)
|
188
|
+
|
189
|
+
/*****************************************************
|
190
|
+
* math2/fftlib/ezfft1.f
|
191
|
+
****************************************************/
|
192
|
+
/*
|
193
|
+
* SUBROUTINE EZFFT1 (N,WA,IFAC)
|
194
|
+
* DIMENSION WA(*) ,IFAC(*) ,NTRYH(4)
|
195
|
+
*/
|
196
|
+
/*
|
197
|
+
SUBROUTINE EZFFT1(
|
198
|
+
INTEGER N
|
199
|
+
REAL (*) o WA
|
200
|
+
INTEGER (*) o IFAC
|
201
|
+
)
|
202
|
+
*/
|
203
|
+
|
204
|
+
/*****************************************************
|
205
|
+
* math2/fftlib/ezfftb.f
|
206
|
+
****************************************************/
|
207
|
+
/*
|
208
|
+
* SUBROUTINE EZFFTB (N,R,AZERO,A,B,WSAVE)
|
209
|
+
* DIMENSION R(*) ,A(*) ,B(*) ,WSAVE(*)
|
210
|
+
*/
|
211
|
+
SUBROUTINE EZFFTB(
|
212
|
+
INTEGER N
|
213
|
+
REAL (N) o R
|
214
|
+
REAL AZERO
|
215
|
+
/* REAL A */
|
216
|
+
REAL (N/2) A
|
217
|
+
REAL (N/2) B
|
218
|
+
REAL (3*N+15) WSAVE
|
219
|
+
)
|
220
|
+
|
221
|
+
/*****************************************************
|
222
|
+
* math2/fftlib/ezfftf.f
|
223
|
+
****************************************************/
|
224
|
+
/*
|
225
|
+
* SUBROUTINE EZFFTF (N,R,AZERO,A,B,WSAVE)
|
226
|
+
* DIMENSION R(*) ,A(*) ,B(*) ,WSAVE(*)
|
227
|
+
*/
|
228
|
+
SUBROUTINE EZFFTF(
|
229
|
+
INTEGER N
|
230
|
+
REAL (N) R
|
231
|
+
REAL o AZERO
|
232
|
+
/* REAL o A */
|
233
|
+
REAL (N/2) o A
|
234
|
+
REAL (N/2) o B
|
235
|
+
REAL (3*N+15) WSAVE
|
236
|
+
)
|
237
|
+
|
238
|
+
/*****************************************************
|
239
|
+
* math2/fftlib/ezffti.f
|
240
|
+
****************************************************/
|
241
|
+
/*
|
242
|
+
* SUBROUTINE EZFFTI (N,WSAVE)
|
243
|
+
* DIMENSION WSAVE(*)
|
244
|
+
*/
|
245
|
+
SUBROUTINE EZFFTI(
|
246
|
+
INTEGER N
|
247
|
+
REAL (3*N+15) o WSAVE
|
248
|
+
)
|
249
|
+
|
250
|
+
/*****************************************************
|
251
|
+
* math2/fftlib/passb.f
|
252
|
+
****************************************************/
|
253
|
+
/*
|
254
|
+
* SUBROUTINE PASSB (NAC,IDO,IP,L1,IDL1,CC,C1,C2,CH,CH2,WA)
|
255
|
+
* DIMENSION CH(IDO,L1,IP) ,CC(IDO,IP,L1) ,
|
256
|
+
*/
|
257
|
+
/*
|
258
|
+
SUBROUTINE PASSB(
|
259
|
+
INTEGER o NAC
|
260
|
+
INTEGER IDO
|
261
|
+
INTEGER IP
|
262
|
+
INTEGER L1
|
263
|
+
INTEGER IDL1
|
264
|
+
REAL (IDO,IP,L1) CC
|
265
|
+
REAL io C1
|
266
|
+
REAL io C2
|
267
|
+
REAL (IDO,L1,IP) io CH
|
268
|
+
REAL io CH2
|
269
|
+
REAL WA
|
270
|
+
)
|
271
|
+
*/
|
272
|
+
|
273
|
+
/*****************************************************
|
274
|
+
* math2/fftlib/passb2.f
|
275
|
+
****************************************************/
|
276
|
+
/*
|
277
|
+
* SUBROUTINE PASSB2 (IDO,L1,CC,CH,WA1)
|
278
|
+
* DIMENSION CC(IDO,2,L1) ,CH(IDO,L1,2) ,
|
279
|
+
*/
|
280
|
+
/*
|
281
|
+
SUBROUTINE PASSB2(
|
282
|
+
INTEGER IDO
|
283
|
+
INTEGER L1
|
284
|
+
REAL (IDO,2,L1) CC
|
285
|
+
REAL (IDO,L1,2) o CH
|
286
|
+
REAL WA1
|
287
|
+
)
|
288
|
+
*/
|
289
|
+
|
290
|
+
/*****************************************************
|
291
|
+
* math2/fftlib/passb3.f
|
292
|
+
****************************************************/
|
293
|
+
/*
|
294
|
+
* SUBROUTINE PASSB3 (IDO,L1,CC,CH,WA1,WA2)
|
295
|
+
* DIMENSION CC(IDO,3,L1) ,CH(IDO,L1,3) ,
|
296
|
+
*/
|
297
|
+
/*
|
298
|
+
SUBROUTINE PASSB3(
|
299
|
+
INTEGER IDO
|
300
|
+
INTEGER L1
|
301
|
+
REAL (IDO,3,L1) CC
|
302
|
+
REAL (IDO,L1,3) o CH
|
303
|
+
REAL WA1
|
304
|
+
REAL WA2
|
305
|
+
)
|
306
|
+
*/
|
307
|
+
|
308
|
+
/*****************************************************
|
309
|
+
* math2/fftlib/passb4.f
|
310
|
+
****************************************************/
|
311
|
+
/*
|
312
|
+
* SUBROUTINE PASSB4 (IDO,L1,CC,CH,WA1,WA2,WA3)
|
313
|
+
* DIMENSION CC(IDO,4,L1) ,CH(IDO,L1,4) ,
|
314
|
+
*/
|
315
|
+
/*
|
316
|
+
SUBROUTINE PASSB4(
|
317
|
+
INTEGER IDO
|
318
|
+
INTEGER L1
|
319
|
+
REAL (IDO,4,L1) CC
|
320
|
+
REAL (IDO,L1,4) o CH
|
321
|
+
REAL WA1
|
322
|
+
REAL WA2
|
323
|
+
REAL WA3
|
324
|
+
)
|
325
|
+
*/
|
326
|
+
|
327
|
+
/*****************************************************
|
328
|
+
* math2/fftlib/passb5.f
|
329
|
+
****************************************************/
|
330
|
+
/*
|
331
|
+
* SUBROUTINE PASSB5 (IDO,L1,CC,CH,WA1,WA2,WA3,WA4)
|
332
|
+
* DIMENSION CC(IDO,5,L1) ,CH(IDO,L1,5) ,
|
333
|
+
*/
|
334
|
+
/*
|
335
|
+
SUBROUTINE PASSB5(
|
336
|
+
INTEGER IDO
|
337
|
+
INTEGER L1
|
338
|
+
REAL (IDO,5,L1) CC
|
339
|
+
REAL (IDO,L1,5) o CH
|
340
|
+
REAL WA1
|
341
|
+
REAL WA2
|
342
|
+
REAL WA3
|
343
|
+
REAL WA4
|
344
|
+
)
|
345
|
+
*/
|
346
|
+
|
347
|
+
/*****************************************************
|
348
|
+
* math2/fftlib/passf.f
|
349
|
+
****************************************************/
|
350
|
+
/*
|
351
|
+
* SUBROUTINE PASSF (NAC,IDO,IP,L1,IDL1,CC,C1,C2,CH,CH2,WA)
|
352
|
+
* DIMENSION CH(IDO,L1,IP) ,CC(IDO,IP,L1) ,
|
353
|
+
*/
|
354
|
+
/*
|
355
|
+
SUBROUTINE PASSF(
|
356
|
+
INTEGER o NAC
|
357
|
+
INTEGER IDO
|
358
|
+
INTEGER IP
|
359
|
+
INTEGER L1
|
360
|
+
INTEGER IDL1
|
361
|
+
REAL (IDO,IP,L1) CC
|
362
|
+
REAL io C1
|
363
|
+
REAL io C2
|
364
|
+
REAL (IDO,L1,IP) io CH
|
365
|
+
REAL io CH2
|
366
|
+
REAL WA
|
367
|
+
)
|
368
|
+
*/
|
369
|
+
|
370
|
+
/*****************************************************
|
371
|
+
* math2/fftlib/passf2.f
|
372
|
+
****************************************************/
|
373
|
+
/*
|
374
|
+
* SUBROUTINE PASSF2 (IDO,L1,CC,CH,WA1)
|
375
|
+
* DIMENSION CC(IDO,2,L1) ,CH(IDO,L1,2) ,
|
376
|
+
*/
|
377
|
+
/*
|
378
|
+
SUBROUTINE PASSF2(
|
379
|
+
INTEGER IDO
|
380
|
+
INTEGER L1
|
381
|
+
REAL (IDO,2,L1) CC
|
382
|
+
REAL (IDO,L1,2) o CH
|
383
|
+
REAL WA1
|
384
|
+
)
|
385
|
+
*/
|
386
|
+
|
387
|
+
/*****************************************************
|
388
|
+
* math2/fftlib/radb2.f
|
389
|
+
****************************************************/
|
390
|
+
/*
|
391
|
+
* SUBROUTINE RADB2 (IDO,L1,CC,CH,WA1)
|
392
|
+
* DIMENSION CC(IDO,2,L1) ,CH(IDO,L1,2) ,
|
393
|
+
*/
|
394
|
+
/*
|
395
|
+
SUBROUTINE RADB2(
|
396
|
+
INTEGER IDO
|
397
|
+
INTEGER L1
|
398
|
+
REAL (IDO,2,L1) CC
|
399
|
+
REAL (IDO,L1,2) o CH
|
400
|
+
REAL WA1
|
401
|
+
)
|
402
|
+
*/
|
403
|
+
|
404
|
+
/*****************************************************
|
405
|
+
* math2/fftlib/passf3.f
|
406
|
+
****************************************************/
|
407
|
+
/*
|
408
|
+
* SUBROUTINE PASSF3 (IDO,L1,CC,CH,WA1,WA2)
|
409
|
+
* DIMENSION CC(IDO,3,L1) ,CH(IDO,L1,3) ,
|
410
|
+
*/
|
411
|
+
/*
|
412
|
+
SUBROUTINE PASSF3(
|
413
|
+
INTEGER IDO
|
414
|
+
INTEGER L1
|
415
|
+
REAL (IDO,3,L1) CC
|
416
|
+
REAL (IDO,L1,3) o CH
|
417
|
+
REAL WA1
|
418
|
+
REAL WA2
|
419
|
+
)
|
420
|
+
*/
|
421
|
+
|
422
|
+
/*****************************************************
|
423
|
+
* math2/fftlib/passf4.f
|
424
|
+
****************************************************/
|
425
|
+
/*
|
426
|
+
* SUBROUTINE PASSF4 (IDO,L1,CC,CH,WA1,WA2,WA3)
|
427
|
+
* DIMENSION CC(IDO,4,L1) ,CH(IDO,L1,4) ,
|
428
|
+
*/
|
429
|
+
/*
|
430
|
+
SUBROUTINE PASSF4(
|
431
|
+
INTEGER IDO
|
432
|
+
INTEGER L1
|
433
|
+
REAL (IDO,4,L1) CC
|
434
|
+
REAL (IDO,L1,4) o CH
|
435
|
+
REAL WA1
|
436
|
+
REAL WA2
|
437
|
+
REAL WA3
|
438
|
+
)
|
439
|
+
*/
|
440
|
+
|
441
|
+
/*****************************************************
|
442
|
+
* math2/fftlib/passf5.f
|
443
|
+
****************************************************/
|
444
|
+
/*
|
445
|
+
* SUBROUTINE PASSF5 (IDO,L1,CC,CH,WA1,WA2,WA3,WA4)
|
446
|
+
* DIMENSION CC(IDO,5,L1) ,CH(IDO,L1,5) ,
|
447
|
+
*/
|
448
|
+
/*
|
449
|
+
SUBROUTINE PASSF5(
|
450
|
+
INTEGER IDO
|
451
|
+
INTEGER L1
|
452
|
+
REAL (IDO,5,L1) CC
|
453
|
+
REAL (IDO,L1,5) o CH
|
454
|
+
REAL WA1
|
455
|
+
REAL WA2
|
456
|
+
REAL WA3
|
457
|
+
REAL WA4
|
458
|
+
)
|
459
|
+
*/
|
460
|
+
|
461
|
+
/*****************************************************
|
462
|
+
* math2/fftlib/radb3.f
|
463
|
+
****************************************************/
|
464
|
+
/*
|
465
|
+
* SUBROUTINE RADB3 (IDO,L1,CC,CH,WA1,WA2)
|
466
|
+
* DIMENSION CC(IDO,3,L1) ,CH(IDO,L1,3) ,
|
467
|
+
*/
|
468
|
+
/*
|
469
|
+
SUBROUTINE RADB3(
|
470
|
+
INTEGER IDO
|
471
|
+
INTEGER L1
|
472
|
+
REAL (IDO,3,L1) CC
|
473
|
+
REAL (IDO,L1,3) o CH
|
474
|
+
REAL WA1
|
475
|
+
REAL WA2
|
476
|
+
)
|
477
|
+
*/
|
478
|
+
|
479
|
+
/*****************************************************
|
480
|
+
* math2/fftlib/radb4.f
|
481
|
+
****************************************************/
|
482
|
+
/*
|
483
|
+
* SUBROUTINE RADB4 (IDO,L1,CC,CH,WA1,WA2,WA3)
|
484
|
+
* DIMENSION CC(IDO,4,L1) ,CH(IDO,L1,4) ,
|
485
|
+
*/
|
486
|
+
/*
|
487
|
+
SUBROUTINE RADB4(
|
488
|
+
INTEGER IDO
|
489
|
+
INTEGER L1
|
490
|
+
REAL (IDO,4,L1) CC
|
491
|
+
REAL (IDO,L1,4) o CH
|
492
|
+
REAL WA1
|
493
|
+
REAL WA2
|
494
|
+
REAL WA3
|
495
|
+
)
|
496
|
+
*/
|
497
|
+
|
498
|
+
/*****************************************************
|
499
|
+
* math2/fftlib/radb5.f
|
500
|
+
****************************************************/
|
501
|
+
/*
|
502
|
+
* SUBROUTINE RADB5 (IDO,L1,CC,CH,WA1,WA2,WA3,WA4)
|
503
|
+
* DIMENSION CC(IDO,5,L1) ,CH(IDO,L1,5) ,
|
504
|
+
*/
|
505
|
+
/*
|
506
|
+
SUBROUTINE RADB5(
|
507
|
+
INTEGER IDO
|
508
|
+
INTEGER L1
|
509
|
+
REAL (IDO,5,L1) CC
|
510
|
+
REAL (IDO,L1,5) o CH
|
511
|
+
REAL WA1
|
512
|
+
REAL WA2
|
513
|
+
REAL WA3
|
514
|
+
REAL WA4
|
515
|
+
)
|
516
|
+
*/
|
517
|
+
|
518
|
+
/*****************************************************
|
519
|
+
* math2/fftlib/radbg.f
|
520
|
+
****************************************************/
|
521
|
+
/*
|
522
|
+
* SUBROUTINE RADBG (IDO,IP,L1,IDL1,CC,C1,C2,CH,CH2,WA)
|
523
|
+
* DIMENSION CH(IDO,L1,IP) ,CC(IDO,IP,L1) ,
|
524
|
+
*/
|
525
|
+
/*
|
526
|
+
SUBROUTINE RADBG(
|
527
|
+
INTEGER IDO
|
528
|
+
INTEGER IP
|
529
|
+
INTEGER L1
|
530
|
+
INTEGER IDL1
|
531
|
+
REAL (IDO,IP,L1) CC
|
532
|
+
REAL io C1
|
533
|
+
REAL io C2
|
534
|
+
REAL (IDO,L1,IP) io CH
|
535
|
+
REAL io CH2
|
536
|
+
REAL WA
|
537
|
+
)
|
538
|
+
*/
|
539
|
+
|
540
|
+
/*****************************************************
|
541
|
+
* math2/fftlib/radf2.f
|
542
|
+
****************************************************/
|
543
|
+
/*
|
544
|
+
* SUBROUTINE RADF2 (IDO,L1,CC,CH,WA1)
|
545
|
+
* DIMENSION CH(IDO,2,L1) ,CC(IDO,L1,2) ,
|
546
|
+
*/
|
547
|
+
/*
|
548
|
+
SUBROUTINE RADF2(
|
549
|
+
INTEGER IDO
|
550
|
+
INTEGER L1
|
551
|
+
REAL (IDO,L1,2) CC
|
552
|
+
REAL (IDO,2,L1) o CH
|
553
|
+
REAL WA1
|
554
|
+
)
|
555
|
+
*/
|
556
|
+
|
557
|
+
/*****************************************************
|
558
|
+
* math2/fftlib/radf3.f
|
559
|
+
****************************************************/
|
560
|
+
/*
|
561
|
+
* SUBROUTINE RADF3 (IDO,L1,CC,CH,WA1,WA2)
|
562
|
+
* DIMENSION CH(IDO,3,L1) ,CC(IDO,L1,3) ,
|
563
|
+
*/
|
564
|
+
/*
|
565
|
+
SUBROUTINE RADF3(
|
566
|
+
INTEGER IDO
|
567
|
+
INTEGER L1
|
568
|
+
REAL (IDO,L1,3) CC
|
569
|
+
REAL (IDO,3,L1) o CH
|
570
|
+
REAL WA1
|
571
|
+
REAL WA2
|
572
|
+
)
|
573
|
+
*/
|
574
|
+
|
575
|
+
/*****************************************************
|
576
|
+
* math2/fftlib/radf4.f
|
577
|
+
****************************************************/
|
578
|
+
/*
|
579
|
+
* SUBROUTINE RADF4 (IDO,L1,CC,CH,WA1,WA2,WA3)
|
580
|
+
* DIMENSION CC(IDO,L1,4) ,CH(IDO,4,L1) ,
|
581
|
+
*/
|
582
|
+
/*
|
583
|
+
SUBROUTINE RADF4(
|
584
|
+
INTEGER IDO
|
585
|
+
INTEGER L1
|
586
|
+
REAL (IDO,L1,4) CC
|
587
|
+
REAL (IDO,4,L1) o CH
|
588
|
+
REAL WA1
|
589
|
+
REAL WA2
|
590
|
+
REAL WA3
|
591
|
+
)
|
592
|
+
*/
|
593
|
+
|
594
|
+
/*****************************************************
|
595
|
+
* math2/fftlib/radf5.f
|
596
|
+
****************************************************/
|
597
|
+
/*
|
598
|
+
* SUBROUTINE RADF5 (IDO,L1,CC,CH,WA1,WA2,WA3,WA4)
|
599
|
+
* DIMENSION CC(IDO,L1,5) ,CH(IDO,5,L1) ,
|
600
|
+
*/
|
601
|
+
/*
|
602
|
+
SUBROUTINE RADF5(
|
603
|
+
INTEGER IDO
|
604
|
+
INTEGER L1
|
605
|
+
REAL (IDO,L1,5) CC
|
606
|
+
REAL (IDO,5,L1) o CH
|
607
|
+
REAL WA1
|
608
|
+
REAL WA2
|
609
|
+
REAL WA3
|
610
|
+
REAL WA4
|
611
|
+
)
|
612
|
+
*/
|
613
|
+
|
614
|
+
/*****************************************************
|
615
|
+
* math2/fftlib/radfg.f
|
616
|
+
****************************************************/
|
617
|
+
/*
|
618
|
+
* SUBROUTINE RADFG (IDO,IP,L1,IDL1,CC,C1,C2,CH,CH2,WA)
|
619
|
+
* DIMENSION CH(IDO,L1,IP) ,CC(IDO,IP,L1) ,
|
620
|
+
*/
|
621
|
+
/*
|
622
|
+
SUBROUTINE RADFG(
|
623
|
+
INTEGER IDO
|
624
|
+
INTEGER IP
|
625
|
+
INTEGER L1
|
626
|
+
INTEGER IDL1
|
627
|
+
REAL (IDO,IP,L1) CC
|
628
|
+
REAL io C1
|
629
|
+
REAL io C2
|
630
|
+
REAL (IDO,L1,IP) io CH
|
631
|
+
REAL io CH2
|
632
|
+
REAL WA
|
633
|
+
)
|
634
|
+
*/
|
635
|
+
|
636
|
+
/*****************************************************
|
637
|
+
* math2/fftlib/rfftb.f
|
638
|
+
****************************************************/
|
639
|
+
/*
|
640
|
+
* SUBROUTINE RFFTB (N,R,WSAVE)
|
641
|
+
* DIMENSION R(*) ,WSAVE(*)
|
642
|
+
*/
|
643
|
+
SUBROUTINE RFFTB(
|
644
|
+
INTEGER N
|
645
|
+
REAL (N) io R
|
646
|
+
REAL (2*N+15) WSAVE
|
647
|
+
)
|
648
|
+
|
649
|
+
/*****************************************************
|
650
|
+
* math2/fftlib/rfftb1.f
|
651
|
+
****************************************************/
|
652
|
+
/*
|
653
|
+
* SUBROUTINE RFFTB1 (N,C,CH,WA,IFAC)
|
654
|
+
* DIMENSION CH(*) ,C(*) ,WA(*) ,IFAC(*)
|
655
|
+
*/
|
656
|
+
/*
|
657
|
+
SUBROUTINE RFFTB1(
|
658
|
+
INTEGER N
|
659
|
+
REAL (*) io C
|
660
|
+
REAL (*) t CH
|
661
|
+
REAL (*) t WA
|
662
|
+
INTEGER (*) t IFAC
|
663
|
+
)
|
664
|
+
*/
|
665
|
+
|
666
|
+
/*****************************************************
|
667
|
+
* math2/fftlib/rfftf.f
|
668
|
+
****************************************************/
|
669
|
+
/*
|
670
|
+
* SUBROUTINE RFFTF (N,R,WSAVE)
|
671
|
+
* DIMENSION R(*) ,WSAVE(*)
|
672
|
+
*/
|
673
|
+
SUBROUTINE RFFTF(
|
674
|
+
INTEGER N
|
675
|
+
REAL (N) io R
|
676
|
+
REAL (2*N+15) WSAVE
|
677
|
+
)
|
678
|
+
|
679
|
+
/*****************************************************
|
680
|
+
* math2/fftlib/rfftf1.f
|
681
|
+
****************************************************/
|
682
|
+
/*
|
683
|
+
* SUBROUTINE RFFTF1 (N,C,CH,WA,IFAC)
|
684
|
+
* DIMENSION CH(*) ,C(*) ,WA(*) ,IFAC(*)
|
685
|
+
*/
|
686
|
+
/*
|
687
|
+
SUBROUTINE RFFTF1(
|
688
|
+
INTEGER N
|
689
|
+
REAL (*) io C
|
690
|
+
REAL (*) t CH
|
691
|
+
REAL (*) t WA
|
692
|
+
INTEGER (*) t IFAC
|
693
|
+
)
|
694
|
+
*/
|
695
|
+
|
696
|
+
/*****************************************************
|
697
|
+
* math2/fftlib/rffti.f
|
698
|
+
****************************************************/
|
699
|
+
/*
|
700
|
+
* SUBROUTINE RFFTI (N,WSAVE)
|
701
|
+
* DIMENSION WSAVE(*)
|
702
|
+
*/
|
703
|
+
SUBROUTINE RFFTI(
|
704
|
+
INTEGER N
|
705
|
+
REAL (2*N+15) o WSAVE
|
706
|
+
)
|
707
|
+
|
708
|
+
/*****************************************************
|
709
|
+
* math2/fftlib/rffti1.f
|
710
|
+
****************************************************/
|
711
|
+
/*
|
712
|
+
* SUBROUTINE RFFTI1 (N,WA,IFAC)
|
713
|
+
* DIMENSION WA(*) ,IFAC(*) ,NTRYH(4)
|
714
|
+
*/
|
715
|
+
/*
|
716
|
+
SUBROUTINE RFFTI1(
|
717
|
+
INTEGER N
|
718
|
+
REAL (*) o WA
|
719
|
+
INTEGER (*) o IFAC
|
720
|
+
)
|
721
|
+
*/
|
722
|
+
|
723
|
+
/*****************************************************
|
724
|
+
* math2/fftlib/sinqb.f
|
725
|
+
****************************************************/
|
726
|
+
/*
|
727
|
+
* SUBROUTINE SINQB (N,X,WSAVE)
|
728
|
+
* DIMENSION X(*) ,WSAVE(*)
|
729
|
+
*/
|
730
|
+
SUBROUTINE SINQB(
|
731
|
+
INTEGER N
|
732
|
+
REAL (N) io X
|
733
|
+
REAL (3*N+15) WSAVE
|
734
|
+
)
|
735
|
+
|
736
|
+
/*****************************************************
|
737
|
+
* math2/fftlib/sinqf.f
|
738
|
+
****************************************************/
|
739
|
+
/*
|
740
|
+
* SUBROUTINE SINQF (N,X,WSAVE)
|
741
|
+
* DIMENSION X(*) ,WSAVE(*)
|
742
|
+
*/
|
743
|
+
SUBROUTINE SINQF(
|
744
|
+
INTEGER N
|
745
|
+
REAL (N) io X
|
746
|
+
REAL (3*N+15) WSAVE
|
747
|
+
)
|
748
|
+
|
749
|
+
/*****************************************************
|
750
|
+
* math2/fftlib/sinqi.f
|
751
|
+
****************************************************/
|
752
|
+
/*
|
753
|
+
* SUBROUTINE SINQI (N,WSAVE)
|
754
|
+
* DIMENSION WSAVE(*)
|
755
|
+
*/
|
756
|
+
SUBROUTINE SINQI(
|
757
|
+
INTEGER N
|
758
|
+
REAL (3*N+15) o WSAVE
|
759
|
+
)
|
760
|
+
|
761
|
+
/*****************************************************
|
762
|
+
* math2/fftlib/sint.f
|
763
|
+
****************************************************/
|
764
|
+
/*
|
765
|
+
* SUBROUTINE SINT (N,X,WSAVE)
|
766
|
+
* DIMENSION X(*) ,WSAVE(*)
|
767
|
+
*/
|
768
|
+
SUBROUTINE SINT(
|
769
|
+
INTEGER N
|
770
|
+
REAL (N) io X
|
771
|
+
REAL (2*N+N/2+15) WSAVE
|
772
|
+
)
|
773
|
+
|
774
|
+
/*****************************************************
|
775
|
+
* math2/fftlib/sint1.f
|
776
|
+
****************************************************/
|
777
|
+
/*
|
778
|
+
* SUBROUTINE SINT1(N,WAR,WAS,XH,X,IFAC)
|
779
|
+
* DIMENSION WAR(*),WAS(*),X(*),XH(*),IFAC(*)
|
780
|
+
*/
|
781
|
+
/*
|
782
|
+
SUBROUTINE SINT1(
|
783
|
+
INTEGER N
|
784
|
+
REAL (*) io WAR
|
785
|
+
REAL (*) t WAS
|
786
|
+
REAL (*) t XH
|
787
|
+
REAL t X
|
788
|
+
INTEGER (*) t IFAC
|
789
|
+
)
|
790
|
+
*/
|
791
|
+
|
792
|
+
/*****************************************************
|
793
|
+
* math2/fftlib/sinti.f
|
794
|
+
****************************************************/
|
795
|
+
/*
|
796
|
+
* SUBROUTINE SINTI (N,WSAVE)
|
797
|
+
* DIMENSION WSAVE(*)
|
798
|
+
*/
|
799
|
+
SUBROUTINE SINTI(
|
800
|
+
INTEGER N
|
801
|
+
REAL (2*N+N/2+15) o WSAVE
|
802
|
+
)
|
803
|
+
|