trunk/models/set_intervals.tcl

# $Revision: 1.1 $
# $Date: 1997/09/26 16:40:22 $
# $Author: ballan $
# Interval list creation functions
# Copyright 1996 Benjamin Allan (ballan@edrc.cmu.edu)
#
# This file defines several functions for creating lists of
# numbers and passing them to the integration interface in
# ASCEND 4. These lists are the sample points at which we
# wish to record values during integration.
#
# We have functions for:
# - evenly spaced intervals starting at 0.0
#   (set_int)
# - evenly spaced intervals starting at 0.0 with
#   intermediate sample points at the 2nd order lagrange locations
#   (set_lagrangeint)
#
# These functions should take a starting point instead of assuming 0 but don't.
#
# We want functions for:
#  logarithmically spaced points, starting at some number and ending
#  at some much larger number, for sampling growing functions.
#
#  points spaced at multiple frequencies
#
#  building lists of sample times
#
#  taking a list, inserting a small initial sample to trick the
#  integrator into thinking the problem stiff, and then generating
#  whatever the desired sampling schedule is after that.
#
# Bugs:
#  Poor design. There should be a suite of commands for building
#  sample lists up, then a single command that takes such a list
#  and sends it to the integrator.
#
# proc set_int {nsteps step units}
#-------------------------------------------------------------------------
# this proc sets intervals incrementally for blsode
#-------------------------------------------------------------------------
proc set_int {nsteps step units} {
    set SetIntervals "integrate_set_samples $units";
    for {set i 0} {$i < $nsteps} {incr i} {
    lappend SetIntervals [expr $i*$step];
    }
    eval $SetIntervals;
}

#
# proc set_lagrangeint {nsteps step units}
#-------------------------------------------------------------------------
# this proc sets intervals incrementally for blsode
# adding the lagrange nodes within each regularly spaced step
#-------------------------------------------------------------------------
proc set_lagrangeint {nsteps step units} {
    set SetIntervals "integrate_set_samples $units";
    for {set i 0} {$i < $nsteps} {incr i} {
    lappend SetIntervals [expr $i*$step];
        lappend SetIntervals [expr $i*$step+0.21132486*$step];
        lappend SetIntervals [expr $i*$step+0.7886751*$step];
    }
    eval $SetIntervals;
}

1	aw0a	1	# $Revision: 1.1 $
2			# $Date: 1997/09/26 16:40:22 $
3			# $Author: ballan $
4			# Interval list creation functions
5			# Copyright 1996 Benjamin Allan (ballan@edrc.cmu.edu)
6			#
7			# This file defines several functions for creating lists of
8			# numbers and passing them to the integration interface in
9			# ASCEND 4. These lists are the sample points at which we
10			# wish to record values during integration.
11			#
12			# We have functions for:
13			# - evenly spaced intervals starting at 0.0
14			# (set_int)
15			# - evenly spaced intervals starting at 0.0 with
16			# intermediate sample points at the 2nd order lagrange locations
17			# (set_lagrangeint)
18			#
19			# These functions should take a starting point instead of assuming 0 but don't.
20			#
21			# We want functions for:
22			# logarithmically spaced points, starting at some number and ending
23			# at some much larger number, for sampling growing functions.
24			#
25			# points spaced at multiple frequencies
26			#
27			# building lists of sample times
28			#
29			# taking a list, inserting a small initial sample to trick the
30			# integrator into thinking the problem stiff, and then generating
31			# whatever the desired sampling schedule is after that.
32			#
33			# Bugs:
34			# Poor design. There should be a suite of commands for building
35			# sample lists up, then a single command that takes such a list
36			# and sends it to the integrator.
37			#
38			# proc set_int {nsteps step units}
39			#-------------------------------------------------------------------------
40			# this proc sets intervals incrementally for blsode
41			#-------------------------------------------------------------------------
42			proc set_int {nsteps step units} {
43			set SetIntervals "integrate_set_samples $units";
44			for {set i 0} {$i < $nsteps} {incr i} {
45			lappend SetIntervals [expr $i*$step];
46			}
47			eval $SetIntervals;
48			}
49
50			#
51			# proc set_lagrangeint {nsteps step units}
52			#-------------------------------------------------------------------------
53			# this proc sets intervals incrementally for blsode
54			# adding the lagrange nodes within each regularly spaced step
55			#-------------------------------------------------------------------------
56			proc set_lagrangeint {nsteps step units} {
57			set SetIntervals "integrate_set_samples $units";
58			for {set i 0} {$i < $nsteps} {incr i} {
59			lappend SetIntervals [expr $i*$step];
60			lappend SetIntervals [expr $i$step+0.21132486$step];
61			lappend SetIntervals [expr $i$step+0.7886751$step];
62			}
63			eval $SetIntervals;
64			}
65