annotate doc/interpreter/eval.txi @ 6641:27ec13d8499b

[project @ 2007-05-21 19:53:58 by jwe]
author jwe
date Mon, 21 May 2007 19:53:58 +0000
parents 5a5a09d7deb8
children 083721ae3dfa
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3294
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
1 @c Copyright (C) 1996, 1997 John W. Eaton
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
2 @c This is part of the Octave manual.
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
3 @c For copying conditions, see the file gpl.texi.
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
4
4167
aae05d51353c [project @ 2002-11-12 02:52:50 by jwe]
jwe
parents: 3371
diff changeset
5 @node Evaluation
3294
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
6 @chapter Evaluation
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
7
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
8 Normally, you evaluate expressions simply by typing them at the Octave
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
9 prompt, or by asking Octave to interpret commands that you have saved in
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
10 a file.
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
11
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
12 Sometimes, you may find it necessary to evaluate an expression that has
6641
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
13 been computed and stored in a string, which is exactly what the
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
14 @code{eval} lets you do.
3294
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
15
3371
86873384cd10 [project @ 1999-11-21 17:31:07 by jwe]
jwe
parents: 3294
diff changeset
16 @DOCSTRING(eval)
3294
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
17
6641
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
18 @menu
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
19 * Calling a Function by its Name::
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
20 * Evaluation in a Different Context::
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
21 @end menu
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
22
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
23 @node Calling a Function by its Name
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
24 @section Calling a Function by its Name
3294
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
25
6641
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
26 The @code{feval} function allow you to call a function from a string
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
27 containing its name. This is useful when writing a function that need to
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
28 call user-supplied functions. The @code{feval} function takes the name
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
29 of the function to call as its first argument, and the remaining
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
30 arguments are given to the function.
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
31
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
32 The following example is a simple-minded function using @code{feval}
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
33 that finds the root of a user-supplied function of one variable using
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
34 Newton's method.
3294
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
35
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
36 @example
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
37 @group
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
38 @cindex Fordyce, A. P.
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
39 @findex newtroot
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
40 function result = newtroot (fname, x)
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
41
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
42 # usage: newtroot (fname, x)
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
43 #
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
44 # fname : a string naming a function f(x).
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
45 # x : initial guess
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
46
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
47 delta = tol = sqrt (eps);
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
48 maxit = 200;
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
49 fx = feval (fname, x);
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
50 for i = 1:maxit
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
51 if (abs (fx) < tol)
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
52 result = x;
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
53 return;
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
54 else
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
55 fx_new = feval (fname, x + delta);
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
56 deriv = (fx_new - fx) / delta;
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
57 x = x - fx / deriv;
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
58 fx = fx_new;
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
59 endif
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
60 endfor
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
61
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
62 result = x;
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
63
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
64 endfunction
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
65 @end group
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
66 @end example
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
67
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
68 Note that this is only meant to be an example of calling user-supplied
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
69 functions and should not be taken too seriously. In addition to using a
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
70 more robust algorithm, any serious code would check the number and type
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
71 of all the arguments, ensure that the supplied function really was a
6641
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
72 function, etc. @xref{Predicates for Numeric Objects}, for example,
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
73 for a list of predicates for numeric objects, and see @ref{Status of
3294
bfe1573bd2ae [project @ 1999-10-19 10:06:07 by jwe]
jwe
parents:
diff changeset
74 Variables}, for a description of the @code{exist} function.
6549
5a5a09d7deb8 [project @ 2007-04-20 06:55:29 by jwe]
jwe
parents: 4463
diff changeset
75
6641
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
76 @DOCSTRING(feval)
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
77
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
78 @node Evaluation in a Different Context
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
79 @section Evaluation in a Different Context
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
80
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
81 Before you evaluate an expression you need to substitute
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
82 the values of the variables used in the expression. These
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
83 are stored in the symbol table. Whenever the interpreter
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
84 starts a new function it saves the current symbol table
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
85 and creates a new one, initializing it with the list of
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
86 function parameters and a couple of predefined variables
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
87 such as @code{nargin}. Expressions inside the function use the
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
88 new symbol table.
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
89
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
90 Sometimes you want to write a function so that when you
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
91 call it, it modifies variables in your own context. This
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
92 allows you to use a pass-by-name style of function,
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
93 which is similar to using a pointer in programming languages such
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
94 as C.
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
95
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
96 Consider how you might write @code{save} and @code{load} as
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
97 m-files. For example,
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
98
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
99 @example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
100 function create_data
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
101 x = linspace (0, 10, 10);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
102 y = sin (x);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
103 save mydata x y
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
104 endfunction
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
105 @end example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
106
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
107 With @code{evalin}, you could write save as follows:
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
108
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
109 @example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
110 function save (file, name1, name2)
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
111 f = open_save_file (file);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
112 save_var(f, name1, evalin ("caller", name1));
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
113 save_var(f, name2, evalin ("caller", name2));
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
114 endfunction
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
115 @end example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
116
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
117 @noindent
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
118 Here, @samp{caller} is the @code{create_data} function and @code{name1}
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
119 is the string @code{"x"}, which evaluates simply as the value of @code{x}.
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
120
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
121 You later want to load the values back from @code{mydata}
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
122 in a different context:
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
123
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
124 @example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
125 function process_data
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
126 load mydata
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
127 @dots{} do work @dots{}
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
128 endfunction
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
129 @end example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
130
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
131 @noindent
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
132 With @code{assignin}, you could write @code{load} as follows:
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
133
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
134 @example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
135 function load (file)
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
136 f = open_load_file (file);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
137 [name, val] = load_var (f);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
138 assignin ("caller", name, val);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
139 [name, val] = load_var (f);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
140 assignin ("caller", name, val);
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
141 endfunction
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
142 @end example
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
143
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
144 @noindent
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
145 Here, @samp{caller} is the @code{process_data} function.
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
146
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
147 You can set and use variables at the command prompt
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
148 using the context @samp{base} rather than @samp{caller}.
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
149
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
150 These functions are rarely used in practice. One
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
151 example is the @code{fail (@samp{code}, @samp{pattern})} function
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
152 which evaluates @samp{code} in the caller's context and
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
153 checks that the error message it produces matches
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
154 the given pattern. Other examples such as @code{save} and @code{load}
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
155 are written in C++ where all octave variables
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
156 are in the @samp{caller} context and @code{evalin} is not needed.
27ec13d8499b [project @ 2007-05-21 19:53:58 by jwe]
jwe
parents: 6549
diff changeset
157
6549
5a5a09d7deb8 [project @ 2007-04-20 06:55:29 by jwe]
jwe
parents: 4463
diff changeset
158 @DOCSTRING(evalin)
5a5a09d7deb8 [project @ 2007-04-20 06:55:29 by jwe]
jwe
parents: 4463
diff changeset
159
5a5a09d7deb8 [project @ 2007-04-20 06:55:29 by jwe]
jwe
parents: 4463
diff changeset
160 @DOCSTRING(assignin)