ruby-dcl 1.6.2
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- 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
|
+
|