|
1
|
1 // builtins.cc -*- C++ -*- |
|
|
2 /* |
|
|
3 |
|
|
4 Copyright (C) 1992, 1993 John W. Eaton |
|
|
5 |
|
|
6 This file is part of Octave. |
|
|
7 |
|
|
8 Octave is free software; you can redistribute it and/or modify it |
|
|
9 under the terms of the GNU General Public License as published by the |
|
|
10 Free Software Foundation; either version 2, or (at your option) any |
|
|
11 later version. |
|
|
12 |
|
|
13 Octave is distributed in the hope that it will be useful, but WITHOUT |
|
|
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
|
|
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
|
|
16 for more details. |
|
|
17 |
|
|
18 You should have received a copy of the GNU General Public License |
|
|
19 along with Octave; see the file COPYING. If not, write to the Free |
|
|
20 Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. |
|
|
21 |
|
|
22 */ |
|
|
23 |
|
|
24 #ifdef __GNUG__ |
|
|
25 #pragma implementation |
|
|
26 #endif |
|
|
27 |
|
|
28 #include <iostream.h> |
|
|
29 #include <math.h> |
|
|
30 #include <float.h> |
|
|
31 |
|
|
32 #include "tree-const.h" |
|
|
33 #include "symtab.h" |
|
|
34 #include "t-builtins.h" |
|
|
35 #include "g-builtins.h" |
|
|
36 #include "builtins.h" |
|
|
37 #include "octave.h" |
|
|
38 #include "utils.h" |
|
|
39 #include "tree.h" |
|
|
40 #include "mappers.h" |
|
|
41 #include "user-prefs.h" |
|
|
42 #include "variables.h" |
|
|
43 |
|
|
44 // NOTE: nargin == 1 means that the function takes no arguments (just |
|
|
45 // like C, the first argument is (or should be, anyway) the function |
|
|
46 // name). Also, -1 is shorthand for infinity. |
|
|
47 |
|
|
48 // The following initializations may eventually need to be reworked |
|
|
49 // like the builtin functions in bash were around version 1.12... |
|
|
50 |
|
|
51 static builtin_mapper_functions mapper_functions[] = |
|
|
52 { |
|
|
53 { "abs", 2, 1, 0, fabs, abs, NULL, |
|
|
54 "compute abs(X) for each element of X\n", }, |
|
|
55 |
|
|
56 { "acos", 2, 1, 0, acos, NULL, acos, |
|
|
57 "compute acos(X) for each element of X\n", }, |
|
|
58 |
|
|
59 { "acosh", 2, 1, 0, acosh, NULL, acosh, |
|
|
60 "compute acosh(X) for each element of X\n", }, |
|
|
61 |
|
|
62 { "angle", 2, 1, 0, arg, arg, NULL, |
|
|
63 "compute arg(X) for each element of X\n", }, |
|
|
64 |
|
|
65 { "arg", 2, 1, 0, arg, arg, NULL, |
|
|
66 "compute arg(X) for each element of X\n", }, |
|
|
67 |
|
|
68 { "asin", 2, 1, 0, asin, NULL, asin, |
|
|
69 "compute asin(X) for each element of X\n", }, |
|
|
70 |
|
|
71 { "asinh", 2, 1, 0, asinh, NULL, asinh, |
|
|
72 "compute asinh(X) for each element of X\n", }, |
|
|
73 |
|
|
74 { "atan", 2, 1, 0, atan, NULL, atan, |
|
|
75 "compute atan(X) for each element of X\n", }, |
|
|
76 |
|
|
77 { "atanh", 2, 1, 0, atanh, NULL, atanh, |
|
|
78 "compute atanh(X) for each element of X\n", }, |
|
|
79 |
|
|
80 { "ceil", 2, 1, 0, ceil, NULL, ceil, |
|
|
81 "ceil(X): round elements of X toward +Inf\n", }, |
|
|
82 |
|
|
83 { "conj", 2, 1, 0, conj, NULL, conj, |
|
|
84 "compute complex conjugate for each element of X\n", }, |
|
|
85 |
|
|
86 { "cos", 2, 1, 0, cos, NULL, cos, |
|
|
87 "compute cos(X) for each element of X\n", }, |
|
|
88 |
|
|
89 { "cosh", 2, 1, 0, cosh, NULL, cosh, |
|
|
90 "compute cosh(X) for each element of X\n", }, |
|
|
91 |
|
|
92 { "exp", 2, 1, 0, exp, NULL, exp, |
|
|
93 "compute exp(X) for each element of X\n", }, |
|
|
94 |
|
|
95 { "finite", 2, 1, 0, xfinite, xfinite, NULL, |
|
|
96 "finite(X): return 1 for finite elements of X\n", }, |
|
|
97 |
|
|
98 { "fix", 2, 1, 0, fix, NULL, fix, |
|
|
99 "fix(X): round elements of X toward zero\n", }, |
|
|
100 |
|
|
101 { "floor", 2, 1, 0, floor, NULL, floor, |
|
|
102 "floor(X): round elements of X toward -Inf\n", }, |
|
|
103 |
|
|
104 { "isinf", 2, 1, 0, xisinf, xisinf, NULL, |
|
|
105 "isinf(X): return 1 for elements of X infinite\n", }, |
|
|
106 |
|
|
107 { "imag", 2, 1, 0, imag, imag, NULL, |
|
|
108 "return imaginary part for each elements of X\n", }, |
|
|
109 |
|
|
110 #ifdef HAVE_ISNAN |
|
|
111 { "isnan", 2, 1, 0, xisnan, xisnan, NULL, |
|
|
112 "isnan(X): return 1 where elements of X are NaNs\n", }, |
|
|
113 #endif |
|
|
114 |
|
|
115 { "log", 2, 1, 1, log, NULL, log, |
|
|
116 "compute log(X) for each element of X\n", }, |
|
|
117 |
|
|
118 { "log10", 2, 1, 1, log10, NULL, log10, |
|
|
119 "compute log10(X) for each element of X\n", }, |
|
|
120 |
|
|
121 { "real", 2, 1, 0, real, real, NULL, |
|
|
122 "return real part for each element of X\n", }, |
|
|
123 |
|
|
124 { "round", 2, 1, 0, round, NULL, round, |
|
|
125 "round(X): round elements of X to nearest integer\n", }, |
|
|
126 |
|
|
127 { "sign", 2, 1, 0, signum, NULL, signum, |
|
|
128 "sign(X): apply signum function to elements of X\n", }, |
|
|
129 |
|
|
130 { "sin", 2, 1, 0, sin, NULL, sin, |
|
|
131 "compute sin(X) for each element of X\n", }, |
|
|
132 |
|
|
133 { "sinh", 2, 1, 0, sinh, NULL, sinh, |
|
|
134 "compute sinh(X) for each element of X\n", }, |
|
|
135 |
|
|
136 { "sqrt", 2, 1, 1, sqrt, NULL, sqrt, |
|
|
137 "compute sqrt(X) for each element of X\n", }, |
|
|
138 |
|
|
139 { "tan", 2, 1, 0, tan, NULL, tan, |
|
|
140 "compute tan(X) for each element of X\n", }, |
|
|
141 |
|
|
142 { "tanh", 2, 1, 0, tanh, NULL, tanh, |
|
|
143 "compute tanh(X) for each element of X\n", }, |
|
|
144 |
|
|
145 { NULL, -1, -1, -1, NULL, NULL, NULL, NULL, }, |
|
|
146 }; |
|
|
147 |
|
|
148 static builtin_text_functions text_functions[] = |
|
|
149 { |
|
|
150 { "casesen", 2, builtin_casesen, |
|
|
151 "print warning if user tries to change case sensitivity\n", }, |
|
|
152 |
|
|
153 { "cd", 2, builtin_cd, |
|
|
154 "change current working directory\n", }, |
|
|
155 |
|
|
156 { "clear", -1, builtin_clear, |
|
|
157 "clear symbol(s) from symbol table\n", }, |
|
|
158 |
|
|
159 { "dir", -1, builtin_ls, |
|
|
160 "print a directory listing\n", }, |
|
|
161 |
|
|
162 { "document", -1, builtin_document, |
|
|
163 "define help string for symbol\n", }, |
|
|
164 |
|
|
165 { "edit_history", -1, builtin_edit_history, |
|
|
166 "usage: edit_history [first] [last]\n", }, |
|
|
167 |
|
|
168 { "format", -1, builtin_format, |
|
|
169 "set output formatting style\n", }, |
|
|
170 |
|
|
171 { "help", -1, builtin_help, |
|
|
172 "print cryptic yet witty messages\n", }, |
|
|
173 |
|
|
174 { "history", -1, builtin_history, |
|
|
175 "print/save/load command history\n", }, |
|
|
176 |
|
|
177 { "load", -1, builtin_load, |
|
|
178 "load variables from a file\n", }, |
|
|
179 |
|
|
180 { "ls", -1, builtin_ls, |
|
|
181 "print a directory listing\n", }, |
|
|
182 |
|
|
183 { "save", -1, builtin_save, |
|
|
184 "save variables to a file\n", }, |
|
|
185 |
|
|
186 { "set", -1, builtin_set, |
|
|
187 "set plotting options\n", }, |
|
|
188 |
|
|
189 { "show", -1, builtin_show, |
|
|
190 "show plotting options\n", }, |
|
|
191 |
|
|
192 { "who", -1, builtin_who, |
|
|
193 "list symbols\n", }, |
|
|
194 |
|
|
195 { NULL, -1, NULL, NULL, }, |
|
|
196 }; |
|
|
197 |
|
|
198 static builtin_general_functions general_functions[] = |
|
|
199 { |
|
|
200 { "all", 2, 1, builtin_all, |
|
|
201 "all(X): are all elements of X nonzero?\n", }, |
|
|
202 |
|
|
203 { "any", 2, 1, builtin_any, |
|
|
204 "any(X): are any elements of X nonzero?\n", }, |
|
|
205 |
|
19
|
206 { "balance", 4, 4, builtin_balance, |
|
|
207 "balancing for algebraic and generalized eigenvalue prblems\n\ |
|
|
208 Usage: \n\ |
|
|
209 algebraic eigenvalues problem:\n\ |
|
|
210 aa = balance(a{,opt}) or [{dd,}aa] = balance(a{,opt})\n\ |
|
|
211 generalized eigenvalue problem:\n\ |
|
|
212 [cc, dd, aa, bb] = balance (a, b {,opt})\n\ |
|
|
213 where 'opt' is an optional single character argument as follows: \n\ |
|
|
214 'N' or 'n': no balancing; arguments copied, transformation(s) \n\ |
|
|
215 set to identity\n\ |
|
|
216 'P' or 'p': permute argument(s) to isolate eigenvalues where possible\n\ |
|
|
217 'S' or 's': scale to improve accuracy of computed eigenvalues\n\ |
|
|
218 'B' or 'b': (default) permute and scale, in that order. Rows/columns of a \n\ |
|
|
219 (and b) that are isolated by permutation are not scaled\n\ |
|
|
220 [dd, aa] = balance (a, opt) returns aa = dd\a*dd,\n\ |
|
|
221 [cc, dd, aa, bb] = balance (a, b, opt) returns aa (bb) = cc*a*dd (cc*b*dd)\n", }, |
|
|
222 |
|
1
|
223 { "clc", 1, 0, builtin_clc, |
|
|
224 "clear screen\n", }, |
|
|
225 |
|
|
226 { "clock", 1, 0, builtin_clock, |
|
|
227 "return current date and time in vector\n", }, |
|
|
228 |
|
|
229 { "closeplot", 1, 0, builtin_closeplot, |
|
|
230 "close the stream to plotter\n", }, |
|
|
231 |
|
|
232 { "colloc", 7, 4, builtin_colloc, |
|
|
233 "[r, A, B, q] = colloc (n [, 'left'] [, 'right']): collocation weights\n", }, |
|
|
234 |
|
|
235 { "cumprod", 2, 1, builtin_cumprod, |
|
|
236 "cumprod (X): cumulative products\n", }, |
|
|
237 |
|
|
238 { "cumsum", 2, 1, builtin_cumsum, |
|
|
239 "cumsum (X): cumulative sums\n", }, |
|
|
240 |
|
|
241 { "dassl", 5, 1, builtin_dassl, |
|
|
242 "solve DAEs using Petzold's DASSL. Usage:\n\ |
|
|
243 \n\ |
|
|
244 dassl ('function_name', x_0, xdot_0, t_out)\n\ |
|
|
245 dassl ('function_name', x_0, xdot_0, t_out, t_crit)\n\ |
|
|
246 \n\ |
|
|
247 The first argument is the name of the function to call to\n\ |
|
|
248 compute the vector of residuals. It must have the form\n\ |
|
|
249 \n\ |
|
|
250 res = f (x, xdot, t)\n\ |
|
|
251 \n\ |
|
|
252 where x, xdot, and res are vectors, and t is a scalar.\n\n", }, |
|
|
253 |
|
|
254 { "date", 1, 0, builtin_date, |
|
|
255 "return current date in a string\n", }, |
|
|
256 |
|
|
257 { "det", 2, 1, builtin_det, |
|
|
258 "det(X): determinant of a square matrix\n", }, |
|
|
259 |
|
|
260 { "diag", 3, 1, builtin_diag, |
|
|
261 "diag(X [,k]): form/extract diagonals\n", }, |
|
|
262 |
|
|
263 { "disp", 3, 1, builtin_disp, |
|
|
264 "disp(X): display value\n", }, |
|
|
265 |
|
|
266 { "eig", 2, 1, builtin_eig, |
|
|
267 "eig (x): compute eigenvalues and eigenvectors of x\n", }, |
|
|
268 |
|
|
269 { "error", 2, 1, builtin_error, |
|
|
270 "print message and jump to top level\n", }, |
|
|
271 |
|
|
272 { "eval", 2, 1, builtin_eval, |
|
|
273 "evaluate text as octave source\n", }, |
|
|
274 |
|
|
275 { "exist", 2, 1, builtin_exist, |
|
|
276 "check if variable or file exists\n", }, |
|
|
277 |
|
|
278 { "exit", 1, 0, builtin_quit, |
|
|
279 "exit Octave gracefully\n", }, |
|
|
280 |
|
|
281 { "expm", 2, 1, builtin_expm, |
|
|
282 "Matrix exponential, e^A\n", }, |
|
|
283 |
|
|
284 { "eye", 3, 1, builtin_eye, |
|
|
285 "create an identity matrix\n", }, |
|
|
286 |
|
|
287 { "fclose", 2, 1, builtin_fclose, |
|
|
288 "fclose('filename' or filenum): close a file\n", }, |
|
|
289 |
|
|
290 { "feval", -1, 1, builtin_feval, |
|
|
291 "evaluate argument as function\n", }, |
|
|
292 |
|
|
293 { "fflush", 2, 1, builtin_fflush, |
|
|
294 "fflush('filename' or filenum): flush buffered data to output file\n", }, |
|
|
295 |
|
|
296 { "fft", 2, 1, builtin_fft, |
|
|
297 "fft(X): fast fourier transform of a vector\n", }, |
|
|
298 |
|
|
299 { "fgets",3, 2, builtin_fgets, |
|
|
300 "[string, length] = fgets('filename' or filenum, length): read a string from a file\n", }, |
|
|
301 |
|
|
302 { "find", -1, 1, builtin_find, |
|
|
303 "find (x): return vector of indices of nonzero elements\n", }, |
|
|
304 |
|
|
305 { "flops", 2, 1, builtin_flops, |
|
|
306 "count floating point operations\n", }, |
|
|
307 |
|
|
308 { "fopen", 3, 1, builtin_fopen, |
|
|
309 "filenum = fopen('filename','mode'): open a file\n", }, |
|
|
310 |
|
|
311 { "fprintf", -1, 1, builtin_fprintf, |
|
|
312 "fprintf ('file', 'fmt', ...)\n", }, |
|
|
313 |
|
|
314 { "freport", 1, 1, builtin_freport, |
|
|
315 "freport: list open files and their status\n", }, |
|
|
316 |
|
|
317 { "frewind", 2, 1, builtin_frewind, |
|
|
318 "frewind('filename' or filenum): set file position at beginning of file\n", }, |
|
|
319 |
|
|
320 { "fscanf", 3, -1, builtin_fscanf, |
|
|
321 "[a,b,c...] = fscanf ('file', 'fmt')\n", }, |
|
|
322 |
|
|
323 { "fseek", 4, 1, builtin_fseek, |
|
|
324 "fseek('filename' or filenum): set file position for reading or writing\n", }, |
|
|
325 |
|
|
326 { "fsolve", 5, 1, builtin_fsolve, |
|
|
327 "Solve nonlinear equations using Minpack. Usage:\n\ |
|
|
328 \n\ |
|
|
329 [x, info] = fsolve ('f', x0)\n\ |
|
|
330 \n\ |
|
|
331 Where the first argument is the name of the function to call to\n\ |
|
|
332 compute the vector of function values. It must have the form\n\ |
|
|
333 \n\ |
|
|
334 y = f (x) |
|
|
335 \n\ |
|
|
336 where y and x are vectors.\n\n", }, |
|
|
337 |
|
|
338 { "fsqp", 11, 3, builtin_fsqp, |
|
|
339 "solve NLPs\n", }, |
|
|
340 |
|
|
341 { "ftell", 2, 1, builtin_ftell, |
|
|
342 "position = ftell ('filename' or filenum): returns the current file position\n", }, |
|
|
343 |
|
|
344 { "getenv", 2, 1, builtin_getenv, |
|
|
345 "get environment variable values\n", }, |
|
|
346 |
|
|
347 { "hess", 2, 2, builtin_hess, |
|
|
348 "Hessenburg decomposition\n",}, |
|
|
349 |
|
|
350 { "home", 1, 0, builtin_clc, |
|
|
351 "clear screen\n", }, |
|
|
352 |
|
|
353 { "input", 3, 1, builtin_input, |
|
|
354 "input('prompt' [,'s']): prompt user for [string] input\n", }, |
|
|
355 |
|
|
356 { "ifft", 2, 1, builtin_ifft, |
|
|
357 "ifft(X): inverse fast fourier transform of a vector\n", }, |
|
|
358 |
|
|
359 { "inv", 2, 1, builtin_inv, |
|
|
360 "inv(X): inverse of a square matrix\n", }, |
|
|
361 |
|
|
362 { "inverse", 2, 1, builtin_inv, |
|
|
363 "inverse(X): inverse of a square matrix\n", }, |
|
|
364 |
|
|
365 { "isstr", 2, 1, builtin_isstr, |
|
|
366 "isstr(X): return 1 if X is a string\n", }, |
|
|
367 |
|
|
368 { "keyboard", 2, 1, builtin_keyboard, |
|
|
369 "maybe help in debugging M-files\n", }, |
|
|
370 |
|
|
371 { "logm", 2, 1, builtin_logm, |
|
|
372 "Matrix logarithm, log (A)\n", }, |
|
|
373 |
|
|
374 { "lpsolve", 11, 3, builtin_lpsolve, |
|
|
375 "Solve linear programs using lp_solve.\n", }, |
|
|
376 |
|
|
377 { "lsode", 6, 1, builtin_lsode, |
|
|
378 "solve ODEs using Hindmarsh's LSODE. Usage:\n\ |
|
|
379 \n\ |
|
|
380 lsode ('function_name', x0, t_out\n\ |
|
|
381 lsode ('function_name', x0, t_out, t_crit)\n\ |
|
|
382 \n\ |
|
|
383 The first argument is the name of the function to call to\n\ |
|
|
384 compute the vector of right hand sides. It must have the form\n\ |
|
|
385 \n\ |
|
|
386 xdot = f (x, t)\n\ |
|
|
387 \n\ |
|
|
388 where xdot and x are vectors and t is a scalar.\n\n", }, |
|
|
389 |
|
|
390 { "lu", 2, 3, builtin_lu, |
|
|
391 "[L, U, P] = lu (A) -- LU factorization\n", }, |
|
|
392 |
|
|
393 { "max", 3, 2, builtin_max, |
|
|
394 "maximum value(s) of a vector (matrix)\n", }, |
|
|
395 |
|
|
396 { "min", 3, 2, builtin_min, |
|
|
397 "minimum value(s) of a vector (matrix)\n", }, |
|
|
398 |
|
|
399 { "npsol", 11, 3, builtin_npsol, |
|
|
400 "Solve nonlinear programs using Gill and Murray's NPSOL. Usage:\n\ |
|
|
401 \n\ |
|
|
402 [x, obj, info, lambda] = npsol (x, 'phi' [, lb, ub] [, lb, A, ub] [, lb, 'g', ub])\n\n\ |
|
|
403 Groups of arguments surrounded in `[]' are optional, but\n\ |
|
|
404 must appear in the same relative order shown above.\n\ |
|
|
405 \n\ |
|
|
406 The second argument is a string containing the name of the objective\n\ |
|
|
407 function to call. The objective function must be of the form\n\ |
|
|
408 \n\ |
|
|
409 y = phi (x)\n\ |
|
|
410 \n\ |
|
|
411 where x is a vector and y is a scalar.\n\n", }, |
|
|
412 |
|
|
413 { "ones", 3, 1, builtin_ones, |
|
|
414 "create a matrix of all ones\n", }, |
|
|
415 |
|
|
416 { "pause", 1, 0, builtin_pause, |
|
|
417 "suspend program execution\n", }, |
|
|
418 |
|
|
419 { "purge_tmp_files", 5, 1, builtin_purge_tmp_files, |
|
|
420 "delete temporary data files used for plotting\n", }, |
|
|
421 |
|
|
422 { "printf", -1, 1, builtin_printf, |
|
|
423 "printf ('fmt', ...)\n", }, |
|
|
424 |
|
|
425 { "prod", 2, 1, builtin_prod, |
|
|
426 "prod (X): products\n", }, |
|
|
427 |
|
|
428 { "pwd", 1, 0, builtin_pwd, |
|
|
429 "print current working directory\n", }, |
|
|
430 |
|
|
431 { "qpsol", 9, 3, builtin_qpsol, |
|
|
432 "Solve nonlinear programs using Gill and Murray's QPSOL. Usage:\n\ |
|
|
433 \n\ |
|
|
434 [x, obj, info, lambda] = qpsol (x, H, c [, lb, ub] [, lb, A, ub])\n\ |
|
|
435 \n\ |
|
|
436 Groups of arguments surrounded in `[]' are optional, but\n\ |
|
|
437 must appear in the same relative order shown above.", }, |
|
|
438 |
|
|
439 { "qr", 2, 2, builtin_qr, |
|
|
440 "[q, r] = qr (X): form QR factorization of X\n", }, |
|
|
441 |
|
|
442 { "quad", 6, 3, builtin_quad, |
|
|
443 "Integrate a nonlinear function of one variable using Quadpack.\n\ |
|
|
444 Usage:\n\ |
|
|
445 \n\ |
|
|
446 [v, ier, nfun] = quad ('f', a, b [, tol] [, sing])\n\ |
|
|
447 \n\ |
|
|
448 Where the first argument is the name of the function to call to\n\ |
|
|
449 compute the value of the integrand. It must have the form\n\ |
|
|
450 \n\ |
|
|
451 y = f (x) |
|
|
452 \n\ |
|
|
453 where y and x are scalars.\n\ |
|
|
454 \n\ |
|
|
455 The second and third arguments are limits of integration. Either or\n\ |
|
|
456 both may be infinite. The optional argument tol specifies the desired\n\ |
|
|
457 accuracy of the result. The optional argument sing is a vector of\n\ |
|
|
458 at which the integrand is singular.\n\n", }, |
|
|
459 |
|
|
460 { "quit", 1, 0, builtin_quit, |
|
|
461 "exit Octave gracefully\n", }, |
|
|
462 |
|
|
463 { "rand", 2, 1, builtin_rand, |
|
|
464 "matrices with random elements\n", }, |
|
|
465 |
|
|
466 { "replot", 1, 0, builtin_replot, |
|
|
467 "redisplay current plot\n", }, |
|
|
468 |
|
|
469 { "scanf", 2, -1, builtin_scanf, |
|
|
470 "[a,b,c...] = scanf ('fmt')\n", }, |
|
|
471 |
|
|
472 { "setstr", 2, 1, builtin_setstr, |
|
|
473 "setstr (v): convert a vector to a string\n", }, |
|
|
474 |
|
|
475 { "shell_cmd", 2, 1, builtin_shell_command, |
|
|
476 "shell_cmd (string [, return_output]): execute shell commands\n", }, |
|
|
477 |
|
|
478 { "schur", 3, 2, builtin_schur, |
|
|
479 "returns the Schur (straight or ordered) decomposition of matrix\n", }, |
|
|
480 |
|
|
481 { "size", 2, 1, builtin_size, |
|
|
482 "size(X): return rows and columns of X\n", }, |
|
|
483 |
|
|
484 { "sort", 2, 2, builtin_sort, |
|
|
485 "[s,i] = sort(x): sort the columns of x, optionally return sort index\n", }, |
|
|
486 |
|
|
487 { "sqrtm", 2, 1, builtin_sqrtm, |
|
|
488 "Matrix sqrt, sqrt (A)\n", }, |
|
|
489 |
|
|
490 { "sprintf", -1, 1, builtin_sprintf, |
|
|
491 "s = sprintf ('fmt', ...)\n", }, |
|
|
492 |
|
|
493 { "sscanf", 3, -1, builtin_sscanf, |
|
|
494 "[a,b,c...] = sscanf (string, 'fmt')\n", }, |
|
|
495 |
|
|
496 { "sum", 2, 1, builtin_sum, |
|
|
497 "sum (X): sum of elements\n", }, |
|
|
498 |
|
|
499 { "sumsq", 2, 1, builtin_sumsq, |
|
|
500 "sumsq (X): sum of squares of elements\n", }, |
|
|
501 |
|
|
502 { "svd", 2, 3, builtin_svd, |
|
|
503 "[U,S,V] = svd(X): return SVD of X\n", }, |
|
|
504 |
|
|
505 { "warranty", 1, 0, builtin_warranty, |
|
|
506 "describe copying conditions\n", }, |
|
|
507 |
|
|
508 { "zeros", 3, 1, builtin_zeros, |
|
|
509 "create a matrix of all zeros\n", }, |
|
|
510 |
|
|
511 { NULL, -1, -1, NULL, NULL, }, |
|
|
512 }; |
|
|
513 |
|
|
514 // This is a lie. Some of these get reassigned to be numeric |
|
|
515 // variables. See below. |
|
|
516 |
|
|
517 static builtin_string_variables string_variables[] = |
|
|
518 { |
|
|
519 { "I", "??", NULL, |
|
|
520 "sqrt (-1)\n", }, |
|
|
521 |
|
|
522 { "Inf", "??", NULL, |
|
|
523 "infinity\n", }, |
|
|
524 |
|
|
525 { "J", "??", NULL, |
|
|
526 "sqrt (-1)\n", }, |
|
|
527 |
|
|
528 #if defined (HAVE_ISNAN) |
|
|
529 { "NaN", "??", NULL, |
|
|
530 "not a number\n", }, |
|
|
531 #endif |
|
|
532 |
|
|
533 { "LOADPATH", "??", sv_loadpath, |
|
|
534 "colon separated list of directories to search for scripts\n", }, |
|
|
535 |
|
|
536 { "PAGER", "??", sv_pager_binary, |
|
|
537 "path to pager binary\n", }, |
|
|
538 |
|
|
539 { "PS1", "\\s:\\#> ", sv_ps1, |
|
|
540 "primary prompt string\n", }, |
|
|
541 |
|
|
542 { "PS2", "> ", sv_ps2, |
|
|
543 "secondary prompt string\n", }, |
|
|
544 |
|
|
545 { "PWD", "??PWD??", sv_pwd, |
|
|
546 "current working directory\n", }, |
|
|
547 |
|
|
548 { "SEEK_SET", "??", NULL, |
|
|
549 "used with fseek to position file relative to the beginning\n", }, |
|
|
550 |
|
|
551 { "SEEK_CUR", "??", NULL, |
|
|
552 "used with fseek to position file relative to the current position\n", }, |
|
|
553 |
|
|
554 { "SEEK_END", "??", NULL, |
|
|
555 "used with fseek to position file relative to the end\n", }, |
|
|
556 |
|
|
557 { "do_fortran_indexing", "false", do_fortran_indexing, |
|
|
558 "allow single indices for matrices\n", }, |
|
|
559 |
|
|
560 { "empty_list_elements_ok", "warn", empty_list_elements_ok, |
|
|
561 "ignore the empty element in expressions like `a = [[], 1]'\n", }, |
|
|
562 |
|
|
563 { "eps", "??", NULL, |
|
|
564 "machine precision\n", }, |
|
|
565 |
|
|
566 { "gnuplot_binary", "gnuplot", sv_gnuplot_binary, |
|
|
567 "path to gnuplot binary\n", }, |
|
|
568 |
|
|
569 { "i", "??", NULL, |
|
|
570 "sqrt (-1)\n", }, |
|
|
571 |
|
|
572 { "implicit_str_to_num_ok", "false", implicit_str_to_num_ok, |
|
|
573 "allow implicit string to number conversion\n", }, |
|
|
574 |
|
|
575 { "inf", "??", NULL, |
|
|
576 "infinity\n", }, |
|
|
577 |
|
|
578 { "j", "??", NULL, |
|
|
579 "sqrt (-1)\n", }, |
|
|
580 |
|
|
581 #if defined (HAVE_ISNAN) |
|
|
582 { "nan", "??", NULL, |
|
|
583 "not a number\n", }, |
|
|
584 #endif |
|
|
585 |
|
|
586 { "ok_to_lose_imaginary_part", "warn", ok_to_lose_imaginary_part, |
|
|
587 "silently convert from complex to real by dropping imaginary part\n", }, |
|
|
588 |
|
|
589 { "output_max_field_width", "??", set_output_max_field_width, |
|
|
590 "maximum width of an output field for numeric output\n", }, |
|
|
591 |
|
|
592 { "output_precision", "??", set_output_precision, |
|
|
593 "number of significant figures to display for numeric output\n", }, |
|
|
594 |
|
|
595 { "page_screen_output", "true", page_screen_output, |
|
|
596 "if possible, send output intended for the screen through the pager\n", }, |
|
|
597 |
|
|
598 { "pi", "??", NULL, |
|
|
599 "ratio of the circumference of a circle to its diameter\n", }, |
|
|
600 |
|
|
601 { "prefer_column_vectors", "true", prefer_column_vectors, |
|
|
602 "prefer column/row vectors\n", }, |
|
|
603 |
|
|
604 { "prefer_zero_one_indexing", "false", prefer_zero_one_indexing, |
|
|
605 "when there is a conflict, prefer zero-one style indexing\n", }, |
|
|
606 |
|
|
607 { "print_answer_id_name", "true", print_answer_id_name, |
|
|
608 "set output style to print `var_name = ...'\n", }, |
|
|
609 |
|
|
610 { "print_empty_dimensions", "true", print_empty_dimensions, |
|
|
611 "also print dimensions of empty matrices\n", }, |
|
|
612 |
|
|
613 { "propagate_empty_matrices", "true", propagate_empty_matrices, |
|
|
614 "operations on empty matrices return an empty matrix, not an error\n", }, |
|
|
615 |
|
|
616 { "resize_on_range_error", "true", resize_on_range_error, |
|
|
617 "enlarge matrices on assignment\n", }, |
|
|
618 |
|
|
619 { "return_last_computed_value", "false", return_last_computed_value, |
|
|
620 "if a function does not return any values explicitly, return the\n\ |
|
|
621 last computed value\n", }, |
|
|
622 |
|
|
623 { "silent_functions", "false", silent_functions, |
|
|
624 "suppress printing results in called functions\n", }, |
|
|
625 |
|
|
626 { "split_long_rows", "true", split_long_rows, |
|
|
627 "split long matrix rows instead of wrapping\n", }, |
|
|
628 |
|
|
629 { "stdin", "??", NULL, |
|
|
630 "file number of the standard input stream\n", }, |
|
|
631 |
|
|
632 { "stdout", "??", NULL, |
|
|
633 "file number of the standard output stream\n", }, |
|
|
634 |
|
|
635 { "stderr", "??", NULL, |
|
|
636 "file number of the standard error stream\n", }, |
|
|
637 |
|
|
638 { "treat_neg_dim_as_zero", "false", treat_neg_dim_as_zero, |
|
|
639 "convert negative dimensions to zero\n", }, |
|
|
640 |
|
|
641 { "warn_assign_as_truth_value", "true", warn_assign_as_truth_value, |
|
|
642 "produce warning for assignments used as truth values\n", }, |
|
|
643 |
|
|
644 { "warn_comma_in_global_decl", "true", warn_comma_in_global_decl, |
|
|
645 "produce warning for commas in global declarations\n", }, |
|
|
646 |
|
|
647 { "warn_divide_by_zero", "true", warn_divide_by_zero, |
|
|
648 "on IEEE machines, allow divide by zero errors to be suppressed\n", }, |
|
|
649 |
|
|
650 { NULL, NULL, NULL, NULL, }, |
|
|
651 }; |
|
|
652 |
|
|
653 void |
|
|
654 make_eternal (char *s) |
|
|
655 { |
|
|
656 symbol_record *sym_rec = curr_sym_tab->lookup (s, 0, 0); |
|
|
657 if (sym_rec != (symbol_record *) NULL) |
|
|
658 sym_rec->make_eternal (); |
|
|
659 } |
|
|
660 |
|
|
661 void |
|
|
662 install_builtins (void) |
|
|
663 { |
|
|
664 symbol_record *sym_rec; |
|
|
665 |
|
|
666 tree_builtin *tb_tmp; |
|
|
667 |
|
|
668 // So that the clear function can't delete other builtin variables and |
|
|
669 // functions, they are given eternal life. |
|
|
670 |
|
|
671 builtin_mapper_functions *mfptr = mapper_functions; |
|
|
672 while (mfptr->name != (char *) NULL) |
|
|
673 { |
|
|
674 sym_rec = curr_sym_tab->lookup (mfptr->name, 1); |
|
|
675 sym_rec->unprotect (); |
|
|
676 |
|
|
677 Mapper_fcn mfcn; |
|
|
678 mfcn.neg_arg_complex = mfptr->neg_arg_complex; |
|
|
679 mfcn.d_d_mapper = mfptr->d_d_mapper; |
|
|
680 mfcn.d_c_mapper = mfptr->d_c_mapper; |
|
|
681 mfcn.c_c_mapper = mfptr->c_c_mapper; |
|
|
682 |
|
|
683 tb_tmp = new tree_builtin (mfptr->nargin_max, mfptr->nargout_max, |
|
|
684 mfcn, sym_rec); |
|
|
685 |
|
|
686 sym_rec->define (tb_tmp); |
|
|
687 sym_rec->document (mfptr->help_string); |
|
|
688 sym_rec->make_eternal (); |
|
|
689 sym_rec->protect (); |
|
|
690 mfptr++; |
|
|
691 } |
|
|
692 |
|
|
693 builtin_text_functions *tfptr = text_functions; |
|
|
694 while (tfptr->name != (char *) NULL) |
|
|
695 { |
|
|
696 sym_rec = curr_sym_tab->lookup (tfptr->name, 1); |
|
|
697 sym_rec->unprotect (); |
|
|
698 |
|
|
699 tb_tmp = new tree_builtin (tfptr->nargin_max, 1, |
|
|
700 tfptr->text_fcn, sym_rec); |
|
|
701 |
|
|
702 sym_rec->define (tb_tmp); |
|
|
703 sym_rec->document (tfptr->help_string); |
|
|
704 sym_rec->make_eternal (); |
|
|
705 sym_rec->protect (); |
|
|
706 tfptr++; |
|
|
707 } |
|
|
708 |
|
|
709 builtin_general_functions *gfptr = general_functions; |
|
|
710 while (gfptr->name != (char *) NULL) |
|
|
711 { |
|
|
712 sym_rec = curr_sym_tab->lookup (gfptr->name, 1); |
|
|
713 sym_rec->unprotect (); |
|
|
714 |
|
|
715 tb_tmp = new tree_builtin (gfptr->nargin_max, gfptr->nargout_max, |
|
|
716 gfptr->general_fcn, sym_rec); |
|
|
717 |
|
|
718 sym_rec->define (tb_tmp); |
|
|
719 sym_rec->document (gfptr->help_string); |
|
|
720 sym_rec->make_eternal (); |
|
|
721 sym_rec->protect (); |
|
|
722 gfptr++; |
|
|
723 } |
|
|
724 |
|
|
725 // Most built-in variables are not protected because the user should |
|
|
726 // be able to redefine them. |
|
|
727 |
|
|
728 builtin_string_variables *svptr = string_variables; |
|
|
729 while (svptr->name != (char *) NULL) |
|
|
730 { |
|
|
731 sym_rec = curr_sym_tab->lookup (svptr->name, 1); |
|
|
732 sym_rec->unprotect (); |
|
|
733 |
|
|
734 tree_constant *tmp = new tree_constant (svptr->value); |
|
|
735 |
|
|
736 sym_rec->set_sv_function (svptr->sv_function); |
|
|
737 sym_rec->define (tmp); |
|
|
738 sym_rec->document (svptr->help_string); |
|
|
739 sym_rec->make_eternal (); |
|
|
740 svptr++; |
|
|
741 } |
|
|
742 |
|
|
743 // XXX FIXME XXX -- Need a convenient way to document these variables. |
|
|
744 |
|
|
745 // IMPORTANT: Always create a new tree_constant for each variable. |
|
|
746 |
|
|
747 tree_constant *tmp = NULL_TREE_CONST; |
|
|
748 bind_variable ("ans", tmp); |
|
|
749 |
|
|
750 Complex ctmp (0.0, 1.0); |
|
|
751 tmp = new tree_constant (ctmp); |
|
|
752 bind_protected_variable ("I", tmp); |
|
|
753 make_eternal ("I"); |
|
|
754 |
|
|
755 tmp = new tree_constant (ctmp); |
|
|
756 bind_protected_variable ("J", tmp); |
|
|
757 make_eternal ("J"); |
|
|
758 |
|
|
759 // Let i and j be functions so they can be redefined without being |
|
|
760 // wiped out. |
|
|
761 |
|
|
762 char *tmp_help; |
|
|
763 |
|
|
764 tmp = new tree_constant (ctmp); |
|
|
765 sym_rec = curr_sym_tab->lookup ("i", 1); |
|
|
766 tmp_help = sym_rec->help (); |
|
|
767 sym_rec->define_as_fcn (tmp); |
|
|
768 sym_rec->document (tmp_help); |
|
|
769 sym_rec->protect (); |
|
|
770 sym_rec->make_eternal (); |
|
|
771 |
|
|
772 tmp = new tree_constant (ctmp); |
|
|
773 sym_rec = curr_sym_tab->lookup ("j", 1); |
|
|
774 tmp_help = sym_rec->help (); |
|
|
775 sym_rec->define_as_fcn (tmp); |
|
|
776 sym_rec->document (tmp_help); |
|
|
777 sym_rec->protect (); |
|
|
778 sym_rec->make_eternal (); |
|
|
779 |
|
|
780 tmp = new tree_constant (get_working_directory ("initialize_globals")); |
|
|
781 bind_protected_variable ("PWD", tmp); |
|
|
782 make_eternal ("PWD"); |
|
|
783 |
|
|
784 tmp = new tree_constant (load_path); |
|
|
785 bind_variable ("LOADPATH", tmp); |
|
|
786 make_eternal ("LOADPATH"); |
|
|
787 |
|
|
788 tmp = new tree_constant (default_pager ()); |
|
|
789 bind_variable ("PAGER", tmp); |
|
|
790 make_eternal ("PAGER"); |
|
|
791 |
|
|
792 tmp = new tree_constant (0.0); |
|
|
793 bind_variable ("SEEK_SET", tmp); |
|
|
794 make_eternal ("SEEK_SET"); |
|
|
795 |
|
|
796 tmp = new tree_constant (1.0); |
|
|
797 bind_variable ("SEEK_CUR", tmp); |
|
|
798 make_eternal ("SEEK_CUR"); |
|
|
799 |
|
|
800 tmp = new tree_constant (2.0); |
|
|
801 bind_variable ("SEEK_END", tmp); |
|
|
802 make_eternal ("SEEK_END"); |
|
|
803 |
|
|
804 tmp = new tree_constant (DBL_EPSILON); |
|
|
805 bind_protected_variable ("eps", tmp); |
|
|
806 make_eternal ("eps"); |
|
|
807 |
|
|
808 tmp = new tree_constant (10.0); |
|
|
809 bind_variable ("output_max_field_width", tmp); |
|
|
810 make_eternal ("output_max_field_width"); |
|
|
811 |
|
|
812 tmp = new tree_constant (5.0); |
|
|
813 bind_variable ("output_precision", tmp); |
|
|
814 make_eternal ("output_precision"); |
|
|
815 |
|
|
816 tmp = new tree_constant (4.0 * atan (1.0)); |
|
|
817 bind_protected_variable ("pi", tmp); |
|
|
818 make_eternal ("pi"); |
|
|
819 |
|
|
820 tmp = new tree_constant (0.0); |
|
|
821 bind_protected_variable ("stdin", tmp); |
|
|
822 make_eternal ("stdin"); |
|
|
823 |
|
|
824 tmp = new tree_constant (1.0); |
|
|
825 bind_protected_variable ("stdout", tmp); |
|
|
826 make_eternal ("stdout"); |
|
|
827 |
|
|
828 tmp = new tree_constant (2.0); |
|
|
829 bind_protected_variable ("stderr", tmp); |
|
|
830 make_eternal ("stderr"); |
|
|
831 |
|
|
832 #if defined (HAVE_ISINF) || defined (HAVE_FINITE) |
|
|
833 #ifdef linux |
|
|
834 tmp = new tree_constant (HUGE_VAL); |
|
|
835 #else |
|
|
836 tmp = new tree_constant (1.0/0.0); |
|
|
837 #endif |
|
|
838 bind_protected_variable ("Inf", tmp); |
|
|
839 make_eternal ("Inf"); |
|
|
840 |
|
|
841 #ifdef linux |
|
|
842 tmp = new tree_constant (HUGE_VAL); |
|
|
843 #else |
|
|
844 tmp = new tree_constant (1.0/0.0); |
|
|
845 #endif |
|
|
846 bind_protected_variable ("inf", tmp); |
|
|
847 make_eternal ("inf"); |
|
|
848 |
|
|
849 #else |
|
|
850 |
|
|
851 // This is sort of cheesy, but what can we do, other than blowing it |
|
|
852 // off completely, or writing an entire IEEE emulation package? |
|
|
853 |
|
|
854 tmp = new tree_constant (DBL_MAX); |
|
|
855 bind_protected_variable ("Inf", tmp); |
|
|
856 make_eternal ("Inf"); |
|
|
857 |
|
|
858 tmp = new tree_constant (DBL_MAX); |
|
|
859 bind_protected_variable ("inf", tmp); |
|
|
860 make_eternal ("inf"); |
|
|
861 #endif |
|
|
862 |
|
|
863 #if defined (HAVE_ISNAN) |
|
|
864 #ifdef linux |
|
|
865 tmp = new tree_constant (NAN); |
|
|
866 #else |
|
|
867 tmp = new tree_constant (0.0/0.0); |
|
|
868 #endif |
|
|
869 bind_protected_variable ("NaN", tmp); |
|
|
870 make_eternal ("NaN"); |
|
|
871 |
|
|
872 #ifdef linux |
|
|
873 tmp = new tree_constant (NAN); |
|
|
874 #else |
|
|
875 tmp = new tree_constant (0.0/0.0); |
|
|
876 #endif |
|
|
877 bind_protected_variable ("nan", tmp); |
|
|
878 make_eternal ("nan"); |
|
|
879 #endif |
|
|
880 } |
|
|
881 |
|
|
882 int |
|
|
883 is_text_function_name (char *s) |
|
|
884 { |
|
|
885 int retval = 0; |
|
|
886 |
|
|
887 builtin_text_functions *tfptr = text_functions; |
|
|
888 while (tfptr->name != (char *) NULL) |
|
|
889 { |
|
|
890 if (strcmp (tfptr->name, s) == 0) |
|
|
891 { |
|
|
892 retval = 1; |
|
|
893 break; |
|
|
894 } |
|
|
895 tfptr++; |
|
|
896 } |
|
|
897 |
|
|
898 return retval; |
|
|
899 } |
|
|
900 |
|
|
901 help_list * |
|
|
902 builtin_mapper_functions_help (void) |
|
|
903 { |
|
|
904 int count = 0; |
|
|
905 builtin_mapper_functions *mfptr; |
|
|
906 |
|
|
907 mfptr = mapper_functions; |
|
|
908 while (mfptr->name != (char *) NULL) |
|
|
909 { |
|
|
910 count++; |
|
|
911 mfptr++; |
|
|
912 } |
|
|
913 |
|
|
914 if (count == 0) |
|
|
915 return (help_list *) NULL; |
|
|
916 |
|
|
917 help_list *hl = new help_list [count+1]; |
|
|
918 |
|
|
919 int i = 0; |
|
|
920 mfptr = mapper_functions; |
|
|
921 while (mfptr->name != (char *) NULL) |
|
|
922 { |
|
|
923 hl[i].name = mfptr->name; |
|
|
924 hl[i].help = mfptr->help_string; |
|
|
925 i++; |
|
|
926 mfptr++; |
|
|
927 } |
|
|
928 |
|
|
929 hl[count].name = (char *) NULL; |
|
|
930 hl[count].help = (char *) NULL; |
|
|
931 |
|
|
932 return hl; |
|
|
933 } |
|
|
934 |
|
|
935 help_list * |
|
|
936 builtin_general_functions_help (void) |
|
|
937 { |
|
|
938 int count = 0; |
|
|
939 builtin_general_functions *gfptr; |
|
|
940 |
|
|
941 gfptr = general_functions; |
|
|
942 while (gfptr->name != (char *) NULL) |
|
|
943 { |
|
|
944 count++; |
|
|
945 gfptr++; |
|
|
946 } |
|
|
947 |
|
|
948 if (count == 0) |
|
|
949 return (help_list *) NULL; |
|
|
950 |
|
|
951 help_list *hl = new help_list [count+1]; |
|
|
952 |
|
|
953 int i = 0; |
|
|
954 gfptr = general_functions; |
|
|
955 while (gfptr->name != (char *) NULL) |
|
|
956 { |
|
|
957 hl[i].name = gfptr->name; |
|
|
958 hl[i].help = gfptr->help_string; |
|
|
959 i++; |
|
|
960 gfptr++; |
|
|
961 } |
|
|
962 |
|
|
963 hl[count].name = (char *) NULL; |
|
|
964 hl[count].help = (char *) NULL; |
|
|
965 |
|
|
966 return hl; |
|
|
967 } |
|
|
968 |
|
|
969 help_list * |
|
|
970 builtin_text_functions_help (void) |
|
|
971 { |
|
|
972 int count = 0; |
|
|
973 builtin_text_functions *tfptr; |
|
|
974 |
|
|
975 tfptr = text_functions; |
|
|
976 while (tfptr->name != (char *) NULL) |
|
|
977 { |
|
|
978 count++; |
|
|
979 tfptr++; |
|
|
980 } |
|
|
981 |
|
|
982 if (count == 0) |
|
|
983 return (help_list *) NULL; |
|
|
984 |
|
|
985 help_list *hl = new help_list [count+1]; |
|
|
986 |
|
|
987 int i = 0; |
|
|
988 tfptr = text_functions; |
|
|
989 while (tfptr->name != (char *) NULL) |
|
|
990 { |
|
|
991 hl[i].name = tfptr->name; |
|
|
992 hl[i].help = tfptr->help_string; |
|
|
993 i++; |
|
|
994 tfptr++; |
|
|
995 } |
|
|
996 |
|
|
997 hl[count].name = (char *) NULL; |
|
|
998 hl[count].help = (char *) NULL; |
|
|
999 |
|
|
1000 return hl; |
|
|
1001 } |
|
|
1002 |
|
|
1003 help_list * |
|
|
1004 builtin_variables_help (void) |
|
|
1005 { |
|
|
1006 int count = 0; |
|
|
1007 builtin_string_variables *svptr; |
|
|
1008 |
|
|
1009 svptr = string_variables; |
|
|
1010 while (svptr->name != (char *) NULL) |
|
|
1011 { |
|
|
1012 count++; |
|
|
1013 svptr++; |
|
|
1014 } |
|
|
1015 |
|
|
1016 if (count == 0) |
|
|
1017 return (help_list *) NULL; |
|
|
1018 |
|
|
1019 help_list *hl = new help_list [count+1]; |
|
|
1020 |
|
|
1021 int i = 0; |
|
|
1022 svptr = string_variables; |
|
|
1023 while (svptr->name != (char *) NULL) |
|
|
1024 { |
|
|
1025 hl[i].name = svptr->name; |
|
|
1026 hl[i].help = svptr->help_string; |
|
|
1027 i++; |
|
|
1028 svptr++; |
|
|
1029 } |
|
|
1030 |
|
|
1031 hl[count].name = (char *) NULL; |
|
|
1032 hl[count].help = (char *) NULL; |
|
|
1033 |
|
|
1034 return hl; |
|
|
1035 } |
|
|
1036 |
|
|
1037 /* |
|
|
1038 ;;; Local Variables: *** |
|
|
1039 ;;; mode: C++ *** |
|
|
1040 ;;; page-delimiter: "^/\\*" *** |
|
|
1041 ;;; End: *** |
|
|
1042 */ |
