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,32 @@
|
|
1
|
+
/*
|
2
|
+
* math2_intrlib.fp
|
3
|
+
*
|
4
|
+
* $Id: math2_intrlib.fp,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
|
5
|
+
*/
|
6
|
+
|
7
|
+
/*****************************************************
|
8
|
+
* math2/intrlib/vrintr.f
|
9
|
+
****************************************************/
|
10
|
+
/*
|
11
|
+
* SUBROUTINE VRINTR(RX,N,JX)
|
12
|
+
* REAL RX(*)
|
13
|
+
*/
|
14
|
+
SUBROUTINE VRINTR(
|
15
|
+
REAL (JX*(N-1)+1) io RX
|
16
|
+
INTEGER N
|
17
|
+
INTEGER JX
|
18
|
+
)
|
19
|
+
|
20
|
+
/*****************************************************
|
21
|
+
* math2/intrlib/vcintr.f
|
22
|
+
****************************************************/
|
23
|
+
/*
|
24
|
+
* SUBROUTINE VCINTR(CX,N,JX)
|
25
|
+
* COMPLEX CX(*),CMISS,CD,C1
|
26
|
+
*/
|
27
|
+
SUBROUTINE VCINTR(
|
28
|
+
COMPLEX (JX*(N-1)+1) io CX
|
29
|
+
INTEGER N
|
30
|
+
INTEGER JX
|
31
|
+
)
|
32
|
+
|
@@ -0,0 +1,568 @@
|
|
1
|
+
/*
|
2
|
+
* math2_odelib.fp
|
3
|
+
*
|
4
|
+
* $Id: math2_odelib.fp_notused,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
|
5
|
+
*/
|
6
|
+
|
7
|
+
/*****************************************************
|
8
|
+
* math2/odelib/odrk1.f
|
9
|
+
****************************************************/
|
10
|
+
/*
|
11
|
+
* SUBROUTINE ODRK1(N, FCN, T, DT, X, DX, XOUT, WORK)
|
12
|
+
* DIMENSION X(N), XOUT(N), DX(N)
|
13
|
+
*/
|
14
|
+
SUBROUTINE ODRK1(
|
15
|
+
INTEGER N
|
16
|
+
REAL FCN
|
17
|
+
REAL T
|
18
|
+
REAL DT
|
19
|
+
REAL (N) X
|
20
|
+
REAL (N) DX
|
21
|
+
REAL (N) o XOUT
|
22
|
+
REAL t WORK
|
23
|
+
)
|
24
|
+
|
25
|
+
/*****************************************************
|
26
|
+
* math2/odelib/odrk2.f
|
27
|
+
****************************************************/
|
28
|
+
/*
|
29
|
+
* SUBROUTINE ODRK2(N, FCN, T, DT, X, DX, XOUT, WORK)
|
30
|
+
* DIMENSION X(N), XOUT(N), DX(N), WORK(N)
|
31
|
+
*/
|
32
|
+
SUBROUTINE ODRK2(
|
33
|
+
INTEGER N
|
34
|
+
REAL FCN
|
35
|
+
REAL T
|
36
|
+
REAL DT
|
37
|
+
REAL (N) X
|
38
|
+
REAL (N) DX
|
39
|
+
REAL (N) o XOUT
|
40
|
+
REAL (N) t WORK
|
41
|
+
)
|
42
|
+
|
43
|
+
/*****************************************************
|
44
|
+
* math2/odelib/odrk4.f
|
45
|
+
****************************************************/
|
46
|
+
/*
|
47
|
+
* SUBROUTINE ODRK4(N, FCN, T, DT, X, DX, XOUT, WORK)
|
48
|
+
* DIMENSION X(N), XOUT(N), DX(N), WORK(N,3)
|
49
|
+
*/
|
50
|
+
SUBROUTINE ODRK4(
|
51
|
+
INTEGER N
|
52
|
+
REAL FCN
|
53
|
+
REAL T
|
54
|
+
REAL DT
|
55
|
+
REAL (N) X
|
56
|
+
REAL (N) DX
|
57
|
+
REAL (N) o XOUT
|
58
|
+
REAL (N,3) t WORK
|
59
|
+
)
|
60
|
+
|
61
|
+
/*****************************************************
|
62
|
+
* math2/odelib/odrk4r.f
|
63
|
+
****************************************************/
|
64
|
+
/*
|
65
|
+
* SUBROUTINE ODRK4R(N, FCN, T, DT, EPSL, X, WORK)
|
66
|
+
* DIMENSION X(N), WORK(N,7)
|
67
|
+
*/
|
68
|
+
SUBROUTINE ODRK4R(
|
69
|
+
INTEGER N
|
70
|
+
REAL FCN
|
71
|
+
REAL io T
|
72
|
+
REAL io DT
|
73
|
+
REAL EPSL
|
74
|
+
REAL (N) io X
|
75
|
+
REAL (N,7) t WORK
|
76
|
+
)
|
77
|
+
|
78
|
+
/*****************************************************
|
79
|
+
* math2/odelib/odrk4s.f
|
80
|
+
****************************************************/
|
81
|
+
/*
|
82
|
+
* SUBROUTINE ODRK4S(N, FCN, T, DT, EPSL, X, WORK)
|
83
|
+
* DIMENSION X(N), WORK(N,5)
|
84
|
+
*/
|
85
|
+
SUBROUTINE ODRK4S(
|
86
|
+
INTEGER N
|
87
|
+
REAL FCN
|
88
|
+
REAL io T
|
89
|
+
REAL io DT
|
90
|
+
REAL EPSL
|
91
|
+
REAL (N) io X
|
92
|
+
REAL (N,5) t WORK
|
93
|
+
)
|
94
|
+
|
95
|
+
/*****************************************************
|
96
|
+
* math2/odelib/odrkdu.f
|
97
|
+
****************************************************/
|
98
|
+
/* NB: (ODRKDU) WORK size is (N,5) (ALGR=ODRK4), (N,3) (otherwise). */
|
99
|
+
/*
|
100
|
+
* SUBROUTINE ODRKDU(N, ALGR, FCN, T, TEND, ISTEP, X, WORK)
|
101
|
+
* DIMENSION X(N), WORK(N,*)
|
102
|
+
*/
|
103
|
+
SUBROUTINE ODRKDU(
|
104
|
+
INTEGER N
|
105
|
+
REAL ALGR
|
106
|
+
REAL FCN
|
107
|
+
REAL io T
|
108
|
+
REAL TEND
|
109
|
+
INTEGER ISTEP
|
110
|
+
REAL (N) io X
|
111
|
+
REAL (N,*) t WORK
|
112
|
+
)
|
113
|
+
|
114
|
+
/*****************************************************
|
115
|
+
* math2/odelib/odrkdv.f
|
116
|
+
****************************************************/
|
117
|
+
/* NB: (ODRKDV) WORK size is (N,7) (ALGR=ODRK4), (N,5) (otherwise). */
|
118
|
+
/*
|
119
|
+
* SUBROUTINE ODRKDV(N, STEPER, FCN, T, TEND, DT, X, WORK)
|
120
|
+
* DIMENSION X(N), WORK(N,*)
|
121
|
+
*/
|
122
|
+
SUBROUTINE ODRKDV(
|
123
|
+
INTEGER N
|
124
|
+
REAL STEPER
|
125
|
+
REAL FCN
|
126
|
+
REAL io T
|
127
|
+
REAL TEND
|
128
|
+
REAL io DT
|
129
|
+
REAL (N) io X
|
130
|
+
REAL (N,*) t WORK
|
131
|
+
)
|
132
|
+
|
133
|
+
/*****************************************************
|
134
|
+
* math2/odelib/odrkg.f
|
135
|
+
****************************************************/
|
136
|
+
/*
|
137
|
+
* SUBROUTINE ODRKG(N, FCN, T, DT, X, DX, XOUT, WORK)
|
138
|
+
* DIMENSION X(N), DX(N), XOUT(N), WORK(N)
|
139
|
+
*/
|
140
|
+
SUBROUTINE ODRKG(
|
141
|
+
INTEGER N
|
142
|
+
REAL FCN
|
143
|
+
REAL T
|
144
|
+
REAL DT
|
145
|
+
REAL X
|
146
|
+
REAL (N) io DX
|
147
|
+
REAL (N) o XOUT
|
148
|
+
REAL (N) t WORK
|
149
|
+
)
|
150
|
+
|
151
|
+
/*****************************************************
|
152
|
+
* math2/odelib/odrkgr.f
|
153
|
+
****************************************************/
|
154
|
+
/*
|
155
|
+
* SUBROUTINE ODRKGR(N, FCN, T, DT, EPSL, X, WORK)
|
156
|
+
* DIMENSION X(N), WORK(N,5)
|
157
|
+
*/
|
158
|
+
SUBROUTINE ODRKGR(
|
159
|
+
INTEGER N
|
160
|
+
REAL FCN
|
161
|
+
REAL io T
|
162
|
+
REAL io DT
|
163
|
+
REAL EPSL
|
164
|
+
REAL (N) io X
|
165
|
+
REAL (N,5) t WORK
|
166
|
+
)
|
167
|
+
|
168
|
+
/*****************************************************
|
169
|
+
* math2/odelib/odrkgs.f
|
170
|
+
****************************************************/
|
171
|
+
/*
|
172
|
+
* SUBROUTINE ODRKGS(N, FCN, T, DT, EPSL, X, WORK)
|
173
|
+
* DIMENSION X(N), WORK(N,3)
|
174
|
+
*/
|
175
|
+
SUBROUTINE ODRKGS(
|
176
|
+
INTEGER N
|
177
|
+
REAL FCN
|
178
|
+
REAL io T
|
179
|
+
REAL io DT
|
180
|
+
REAL EPSL
|
181
|
+
REAL (N) io X
|
182
|
+
REAL (N,3) t WORK
|
183
|
+
)
|
184
|
+
|
185
|
+
/*****************************************************
|
186
|
+
* math2/odelib/odiget.f
|
187
|
+
****************************************************/
|
188
|
+
/*
|
189
|
+
* SUBROUTINE ODIGET(CP, IPARA)
|
190
|
+
* CHARACTER CP*(*)
|
191
|
+
*/
|
192
|
+
SUBROUTINE ODIGET(
|
193
|
+
CHARACTER*(*) CP
|
194
|
+
INTEGER o IPARA
|
195
|
+
)
|
196
|
+
|
197
|
+
/*
|
198
|
+
* ENTRY ODISET(CP, IPARA)
|
199
|
+
* CHARACTER CP*(*)
|
200
|
+
*/
|
201
|
+
SUBROUTINE ODISET(
|
202
|
+
CHARACTER*(*) CP
|
203
|
+
INTEGER IPARA
|
204
|
+
)
|
205
|
+
|
206
|
+
/*
|
207
|
+
* ENTRY ODISTX(CP, IPARA)
|
208
|
+
* CHARACTER CP*(*)
|
209
|
+
*/
|
210
|
+
SUBROUTINE ODISTX(
|
211
|
+
CHARACTER*(*) CP
|
212
|
+
INTEGER IPARA
|
213
|
+
)
|
214
|
+
|
215
|
+
/*****************************************************
|
216
|
+
* math2/odelib/odiqnp.f
|
217
|
+
****************************************************/
|
218
|
+
/*
|
219
|
+
* SUBROUTINE ODIQNP(NCP)
|
220
|
+
*/
|
221
|
+
SUBROUTINE ODIQNP(
|
222
|
+
INTEGER o NCP
|
223
|
+
)
|
224
|
+
|
225
|
+
/*
|
226
|
+
* ENTRY ODIQID(CP,IDX)
|
227
|
+
* CHARACTER CP*(*)
|
228
|
+
*/
|
229
|
+
SUBROUTINE ODIQID(
|
230
|
+
CHARACTER*(*) CP
|
231
|
+
INTEGER o IDX
|
232
|
+
)
|
233
|
+
|
234
|
+
/*
|
235
|
+
* ENTRY ODIQCP(IDX,CP)
|
236
|
+
* CHARACTER CP*(*)
|
237
|
+
*/
|
238
|
+
SUBROUTINE ODIQCP(
|
239
|
+
INTEGER IDX
|
240
|
+
CHARACTER*(*) o CP
|
241
|
+
)
|
242
|
+
|
243
|
+
/*
|
244
|
+
* ENTRY ODIQCL(IDX, CP)
|
245
|
+
* CHARACTER CP*(*)
|
246
|
+
*/
|
247
|
+
SUBROUTINE ODIQCL(
|
248
|
+
INTEGER IDX
|
249
|
+
CHARACTER*(*) o CP
|
250
|
+
)
|
251
|
+
|
252
|
+
/*
|
253
|
+
* ENTRY ODIQVL(IDX,IPARA)
|
254
|
+
*/
|
255
|
+
SUBROUTINE ODIQVL(
|
256
|
+
INTEGER IDX
|
257
|
+
INTEGER o IPARA
|
258
|
+
)
|
259
|
+
|
260
|
+
/*
|
261
|
+
* ENTRY ODISVL(IDX,IPARA)
|
262
|
+
*/
|
263
|
+
SUBROUTINE ODISVL(
|
264
|
+
INTEGER IDX
|
265
|
+
INTEGER IPARA
|
266
|
+
)
|
267
|
+
|
268
|
+
/*
|
269
|
+
* ENTRY ODIQIN(CP, IN)
|
270
|
+
* CHARACTER CP*(*)
|
271
|
+
*/
|
272
|
+
SUBROUTINE ODIQIN(
|
273
|
+
CHARACTER*(*) CP
|
274
|
+
INTEGER o IN
|
275
|
+
)
|
276
|
+
|
277
|
+
/*****************************************************
|
278
|
+
* math2/odelib/odlget.f
|
279
|
+
****************************************************/
|
280
|
+
/*
|
281
|
+
* SUBROUTINE ODLGET(CP, LPARA)
|
282
|
+
* CHARACTER CP*(*)
|
283
|
+
* LOGICAL LPARA
|
284
|
+
*/
|
285
|
+
SUBROUTINE ODLGET(
|
286
|
+
CHARACTER*(*) CP
|
287
|
+
LOGICAL o LPARA
|
288
|
+
)
|
289
|
+
|
290
|
+
/*
|
291
|
+
* ENTRY ODLSET(CP, LPARA)
|
292
|
+
* CHARACTER CP*(*)
|
293
|
+
* LOGICAL LPARA
|
294
|
+
*/
|
295
|
+
SUBROUTINE ODLSET(
|
296
|
+
CHARACTER*(*) CP
|
297
|
+
LOGICAL LPARA
|
298
|
+
)
|
299
|
+
|
300
|
+
/*
|
301
|
+
* ENTRY ODLSTX(CP, LPARA)
|
302
|
+
* CHARACTER CP*(*)
|
303
|
+
* LOGICAL LPARA
|
304
|
+
*/
|
305
|
+
SUBROUTINE ODLSTX(
|
306
|
+
CHARACTER*(*) CP
|
307
|
+
LOGICAL LPARA
|
308
|
+
)
|
309
|
+
|
310
|
+
/*****************************************************
|
311
|
+
* math2/odelib/odlqnp.f
|
312
|
+
****************************************************/
|
313
|
+
/*
|
314
|
+
* SUBROUTINE ODLQNP(NCP)
|
315
|
+
*/
|
316
|
+
SUBROUTINE ODLQNP(
|
317
|
+
INTEGER o NCP
|
318
|
+
)
|
319
|
+
|
320
|
+
/*
|
321
|
+
* ENTRY ODLQID(CP,IDX)
|
322
|
+
* CHARACTER CP*(*)
|
323
|
+
*/
|
324
|
+
SUBROUTINE ODLQID(
|
325
|
+
CHARACTER*(*) CP
|
326
|
+
INTEGER o IDX
|
327
|
+
)
|
328
|
+
|
329
|
+
/*
|
330
|
+
* ENTRY ODLQCP(IDX,CP)
|
331
|
+
* CHARACTER CP*(*)
|
332
|
+
*/
|
333
|
+
SUBROUTINE ODLQCP(
|
334
|
+
INTEGER IDX
|
335
|
+
CHARACTER*(*) o CP
|
336
|
+
)
|
337
|
+
|
338
|
+
/*
|
339
|
+
* ENTRY ODLQCL(IDX, CP)
|
340
|
+
* CHARACTER CP*(*)
|
341
|
+
*/
|
342
|
+
SUBROUTINE ODLQCL(
|
343
|
+
INTEGER IDX
|
344
|
+
CHARACTER*(*) o CP
|
345
|
+
)
|
346
|
+
|
347
|
+
/*
|
348
|
+
* ENTRY ODLQVL(IDX, LPARA)
|
349
|
+
* LOGICAL LPARA
|
350
|
+
*/
|
351
|
+
SUBROUTINE ODLQVL(
|
352
|
+
INTEGER IDX
|
353
|
+
LOGICAL o LPARA
|
354
|
+
)
|
355
|
+
|
356
|
+
/*
|
357
|
+
* ENTRY ODLSVL(IDX,LPARA)
|
358
|
+
* LOGICAL LPARA
|
359
|
+
*/
|
360
|
+
SUBROUTINE ODLSVL(
|
361
|
+
INTEGER IDX
|
362
|
+
LOGICAL LPARA
|
363
|
+
)
|
364
|
+
|
365
|
+
/*
|
366
|
+
* ENTRY ODLQIN(CP, IN)
|
367
|
+
* CHARACTER CP*(*)
|
368
|
+
*/
|
369
|
+
SUBROUTINE ODLQIN(
|
370
|
+
CHARACTER*(*) CP
|
371
|
+
INTEGER o IN
|
372
|
+
)
|
373
|
+
|
374
|
+
/*****************************************************
|
375
|
+
* math2/odelib/odpget.f
|
376
|
+
****************************************************/
|
377
|
+
/*
|
378
|
+
* SUBROUTINE ODPGET(CP, IPARA)
|
379
|
+
* CHARACTER CP*(*)
|
380
|
+
*/
|
381
|
+
/* SUBROUTINE ODPGET(
|
382
|
+
* CHARACTER*(*) CP
|
383
|
+
* INTEGER o IPARA
|
384
|
+
* )
|
385
|
+
*/
|
386
|
+
|
387
|
+
/*
|
388
|
+
* ENTRY ODPSET(CP, IPARA)
|
389
|
+
* CHARACTER CP*(*)
|
390
|
+
*/
|
391
|
+
/* SUBROUTINE ODPSET(
|
392
|
+
* CHARACTER*(*) CP
|
393
|
+
* INTEGER IPARA
|
394
|
+
* )
|
395
|
+
*/
|
396
|
+
|
397
|
+
/*
|
398
|
+
* ENTRY ODPSTX(CP, IPARA)
|
399
|
+
* CHARACTER CP*(*)
|
400
|
+
*/
|
401
|
+
/* SUBROUTINE ODPSTX(
|
402
|
+
* CHARACTER*(*) CP
|
403
|
+
* INTEGER IPARA
|
404
|
+
* )
|
405
|
+
*/
|
406
|
+
|
407
|
+
/*****************************************************
|
408
|
+
* math2/odelib/odpqnp.f
|
409
|
+
****************************************************/
|
410
|
+
/*
|
411
|
+
* SUBROUTINE ODPQNP(NCP)
|
412
|
+
*/
|
413
|
+
SUBROUTINE ODPQNP(
|
414
|
+
INTEGER o NCP
|
415
|
+
)
|
416
|
+
|
417
|
+
/*
|
418
|
+
* ENTRY ODPQID(CP, IDX)
|
419
|
+
* CHARACTER CP*(*)
|
420
|
+
*/
|
421
|
+
SUBROUTINE ODPQID(
|
422
|
+
CHARACTER*(*) CP
|
423
|
+
INTEGER o IDX
|
424
|
+
)
|
425
|
+
|
426
|
+
/*
|
427
|
+
* ENTRY ODPQCP(IDX, CP)
|
428
|
+
* CHARACTER CP*(*)
|
429
|
+
*/
|
430
|
+
SUBROUTINE ODPQCP(
|
431
|
+
INTEGER IDX
|
432
|
+
CHARACTER*(*) o CP
|
433
|
+
)
|
434
|
+
|
435
|
+
/*
|
436
|
+
* ENTRY ODPQCL(IDX, CP)
|
437
|
+
* CHARACTER CP*(*)
|
438
|
+
*/
|
439
|
+
SUBROUTINE ODPQCL(
|
440
|
+
INTEGER IDX
|
441
|
+
CHARACTER*(*) o CP
|
442
|
+
)
|
443
|
+
|
444
|
+
/*
|
445
|
+
* ENTRY ODPQIT(IDX, ITP)
|
446
|
+
*/
|
447
|
+
SUBROUTINE ODPQIT(
|
448
|
+
INTEGER IDX
|
449
|
+
INTEGER o ITP
|
450
|
+
)
|
451
|
+
|
452
|
+
/*
|
453
|
+
* ENTRY ODPQVL(IDX, IPARA)
|
454
|
+
*/
|
455
|
+
SUBROUTINE ODPQVL(
|
456
|
+
INTEGER IDX
|
457
|
+
INTEGER o IPARA
|
458
|
+
)
|
459
|
+
|
460
|
+
/*
|
461
|
+
* ENTRY ODPSVL(IDX, IPARA)
|
462
|
+
*/
|
463
|
+
SUBROUTINE ODPSVL(
|
464
|
+
INTEGER IDX
|
465
|
+
INTEGER IPARA
|
466
|
+
)
|
467
|
+
|
468
|
+
/*
|
469
|
+
* ENTRY ODPQIN(CP, IN)
|
470
|
+
* CHARACTER CP*(*)
|
471
|
+
*/
|
472
|
+
SUBROUTINE ODPQIN(
|
473
|
+
CHARACTER*(*) CP
|
474
|
+
INTEGER o IN
|
475
|
+
)
|
476
|
+
|
477
|
+
/*****************************************************
|
478
|
+
* math2/odelib/odrget.f
|
479
|
+
****************************************************/
|
480
|
+
/*
|
481
|
+
* SUBROUTINE ODRGET(CP, RPARA)
|
482
|
+
* CHARACTER CP*(*)
|
483
|
+
*/
|
484
|
+
SUBROUTINE ODRGET(
|
485
|
+
CHARACTER*(*) CP
|
486
|
+
REAL o RPARA
|
487
|
+
)
|
488
|
+
|
489
|
+
/*
|
490
|
+
* ENTRY ODRSET(CP, RPARA)
|
491
|
+
* CHARACTER CP*(*)
|
492
|
+
*/
|
493
|
+
SUBROUTINE ODRSET(
|
494
|
+
CHARACTER*(*) CP
|
495
|
+
REAL RPARA
|
496
|
+
)
|
497
|
+
|
498
|
+
/*
|
499
|
+
* ENTRY ODRSTX(CP, RPARA)
|
500
|
+
* CHARACTER CP*(*)
|
501
|
+
*/
|
502
|
+
SUBROUTINE ODRSTX(
|
503
|
+
CHARACTER*(*) CP
|
504
|
+
REAL RPARA
|
505
|
+
)
|
506
|
+
|
507
|
+
/*****************************************************
|
508
|
+
* math2/odelib/odrqnp.f
|
509
|
+
****************************************************/
|
510
|
+
/*
|
511
|
+
* SUBROUTINE ODRQNP(NCP)
|
512
|
+
*/
|
513
|
+
SUBROUTINE ODRQNP(
|
514
|
+
INTEGER o NCP
|
515
|
+
)
|
516
|
+
|
517
|
+
/*
|
518
|
+
* ENTRY ODRQID(CP,IDX)
|
519
|
+
* CHARACTER CP*(*)
|
520
|
+
*/
|
521
|
+
SUBROUTINE ODRQID(
|
522
|
+
CHARACTER*(*) CP
|
523
|
+
INTEGER o IDX
|
524
|
+
)
|
525
|
+
|
526
|
+
/*
|
527
|
+
* ENTRY ODRQCP(IDX,CP)
|
528
|
+
* CHARACTER CP*(*)
|
529
|
+
*/
|
530
|
+
SUBROUTINE ODRQCP(
|
531
|
+
INTEGER IDX
|
532
|
+
CHARACTER*(*) o CP
|
533
|
+
)
|
534
|
+
|
535
|
+
/*
|
536
|
+
* ENTRY ODRQCL(IDX, CP)
|
537
|
+
* CHARACTER CP*(*)
|
538
|
+
*/
|
539
|
+
SUBROUTINE ODRQCL(
|
540
|
+
INTEGER IDX
|
541
|
+
CHARACTER*(*) o CP
|
542
|
+
)
|
543
|
+
|
544
|
+
/*
|
545
|
+
* ENTRY ODLQVR(IDX, RPARA)
|
546
|
+
*/
|
547
|
+
SUBROUTINE ODLQVR(
|
548
|
+
INTEGER IDX
|
549
|
+
REAL o RPARA
|
550
|
+
)
|
551
|
+
|
552
|
+
/*
|
553
|
+
* ENTRY ODRSVL(IDX,RPARA)
|
554
|
+
*/
|
555
|
+
SUBROUTINE ODRSVL(
|
556
|
+
INTEGER IDX
|
557
|
+
REAL RPARA
|
558
|
+
)
|
559
|
+
|
560
|
+
/*
|
561
|
+
* ENTRY ODRQIN(CP, IN)
|
562
|
+
* CHARACTER CP*(*)
|
563
|
+
*/
|
564
|
+
SUBROUTINE ODRQIN(
|
565
|
+
CHARACTER*(*) CP
|
566
|
+
INTEGER o IN
|
567
|
+
)
|
568
|
+
|