UserGuide10 PDF
UserGuide10 PDF
GAUSS
User Guide
c
Copyright Aptech Systems, Inc. Black Diamond WA 1984-2010
All Rights Reserved Worldwide.
c
SuperLU.
Copyright 2003, The Regents of the University of California, through
Lawrence Berkeley National Laboratory (subject to receipt of any required approvals
from U.S. Dept. of Energy). All Rights Reserved. See GAUSS Software Product
License for additional terms and conditions.
c
TAUCS Version 2.0, November 29, 2001.
Copyright 2001, 2002, 2003 by Sivan
Toledo, Tel-Aviv University, [email protected]. All Rights Reserved. See GAUSS
Software License for additional terms and conditions.
GAUSS, GAUSS Engine and GAUSS Light are trademarks of Aptech Systems, Inc.
GEM is a trademark of Digital Research, Inc.
Lotus is a trademark of Lotus Development Corp.
HP LaserJet and HP-GL are trademarks of Hewlett-Packard Corp.
PostScript is a trademark of Adobe Systems Inc.
IBM is a trademark of International Business Machines Corporation
Hercules is a trademark of Hercules Computer Technology, Inc.
GraphiC is a trademark of Scientific Endeavors Corporation
Tektronix is a trademark of Tektronix, Inc.
Windows is a registered trademark of Microsoft Corporation.
Other trademarks are the property of their respective owners.
Part Number: 006903
Version 10
Documentation Revision: 992 May 26, 2010
Contents
Contents
1 Introduction
2 Getting Started
v
GAUSS User Guide
vi
Contents
6 Matrix Editor
7 Library Tool
vii
GAUSS User Guide
9 GAUSS Help
10 Language Fundamentals
viii
Contents
11 Operators
ix
GAUSS User Guide
13 Sparse Matrices
14 N-Dimensional Arrays
x
Contents
16 Structures
xi
GAUSS User Guide
xii
Contents
19 Libraries
20 Compiler
21 File I/O
xiii
GAUSS User Guide
23 Data Transformations
xiv
Contents
26 Graphics Editor
xv
GAUSS User Guide
28 ATOG
29 Error Messages
30 Maximizing Performance
xvi
Contents
A Fonts
32 Commands by Category
xvii
GAUSS User Guide
33 Command Reference
D Obsolete Commands
E Colors
Index
xviii
List of Figures
List of Figures
4.1 GAUSS Graphical User Interface . . . . . . . . . . . . . . . . . . . . . . . . 4-2
4.2 Main Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-11
4.3 Working Directory Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-12
4.4 Debug Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-13
4.5 Window Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-14
4.6 Status Bar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-15
6.1 Matrix Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-2
7.1 Library Tool . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7-2
8.1 Source Browser . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8-3
16.1 Structure tree for e1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16-7
26.1 Graphics Editor Workspace . . . . . . . . . . . . . . . . . . . . . . . . . . . 26-2
26.2 Graphics Editor Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26-3
26.3 Graphics Editor Status Bar . . . . . . . . . . . . . . . . . . . . . . . . . . . 26-3
xix
Introduction
Introduction 1
TM
GAUSS is a complete analysis environment suitable for performing quick calculations, complex
analysis of millions of data points, or anything in between. Whether you are new to computerized
analysis or a seasoned programmer, the GAUSS family of products combine to offer you an easy
to learn environment that is powerful and versatile enough for virtually any numerical task.
Since its introduction in 1984, GAUSS has been the standard for serious number crunching and
complex modeling of large-scale data. Worldwide acceptance and use in government, industry,
and the academic community is a firm testament to its power and versatility.
The GAUSS System can be described several ways: It is an exceptionally efficient number
cruncher, a comprehensive programming language, and an interactive analysis environment.
GAUSS may be the only numerical tool you will ever need.
1-1
GAUSS User Guide
The following table describes how text formatting is used to identify GAUSS programming
elements:
1-2
Started
Getting
Getting Started 2
2.1 Installation Under UNIX/Linux
2. cd to that directory.
8. Install the license. (To receive a license and license installation instructions, email
[email protected].)
2-1
GAUSS User Guide
• Minimum of 100 MB free hard disk space, more may be needed depending on the size of
matrices and the complexity of the program.
• Monthly defragmenting is recommended.
Insert the GAUSS compact disc into the CD-ROM drive, and setup should start automatically. If
setup does not start automatically, click Start, then click Run. Type D:\setup.exe in the dialog
box (where D is the drive letter of the CD-ROM drive).
You can use this procedure for the initial installation of GAUSS, and for additions or
modifications to GAUSS components.
2-2
Using the Command Line
Interface 3
Line
Command
TGAUSS is the command line version of GAUSS. The executable file, tgauss is located in the
GAUSS installation directory.
-b Execute file in batch mode and then exit. You can execute multiple files by
separating file names with spaces.
-l logfile Set the name of the batch mode log file when using the -b argument. The
default is tmp/gauss.log###, where ### is the process ID.
-e expression Execute a GAUSS expression. This command is not logged when GAUSS is in
batch mode.
-o Suppress the sign-on banner (output only).
-T Turn the dataloop translator on.
-t Turn the dataloop translator off.
3-1
GAUSS User Guide
GAUSS generates .tkf files for graphical output. The default output for graphics is
graphic.tkf. On Windows, you can use vwr.exe to view the graphics file; on
UNIX/Linux/Mac, you can use vwrmp. Two functions are available to convert .tkf files to
PostScript for printing and viewing with external viewers: the tkf2ps function will convert .tkf
files to PostScript (.ps) files, and the tkf2eps function will convert .tkf files to encapsulated
PostScript (.eps) files. For example, to convert the file graphic.tkf to a postscript file named
graphic.ps use:
When you run a command at the TGAUSS prompt, it is added to your command line history,
which is stored in a file called .gauss_prompt_history in your $(HOME) directory on
UNIX/Linux or in your $(HOMEDRIVE)\$(HOMEPATH) directory on Windows. A separate history
for commands entered in the command line debugger is stored in a file called
.gauss_debug_prompt_history in the same directory. By default, the last 500 commands
executed at the TGAUSS and debugger command lines are stored in these files. You can change
this number by changing prompt_hist_num in your gauss.cfg file. The following keystrokes
are supported for movement and editing at the command line and for retrieving the command line
history:
3.2.1 Movement
3-2
Using the Command Line Interface
Line
Command
CTRL+RIGHT ARROW
3.2.2 Editing
3-3
GAUSS User Guide
PAGE DOWN or CTRL+S Retrieves next line in history that matches text to left
of cursor
Note that some of these keystrokes are mapped differently on different computers. For example,
on some computers, SHIFT+RIGHT ARROW behaves the same as RIGHT ARROW, while
ALT+RIGHT ARROW moves the cursor right one word. Therefore, multiple keystroke mappings
have been supported to maximize the availability of these commands on any given machine.
3-4
Using the Command Line Interface
3.3.1 quit
quit
Line
Command
You can also use the system command to exit TGAUSS from either the command line or a
program (see system in the GAUSS L R).
system
3.3.2 ed
The ed command will open an input file in an external text editor (see ed in the GAUSS L
R).
ed filename
3.3.3 browse
The browse command allows you to search for specific symbols in a file and open the file in the
default editor. You can use wildcards to extend search capabilities of the browse command.
3-5
GAUSS User Guide
browse symbol
3.3.4 config
The config command gives you access to the configuration menu allowing you to change the way
GAUSS runs and compiles files.
config
Run Menu
Translator Toggles on/off the translation of a file using dataloop. The translator
is not necessary for GAUSS program files not using dataloop.
Translator line Toggles on/off execution time line number tracking of the original
number tracking file before translation.
Line number Toggles on/off the execution time line number tracking. If the
tracking translator is on, the line numbers refer to the translated file.
Compile Menu
3-6
Using the Command Line Interface
3.4 Debugging
The debug command runs a program under the source level debugger.
Line
Command
The format for debug is:
debug filename
l number Displays a specified number of lines of source code in the current file.
lc Displays source code in the current file starting with the current line.
ll file line Displays source code in the named file starting with the specified line.
ll file Displays source code in the named file starting with the first line.
ll line Displays source code starting with the specified line. File does not change.
ll Displays the next page of source code.
lp Displays the previous page of source code.
3-7
GAUSS User Guide
x number Executes code from the beginning of the program to the specified line count, or
until a breakpoint is hit.
g [[args]] Executes from the current line to the end of the program, stopping at
breakpoints. The optional arguments specify other stopping points. The syntax
for each optional argument is:
filename line period The debugger will stop every period times it reaches the
specified line in the named file.
filename line The debugger will stop when it reaches the specified line
in the named file.
filename ,, period The debugger will stop every period times it reaches any
line in the named file.
line period The debugger will stop every period times it reaches the
specified line in the current file.
filename The debugger will stop at every line in the named file.
line The debugger will stop when it reaches the specified line
in the current file.
procedure period The debugger will stop every period times it reaches the
first line in a called procedure.
procedure The debugger will stop every time it reaches the first line
in a called procedure.
j [[args]] Executes code to a specified line, procedure, or period in the file without
stopping at breakpoints. The optional arguments are the same as g, listed above.
jx number Executes code to the execution count specified (number) without stopping at
breakpoints.
o Executes the remainder of the current procedure (or to a breakpoint) and stops
at the next line in the calling procedure.
3-8
Using the Command Line Interface
v [[vars]] Searches for (a local variable, then a global variable) and displays the value of a
specified variable.
v$ [[vars]] Searches for (a local variable, then a global variable) and displays the specified
character matrix.
The display properties of matrices and string arrays can be set using the following commands.
Line
Command
r Specifies the number of rows to be shown.
num,num Specifies the indices of the upper left corner of the block to be shown.
b [[args]] Sets a breakpoint in the code. The syntax for each optional argument is:
filename line period The debugger will stop every period times it reaches the
specified line in the named file.
filename line The debugger will stop when it reaches the specified line
in the named file.
filename ,, period The debugger will stop every period times it reaches any
line in the named file.
3-9
GAUSS User Guide
line period The debugger will stop every period times it reaches the
specified line in the current file.
filename The debugger will stop at every line in the named file.
line The debugger will stop when it reaches the specified line
in the current file.
procedure period The debugger will stop every period times it reaches the
first line in a called procedure.
procedure The debugger will stop every time it reaches the first line
in a called procedure.
d [[args]] Removes a previously specified breakpoint. The optional arguments are the
same arguments as b, listed above.
3-10
Introduction to the Windows
Interface 4
GUI Intro
Windows
The GAUSS graphical user interface is a multiple document interface. The interface consists of
the Menu Bar, the Toolbar, edit windows, the Command Input-Output window, and the Status bar
(see Figure 4.1).
You can view the commands on a menu by either clicking the menu name or pressing ALT+n,
where n is the underlined letter in the menu name. For example, to display the File menu, you can
either click File or press ALT+F.
The File menu lets you access the file, printer setup, and exit commands. Some of these actions
can also be executed from the toolbar. The File menu contains the following commands:
4-1
GAUSS User Guide
Toolbar
Source
View
Status Bar
4-2
Introduction to the Windows Interface
GUI Intro
Windows
Insert File Opens an existing text file and copies the contents into the active document.
This is similar to pasting text from the Windows clipboard.
Print Prints the active file or selected text from the active window.
Print Setup Specifies the printer you want to use. Other printer options, such as page
orientation and paper tray, are also accessed with this command.
Properties Displays information about the active file.
Change Working Changes the directory where GAUSS looks for the files it uses for normal
Directory operation. This command does not affect the Open or Save As paths.
Clear Working Clears the working directory list.
Directory List
Exit Closes all open files and exits GAUSS. You are prompted to save any file
that has been modified since it was last saved.
Recent Files GAUSS maintains a list of the ten most recent files you opened, at the end of
the File menu. If the file you want to open is on this list, click it and GAUSS
opens it in an Edit window.
4-3
GAUSS User Guide
The Edit menu lets you access the set of editing commands. Some of these actions can also be
executed from the toolbar. The Edit menu contains the following commands:
Redo Restores changes in the active window that you removed using the Undo
Edit command.
Cut Removes selected text from the active window and places it on the Windows
clipboard.
Copy Copies selected text from the active window to the Windows clipboard.
Paste Copies text from the Windows clipboard to the active window at the cursor
position.
Find Finds the specified text in the active window. The search starts at the cursor
position and continues to the end of the text in the active window. The
search can be case sensitive or case insensitive. You may also limit the
search to regular expressions.
Find Again Resumes the search for the next occurrence of the text you specified in the
previous Find action. Subsequent searches for the same text can also be
performed by pressing F3.
Replace Locates the specified text in the active window and replaces it with the text
you entered in the “Replace with” field in the Search dialog box. The search
starts at the cursor position and continues to the end of the text in the active
window. The search can be case sensitive or case insensitive, and the
replacement can be unique or global.
Insert Time/Date Inserts the current time and date at the cursor position. GAUSS uses the
time and date that appears in the Microsoft Windows Date/Time Properties
window.
4-4
Introduction to the Windows Interface
GUI Intro
Windows
The View menu lets you toggle the Main Toolbar, the Status Bar, the Working Directory Toolbar,
and the Debug Toolbar on or off.
Main Toolbar Toggles the Main toolbar on or off. For more information about the Main
toolbar, see M T, Section 4.1.11.
Status Bar The Status Bar is located along the bottom of the GAUSS window. For more
information about the status bar, see S B, Section 4.1.15.
Working Toggles the Working Directory toolbar on or off. For more information
Directory about the working directory toolbar, see W D T,
Section 4.1.12.
Debug Toolbar Toggles the Debug toolbar on or off. For more information about the Debug
toolbar, see D T, Section 4.1.13.
Window Toolbar Toggles the Window toolbar on or off. For more information about the
Window toolbar, see W T, Section 4.1.14.
Error Output Opens or closes the Error Output window.
Source View Displays or undisplays the Source/Symbols window.
4-5
GAUSS User Guide
Preferences Opens the General Preferences window. From the General Preferences
window you can define Run options, Compile options, DOS window
options, and Autosave options. For more information on configuring
GAUSS General Preferences, see P D B, Section 5.6.1.
Editor Opens the Editor Properties window. From the Editor Properties window
Properties you can define colors and fonts, the language syntax, tabs, or general editor
properties. For more information on configuring editor properties, see
E P, Section 5.6.2.
The Run menu lets you run the code you have entered, a block of code you selected, or the active
file, depending on the operating mode.
Insert GAUSS Manually adds the GAUSS prompt (>>) at the cursor position. The GAUSS
Prompt prompt is automatically displayed following the execution of GAUSS code.
Insert Last Cmd Re-enters the last command written to the Input buffer.
Run Selected Runs any text selected from the editor or the Command Input-Output
Text window.
Run Active File Runs the active file. The file then becomes the main file.
Test Compile Compiles the currently selected file. During compilation, any errors are
Active File displayed in the Output window.
Note: This command is different than the GAUSS compile command, which
compiles a program and saves the pseudocode as a file.
Run Main File Runs the file specified in the Main File list.
Test Compile Test compiles the main file. During compilation, any errors are displayed in
4-6
Introduction to the Windows Interface
Edit Main File Opens the specified main file in an edit window.
Stop Program Stops the program currently running and returns control to the editor.
Build GCG File Creates GAUSS pseudocode file that can be run over and over with no
from Main compile time.
Set Main File Makes the active file the main file.
Clear Main Removes all entries in the Main File list on the Main toolbar.
File List
Translate Toggles translate dataloop command on and off. For more information see
Dataloop Cmds D T, Chapter 23.
GUI Intro
Windows
4.1.6 Debug Menu
The Debug menu lets you access the commands used to debug your active file or main file.
Edit Opens a list of all breakpoints in your program. The breakpoints are listed
Breakpoints by line number. Any procedure breakpoints are also listed.
Clear All Removes all line and procedure breakpoints from the active file.
Breakpoints
4-7
GAUSS User Guide
Run to Cursor Runs the program until it reaches the cursor position.
Step Into Runs the next executable line of code in the application and steps into
procedures.
Step Over Runs the next executable line of code in the application but does not step
into procedures.
Step Out Runs the remainder of the current procedure and stops at the next line in the
calling procedure. Step Out returns if a breakpoint is encountered.
Set Watch Opens the Matrix Editor for watching changing variable data. For more
information about viewing variables see V V, Section 5.5.5.
The Tools menu lets you open GAUSS tools windows. The following commands can be used:
Matrix Editor Lets you create or edit data in a matrix (or grid). A cell can be edited by
typing in a new value and pressing ENTER. For more information see
M E, Chapter 6.
Graphics Editor Opens the Graphics Editor, which is an interactive TKF file editor. This
menu item will be inactive if you have not purchased the Graphics Editor.
Source Browser Searches source files for string patterns. For more information see GAUSS
S B, Chapter 8.
Lib Tool Lets you manage the contents of libraries. For more information see L
T, Chapter 7.
4-8
Introduction to the Windows Interface
The Window menu commands let you manage your workspace. You can toggle the focus between
all open windows using CTRL+TAB, or clicking in the window you want active. All open
windows are listed at the end of the Window menu. The following commands can be used:
Cmd Window Makes the Command Input - Output window the active window.
Output Window Splits the output from the Command Input - Output window.
Debug Window Starts the debugger on the current file.
Re-use Window If checked, the next file browsed in the Source Browser will be displayed in
the same window.
Command Log Loads the command log window into an editor.
Close All Closes all TKF File Viewer windows.
Graphics
GUI Intro
Windows
Dual Horizontal Horizontally tiles the program source and execution windows within the
main window, and minimizes all other windows.
Dual Vertical Vertically tiles the program source and execution windows within the main
window, and minimizes all other windows.
Cascade Arranges all open windows on the screen, overlapping each, with the active
window on top.
Tile Horizontal Arranges all open windows horizontally on the screen without any overlap.
Tile Vertical Arranges all open windows vertically on the screen without any overlap.
Arrange Icons Arranges all minimized windows across the bottom of the main GAUSS
window.
Split Splits the active window into two horizontal panes. This allows you to view
Horizontally two different areas of the same document to facilitate split-window editing.
Note: You can move the splitter bar by dragging it with the mouse. You can
remove the splitter bar from the window by dragging it to the end of the
window.
4-9
GAUSS User Guide
Split Vertically Splits the active window into two vertical panes. This allows you to view
two different areas of the same document to facilitate split-window editing.
Note: You can move the splitter bar by dragging it with the mouse. You can
remove the splitter bar from the window by dragging it to the end of the
window.
Open Window GAUSS maintains a list of all the windows you have opened at the end of
List the Window menu. If the window you want to view is on this list, click it and
it becomes the active window.
The Help menu lets you access information in the GAUSS Help system. The GAUSS Help menu
contains the following Commands:
Keyboard Accesses the list of keystrokes you can use for cursor movement, editing,
and text selection.
Reference Accesses the online GAUSS L R, which contains the syntax
for each GAUSS command.
Tip of the Day Displays a tip to help you make better use of the features available in the
GAUSS Windows Interface.
About GAUSS... Provides information about your version of GAUSS, your license type and
ID, as well as copyright information.
The toolbar buttons let you have fast access to the most commonly used commands. Place the
mouse pointer over the button to display a description of the command.
4-10
Introduction to the Windows Interface
GUI Intro
Windows
Save Saves your changes to the file in the active window. If the file is untitled, you
are prompted for a path and filename.
Cut Removes selected text from the active window and places it on the Windows
clipboard.
Copy Copies selected text from the active window to the Windows clipboard.
Paste Copies text from the Windows clipboard to the active window at the cursor
position.
Print Prints the active file or selected text from the active window.
Source Opens the GAUSS Source Browser, which allows you to search for
Browser symbols in a specified file or directory. For more information, see GAUSS
S B, Chapter 8.
Graphics Editor Opens the Graphics Editor, which is an interactive TKF file editor. This
menu item will be inactive if you have not purchased the Graphics Editor.
Run Selected Runs any text selected from the editor or the Command Input-Output
4-11
GAUSS User Guide
Text window.
Run Active Runs the active file. The file then becomes the main file.
File
Main File Displays the name of the main file and lets you quickly change the main file
List to one of the files listed.
Run Main Runs the file specified in the Main File list.
File
Stop Program Stops the program currently running and returns control to the editor.
Test Compile Compiles the main file. During compilation, any errors are displayed in the
Main File Output window.
Note: This command is different than the GAUSS compile command, which
compiles a program and saves the pseudocode as a file.
You can use the Working Directory toolbar to quickly change your working directory.
Current Displays the name of the current working directory and lets you quickly
Working change the working directory to one of the directories listed.
4-12
Introduction to the Windows Interface
Directory List
You can use the Debug toolbar for quick access to commands while debugging a file.
Step Into
Clear All Step Over
Breakpoints Step Out
GUI Intro
Windows
Figure 4.4: Debug Toolbar
Clear All Removes all line and procedure breakpoints from the active file.
Breakpoints
Set Watch Opens the Matrix Editor for watching changing variable data. For more
information about viewing variables see V V, Section 5.5.5.
Step Into Runs the next executable line of code in the application and steps into
procedures.
Step Over Runs the next executable line of code in the application but does not step
into procedures.
4-13
GAUSS User Guide
Step Out Runs the remainder of the current procedure and stops at the next line in the
calling procedure. Step Out returns if a breakpoint is encountered.
Run to Cursor Runs the program until it reaches the cursor position.
You can use the Window toolbar for quick access to window commands.
Activate/ Splits the output from the Command Input-Output window, or deactivates
Deactivate Output window.
Output Window
Tile windows Tiles the active window and the Output or Command Input - Output
horizontally window horizontally.
Tile windows Tiles the active window and the Output or Command Input - Output
4-14
Introduction to the Windows Interface
Src Browser If selected, the next file browsed in the Source Browser will be displayed
replaces window in the same window.
with new file
contents
The status bar is located along the bottom of the GAUSS window. The status of the windows and
processes are shown on the status bar.
The first section of the status bar shows the current GAUSS status. From time to time you are
alerted to the task GAUSS is performing by new messages appearing in the status bar.
GUI Intro
Windows
Dataloop CAP
Cursor Location The line number and column number where the cursor is located appear on
the status bar for the active window. When a block of text is selected, the
values indicate the first position of the selected text.
DATALOOP DATALOOP appears on the status bar to indicate the Dataloop Tranlator is
turned on.
OVR OVR appears on the status bar when typing replaces the existing text with
text you enter. When OVR does not appear on the status bar, typing inserts
4-15
GAUSS User Guide
text without deleting the existing text. Press the INSERT key to toggle
between the two conditions.
CAP CAP appears on the status bar to indicate the Caps Lock key has been
pressed and all text you enter will appear in upper case.
NUM NUM appears on the status bar to indicate the Num Lock key has been
pressed and the keypad numbers are active.
4-16
Using the Windows Interface 5
The GAUSS graphical user interface is a multiple document interface. The interface consists of
edit windows and the Command Input - Output window. Integrated into GAUSS is a full debugger
with breakpoints and watch variables. The GAUSS graphical user interface also incorporates the
Matrix Editor (see Chapter 6), Library Tool (see Chapter 7), and GAUSS Source Browser (see
Chapter 8), as well as a context-sensitive HTML Help system (see Chapter 9).
Windows GUI
5.1 Using the GAUSS Edit Windows
The GAUSS edit windows provide syntax color coding and auto-formatting as well as easy access
to the Matrix Editor and Library Tool, and include an integrated context-sensitive help system
accessible through the F1 key.
The edit windows provide standard text editing features like drag and drop text editing, and find
and replace. The editor also lets you set bookmarks, define keystroke macros, find and replace
using regular expressions, and run selected text from the editor.
5-1
GAUSS User Guide
To begin editing, open an edit window by browsing to the source file, or by typing edit and the
filename in the Command Input - Output window. If more than one file is open, the last file opened
or run becomes the active window.
Bookmarks are efficient placeholders used to identify particular sections or lines of code. To add
or remove bookmarks, place the cursor in the line you want to bookmark and then press
CTRL+F2, or click Toggle Bookmark on the Edit menu. You can jump to the next bookmark by
pressing F2, or go to the previous bookmark by pressing SHIFT+F2.
To edit a list of all currently defined bookmarks, click Edit Bookmarks on the Edit menu. The Edit
Bookmarks window allows you to add, remove, name or select the bookmark to which you wish to
jump.
You can customize the formatting of your code and text by changing font colors, fonts, adding line
indentations, and adding line numbering to your programs. To access these properties, on the
Configure menu click Editor Properties, or right-click on an edit window and click Properties on
the context menu.For more information about the Editor Properties see E P, Section
5.6.2.
To record a keystroke macro, press CTRL+SHIFT+R, or click Record Macro on the Edit menu.
When you start recording the macro, a stop button will appear in the GAUSS window.
5-2
Using the Windows Interface
You create a macro by clicking Record Macro and pressing the keystrokes you want recorded.
Once you have completed recording the macro, you can stop recording with the stop button. Once
you have finished recording the macro, you can select one of ten macro names for it.
Use the following guidelines when creating and using your macro:
• Only keystrokes in the active window are recorded, not keystrokes in a dialog box.
• Only keystrokes are recorded, not mouse movements.
If your macro is lengthy, consider creating a separate file and copying the information from the file
into the active window, rather than using a macro to enter the information.
The margin of the edit window can be used to show currently set bookmarks, currently set
breakpoints, and line numbers. You can also select an entire line of text with a single click in the
Selection Margin.
You can turn on or off the margin in the Misc tab of the Editor Properties dialog box.
Windows GUI
5.1.6 Editing with Split Views
Using split views, you can edit two parts of the same program in the same buffer. To open split
views, click Split Horizontally or Split Vertically on the Window menu.
Along with a standard find and replace function, you can use the edit window to find and replace
regular expressions. To find regular expressions, open the Find dialog box and select the checkbox
for regular expressions.
5-3
GAUSS User Guide
There are three ways you can run selected text. First, highlight the text you want to run, then either
press CTRL+R, drag and drop the selected text into the Command Input - Output window, or click
“Run Selected Text” on the Run menu.
The Command Input - Output window lets you input interactive commands and view the results.
The Command Input - Output window can be split into two separate windows, one for input and
one for output, by clicking Output Window on the Window menu.
Output will be written at the insertion point in the Command Input - Output window or the Output
window, when it is a separate window. GAUSS commands cannot be executed from this window.
From the Command Input - Output window, you can run saved programs. You can view or edit the
data of any variable in the active workspace with the Matrix Editor. You can also open files for
editing or to debug.
The GAUSS Command Input - Output window has many of the same features that the GAUSS text
editor has. You can cut and paste text. You can search the buffer of the Command Input - Output
window. You can also save the contents of the Command Input - Output window to a text file.
The GAUSS interface allows you to run programs that consist of single commands or blocks of
commands executed interactively, as well as large-scale programs that may consist of commands
in one or more files. The file that is run to execute the command is the main file (the file name
displayed in the Main File list).
When you run commands interactively, the actual code being processed is called the “active
block.” The active block is all code between the GAUSS prompt (>>) and the end of the current
line. Thus, the active block can be one or more lines of code.
5-4
Using the Windows Interface
Interactive commands can be entered at the “>>” prompt in the Command Input - Output window
or selected using the mouse and clicking the Run Selected Text button on the Main toolbar.
A block of code can be executed by selecting the block with the mouse and then running that
block using the Run Selected Text function.
Note: The GAUSS prompt (>>) at the beginning of the selected text is ignored.
You can enter multi-line commands into the Command Input - Output window by pressing
CTRL+ENTER at the end of each line. At the end of the final line in a multi-line command, press
ENTER. The Command Input - Output window will automatically place a semicolon at the end of
a single-line command before it is interpreted. For multi-line commands, you must enter a
semicolon at the end of each line.
You can also run multi-line commands by pasting the text of a file at the GAUSS prompt, or
selecting multiple lines of code from the Command Input - Output window and pressing CTRL+R.
You can repeat any of the last 20 lines entered into the command buffer by pressing CTRL+L to
cycle through the last command buffer.
Windows GUI
You can execute the active file by clicking Run Active File on the Run menu, or by clicking the
Run Currently Active File button on the Main toolbar.
You can execute the file displayed in the Main File list (the main file) by clicking Run Main file on
the Run menu, or by clicking the Run Main File button on the Main toolbar.
Source View is a dockable dialog bar with two tabs that provide easy access to source files and
symbols associated with your current GAUSS workspace.
5-5
GAUSS User Guide
The Source tab is a tree view that displays a list of active libraries and the source files they contain.
Under each source file is a list of the symbols and procedures which they define. By using the
right mouse button, you can search for symbols, open source files or view source file properties.
To open a source file, double click the file name or right click the file and click Edit.
To search the source files right click any file name in the source tab and click Find. In the Find
dialog enter a keyword and click OK.
The Symbols tab contains a tree view of the GAUSS workspace global symbols organized by
symbol type: Matrices, Arrays, Strings, String Arrays, and Structures.
To edit or view a symbol, double-clicking on it or right-clicking and selecting Edit from the menu.
To search the source files right click any file name in the source tab and click Find. In the Find
dialog enter a keyword and click OK.
5-6
Using the Windows Interface
The Error Output window allows errors messages to be output to a separate window, instead of the
GAUSS Input - Output window. When an error occurrs, you can open to program of source file
directly from the Error Output window.
To open the program or source file, press F4 or double click the error message. The file will open
at the line the error occurred.
The debugger greatly simplifies program development. With all of the features of a dedicated
debugging system, the debugger can help you to quickly identify and solve logic errors at run-time.
The debugger is integrated into the multiple document interface of GAUSS; it uses the interface
tools, such as the edit windows, the Matrix Editor, and the Command Input - Output window for
debugging. So while using the debugger, you still have all the features of the edit windows and
Matrix Editor, along with GAUSS’s suite of debugging tools.
You use the debugger to watch the program code as it runs. Prior to running the debugger,
breakpoints and watch variables can be set to stop the program at points you set and provide
additional data as the code is run.
Windows GUI
5.5.1 Starting and Stopping the Debugger
You can start the debugger by clicking Go on the Debug menu or the Debug toolbar.
When starting the debugger, you can choose to debug the active file or to debug the main file of a
program. If you are debugging a single file and already have the file open, you can use the menu or
toolbar to start the debugger on the file, or simply type debug and the filename in the Command
Input - Output window.
When you start the debugger, the debugger automatically highlights the first line of code to be run.
Any breakpoints are shown in the left margin of the window.
5-7
GAUSS User Guide
You can stop the debugger at any time by clicking Stop on the Debug menu or the Debug toolbar.
Breakpoints stop code execution where you have inserted them. Breakpoints are normally set prior
to running the debugger, but can also be set or cleared during debugging by clicking the Set/Clear
Breakpoint command on the Debug menu.
The debugger supports two types of breakpoints: procedure breakpoints and line number
breakpoints. Procedure breakpoints pause execution when the specified procedure or function is
reached. Line number breakpoints pause execution when the specified line is reached. In either
case, the break occurs before any of the GAUSS code for the procedure or line is executed. The
debugger also allows you to specify a certain cycle of execution for a line number or procedure
where you want the execution to be paused. The cycle count is for the occurrence of the line
number or procedure, not the number of times a line is to be skipped.
You can set or clear a line breakpoint in the highlighted line of code by clicking Set/Clear
Breakpoint on the Debug menu or by pressing the F9 key.
To set breakpoints in any part of the file not currently being executed, just click the line where you
want the breakpoint to be, then click Toggle Breakpoint.
To clear breakpoints in the file, click a line of code that has a breakpoint set and then click
Set/Clear Breakpoint. You can also clear all breakpoints from the active file by clicking Clear All
Breakpoints.
The Breakpoint Editor allows you to set or clear both line and procedure breakpoints. It also lets
you specify cycles of execution for breakpoints. With the Breakpoint Editor, you can set or clear
breakpoints in any program currently in your working directory.
5-8
Using the Windows Interface
GAUSS’s debugger includes the ability to step into, step out of, and step over code during
debugging.
Use Step Into to execute the line of code currently highlighted by the debugger.
Use Step Out to execute to the end of the current function without pause and return to the calling
function.
Use Step Over to execute the line of code currently highlighted by the debugger without entering
the functions that are called.
GAUSS allows you to view and edit the values of variables during debugging.
Once the debugger is started, the editor window uses floatover variable windows for viewing
variable data. Floatover variable windows give a quick view of the value a variable currently holds
Windows GUI
by simply moving your mouse over the variable name in the edit window.
The floatover variable window is only intended to give a quick view of the data, so it may not
show all data held by the variable. If the variable data is incomplete, the floatover variable window
will display an arrow to show that there is more data. If you need to view more data, open the
Matrix Editor by highlighting the variable name and pressing CTRL+E.
The debugger integrates the Matrix Editor to edit values of loaded variables, or to use as a watch
window to view the changing values of variables as you step through a program.
5-9
GAUSS User Guide
To edit a variable value, highlight the variable in the edit window, or the Command Input - Output
window and then open the Matrix Editor. You can use the menu or toolbar to start the Matrix
Editor, or simply type CTRL+E.
You can make the Matrix Editor a Watch window, allowing you to watch the changing value of a
variable as the lines of the program are executed. You can activate the Watch window by clicking
Set Watch on the Debug menu, or by highlighting a variable name in the debugger window and
pressing CTRL+E.
You use a Watch window to see how variables change in value during debugging. Watch variables
can be specified prior to running the debugger or during a debugging session.
The debugger searches for a watch variable using the following order:
2. A global variable.
A watch variable can be the name of a matrix, a scalar, a string array, or a string. For a matrix or a
string array, the first element is displayed. If a matrix element is clicked, the Matrix Editor is
loaded with the matrix. The matrix elements can be changed during the debugging session.
The Preferences dialog box lets you specify how GAUSS operates. To open the Preferences dialog
box, click Preferences... on the Configure menu. The changes you make in the Preferences dialog
box remain set between sessions.
5-10
Using the Windows Interface
Run Options
Dataloop Specifies whether or not GAUSS will translate data loops into
Translator procedures.
Translate Line Specifies whether or not GAUSS will preserve the line numbers of data
Number loops after being translated to procedures.
Tracking
Line Number Specifies whether or not GAUSS will preserve line numbers of a file
Tracking being compiled for the interpreter.
Sound at End Determines whether or not a sound is played at the end of the execution
of Job of GAUSS code. The sound can be selected using the Select button
and played using the Test button.
The default is OFF.
Compile Options
The Compile tab contains options that let you control how GAUSS compiles a program before it is
run.
Windows GUI
Autoload Specifies whether the autoloader will automatically resolve references
in your code. If Autoload is off, you must define all symbols used in
your program.
GAUSS Library Specifies whether the autoloader will use the standard GAUSS library
in compiling your code.
User Library Specifies whether the autoloader will use the User Libraries in
compiling your code.
Declare Specifies whether the GAUSS compiler will display declare warnings
5-11
GAUSS User Guide
Compiler Trace Specifies whether you would like to trace the file compilation by file
opening and closing, specific lines, or whether you would like to trace
by local and global symbols.
Cmd Window
The Cmd Window tab contains options that let you control how the GAUSS Command Window
operates.
Action on Specifies whether pressing ENTER executes the current whole line
Enter always or only when the cursor is at the end of a line. Also, specifies
whether placing a semi-colon at the end of a line causes GAUSS to
enter multi-line mode.
Performance Specifies whether or not output is buffered, and sets the buffer size in
kilobytes.
Cmd Prompt Specifies whether new GAUSS prompts are inserted at the current
cursor location, appended to the text in the Command window, or
relocated on the line following the cursor.
DOS Compatibility
The DOS Compatibility tab lets you control the appearance of the DOS Compatibility window.
Change Font Specifies what font the DOS Compatibility window will use.
5-12
Using the Windows Interface
Stay on top of Specifies whether the DOS Compatibility window will stay on top of
GAUSS GAUSS.
File
The File tab contains options that let you control how GAUSS auto-saves your work.
Misc
Windows GUI
The Misc tab contains several general options to control GAUSS.
Show Tip of the Turns on/off the Tip of the Day at startup.
Day at startup
Keep Help Specifies whether the Help window always stays on top of GAUSS
Window On Top when opened.
of GAUSS
Set Initial Specifies initial position of the Help
Window Position window.
Graphics Editor If you have purchased the Graphics Editor,
License this is where you will enter or change the Graphics Editor license key.
5-13
GAUSS User Guide
You can customize the formatting of your code and text by changing font colors, fonts, adding line
indentations, and adding line numbering to your programs. To access these properties, on the
Configure menu click Editor Properties.
Color/Font
Language/Tabs
Auto Indentation Specifies how the autoindenter will indent your code.
Style
Language Specifies what syntax the GAUSS editor will recognize for syntax
coloring.
Fixup Text Case Specifies whether the editor will automatically change the case of
While Typing GAUSS keywords when they use the wrong case.
Language
Keywords
Misc
Show Left Enables or disables the editor’s margin. The margin is used for
Margin showing breakpoints, bookmarks, or line numbers.
Line Tooltips Shows the first line number on screen as a tooltip as you scroll up and
on Scroll down the file.
5-14
Using the Windows Interface
Confine Caret Tells the GAUSS editor to interpret carets as text only rather than as
to Text substitution symbols or text.
Allow Vertical Toggles on or off the ability to split editor panes vertically.
Splitting
Line Numbering Specifies the style and starting digit for line numbering.
Max Undoable Sets the number of actions that you can undo.
Windows GUI
Actions
5-15
GAUSS User Guide
SHIFT+DEL or CTRL+X Delete selected text and place it onto Windows clipboard
SHIFT+INS or CTRL+V Paste text from Windows clipboard at the cursor position
5-16
Using the Windows Interface
Windows GUI
CTRL+A Redo
5-17
GAUSS User Guide
CTRL+Z Undo
F2 Go to next bookmark
F3 Find again
F8 Step Into
F9 Set/Clear breakpoint
5-18
Using the Windows Interface
CTRL+F1 Searches the active libraries for the source code of a function.
Windows GUI
ALT+F File menu
5-19
Matrix Editor 6
6.1 Using the Matrix Editor
The Matrix Editor lets you view and edit matrix data in your current workspace. You can open the
Matrix Editor from either the Command Input - Output window or a GAUSS edit window by
highlighting a matrix variable name and typing CTRL+E. You can view multiple matrices at the
same time by opening more than one Matrix Editor.
The Matrix Editor will allow you to format matrices in decimal, scientific, Hexadecimal, or as text
Matrix Editor
characters.
Just like a spreadsheet, when using the Matrix Editor, you can use your keyboard’s arrow keys to
quickly move between matrix positions. To edit a scalar value, select a cell and press Enter. You
can use the Home and End keys to move to the beginning or end of a scalar. When finished
editing, press Enter again.
6-1
GAUSS User Guide
6-2
Matrix Editor
All variables are treated as matrices in GAUSS. A scalar is simply a 1×1 matrix. A vector is a
(N×1) or (1×N) matrix. So you can use the Matrix Editor to view and monitor the value of any
variable. You can update the value of a variable at any time by using the Reload function. When
using the Matrix Editor to view, edit or monitor smaller matrices, you can minimize space it
occupies on the screen by selecting Minimal View from the View menu.
By using the Auto-reload function, GAUSS will automatically update the values of variables in
the Matrix Editor. Using Auto-reload you can create a watch window.
Watch Variables allow you to see how variables change in value while debugging a program. A
watch variable can be the name of a matrix, a scalar, an array, a string array, or a string.
• a global variable
Matrix Menu
Matrix Editor
The Matrix menu lets you control the data of the Matrix in the Matrix Editor as an entire set.
Load Clears any existing grid and loads any named matrix from the GAUSS
workspace to the grid.
Reload Reloads the existing matrix with the name shown on the Title bar.
6-3
GAUSS User Guide
Auto-Reload Automatically updates the data shown in the Matrix Editor, creating a
watch window.
Save Saves the grid as a matrix in the GAUSS workspace. If a matrix of the
same name already exists in the workspace, it is overwritten.
Format Menu
The Format menu lets you control the way the data is presented in the Matrix Editor.
Edit Menu
The Edit menu gives you tools to control the data in the Matrix Editor.
Clear All Clears the grid of all values but keep the row and column order.
Preferences Sets several matrix options, including the number of digits to the right
of the decimal point, cell height and width, and whether pressing the
Enter key moves the cursor down or over one cell. These options, along
with screen position and window state, are saved between sessions.
6-4
Matrix Editor
View Menu
The View menu lets you control the Matrix Editor window. The View menu also lets you control
your view of imaginary numbers.
Real Parts Specifies that you want the real parts of imaginary numbers to be
displayed in the Matrix Editor.
Imaginary Parts Specifies that you want the imaginary parts of numbers to be displayed
in the Matrix Editor.
Minimal View Minimizes the amount of screen space occupied by the Matrix Editor.
This is especially useful for creating watch windows for single
variables.
Stay on Top Forces the Matrix Editor window to remain visible on the screen even
when the interface focus has shifted to another window.
Matrix Editor
6-5
Library Tool
Library Tool 7
7.1 Using the Library Tool
The Library Tool lets you quickly manage your libraries. You can add and remove libraries and
you can add and remove files within the libraries.
Using the New Library button, you can create a new library for organizing your code. You can
remove a library by selecting the Delete Library button.
To add absolute path names to the library index, use the Add Paths button. To only use file names
for searching libraries, use the Strip Paths button. Use Rebuild to recompile all the files used in the
library, and rebuild the library index file. Use the Revert to Original button to revert to the
configuration the library was in when the Library Tool was opened.
7-1
GAUSS User Guide
7-2
Library Tool
Library Tool
7.1.3 Managing Library Files
You can add files to a library with the Add button. You can remove files from a library with the
Remove button. After changing source files referred to in a library, select the files in the file list
and update the library index with the Update button. To remove multiple files from a library, select
the files in the file selection window, and use the Clear Selection button.
7-3
Browser
Source
GAUSS Source Browser 8
The GAUSS Source Browser lets users quickly find, view, and if necessary, modify source code.
Both the TGAUSS and GAUSS Source Browsers can be used to search for external symbols in
active libraries. The GAUSS Source Browser can also be used to search for symbols in any
directory or source file.
To start the Source Browser in TGAUSS, type BROWSE followed by a symbol name. When the
Source Browser is active, the prompt displays Browse:. GAUSS searches through all active
libraries for the file in which the symbol is defined. If found, the file containing the source code is
opened in the default editor.
Wildcard (*) searches can also be used. When using wildcard searches, each symbol that the string
matches will be displayed on-screen in a numbered list. To select a specific command to view in
the default editor, select the number from the list.
The Source Browser will remain active until you type CTRL-C to return to the GAUSS prompt.
8-1
GAUSS User Guide
To open the Source Browser in GAUSS, from the Tools menu select Source Browser.
Using the Source Browser you can search a file for a specified symbol or search across all the
source files in a directory. Using a comma separated list of files and directories in the Look in: list,
you can search multiple locations in a single search. When searching for symbols, you can use
wildcards (*) to further modify the scope of the search.
Note: The Source Browser does not search recursively through sub-folders. To search sub-folders
during a search, add the sub-folder names to the Look in: list.
Once the search is complete, the Source Browser lists where the specified symbol was found. The
Filename column of the Results List shows the file in which the symbol was found. The Line
column shows the line number where symbol was found and the Line Description column shows
the text of the line where the symbol was found.
Search locations typed into the Look in: text box will persist between Source Browser sessions.
Look in: Limits the scope of the search to specific files or directories. Using a comma
separated list, searches multiple files and directories in a single search.
Stay on top Keeps the Source Browser on top even when another window is active.
Browse Lets you limit the scope of the search to specific files or directories.
8-2
GAUSS Source Browser
Browser
Source
Figure 8.1: Source Browser
8-3
GAUSS User Guide
Status Lists how many occurrences there were, and how many files the symbol
occurred in.
Close Closes the Source Browser.
Double-click the file name to open a file in its own editor window. When opened, the cursor is
placed at the beginning of the line selected in the Results List. By double-clicking different files in
the Source Browser, you can open each file in its own separate editor window.
Use the F4 key to quickly view or edit the next file in the Results List using the active editor
window. Using the F4 key opens the file in the active editor window and places the cursor at the
beginning of the line in which the symbol was found. The F4 key uses the active editor window to
display the source file; it will not open an editor window to display files. You can use the F4 key
from either the Source Browser or from the active editor window to move to the next occurrence
of the symbol shown in the Results List.
Use SHIFT+F4 to quickly view or edit the previous file in the Results List using the active editor
window. Using the F4 key opens the file in the active editor window and places the cursor at the
beginning of the line in which the symbol was found.
8-4
GAUSS Help 9
Help
9.1 Help Menu
From the Help menu, you can directly access the online U G, Keyboard Assignments list,
and L R. Pressing F1 also accesses the Help system, displaying either the
I to the U G or, if an object has focus and help can be directly accessed, help
for that object.
GAUSS integrates a context-sensitive Help system to help you use the GAUSS environment and
the GAUSS language. Context-sensitive means that Help for the object with focus is displayed
without navigating through the Help system. For example, to display Help on a keyword in the
GAUSS language in a GAUSS edit window or the Command Input - Output window, place the
insertion point on the keyword and press F1.
9-1
GAUSS User Guide
• GAUSS windows
• Toolbar buttons
• GAUSS menus
• The GAUSS language
For intrinsic commands and functions, the GAUSS L R for the command is
displayed. For other external procedures in active libraries, a window displays a source code file,
allowing you to scroll to the desired symbol.
If you press SHIFT+F1 or click on the Help toolbar button (an arrow with a question mark), the
pointer changes to a Help pointer (arrow + ?). Click on an object to display the Help system or, if
available, context-sensitive Help for that object.
You can search through all active libraries for any global symbol by placing the cursor on the
symbol name and pressing CTRL+F1.
GAUSS searches through all active libraries for the file that the symbol is defined in. If found, the
file containing the source code is opened in an edit window. If the file contains the code string
“**> symbol name” (without quotes) at the beginning of a line of commented code, the cursor will
be placed at the beginning of that line. If the string is not found in the file, the cursor will be
placed at the beginning of the file.
9-2
GAUSS Help
9.5 ToolTips
A ToolTip is a small label that is displayed when the mouse pointer is held over a GAUSS button.
The ToolTip will give a brief description of the button’s function.
GAUSS includes full online versions of the GAUSS L R and GAUSS U G
in PDF format. These manuals are available for download; contact Aptech Systems for download
Help
instructions.
The Gaussians mail list is an e-mail list providing users of GAUSS an easy way to reach other
GAUSS users. Gaussians provides a forum for information exchange, tips and experiences using
GAUSS. For more information about the Gaussians mail list, see the Resource Library page at
http://www.aptech.com. You can also e-mail [email protected] when you have current
Premier Support.
9-3
Language Fundamentals 10
Fundamentals
Language
GAUSS is a compiled language. GAUSS is also an interpreter. A compiled language, because
GAUSS scans the entire program once and translates it into a binary code before it starts to execute
the program. An interpreter, because the binary code is not the native code of the CPU. When
GAUSS executes the binary pseudocode it must “interpret” each instruction for the computer.
How can GAUSS be so fast if it is an interpreter? Two reasons. First, GAUSS has a fast
interpreter, and the binary compiled code is compact and efficient. Second, and most significantly,
GAUSS is a matrix language. It is designed to tackle problems that can be solved in terms of
matrix or vector equations. Much of the time lost in interpreting the pseudocode is made up in the
matrix or vector operations.
This chapter will enable you to understand the distinction between “compile time” and “execution
time”, two very different stages in the life of a GAUSS program.
10.1 Expressions
10-1
GAUSS User Guide
10.2 Statements
y = x*3;
print x*3;
or
x*3;
output on;
There can be multiple statements on the same line as long as each statement is terminated with a
semicolon.
10-2
Language Fundamentals
Executable statements are statements that can be “executed” over and over during the execution
phase of a GAUSS program (execution time). As an executable statement is compiled, binary code
is added to the program being compiled at the current location of the instruction pointer. This
binary code will be executed whenever the interpreter passes through this section of the program.
If this code is in a loop, it will be executed each iteration of the loop.
y = 34.25;
print y;
x = 1 3 7 2 9 4 0 3 ;
Fundamentals
Language
10.2.2 Nonexecutable Statements
Nonexecutable statements are statements that have an effect only when the program is compiled
(compile time). They generate no executable code at the current location of the instruction pointer.
declare matrix x = 1 2 3 4 ;
Procedure definitions are nonexecutable. They do not generate executable code at the current
location of the instruction pointer.
Here is an example:
10-3
GAUSS User Guide
zed = rndn(3,3);
proc sqrtinv(x);
local y;
y = sqrt(x);
retp(y+inv(x));
endp;
zsi = sqrtinv(zed);
There are two executable statements in the example above: the first line and the last line. In the
binary code that is generated, the last line will follow immediately after the first line. The last line
is the call to the procedure. This generates executable code. The procedure definition generates
no code at the current location of the instruction pointer.
There is code generated in the procedure definition, but it is isolated from the rest of the program.
It is executable only within the scope of the procedure and can be reached only by calling the
procedure.
10.3 Programs
A program is any set of statements that are run together at one time. There are two sections within
a program.
The main section of the program is all of the code that is compiled together WITHOUT relying on
the autoloader. This means code that is in the main file or is included in the compilation of the
main file with an #include statement. ALL executable code should be in the main section.
There must always be a main section even if it consists only of a call to the one and only procedure
called in the program.
10-4
Language Fundamentals
Secondary sections of the program are files that are neither run directly nor included in the main
section with #include statements.
The secondary sections of the program can be left to the autoloader to locate and compile when
they are needed. Secondary sections must have only procedure definitions and other
nonexecutable statements.
#include statements are allowed in secondary sections as long as the file being included does not
violate the above criteria.
proc feq(a,b);
Fundamentals
Language
retp(abs(a-b) <= tol);
endp;
Compiler directives are commands that tell GAUSS how to process a program during compilation.
Directives determine what the final compiled form of a program will be. They can affect part or all
of the source code for a program. Directives are not executable statements and have no effect at
run-time. They do not take a semicolon at the end of the line.
The #include statement mentioned earlier is actually a compiler directive. It tells GAUSS to
compile code from a separate file as though it were actually part of the file being compiled. This
code is compiled in at the position of the #include statement.
10-5
GAUSS User Guide
#lineson Compile program with line number and file name records.
#linesoff Compile program without line number and file name records.
#srcfile Insert source file name record at this point (currently used when
doing data loop translation).
#srcline Insert source file line number record at this point (currently used
when doing data loop translation).
The #define statement can be used to define abstract constants. For example, you could define
the default graphics page size as:
10-6
Language Fundamentals
and then write your program using hpage and vpage. GAUSS will replace them with 9.0 and
6.855 when it compiles the program. This makes a program much more readable.
#ifdef log_10
y = log(x);
#else
y = ln(x);
#endif
This allows the same program to calculate answers using different base logarithms, depending on
whether or not the program has a #define log_10 statement at the top.
#undef allows you to undefine text-replacement or flag variables so they no longer affect a
Fundamentals
Language
program, or so you can #define them again with a different value for a different section of the
program. If you use #definecs to define a case-sensitive variable, you must use the right case
when #undef’ing it.
With #lineson, #linesoff, #srcline, and #srcfile you can include line number and file
name records in your compiled code, so that run-time errors will be easier to track down.
#srcline and #srcfile are currently used by GAUSS when doing data loop translation.
For more information on line number tracking, see D, Section 20.3 and see D
D L, Section 23.3. See also #lineson in the GAUSS L R.
The syntax for #srcfile and #srcline is different than for the other directives that take
arguments. Typically, directives do not take arguments in parentheses; that is, they look like
keywords:
#define red 4
#srcfile and #srcline, however, do take their arguments in parentheses (like procedures):
10-7
GAUSS User Guide
#srcline(12)
This allows you to place #srcline statements in the middle of GAUSS commands, so that line
numbers are reported precisely as you want them. For example:
#srcfile(/gauss/test.e)
10.5 Procedures
A procedure allows you to define a new function which you can then use as if it were an intrinsic
function. It is called in the same way as an intrinsic function.
y = myproc(a,b,c);
Procedures are isolated from the rest of your program and cannot be entered except by calling
them. Some or all of the variables inside a procedure can be local variables . local variables
exist only when the procedure is actually executing and then disappear. Local variables cannot get
mixed up with other variables of the same name in your main program or in other procedures.
For details on defining and calling procedures, see P K, chapter 12.
10-8
Language Fundamentals
There are four basic data types in GAUSS, matrices, N-dimensional arrays, strings and string
arrays. It is not necessary to declare the type of a variable, but it is good programming practice to
respect the types of variables whenever possible. The data type and size can change in the course
of a program.
The declare statement, used for compile-time initialization, enforces type checking.
Short strings of up to 8 bytes can be entered into elements of matrices, to form character matrices
(For details, see C M, Section 10.6.7).
10.6.1 Constants
Fundamentals
Language
Decimal
1.34e-10
1.34e123
-1.34e+10
-1.34d-10
1.34d10
10-9
GAUSS User Guide
1.34d+10
123.456789345
Up to 18 consecutive digits before and after the decimal point(depending on the platform) are
significant, but the final result will be rounded to double precision if necessary. The range is the
same as for matrices (For details, see M, Section 10.6.2.
String
“This is a string.”
Hexadecimal Integer
0x0ab53def2
Hexadecimal floating point constants are prefixed with 0v. This allows you to input a double
precision value exactly as you want using 16 hexadecimal digits. The highest order byte is to the
left:
0vfff8000000000000
10-10
Language Fundamentals
10.6.2 Matrices
Matrices are 2-dimensional arrays of double precision numbers. All matrices are implicitly
complex, although if it consists only of zeros, the imaginary part may take up no space. Matrices
are stored in row major order. A 2×3 real matrix will be stored in the following way from the
lowest addressed element to the highest addressed element:
A 2×3 complex matrix will be stored in the following way from the lowest addressed element to
the highest addressed element:
Conversion between complex and real matrices occurs automatically and is transparent to the user
Fundamentals
Language
in most cases. Functions are provided to provide explicit control when necessary.
All elements of a GAUSS matrix are stored in double precision floating point format, and each
takes up 8 bytes of memory. This is the IEEE 754 format:
Matrices with only one element (1×1 matrices) are referred to as scalars, and matrices with only
one row or column (1×N or N×1 matrices) are referred to as vectors.
Any matrix or vector can be indexed with two indices. Vectors can be indexed with one index.
Scalars can be indexed with one or two indices also, because scalars, vectors, and matrices are the
same data type to GAUSS
The majority of functions and operators in GAUSS take matrices as arguments. The following
functions and operators are used for defining, saving, and loading matrices:
10-11
GAUSS User Guide
[ ] Indexing matrices.
= Assignment operator.
| Vertical concatenation.
∼ Horizontal concatenation.
An assignment statement followed by data enclosed in braces is an implicit let statement. Only
constants are allowed in let statements; operators are illegal. When braces are used in let
statements, commas are used to separate rows. The statement
let x = 1 2 3, 4 5 6, 7 8 9 ;
or
10-12
Language Fundamentals
x = 1 2 3, 4 5 6, 7 8 9 ;
will result in
1 2 3
x= 4 5 6
7 8 9
The statement
let x[3,3] = 1 2 3 4 5 6 7 8 9;
will result in
1 2 3
Fundamentals
Language
x= 4 5 6
7 8 9
The statement
let x[3,3] = 1;
will result in
1 1 1
x= 1 1 1
1 1 1
The statement
let x[3,3];
10-13
GAUSS User Guide
will result in
0 0 0
x= 0 0 0
0 0 0
The statement
let x = 1 2 3 4 5 6 7 8 9;
will result in
1
2
3
4
x= 5
6
7
8
9
Complex constants can be entered in a let statement. In the following example, the + or - is not a
mathematical operator, but connects the two parts of a complex number. There should be no
spaces between the + or - and the parts of the number. If a number has both real and imaginary
parts, the trailing ‘i’ is not necessary. If a number has no real part, you can indicate that it is
imaginary by appending the ‘i’. The statement
will result in
1 + 2i 3 − 4i
x=
5 0 + 6i
10-14
Language Fundamentals
Complex constants can also be used with the declare, con and stof statements.
An “empty matrix” is a matrix that contains no data. Empty matrices are created with the let
statement and braces:
x = {};
Empty matrices are supported by several functions, including rows and cols and the
concatenation (∼,|) operators.
x = {};
hsec0 = hsec;
do until hsec-hsec0 > 6000;
x = x ˜ data_in(hsec-hsec0);
endo;
Fundamentals
Language
You can test whether a matrix is empty by entering rows(x), cols(x) and scalerr(x). If the
matrix is empty rows and cols will return a 0, and scalerr will return 65535.
The ∼ is the horizontal concatenation operator and the | is the vertical concatenation operator. The
statement
y = 1∼2|3∼4;
will be evaluated as
y = (1 ∼ 2) | (3 ∼ 4);
and will result in a 2×2 matrix because horizontal concatenation has precedence over vertical
concatenation:
1 2
3 4
10-15
GAUSS User Guide
The statement
y = 1+1∼2*2|3-2∼6/2;
will be evaluated as
and will result in a 2×2 matrix because the arithmetic operators have precedence over
concatenation:
2 4
1 3
let x[2,2] = 1 2 3 4;
Unlike the concatenation operators, it cannot be used to define matrices in terms of expressions
such as
y = x1-x2∼x2|x3*3∼x4;
The statement
y = x[1:3,5:8];
10-16
Language Fundamentals
will put the intersection of the first three rows and the fifth through eighth columns of x into the
matrix y.
The statement
will create a 3×3 matrix y with the intersection of the specified rows and columns pulled from x
(in the indicated order).
will have the same effect as the previous example, but is more general.
Fundamentals
Language
The statement
y[2,4] = 3;
will set the 2,4 element of the existing matrix y to 3. This statement is illegal if y does not have at
least 2 rows and 4 columns.
The statement
x = con(3,2);
- (1,1)
10-17
GAUSS User Guide
indicating that the user should enter the [1,1] element of the matrix. Entering a number and then
pressing ENTER will cause a prompt for the next element of the matrix to appear. Pressing ? will
display a help screen, and pressing x will exit.
The statement
will load data contained in an ASCII file into an N×1 vector x. (Use rows(x) to find out how
many numbers were loaded, and use reshape(x,N,K) to reshape it to an N×K matrix).
The statement
load x;
will load the matrix x.fmt from disk (using the current load path) into the matrix x in memory.
The statement
open d1 = dat1;
x = readr(d1,100);
will read the first 100 rows of the GAUSS data set dat1.dat.
Many GAUSS operators and commands support the sparse matrix data type. You may use any of
the following commands to create a sparse matrix:
10-18
Language Fundamentals
packedToSp Creates a sparse matrix from a packed matrix of non-zero values and
row and column indices.
Fundamentals
Language
Many GAUSS commands support arrays of N dimensions. The following commands may be used
to create and manipulate an N-dimensional array:
10-19
GAUSS User Guide
10.6.5 Strings
Strings can be used to store the names of files to be opened, messages to be printed, entire files, or
whatever else you might need. Any byte value is legal in a string from 0–255. The buffer where a
string is stored always contains a terminating byte of ASCII 0. This allows passing strings as
arguments to C functions through the Foreign Language Interface.
strsplit Split an N×1 string vector into an N×K string array of the individual
tokens.
10-20
Language Fundamentals
strsplitPad Split a string vector into a string array of the individual tokens. Pads
on the right with null strings.
x = "example string";
Fundamentals
Language
or
or
print x;
print $x;
10-21
GAUSS User Guide
A string can be saved to disk with the save command in a file with a .fst extension and then
loaded with the load command:
save x;
loads x;
or
loads x=x.fst;
The backslash is used as the escape character inside double quotes to enter special characters:
When entering DOS pathnames in double quotes, two backslashes must be used to insert one
backslash:
st = "c:\\gauss\\myprog.prg";
An important use of strings and character elements of matrices is with the substitution operator (ˆ).
In the command
10-22
Language Fundamentals
by default, GAUSS will interpret the olsdat as a literal; that is, the literal name of the GAUSS
data file you want to create. It will also interpret the x as the literal prefix string for the variable
names: x1 x2 x3 x4.
If you want to get the data set name from a string variable, the substitution operator (ˆ) could be
used as:
dataset="olsdat";
create f1=ˆdataset with x,4,2;
If you want to get the data set name from a string variable and the variable names from a character
vector, use
dataset="olsdat";
let vnames=age pay sex;
Fundamentals
Language
create f1=ˆdataset with ˆvnames,0,2;
The substitution operator (ˆ) works with load and save also:
lpath="/gauss/procs";
name="mydata";
load path=ˆlpath x=ˆname;
command="dir *.fmt";
ˆvariable name
Expressions are not allowed. The following commands are supported with the substitution
operator (ˆ):
10-23
GAUSS User Guide
String arrays are N×K matrices of strings. Here is a partial list of the functions for manipulating
string arrays:
0 Transpose operator.
10-24
Language Fundamentals
format Define output format for matrices, string arrays, and strings.
disk and gives the disk file either a .fmt, .fst or .fcg extension.
Fundamentals
Language
show Display global symbol table.
string array.
String arrays are created through the use of the string array concatenation operators. Below is a
contrast of the horizontal string and horizontal string array concatenation operators.
10-25
GAUSS User Guide
x = "age";
y = "pay";
n = "sex";
s = x $+ y $+ n;
sa = x $∼ y $∼ n;
s = agepaysex
Matrices can have either numeric or character elements. For convenience, a matrix containing
character elements is referred to as a character matrix.
A character matrix is not a separate data type, but gives you the ability to store and manipulate
data elements that are composed of ASCII characters as well as floating point numbers. For
example, you may want to concatenate a column vector containing the names of the variables in an
analysis onto a matrix containing the coefficients, standard errors, t-statistic, and p-value. You can
then print out the entire matrix with a separate format for each column with one call to the
function printfm.
The logic of the programs will dictate the type of data assigned to a matrix, and the increased
flexibility allowed by being able to bundle both types of data together in a single matrix can be
very powerful. You could, for instance, create a moment matrix from your data, concatenate a new
row onto it containing the names of the variables and save it to disk with the save command.
Numeric matrices are double precision, which means that each element is stored in 8 bytes. A
character matrix can thus have elements of up to 8 characters.
GAUSS does not automatically keep track of whether a matrix contains character or numeric
information. The ASCII to GAUSS conversion program ATOG will record the types of variables
in a data set when it creates it. The create command will, also. The function vartypef gets a
vector of variable type information from a data set. This vector of ones and zeros can be used by
printfm when printing your data. Since GAUSS does not know whether a matrix has character or
10-26
Language Fundamentals
numeric information, it is up to you to specify which type of data it contains when printing the
contents of the matrix. (For details, see print and printfm in the GAUSS L R.)
Most functions that take a string argument will take an element of a character matrix also,
interpreting it as a string of up to 8 characters.
DT Scalar Format
The DT scalar format is a double precision representation of the date and time. In the DT scalar
format, the number
20010421183207
Fundamentals
Language
DTV Vector Format
The DTV vector is a 1×8 vector. The format for the DTV vector is:
[1] Year
[2] Month, 1-12
[3] Day of month, 1-31
[4] Hour of day, 0-23
[5] Minute of hour, 0-59
[6] Second of minute, 0-59
[7] Day of week, 0-6 where 0 is Sunday
[8] Day since beginning of year, 0-365
The UTC scalar format is the number of seconds since January 1, 1970, Greenwich Mean Time.
10-27
GAUSS User Guide
The IEEE floating point format has many encodings that have special meaning. The print
command will print them accurately so that you can tell if your calculation is producing
meaningful results.
NaN
There are many floating point encodings which do not correspond to a real number. These
encodings are referred to as NaN’s. NaN stands for Not A Number.
Certain numerical errors will cause the math coprocessor to create a NaN called an “indefinite”.
This will be printed as a -NaN when using the print command. These values are created by the
following operations:
• +∞ plus −∞
• +∞ minus +∞
• −∞ minus −∞
• 0∗∞
• ∞/∞
• 0/0
INF
When the math coprocessor overflows, the result will be a properly signed infinity. Subsequent
calculations will not deal well with an infinity; it usually signals an error in your program. The
result of an operation involving an infinity is most often a NaN.
10-28
Language Fundamentals
DEN, UNN
When some math coprocessors underflow, they may do so gradually by shifting the significand of
the number as necessary to keep the exponent in range. The result of this is a denormal (DEN).
When denormals are used in calculations, they are usually handled automatically in an appropriate
way. The result will either be an unnormal (UNN), which like the denormal represents a number
very close to zero, or a normal, depending on how significant the effect of the denormal was in the
calculation. In some cases the result will be a NaN.
Following are some procedures for dealing with these values. These procedures are not defined in
the Run-Time Library. If you want to use them, you will need to define them yourself.
The procedure isindef will return 1 (true) if the matrix passed to it contains any NaN’s that are
the indefinite mentioned earlier. The GAUSS missing value code as well as GAUSS scalar error
codes are NaN’s, but this procedure tests only for indefinite:
proc isindef(x);
Fundamentals
Language
retp(not x $/= __INDEFn);
endp;
Be sure to call gausset before calling isindef. gausset will initialize the value of the global
__INDEFn to a platform-specific encoding.
The procedure normal will return a matrix with all denormals and unnormals set to zero.
proc normal(x);
retp(x .* (abs(x) .> 4.19e-307));
endp;
The procedure isinf, will return 1 (true) if the matrix passed to it contains any infinities:
proc isinf(x);
local plus,minus;
plus = __INFp;
10-29
GAUSS User Guide
minus = __INFn;
retp(not x /= plus or not x /= minus);
endp;
Be sure to call gausset before calling isinf. gausset will initialize the values of the globals
__INFn and __INFp to platform specific encodings.
The order in which an expression is evaluated is determined by the precedence of the operators
involved and the order in which they are used. For example, the * and / operators have a higher
precedence than the + and - operators. In expressions that contain these operators, the operand
pairs associated with the * or / operator are evaluated first. Whether * or / is evaluated first
depends on which comes first in the particular expression. For a listing of the precedence of all
operators, see O P, Section 11.7.
The expression
-5+3/4+6*3
is evaluated as
(−5) + (3/4) + (6 ∗ 3)
Within a term, operators of equal precedence are evaluated from left to right.
The term
2ˆ3ˆ7
10-30
Language Fundamentals
is evaluated
(23 )7
In the expression
f1(x)*f2(y)
Expression Evaluation
a+b*c+d (a + (b ∗ c)) + d
-2+4-6*inv(8)/9 ((−2) + 4) − ((6 ∗ inv(8))/9)
Fundamentals
Language
3.14ˆ5*6/(2+sqrt(3)/4) ((3.145 ) ∗ 6)/(2 + (sqrt(3)/4))
-a+b*cˆ2 (−a) + (b ∗ (c2 ))
a+b-c+d-e (((a + b) − c) + d) − e
aˆbˆc*d ((ab )c ) ∗ d
a*b/d*c ((a ∗ b)/d) ∗ c
aˆb+c*d (ab ) + (c ∗ d)
2ˆ4! 2(4!)
2*3! 2 ∗ (3!)
A computer language needs facilities for decision making and looping to control the order in
which computations are done. GAUSS has several kinds of flow control statements.
10-31
GAUSS User Guide
10.8.1 Looping
do loop
also
The scalar expression is any expression that returns a scalar result. The expression will be
evaluated as TRUE if its real part is nonzero and FALSE if it is zero. There is no counter variable
that is automatically incremented in a do loop. If one is used, it must be set to its initial value
before the loop is entered and explicitly incremented or decremented inside the loop.
The following example illustrates nested do loops that use counter variables.
10-32
Language Fundamentals
do while j <= 3;
print space i comma j;;
j = j+1;
endo;
i = i+1;
print;
endo;
1, 1 1, 2 1, 3
2, 1 2, 2 2, 3
3, 1 3, 2 3, 3
4, 1 4, 2 4, 3
Use the relational and logical operators without the dot ‘.’ in the expression that controls a do
loop. These operators always return a scalar result.
Fundamentals
Language
break and continue are used within do loops to control execution flow. When break is
encountered, the program will jump to the statement following the endo. This terminates the loop.
When continue is encountered, the program will jump up to the top of the loop and reevaluate
the while or until expression. This allows you to reiterate the loop without executing any more
of the statements inside the loop:
10-33
GAUSS User Guide
for loop
counter is the literal name of the counter variable. start, stop and step are scalar expressions. start
is the initial value, stop is the final value and step is the increment.
break and continue are also supported by for loops. (For more information, see for in the
GAUSS L R.)
if scalar expression;
.
.
statements
.
.
elseif scalar expression;
.
.
statements
.
.
10-34
Language Fundamentals
else;
.
.
statements
.
.
endif;
The scalar expression is any expression that returns a scalar result. The expression will be
evaluated as TRUE if its real part is nonzero and FALSE if it is zero.
GAUSS will test the expression after the if statement. If it is TRUE, then the first list of
statements is executed. If it is FALSE, then GAUSS will move to the expression after the first
elseif statement, if there is one, and test it. It will keep testing expressions and will execute the
first list of statements that corresponds to a TRUE expression. If no expression is TRUE, then the
list of statements following the else statement is executed. After the appropriate list of statements
is executed, the program will go to the statement following the endif and continue on.
Use the relational and logical operators without the dot ‘.’ in the expression that controls an if or
Fundamentals
Language
elseif statement. These operators always return a scalar result.
One endif is required per if clause. If an else statement is used, there may be only one per if
clause. There may be as many elseif’s as are required. There need not be any elseif’s or any
else statement within an if clause.
The goto and gosub statements control unconditional branching. The target of both a goto and a
gosub is a label.
goto
10-35
GAUSS User Guide
label:
.
.
goto label;
Parameters can be passed with a goto. The number of parameters is limited by available stack
space. This is helpful for common exit routines:
.
.
goto errout("Matrix singular");
.
.
goto errout("File not found");
.
.
errout:
pop errmsg;
errorlog errmsg;
end;
gosub
With a gosub, the address of the gosub statement is remembered and when a return statement is
encountered, the program will resume executing at the statement following the gosub.
Parameters can be passed with a gosub in the same way as a goto. With a gosub it is also
possible to return parameters with the return statement.
Subroutines are not isolated from the rest of your program and the variables referred to between
the label and the return statement can be accessed from other places in your program.
Since a subroutine is only an address marked by a label, there can be subroutines inside of
procedures. The variables used in these subroutines are the same variables that are known inside
the procedure. They will not be unique to the subroutine, but they may be locals that are unique to
10-36
Language Fundamentals
the procedure that the subroutine is in. (For details, see gosub in the GAUSS L
R.)
10.9 Functions
Single line functions that return one item can be defined with the fn statement.
fn area(r) = pi * r * r;
These functions can be called in the same way as intrinsic functions. The above function could be
used in the following program sequence.
diameter = 3;
Fundamentals
Language
radius = 3 / 2;
a = area(radius);
This section lists the general rules of syntax for GAUSS programs.
10.10.1 Statements
10-37
GAUSS User Guide
Column position is not significant. Blank lines are allowed. Inside a statement and outside of
double quotes, the carriage return/line feed at the end of a physical line will be converted to a
space character as the program is compiled.
A statement containing a quoted string can be continued across several lines with a backslash as
follows.
10.10.2 Case
GAUSS does not distinguish between uppercase and lowercase except inside double quotes.
10.10.3 Comments
// This comments out all text between the ’//’ and the end of
// the line
Extraneous spaces are significant in print and lprint statements where the space is a delimiter
between expressions:
print x y z;
10-38
Language Fundamentals
In print and lprint statements, spaces can be used in expressions that are in parentheses:
print (x * y) (x + y);
The names of matrices, strings, procedures, and functions can be up to 32 characters long. The
characters must be alphanumeric or an underscore. The first character must be alphabetic or an
underscore.
10.10.6 Labels
A label is used as the target of a goto or a gosub. The rules for naming labels are the same as for
matrices, strings, procedures, and functions. A label is followed immediately by a colon:
Fundamentals
Language
here:
goto here;
y = x + z;
10-39
GAUSS User Guide
mant,pow = base10(x);
if x =\,= y;
print "x is equal to y";
endif;
y = sqrt(x);
x = { 1 2 3,
3 7 5,
3 7 4,
8 9 5,
6 1 8 };
y = x[3,3];
z = x[1 2:4,1 3];
v = 1 2 3 4 5 6 7 8 9 ;
k = v[3];
j = v[1,6:9];
10-40
Language Fundamentals
x[2,3] returns the element in the second row and the third column of x.
x[1 3 5,4 7] returns the submatrix that is the intersection of rows 1, 3, and 5 and columns 4 and
7.
x[3:5,.] returns the submatrix containing the third through the fifth rows of x.
The indexing operator will take vector arguments for submatrix extraction or submatrix
assignments:
y = x[rv,cv];
y[rv,cv] = x;
Fundamentals
Language
rv and cv can be any expressions returning vectors or matrices. The elements of rv will be used
as the row indices and the elements of cv will be used as the column indices. If rv is a scalar 0, all
rows will be used; if cv is a scalar 0, all columns will be used. If a vector is used in an index
expression, it is illegal to use the space operator or the colon operator on the same side of the
comma as the vector.
In this example, a set of matrix names is assigned to mvec. The name y is indexed from mvec and
passed to varget which will return the global matrix y. The returned matrix is inverted and
assigned to g:
mvec = { x y z a };
i = 2;
g = inv(varget(mvec[i]));
10-41
GAUSS User Guide
The following procedure can be used to index the matrices in mvec more directly:
proc imvec(i);
retp(varget(mvec[i]));
endp;
Then imvec(i) will equal the matrix whose name is in the ith element of mvec.
In the example above, the procedure imvec() was written so that it always operates on the vector
mvec. The following procedure makes it possible to pass in the vector of names being used:
proc get(array,i);
retp(varget(array[i]));
endp;
proc put(x,array,i);
retp(varput(x,array[i]));
endp;
And put(x,mvec,3) will assign x to the 3rd matrix listed in mvec and return a 1 if successful or a
0 if it fails.
It is also possible to index procedures. The ampersand operator (&) is used to return a pointer to a
procedure.
Assume that f1, f2, and f3 are procedures that take a single argument. The following code
defines a procedure fi that will return the value of the ith procedure, evaluated at x.
10-42
Language Fundamentals
proc fi(x,i);
local f;
f = nms[i];
local f:proc;
retp( f(x) );
endp;
fi(x,2) will return f2(x). The ampersand is used to return the pointers to the procedures. nms is
a numeric vector that contains a set of pointers. The local statement is used twice. The first tells
the compiler that f is a local matrix. The ith pointer, which is just a number, is assigned to f. Then
the second local statement tells the compiler to treat f as a procedure from this point on; thus the
subsequent statement f(x) is interpreted as a procedure call.
Fundamentals
Language
10-43
Operators 11
11.1 Element-by-Element Operators
Element-by-element operators share common rules of conformability. Some functions that have
two arguments also operate according to the same rules.
Element-by-element operators handle those situations in which matrices are not conformable
Operators
according to standard rules of matrix algebra. When a matrix is said to be E×E conformable, it
refers to this element-by-element conformability . The following cases are supported:
matrix op matrix
matrix op scalar
scalar op matrix
matrix op vector
vector op matrix
vector op vector
11-1
GAUSS User Guide
z = x + y;
• If x and y are the same size, the operations are carried out corresponding element by
corresponding element:
1 3 2
x= 4 5 1
3 7 4
2 4 3
y= 3 1 4
6 1 2
3 7 5
z= 7 6 5
9 8 6
• If x is a matrix and y is a scalar, or vice versa, then the scalar is operated on with respect to
every element in the matrix. For example, x + 2 will add 2 to every element of x:
1 3 2
x= 4 5 1
3 7 4
y= 2
3 5 4
z= 6 7 3
5 9 6
11-2
Operators
• If x is an N×1 column vector and y is an N×K matrix, or vice versa, the vector is swept
“across” the matrix:
vector matrix
1 −→ 2 4 3
4 −→ 3 1 4
3 −→ 6 1 2
result
3 5 4
7 5 8
9 4 5
• If x is an 1×K column vector and y is an N×K matrix, or vice versa, then the vector is swept
“down” the matrix:
vector 2 4 3
↓ ↓ ↓
Operators
2 4 3
matrix 3 1 4
6 1 2
4 8 6
result 5 5 7
8 5 5
• When one argument is a row vector and the other is a column vector, the result of an
element-by-element operation will be the “table” of the two:
11-3
GAUSS User Guide
row vector 2 4 3 1
3 5 7 6 4
column vector 2 4 6 5 3
5 7 9 8 6
If x and y are such that none of these conditions apply, the matrices are not conformable to these
operations and an error message will be generated.
The following operators work on matrices. Some assume numeric data and others will work on
either character or numeric data.
For details on how matrix conformability is defined for element-by-element operators, see
E--E O, Section 11.1.
+ Addition
y = x + z;
− Subtraction or negation
y = x - z;
y = -k;
11-4
Operators
y = x * z;
When z has the same number of rows as x has columns, this will perform matrix
multiplication (inner product). If x or z are scalar, this performs standard
element-by-element multiplication.
x = b / A;
If A and b are scalars, this performs standard division. If one of the operands is a matrix and
the other is scalar, the result is a matrix the same size with the results of the divisions
between the scalar and the corresponding elements of the matrix. Use ./ for
element-by-element division of matrices.
If b and A are conformable, this operator solves the linear matrix equation
Ax = b
Operators
• If A is a square matrix and has the same number of rows as b, this statement will solve
the system of linear equations using an LU decomposition.
• If A is rectangular with the same number of rows as b, this statement will produce the
least squares solutions by forming the normal equations and using the Cholesky
decomposition to get the solution:
A0 b
x=
A0 A
If trap 2 is set, missing values will be handled with pairwise deletion.
% Modulo division
11-5
GAUSS User Guide
y = x %z;
For integers, this returns the integer value that is the remainder of the integer division of x
by z. If x or z is noninteger, it will first be rounded to the nearest integer. This is an
element-by-element operator.
! Factorial
y = x!;
Computes the factorial of every element in the matrix x. Nonintegers are rounded to the
nearest integer before the factorial operator is applied. This will not work with complex
matrices. If x is complex, a fatal error will be generated.
.* Element-by-element multiplication
y = x .* z;
If x is a column vector, and z is a row vector (or vice versa), the “outer product” or “table” of
the two will be computed. (For comformability rules, see E--E O,
Section 11.1.)
./ Element-by-element division
y = x ./ z;
ˆ Element-by-element exponentiation
y = xˆz;
.ˆ Same as ˆ
11-6
Operators
y = x .*. z;
This results in a matrix in which every element in x has been multiplied (scalar
multiplication) by the matrix z. For example:
x = { 1 2,
3 4 };
z = { 4 5 6,
7 8 9 };
y = x .*. z;
1 2
x=
3 4
4 5 6
z=
7 8 9
4 5 6 8 10 12
7 8 9 14 16 18
y=
12 15 18 16 20 24
21 24 27 28 32 36
Operators
z = x *∼ y;
1 2
x=
3 4
5 6
y=
7 8
5 6 10 12
z=
21 24 28 32
The input matrices x and y must have the same number of rows. The result will have
cols(x) * cols(y) columns.
11-7
GAUSS User Guide
0 Transpose operator
y = x0 ;
The columns of y will contain the same values as the rows of x and the rows of y will
contain the same values as the columns of x. For complex matrices this computes the
complex conjugate transpose.
If an operand immediately follows the transpose operator, the 0 will be interpreted as 0 *.
Thus y = x0 x is equivalent to y = x0 *x.
y = x.0 ;
This is provided primarily as a matrix handling tool for complex matrices. For all matrices,
the columns of y will contain the same values as the rows of x and the rows of y will contain
the same values as the columns of x. The complex conjugate transpose is NOT computed
when you use .0 .
If an operand immediately follows the bookkeeping transpose operator, the .0 will be
interpreted as .0 *. Thus y = x.0 x is equivalent to y = x.0 *x.
| Vertical concatenation
z = x|y;
1 2 3
x=
3 4 5
y= 7 8 9
1 2 3
z= 3 4 5
7 8 9
11-8
Operators
∼ Horizontal concatenation
z = x∼y;
1 2
x=
3 4
5 6
y=
7 8
1 2 5 6
z=
3 4 7 8
For details on how matrix conformability is defined for element-by-element operators, see
E--E O, Section 11.1
Each of these operators has two equivalent representations. Either can be used (for example, < or
lt), depending only upon preference. The alphabetic form should be surrounded by spaces.
Operators
A third form of these operators has a ‘$’ and is used for comparisons between character data and
for comparisons between strings or string arrays. The comparisons are done byte by byte starting
with the lowest addressed byte of the elements being compared.
The equality comparison operators (<=, = =, >=, /=) and their dot equivalents can be used to test
for missing values and the NaN that is created by floating point exceptions. Less than and greater
than comparisons are not meaningful with missings or NaN’s, but equal and not equal are valid.
These operators are sign-insensitive for missings, NaN’s, and zeros.
The string ‘$’ versions of these operators can also be used to test missings, NaN’s and zeros.
Because they do a strict byte-to-byte comparison, they are sensitive to the sign bit. Missings,
NaN’s, and zeros can all have the sign bit set to 0 or 1, depending on how they were generated and
have been used in a program.
11-9
GAUSS User Guide
If the relational operator is NOT preceded by a dot ‘.’, then the result is always a scalar 1 or 0,
based upon a comparison of all elements of x and y. All comparisons must be true for the
relational operator to return TRUE.
if x /= y;
is interpreted as: “if every element of x is not equal to the corresponding element of y”. To check
if two matrices are not identical, use
if not x = = y;
For complex matrices, the = =, /=, .= = and ./= operators compare both the real and imaginary
parts of the matrices; all other relational operators compare only the real parts.
• Less than
z = x < y;
z = x lt y;
z = x $< y;
z = x <= y;
z = x le y;
z = x $<= y;
• Equal to
z = x = = y;
11-10
Operators
z = x eq y;
z = x $= = y;
• Not equal
z = x /= y;
z = x ne y;
z = x $/= y;
z = x >= y;
z = x ge y;
z = x $>= y;
• Greater than
z = x > y;
z = x gt y;
Operators
z = x $> y;
If the relational operator IS preceded by a dot ‘.’, then the result will be a matrix of 1’s and 0’s,
based upon an element-by-element comparison of x and y.
z = x .< y;
z = x .lt y;
z = x .$< y;
11-11
GAUSS User Guide
z = x .<= y;
z = x .le y;
z = x .$<= y;
• Element-by-element equal to
z = x .= = y;
z = x .eq y;
z = x .$= = y;
z = x ./= y;
z = x .ne y;
z = x .$/= y;
z = x .>= y;
z = x .ge y;
z = x .$>= y;
z = x .> y;
z = x .gt y;
z = x .$> y;
11-12
Operators
The logical operators perform logical or Boolean operations on numeric values. On input a
nonzero value is considered TRUE and a zero value is considered FALSE. The logical operators
return a 1 if TRUE and a 0 if FALSE. Decisions are based on the following truth tables:
Complement
X not X
T F
F T
Conjunction
X Y X and Y
T T T
T F F
Operators
F T F
F F F
Disjunction
X Y X or Y
T T T
T F T
F T T
F F F
11-13
GAUSS User Guide
Exclusive Or
X Y X xor Y
T T F
T F T
F T T
F F F
Equivalence
X Y X eqv Y
T T T
T F F
F T F
F F T
For complex matrices, the logical operators consider only the real part of the matrices.
The following operators require scalar arguments. These are the ones to use in if and do
statements:
• Complement
z = not x;
• Conjunction
z = x and y;
• Disjunction
z = x or y;
• Exclusive or
z = x xor y;
11-14
Operators
• Equivalence
z = x eqv y;
If the logical operator is preceded by a dot ‘.’, the result will be a matrix of 1’s and 0’s based upon
an element-by-element logical comparison of x and y:
z = .not x;
• Element-by-element conjunction
z = x .and y;
• Element-by-element disjunction
z = x .or y;
• Element-by-element exclusive or
z = x .xor y;
• Element-by-element equivalence
Operators
z = x .eqv y;
Assignment Operator
y = 3;
11-15
GAUSS User Guide
Comma
clear x,y,z;
y = x[3,5];
y = momentd(x,d);
Period
y = x[.,5];
Space
No extraneous spaces are allowed immediately before or after the comma, or immediately after the
left bracket or before the right bracket.
Spaces are also used in print and lprint statements to separate the separate expressions to be
printed:
11-16
Operators
No extraneous spaces are allowed within expressions in print or lprint statements unless the
expression is enclosed in parentheses:
print (x / 2) (2 * sqrt(x));
Colon
y = x[1:5,.];
Ampersand
The (&) ampersand operator will return a pointer to a procedure (proc), function (fn), or structure
(struct). It is used when passing procedures or functions to other functions, when indexing
procedures, and when initializing structure pointers. (For more information, see I
P, Section 12.5 or S P, Section 16.2.)
Operators
String Concatenation
x = "dog";
y = "cat";
z = x $+ y;
print z;
dogcat
If the first argument is of type string, the result will be of type string. If the first argument is of
type matrix, the result will be of type matrix. Here are some examples:
11-17
GAUSS User Guide
y = 0 $+ "caterpillar";
y = zeros(3,1) $+ "cat";
k = y $+ "fish";
The result will be a 3×1 matrix with each element containing ‘catfish’.
t = "" $+ k[1,1];
z = "dog" $+ k[1,1];
11-18
Operators
x = "dog";
y = "fish";
k = x $| y;
print k;
dog
fish
x = "dog";
y = "fish";
k = x $˜ y;
print k;
dog fish
Operators
by default GAUSS will interpret olsdat as the literal name of the GAUSS data file you want to
create. It will also interpret x as the literal prefix string for the variable names x1 x2 x3 x4.
To get the data set name from a string variable, the substitution operator (ˆ) could be used as
follows:
dataset = "olsdat";
create f1 = ˆdataset with x,4,2;
To get the data set name from a string variable and the variable names from a character vector, use
the following:
11-19
GAUSS User Guide
dataset = "olsdat";
vnames = { age, pay, sex };
create f1 = ˆdataset with ˆvnames,0,2;
ˆvariable name
The following commands are currently supported with the substitution operator (ˆ) in the current
version.
When you use those operators preceded by a ‘.’ (dot operators) with a scalar integer constant,
insert a space between the constant and any following dot operator. Otherwise, the dot will be
interpreted as part of the scalar; that is, the decimal point. For example:
11-20
Operators
let y = 1 2 3;
x = 2.<y;
x = 2.<y;
is interpreted as
x = 2. < y;
and not as
x = 2 .< y;
Be careful when using the dot relational operators (.<, .<=, .= =, ./=, .>, .>=). The same
problem can occur with other dot operators, also. For example:
Operators
let x = 1 1 1;
y = x./2./x;
y = x./2./x;
is interpreted as
y = (x ./ 2.) / x;
11-21
GAUSS User Guide
not
y = (x ./ 2) ./ x;
The second division, then, is handled as a matrix division rather than an element-by-element
division.
The order in which an expression is evaluated is determined by the precedence of the operators
involved and the order in which they are used. For example, the * and / operators have a higher
precedence than the + and − operators. In expressions that contain the above operators, the
operand pairs associated with the * or / operator are evaluated first. Whether * or / is evaluated
first depends on which comes first in the particular expression.
The expression
-5+3/4+6*3
is evaluated as
(-5)+(3/4)+(6*3)
Within a term, operators of equal precedence are evaluated from left to right. The precedence of
all operators, from the highest to the lowest, is listed in the following table:
11-22
Operators
Operators
11-23
Procedures and Keywords 12
Procedures are multiple-line, recursive functions that can have either local or global variables.
Procedures allow a large computing task to be written as a collection of smaller tasks. These
smaller tasks are easier to work with and keep the details of their operation from the other parts of
the program that do not need to know them. This makes programs easier to understand and easier
to maintain.
Any intrinsic command or function may be used in a procedure, as well as any user-defined
function or other procedure. Procedures can refer to any global variable; that is, any variable in the
Procedures
global symbol table that can be shown with the show command. It is also possible to declare local
variables within a procedure. These variables are known only inside the procedure they are defined
in and cannot be accessed from other procedures or from the main level program code.
All labels and subroutines inside a procedure are local to that procedure and will not be confused
with labels of the same name in other procedures.
12-1
GAUSS User Guide
A procedure definition consists of five parts, four of which are denoted by explicit GAUSS
commands:
There is always one proc statement and one endp statement in a procedure definition. Any
statements that come between these two statements are part of the procedure. Procedure
definitions cannot be nested. local and retp statements are optional. There can be multiple
local and retp statements in a procedure definition. Here is an example:
This could be used as a function that takes two matrix arguments and returns three matrices as a
result. For example: is:
{ b,sd,t } = regress(x,y);
12-2
Procedures and Keywords
The proc statement is the procedure declaration statement. The format is:
rets Optional constant, number of values returned by the procedure. Acceptable values here
are 0-1023; the default is 1.
arg# Names that will be used inside the procedure for the arguments that are passed to the
procedure when it is called. There can be 0-1023 arguments. These names will be
known only in the procedure being defined. Other procedures can use the same names,
but they will be separate entities.
The local statement is used to declare local variables. Local variables are variables known only
to the procedure being defined. The names used in the argument list of the proc statement are
always local. The format of the local statement is:
Local variables can be matrices or strings. If :proc, :fn, or :keyword follows the variable name
in the local statement, the compiler will treat the symbol as if it were a procedure, function, or
Procedures
keyword respectively. This allows passing procedures, functions, and keywords to other
procedures. (For more information, see P P P, Section 12.4.
Variables that are global to the system (that is, variables listed in the global symbol table that can
be shown with the show command) can be accessed by any procedure without any redundant
declaration inside the procedure. If you want to create variables known only to the procedure
12-3
GAUSS User Guide
being defined, the names of these local variables must be listed in a local statement. Once a
variable name is encountered in a local statement, further references to that name inside the
procedure will be to the local rather than to a global having the same name. (See clearg, varget,
and varput in the GAUSS L R for ways of accessing globals from within
procedures that have locals with the same name.)
The local statement does not initialize (set to a value) the local variables. If they are not passed
in as parameters, they must be assigned some value before they are accessed or the program will
terminate with a Variable not initialized error message.
All local and global variables are dynamically allocated and sized automatically during execution.
Local variables, including those that were passed as parameters, can change in size during the
execution of the procedure.
Local variables exist only when the procedure is executing and then disappear. Local variables
cannot be listed with the show command.
The maximum number of locals is limited by stack space and the size of workspace memory. The
limiting factor applies to the total number of active local symbols at any one time during
execution. If cat has 10 locals and it calls dog which has 20 locals, there are 30 active locals
whenever cat is called.
There can be multiple local statements in a procedure. They will affect only the code in the
procedure that follows. Therefore, for example, it is possible to refer to a global x in a procedure
and follow that with a local statement that declares a local x. All subsequent references to x
would be to the local x. (This is not good programming practice, but it demonstrates the principle
that the local statement affects only the code that is physically below it in the procedure
definition.) Another example is a symbol that is declared as a local and then declared as a local
procedure or function later in the same procedure definition. This allows doing arithmetic on local
function pointers before calling them. (For more information, see I P, Section
12.5.
The body of the procedure can have any GAUSS statements necessary to perform the task the
procedure is being written for. Other user-defined functions and other procedures can be
referenced as well as any global matrices and strings.
12-4
Procedures and Keywords
GAUSS procedures are recursive, so the procedure can call itself as long as there is logic in the
procedure to prevent an infinite recursion. The process would otherwise terminate with either an
Insufficient workspace memory message or a Procedure calls too deep message,
depending on the space necessary to store the locals for each separate invocation of the procedure.
The return from the procedure is accomplished with the retp statement:
retp;
retp(expression1,expression2,. . .,expressionN);
The retp statement can have multiple arguments. The number of items returned must coincide
with the number of rets in the proc statement.
If the procedure was defined with no items returned, the retp statement is optional. The endp
statement that ends the procedure will generate an implicit retp with no objects returned. If the
procedure returns one or more objects, there must be an explicit retp statement.
There can be multiple retp statements in a procedure, and they can be anywhere inside the body
of the procedure.
Procedures
endp;
An implicit retp statement that returns nothing is always generated here so it is impossible to run
off the end of a procedure without returning. If the procedure was defined to return one or more
objects, executing this implicit return will result in a Wrong number of returns error message
and the program will terminate.
12-5
GAUSS User Guide
dog(i,j,k); /* no returns */
Procedures are called in the same way that intrinsic functions are called. The procedure name is
followed by a list of arguments in parentheses. The arguments must be separated by commas.
dog(ak,4,3);
or
call dog(ak,4,3);
The arguments passed to procedures can be complicated expressions involving calls to other
functions and procedures. This calling mechanism is completely general. For example,
12-6
Procedures and Keywords
y = dog(cat(3*x,bird(x,y))-2,2,1);
is legal.
12.3 Keywords
A keyword, like a procedure, is a subroutine that can be called interactively or from within a
GAUSS program. A keyword differs from a procedure in that a keyword accepts exactly one
string argument, and returns nothing. Keywords can perform many tasks not as easily
accomplished with procedures.
A keyword definition is much like a procedure definition. Keywords always are defined with 0
returns and 1 argument. The beginning of a keyword definition is the keyword statement:
keyword name(strarg);
strarg Name that will be used inside of the keyword for the argument that is passed to the
keyword when it is called. There is always one argument. The name is known only in
the keyword being defined. Other keywords can use the same name, but they will be
separate entities. This will always be a string. If the keyword is called with no
Procedures
characters following the name of the keyword, this will be a null string.
The rest of the keyword definition is the same as a procedure definition. (For more information,
see D P, Section 12.1. Keywords always return nothing. Any retp statements,
if used, should be empty. For example:
12-7
GAUSS User Guide
keyword add(s);
local tok, sum;
if s $=\,= "";
print "The argument is a null string";
retp;
endif;
sum = 0;
do until s $=\,= "";
{ tok, s } = token(s);
sum = sum + stof(tok);
endo;
format /rd 1,2;
print "The sum is: " sum;
endp;
The keyword defined above will print the string argument passed to it. The argument will be
printed enclosed in single quotes.
When a keyword is called, every character up to the end of the statement, excluding the leading
spaces, is passed to the keyword as one string argument. For example, if you type
add 1 2 3 4 5;
12-8
Procedures and Keywords
add;
Procedures
endp;
12-9
GAUSS User Guide
proc myproc(&f1,&f2,x,y);
local f1:proc, f2:fn, z;
z = f1(x,y);
retp(f2(z));
endp;
The procedure myproc takes four arguments. The first is a procedure f1 that has two arguments.
The second is a function f2 that has one argument. It also has two other arguments that must be
matrices or scalars. In the local statement, f1 is declared to be a procedure and f2 is declared to
be a function. They can be used inside the procedure in the usual way. f1 will be interpreted as a
procedure inside myproc, and f2 will be interpreted as a function. The call to myproc is made as
follows:
The ampersand (&) in front of the function or procedure name in the call to myproc causes a
pointer to the function or procedure to be passed. No argument list should follow the name when it
is preceded by the ampersand.
Inside myproc, the symbol that is declared as a procedure in the local statement is assumed to
contain a pointer to a procedure. It can be called exactly like a procedure is called. It cannot be
save’d but it can be passed on to another procedure. If it is to be passed on to another procedure,
use the ampersand in the same way.
This example assumes there are a set of procedures named f1-f5 that are already defined. A 1×5
vector procvec is defined by horizontally concatenating pointers to these procedures. A new
procedure, g(x,i) is then defined to return the value of the ith procedure evaluated at x:
12-10
Procedures and Keywords
proc g(x,i);
local f;
f = procvec[i];
local f:proc;
retp( f(x) );
endp;
The local statement is used twice. The first time, f is declared to be a local matrix. After f has
been set equal to the ith pointer, f is declared to be a procedure and is called as a procedure in the
retp statement.
Procedures can return multiple items, up to 1023. The procedure is defined like this example of a
complex inverse:
Procedures
It can then be called like this:
{ zr,zi } = cminv(xr,xi);
12-11
GAUSS User Guide
Also, a procedure that returns more than one argument can be used as input to another procedure
or function that takes more than one argument:
The two returned matrices from cmmult() are passed directly to cminv() in the statement above.
This is equivalent to the following statements:
{ tr,ti } = cmmult(xr,xi,yr,yi);
{ zr,zi } = cminv(tr,ti);
12-12
Procedures and Keywords
{ xr,xi } = cmcplx(rndn(3,3));
{ yr,yi } = cmcplx(rndn(3,3));
{ wr,wi } = cmmult(yr,yi,cminv(cmmult(cminv(xr,xi),yr,yi)));
When a file containing a procedure definition is run, the procedure is compiled and is then resident
in memory. The procedure can be called as if it were an intrinsic function. If the new command is
executed or you quit GAUSS and exit to the operating system, the compiled image of the
Procedures
procedure disappears and the file containing the procedure definition will have to be compiled
again.
If a procedure contains no global references, that is, if it does not reference any global matrices or
strings and it does not call any user-defined functions or procedures, it can be saved to disk in
compiled form in a .fcg file with the save command, and loaded later with the loadp command
12-13
GAUSS User Guide
whenever it is needed. This will usually be faster than recompiling. For example:
The name of the file will be the same as the name of the procedure, with a .fcg extension. (For
details, see loadp and save in the GAUSS L R.)
All compiled procedures should be saved in the same subdirectory, so there is no question where
they are located when it is necessary to reload them. The loadp path can be set in your startup file
to reflect this. Then, to load in procedures, use
loadp proc1,proc2,proc3;
Procedures that are saved in .fcg files will NOT be automatically loaded. It is necessary to
explicitly load them with loadp. This feature should be used only when the time necessary for the
autoloader to compile the source is too great. Also, unless these procedures have been compiled
with #lineson, debugging will be more complicated.
12-14
Matrices
Sparse
Sparse Matrices 13
The sparse matrix data type stores only the non-zero values of a 2-dimensional sparse matrix,
which makes working with sparse matrices faster and more efficient.
The sparse matrix data type is strongly typed in GAUSS, which means that a variable must be
defined as a sparse matrix variable before it may be used as such. Once a variable has been defined
as a sparse matrix, it may not be used as another data type. Similarly, once a variable has been
used as a matrix, array, or other non-sparse data type, it may not be redefined as a sparse matrix.
To define a global sparse matrix, you may use either the declare or the let command:
13-1
GAUSS User Guide
To define a local sparse matrix inside of a procedure, use an implicit let statement:
As neither let nor declare support the initialization of a sparse matrix at this time, you must
initialize a sparse matrix with an assignment after defining it.
Several new functions have been added to allow you to create and manipulate sparse matrices.
These functions are:
packedToSp Creates a sparse matrix from a packed matrix of non-zero values and
row and column indices.
spBiconjGradSol Solves the system of linear equations Ax=b using the biconjugate
gradient method.
spConjGradSol Solves the system of linear equations Ax=b for symmetric matrices
using the conjugate gradient method.
13-2
Sparse Matrices
Matrices
Sparse
spCreate Creates a sparse matrix from vectors of non-zero values, row
indices, and column indices.
spDenseSubmat Returns a dense submatrix of sparse matrix.
spDiagRvMat Inserts submatrices along the diagonal of a sparse matrix.
spEigv Computes a specified number of eigenvalues and eigenvectors of a
square, sparse matrix.
spEye Creates a sparse identity matrix.
spGetNZE Returns the non-zero values in a sparse matrix, as well as their
corresponding row and column indices.
spLDL Computes the LDL decomposition of a symmetric sparse matrix.
spLU Computes the LU decomposition of a sparse matrix with partial
pivoting.
spNumNZE Returns the number of non-zero elements in a sparse matrix.
spOnes Generates a sparse matrix containing only ones and zeros
spSubmat Returns a sparse submatrix of sparse matrix.
spToDense Converts a sparse matrix to a dense matrix.
spTrTDense Multiplies a sparse matrix transposed by a dense matrix.
spTScalar Multiplies a sparse matrix by a scalar.
spZeros Creates a sparse matrix containing no non-zero values.
See C R, Chapter 33, for detailed information on each command.
Support for the sparse matrix data type has also been added to many matrix functions and
operators. The following is a complete list of the matrix functions and operators that currently
support the new sparse matrix type:
13-3
GAUSS User Guide
13-4
Sparse Matrices
Matrices
Sparse
0
./= .<=
∼ == abs
| .= = cols
* > maxc
.* .> minc
+ >= print
- .>= rows
/ < scalerr
./ .< show
/= <= type
Indexing is also supported for sparse matrices, using the same syntax as matrix indexing.
Note that printing a sparse matrix results in a table of the non-zero values contained in the sparse
matrix, followed by their corresponding row and column indices, respectively.
The types of the returns for the dyadic operators were decided on a case-by-case basis, using the
following general principles:
1. The return type for dyadic operations on two dense arguments is always dense.
2. The return type for dyadic operations on two sparse arguments is always sparse unless the
result is likely to be significantly less sparse than the sparse arguments.
3. The return type for dyadic operations on a dense argument and a sparse argument (regardless
of order) is dense unless the return is likely to be at least as sparse as the sparse argument.
These general principles have led to the following decisions regarding return types (note that only
the cases that are displayed in these tables have been implemented at this point):
13-5
GAUSS User Guide
Element-by-Element Addition
Result = Left Operator Right
dense = sparse + dense
dense = dense + dense
sparse = sparse + sparse
dense = dense + sparse
Element-by-Element Subtraction
Result = Left Operator Right
dense = sparse - dense
dense = dense - dense
sparse = sparse - sparse
dense = dense - sparse
Element-by-Element Multiplication
Result = Left Operator Right
sparse = sparse .* dense
dense = dense .* dense
sparse = sparse .* sparse
sparse = dense .* sparse
Element-by-Element Division
Result = Left Operator Right
sparse = sparse ./ dense
dense = dense ./ dense
dense = sparse ./ sparse
dense = dense ./ sparse
13-6
Sparse Matrices
Matrices
Sparse
Other Numeric Operators
Matrix Multiplication
Result = Left Operator Right
dense = sparse * dense
dense = dense * dense
sparse = sparse * sparse
Linear Solve
Result = Left Operator Right
dense = dense / dense
dense = dense / sparse
Note that at this time, the dense = dense / sparse case is defined only for real data.
When either of its arguments are sparse, the / operator uses a tolerance to determine the result,
which may be read or set using the sysstate function, case 39. The default tolerance is 1e-14.
Relational Operators
Since the results of element-by-element ’dot’ comparison operators depend largely on the kind of
data inputted, there are both both dense-returning and sparse-returning versions of the dot
comparison operators when one or both arguments is a sparse matrix. The regular dot comparison
operators and their alphabetic counterparts always return dense matrices, and there is a new set of
alphabetic dot comparison operators that all return sparse matrices:
13-7
GAUSS User Guide
Since the element-by-element ’non-dot’ comparison operators (= =, /=, <, <=, >, >=) and their
alphabetic counterparts (eq, ne, lt, le, gt, ge) all return scalars, there are no sparse-returning
versions of them.
Horizontal Concatenation
Result = Left Operator Right
dense = dense ∼ dense
sparse = sparse ∼ sparse
Vertical Concatenation
Result = Left Operator Right
dense = dense | dense
sparse = sparse | sparse
13-8
Arrays
N-Dimensional Arrays 14
In GAUSS, internally, matrices and arrays are separate data types. Matrices, which are
2-dimensional objects, are stored in memory in row major order. Therefore, a 3×2 matrix is stored
as follows:
The slowest moving dimension in memory is indexed on the right, and the fastest moving
dimension is indexed on the left. This is true of N-dimensional arrays as well. A 4×3×2 array is
stored in the following way:
A complex N-dimensional array is stored in memory in the same way. Like complex matrices,
complex arrays are stored with the entire real part first, followed by the entire imaginary part.
14-1
GAUSS User Guide
Every N-dimensional array has a corresponding N×1 vector of orders that contains the sizes of
each dimension of the array. This is stored with the array and can be accessed with getorders.
The first element of the vector of orders corresponds to the slowest moving dimension, and the last
element corresponds to the fastest moving dimension (refer to the sectionnameGlossary of Terms
at the end of the chapter for clear definitions of these terms). The vector of orders for a
6×5×4×3×2 array, which has 5 dimensions, is the following 5×1 vector:
6
5
4
3
2
Two terms that are important in working with N-dimensional arrays are “dimension index” and
“dimension number.” A dimension index specifies a dimension based on indexing the vector of
orders. It is a scalar, 1-to-N, where 1 corresponds to the dimension indicated by the first element
of the vector of orders of the array (the slowest moving dimension) and N corresponds to the
dimension indicated by the last element of the vector of orders (the fastest moving dimension).
A dimension number specifies dimensions by numbering them in the same order that one would
add dimensions to an array. In other words, the dimensions of an N-dimensional array are
numbered such that the fastest moving dimension has a dimension number of 1, and the slowest
moving dimension has a dimension number of N.
A 6×5×4×3×2 array has 5 dimensions, so the first element of the vector of orders (in this case, 6)
refers to the size of dimension number 5. Since the index of this element in the vector of orders is
1, the dimension index of the corresponding dimension (dimension number 5) is also 1.
You will find references to both dimension index and dimension number in the documentation for
the functions that manipulate arrays.
There are a number of functions that have been designed to manipulate arrays. These functions
allow you to manipulate a subarray within the array by passing in a locator vector to index any
subarray that comprises a contiguous block of memory within the larger block. A vector of indices
of an N-dimensional array is a [1-to-N]×1 vector of base 1 indices into the array, where the first
element corresponds to the first element in a vector of orders. An N×1 vector of indices locates the
14-2
N-Dimensional Arrays
scalar whose position is indicated by the indices. For a 4×3×2 array x, the 3×1 vector of indices:
3
2
Arrays
1
indexes the [3,2,1] element of x. A 2×1 vector of indices for this 3-dimensional example,
references the 1-dimensional array whose starting location is given by the indices.
Because the elements of the vector of indices are always in the same order (the first element of the
vector of indices corresponds to the slowest moving dimension of the array, the second element to
the second slowest moving dimension, and so on), each unique vector of indices locates a unique
subarray.
In general, an [N-K]×1 vector of indices locates a K-dimensional subarray that begins at the
position indicated by the indices. The sizes of the dimensions of the K-dimensional subarray
correspond to the last K elements of the vector of orders of the N-dimensional array. For a
6×5×4×3×2 array y, the 2×1 vector of indices:
2
5
locates the 4×3×2 subarray in y that begins at [2,5,1,1,1] and ends at [2,5,4,3,2].
Brackets ‘[ ]’ can be used to index N-dimensional arrays in virtually the same way that they are
used to index matrices. Bracketed indexing is slower than the convenience array functions, such as
getarray and setarray; however, it can be used to index non-contiguous elements. In order to
index an N-dimensional array with brackets, there must be N indices located within the brackets,
where the first index corresponds to the slowest moving dimension of the array and the last index
corresponds to the fastest moving dimension.
14-3
GAUSS User Guide
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16
17 18 19 20
21 22 23 24
10
22
13 16
17 20
21 24
14-4
N-Dimensional Arrays
5 6 7
Arrays
17 18 19
The following describes rules for E×E conformability of arrays for operators and functions with
two or more arguments.
• An array is conformable to a matrix only if the array has fewer than 3 dimensions, and the
array and matrix follow the standard rules of E×E conformability.
• Two arrays are E×E conformable if they comply with one of the following requirements:
– The two arrays have the same number of dimensions, and each dimension has the same
size.
– The two arrays have the same number of dimensions, and each of the N-2 slowest
moving dimensions has the same size. In this case, the 2 fastest moving dimensions of
the arrays must follow the E×E comformability rules that apply to matrices.
– Both of the arrays have fewer than 3 dimensions, and they follow the E×E
conformability rules that apply to matrices.
vector of orders N×1 vector of the sizes of the dimensions of an object, where N is the
number of dimensions, and the first element corresponds to the slowest moving
dimension.
14-5
GAUSS User Guide
vector of indices [1-to-N]×1 vector of indices into an array, where the first element
corresponds to the first element in a vector of orders.
dimension number Scalar [1-to-N], where 1 corresponds to the fastest moving dimension
and N to the slowest moving dimension.
dimension index Scalar [1-to-N], where 1 corresponds to the first element of the vector of
orders or vector of indices.
locator [1-to-N]×1 vector of indices into an array used by array functions to locate a
contiguous block of the array.
14-6
Working with Arrays 15
with Arrays
Working
15.1 Initializing Arrays
The use of N-dimensional arrays in GAUSS is an additional tool for reducing development time
and increasing execution speed of programs. There are multiple ways of handling N-dimensional
arrays and using them to solve problems, and these ways sometimes have implications for a
trade-off between speed of execution and development time. We will try to make this clear in this
chapter.
The term “arrays” specifically refers to N-dimensional arrays and must not be confused with
matrices. Matrices and arrays are distinct types even if in fact they contain identical information.
Functions for conversion from one to the other are described below.
There are five basic ways of creating an array depending on how the contents are specified:
15-1
GAUSS User Guide
15.1.1 areshape
areshape is a method for creating an array with specified contents. arrayinit creates an array
filled with a selected scalar value: areshape will do the same, but with a matrix. For example,
given a matrix, areshape will create an array containing multiple copies of that matrix:
x = reshape(seqa(1,1,4),2,2);
ord = 3 | 2 | 2;
a = areshape(x,ord);
print a;
Plane [1,.,.]
1.0000 2.0000
3.0000 4.0000
Plane [2,.,.]
1.0000 2.0000
3.0000 4.0000
Plane [3,.,.]
1.0000 2.0000
3.0000 4.0000
areshape is a fast way to re-dimension a matrix or array already in memory. For example,
suppose we have a GAUSS data set containing panel data and that it’s small enough to be read in
all at once:
15-2
Working with Arrays
panel = areshape(loadd("panel"),5|100|10);
mn = amean(panel,2); /* 5x1x10 array of means */
/*of each panel */
mm = moment(panel,0); /* 5x10x10 array of moments */
/* of each panel */
/*
** vc is a 5x10x10 array of
** covariance matrices
*/
vc = mm / 100 - amult(atranspose(mn,1|3|2),mn);
with Arrays
Working
panel is a 5×100×10 array, and in this context is 5 panels of 100 cases measured on 10 variables.
A random array of any dimension or size can be quickly created using areshape. Thus, for a
10×10×5×3 array:
The quick and dirty method above uses the linear congruential generator, which is fast but doesn’t
have the properties required for serious Monte Carlo work. For series simulation you will need to
use the KM generator:
sd0 = 345678;
ord = { 10, 10, 5, 3 };
{ z,sd0 } = rndKMu(prodc(ord),1,sd0);
y = areshape(z,ord);
15-3
GAUSS User Guide
For computing the log-likelihood of a variance components model of panel data, it is necessary to
expand a T×T matrix into an NT×T array of these matrices. This is easily accomplished using
areshape. For example:
r = areshape(m,3|3|3);
print r;
Plane [1,.,.]
Plane [2,.,.]
Plane [3,.,.]
15.1.2 aconcat
aconcat creates arrays from conformable sets of matrices or arrays. With this function, contents
are completely specified by the user. This example tries three concatenations, one along each
15-4
Working with Arrays
dimension:
rndseed 345678;
x1 = rndn(2,2);
x2 = arrayinit(2|2,1);
/*
** along the first dimension or rows
*/
a = aconcat(x1,x2,1);
with Arrays
Working
print a;
/*
** along the second dimension or columns
*/
a = aconcat(x1,x2,2);
print a;
-0.4300 -0.2878
-0.1327 -0.0573
1.0000 1.0000
1.0000 1.0000
/*
** along the third dimension
*/
a = aconcat(x1,x2,3);
print a;
Plane [1,.,.]
15-5
GAUSS User Guide
-0.4300 -0.2878
-0.1327 -0.0573
Plane [2,.,.]
1.0000 1.0000
1.0000 1.0000
15.1.3 aeye
aeye creates an array in which the principal diagonal of the two trailing dimensions is set to one.
For example:
ord = 2 | 3 | 3;
a = aeye(ord);
print a;
Plane [1,.,.]
Plane [2,.,.]
15.1.4 arrayinit
arrayinit creates an array with all elements set to a specified value. For example:
ord = 3 | 2 | 3;
15-6
Working with Arrays
a = arrayinit(ord,1);
print a;
Plane [1,.,.]
Plane [2,.,.]
with Arrays
Working
1.0000 1.0000 1.0000
Plane [3,.,.]
15.1.5 arrayalloc
arrayalloc creates an array with specified number and size of dimensions without setting
elements to any values. This requires a vector specifying the order of the array. The length of the
vector determines the number of dimensions, and each element determines the size of the
corresponding dimensions. The array will then have to be filled using any of several methods
described later in this chapter.
rndseed 345678;
ord = 3 | 2 | 2;
a = arrayalloc(ord,0);
for i(1,ord[1],1);
a[i,.,.] = rndn(2,3);
endfor;
15-7
GAUSS User Guide
print a;
Plane [1,.,.]
Plane [2,.,.]
Plane [3,.,.]
The second argument in the call to arrayalloc specifies whether the created array is real or
complex. arrayinit creates only real arrays.
15-8
Working with Arrays
The index operator is the slowest way to extract parts of arrays. The specialized functions are the
fastest when the circumstances are appropriate for their use.
with Arrays
Working
15.2.1 index operator
The index operator will put a subarray into an array in a manner analogous to the use of index
operators on matrices:
a = arrayinit(3|2|2,0);
b = arrayinit(3|1|2,1);
a[.,2,.] = b;
print a;
Plane [1,.,.]
0.00000 0.00000
1.0000 1.0000
Plane [2,.,.]
0.00000 0.00000
1.0000 1.0000
Plane [3,.,.]
0.00000 0.00000
1.0000 1.0000
15-9
GAUSS User Guide
As this example illustrates, the assignment doesn’t have to be contiguous. putMatrix and
setMatrix require a contiguous assignment, but for that reason they are faster.
a[1,.,.] = rndn(2,2);
print a;
Plane [1,.,.]
-1.7906502 -0.61038103
1.2586160 -0.47736360
Plane [2,.,.]
0.00000 0.00000
1.00000 1.00000
Plane [3,.,.]
0.00000 0.00000
1.00000 1.00000
The index operator will extract an array from a subarray in a manner analogous to the use of index
operators on matrices:
a = areshape(seqa(1,1,12),3|2|2);
b = a[.,1,.];
print a;
Plane [1,.,.]
1.0000 2.0000
3.0000 4.0000
15-10
Working with Arrays
Plane [2,.,.]
5.0000 6.0000
7.0000 8.0000
Plane [3,.,.]
9.0000 10.000
11.000 12.000
with Arrays
Working
print b;
Plane [1,.,.]
1.0000 2.0000
Plane [2,.,.]
5.0000 6.0000
Plane [3,.,.]
9.0000 10.000
It is important to note that the result is always an array even if it’s a scalar value:
c = a[1,1,1];
print c;
Plane [1,.,.]
1.0000
15-11
GAUSS User Guide
If you require a matrix result, and if the result has one or two dimensions, use arraytomat to
convert to a matrix, or use getMatrix, or getMatrix4D. Or, if the result is a scalar, use
getScalar3D or getScalar4D.
15.2.2 getArray
a = areshape(seqa(1,1,12),3|2|2);
b = getarray(a,2|1);
print a;
Plane [1,.,.]
1.0000 2.0000
3.0000 4.0000
Plane [2,.,.]
5.0000 6.0000
7.0000 8.0000
Plane [3,.,.]
9.0000 10.000
11.000 12.000
print b;
5.0000 6.0000
getArray can only extract a contiguous part of an array. To get non-contiguous parts you must
use the index operator.
15-12
Working with Arrays
15.2.3 getMatrix
If the result is one or two dimensions, getMatrix returns a portion of an array converted to a
matrix. getMatrix is about 20 percent faster than the index operator:
a = areshape(seqa(1,1,12),3|2|2);
b = getMatrix(a,2);
print b;
5.0000 6.0000
with Arrays
Working
7.0000 8.0000
15.2.4 getMatrix4D
This is a specialized version of getMatrix for 4-dimensional arrays. It behaves just like
getMatrix but is dramatically faster for that type of array. The following illustrates the difference
in timing:
a = arrayinit(100|100|10|10,1);
t0 = date;
for i(1,100,1);
for j(1,100,1);
b = a[i,j,.,.];
endfor;
endfor;
t1 = date;
e1 = ethsec(t0,t1);
print e1;
print;
t2=date;
for i(1,100,1);
15-13
GAUSS User Guide
for j(1,100,1);
b = getMatrix4d(a,i,j);
endfor;
endfor;
t3 = date;
e2 = ethsec(t2,t3);
print e2;
print;
print ftostrC(100*((e1-e2)/e1),
"percent difference - %6.2lf%%");
13.000000
5.0000000
These are specialized versions of getMatrix for retrieving scalar elements of 3-dimensional and
4-dimensional arrays, respectively. They behave just like getMatrix, with scalar results, but are
much faster. For example:
a = arrayinit(100|10|10,1);
t0 = date;
for i(1,100,1);
for j(1,10,1);
for k(1,10,1);
b = a[i,j,k];
endfor;
endfor;
endfor;
15-14
Working with Arrays
t1 = date;
e1 = ethsec(t0,t1);
print e1;
print;
t2=date;
for i(1,100,1);
for j(1,10,1);
for k(1,10,1);
b = getscalar3d(a,i,j,k);
endfor;
with Arrays
Working
endfor;
endfor;
t3 = date;
e2 = ethsec(t2,t3);
print e2;
print;
print ftostrC(100*((e1-e2)/e1),
"percent difference - %6.2lf%%");
7.0000000
2.0000000
15.2.6 putArray
putArray enters a subarray, matrix, or scalar into an N-dimensional array and returns the result in
an array. This function is much faster than the index operator, but it requires the part of the array
being assigned to be contiguous:
a = arrayinit(3|2|2,3);
b = putarray(a,2,eye(2));
15-15
GAUSS User Guide
print b;
Plane [1,.,.]
3.0000 3.0000
3.0000 3.0000
Plane [2,.,.]
1.0000 0.00000
0.00000 1.0000
Plane [3,.,.]
3.0000 3.0000
3.0000 3.0000
15.2.7 setArray
a = arrayinit(3|2|2,3);
setarray a,2,eye(2);
print b;
Plane [1,.,.]
3.0000 3.0000
3.0000 3.0000
Plane [2,.,.]
1.0000 0.0000
0.0000 1.0000
15-16
Working with Arrays
Plane [3,.,.]
3.0000 3.0000
3.0000 3.0000
When working with arrays, for loops and do loops may be used in the usual way. In the
following, let Y be an N×1×L array of L time series, X an N×1×K array of K independent
with Arrays
Working
variables, B a K×L matrix of regression coefficients, phi a P×L×L array of garch coefficients,
theta a Q×L×L array of arch coefficients, and omega a L×L symmetric matrix of constants. The
log-likelihood for a multivariate garch BEKK model can be computed using the index operator:
yord = getOrders(Y);
xord = getOrders(X);
gord = getOrders(phi);
aord = getOrders(theta);
r = maxc(P|Q);
E = Y - amult(X,areshape(B,N|K|L));
sigma = areshape(omega,N|L|L);
for i(r+1,N,1);
for j(1,Q,1);
W = amult(theta[j,.,.],
atranspose(E[i-j,.,.],1|3|2));
sigma[i,.,.] = sigma[i,.,.] + amult(W,atranspose(W,1|3|2));
15-17
GAUSS User Guide
endfor;
for j(1,P,1);
sigma[i,.,.] = sigma[i,.,.] + amult(amult(phi[j,.,.],
sigma[i-j,.,.]),phi[j,.,.]);
endfor;
endfor;
sigmai = invpd(sigma);
lndet = ln(det(sigma));
lnl = -0.5*( L*(N-R)*asum(ln(det(sigmai)),1) +
asum(amult(amult(E,sigmai),atranspose(E,1|3|2)),3);
Instead of index operators, the above computation can be done using getArray and setArray:
yord = getOrders(Y);
xord = getOrders(X);
gord = getOrders(phi);
aord = getOrders(theta);
r = maxc(P|Q);
E = Y - amult(X,areshape(B,N|K|L));
sigma = areshape(omega,N|L|L);
for i(r+1,N,1);
for j(1,Q,1);
W = amult(getArray(theta,j),
atranspose(getArray(E,i-j),2|1));
setarray sigma,i,getArray(sigma,i)+
15-18
Working with Arrays
amult(W,atranspose(W,2|1));
endfor;
for j(1,P,1);
setarray sigma,i,getArray(sigma,i)+
areshape(amult(amult(getArray(phi,j),
getArray(sigma,i-j)),getArray(phi,j)),3|3);
endfor;
endfor;
sigmai = invpd(sigma);
with Arrays
Working
lndet = ln(det(sigma));
lnl = -0.5*( L*(N-R)*asum(ln(det(sigmai)),1)+
asum(amult(amult(E,sigmai),atranspose(E,1|3|2)),3)
Putting the two code fragments above into loops that called them a hundred times and measuring
the time, produced the following results:
Thus, the getArray and setArray methods are more than twice as fast.
15.3.1 loopnextindex
Several keyword functions are available in GAUSS for looping with arrays. The problem in the
previous section, for example, can be written using these functions rather than with for loops:
sigind = r + 1;
sigloop:
15-19
GAUSS User Guide
sig0ind = sigind[1];
thetaind = 1;
thetaloop:
sig0ind = sig0ind - 1;
W = amult(getArray(theta,thetaind),
atranspose(getArray(E,sig0ind),2|1));
setarray sigma,sigind,getArray(sigma,sigind)+
amult(W,atranspose(W,2|1));
loopnextindex thetaloop,thetaind,aord;
sig0ind = sigind;
phiind = 1;
philoop:
sig0ind[1] = sig0ind[1] - 1;
setarray sigma,sigind,getArray(sigma,sigind)+
areshape(amult(amult(getArray(phi,phiind),
getArray(sigma,sig0ind)),
getArray(phi,phiind)),3|3);
loopnextindex philoop,phiind,gord;
loopnextindex sigloop,sigind,sigord;
The loopnextindex function in this example isn’t faster than the for loop used in the previous
section primarily because the code is looping only through the first dimension in each loop. The
advantages of loopnextindex, previousindex, nextindex, and walkindex are when the
code is looping through the higher dimensions of a highly dimensioned array. In this case, looping
through an array can be very complicated and difficult to manage using for loops.
loopnextindex can be faster and more useful.
The next example compares two ways of extracting a subarray from a 5-dimensional array:
ord = 3|3|3|3|3;
a = areshape(seqa(1,1,prodc(ord)),ord);
b = eye(3);
for i(1,3,1);
15-20
Working with Arrays
for j(1,3,1);
for k(1,3,1);
setarray a,i|j|k,b;
endfor;
endfor;
endfor;
ind = { 1,1,1 };
loopi:
setarray a,ind,b;
loopnextindex loopi,ind,ord;
with Arrays
Working
Calling each loop 10,000 times and measuring the time each takes, we get
In other words, loopnextindex is about four times faster, a very significant difference.
15.4.1 atranspose
a = areshape(seqa(1,1,12),2|3|2);
print a;
Plane [1,.,.]
15-21
GAUSS User Guide
1.0000 2.0000
3.0000 4.0000
5.0000 6.0000
Plane [2,.,.]
7.0000 8.0000
9.0000 10.000
11.000 12.000
/*
** swap 2nd and 3rd dimension
*/
print atranspose(a,1|3|2);
Plane [1,.,.]
Plane [2,.,.]
/*
** swap 1st and 3rd dimension
*/
print atranspose(a,3|2|1);
Plane [1,.,.]
1.0000 7.0000
3.0000 9.0000
5.0000 11.000
15-22
Working with Arrays
Plane [2,.,.]
2.0000 8.0000
4.0000 10.000
6.0000 12.000
/*
** move 3rd into the front
*/
with Arrays
Working
print atranspose(a,3|1|2);
Plane [1,.,.]
Plane [2,.,.]
15.4.2 amult
This function performs a matrix multiplication on the last two trailing dimensions of an array. The
leading dimensions must be strictly conformable, and the last two trailing dimensions must be
conformable in the matrix product sense. For example:
a = areshape(seqa(1,1,12),2|3|2);
b = areshape(seqa(1,1,16),2|2|4);
c = amult(a,b);
print a;
Plane [1,.,.]
15-23
GAUSS User Guide
1.0000 2.0000
3.0000 4.0000
5.0000 6.0000
Plane [2,.,.]
7.0000 8.0000
9.0000 10.000
11.000 12.000
print b;
Plane [1,.,.]
Plane [2,.,.]
print c;
Plane [1,.,.]
Plane [2,.,.]
15-24
Working with Arrays
Suppose we have a matrix of data sets, a 2×2 matrix of 100×5 data sets that we’ve stored in a
2×2×100×5 array called x. The moment matrices of these data sets can easily and quickly be
computed using atranspose and amult:
vc = amult(atranspose(x,1|2|4|3),x);
These functions compute the means, minimums, and maximums, respectively, across a dimension
with Arrays
Working
of an array. The size of the selected dimension of the resulting array is shrunk to one and contains
the means, minimums, or maximums depending on the function called. For example:
a = areshape(seqa(1,1,12),2|3|2);
print a;
Plane [1,.,.]
1.0000 2.0000
3.0000 4.0000
5.0000 6.0000
Plane [2,.,.]
7.0000 8.0000
9.0000 10.000
11.000 12.000
/*
** compute means along third dimension
*/
print amean(a,3);
Plane [1,.,.]
15-25
GAUSS User Guide
4.0000 5.0000
6.0000 7.0000
8.0000 9.0000
/*
** print means along the second dimension, i.e.,
** down the columns
*/
print amean(a,2);
Plane [1,.,.]
3.0000 4.0000
Plane [2,.,.]
9.0000 10.000
/*
** print the minimums down the columns
*/
print amin(a,2);
Plane [1,.,.]
1.0000 2.0000
Plane [2,.,.]
7.0000 8.0000
/*
** print the maximums along the third dimension
*/
15-26
Working with Arrays
print amax(a,3);
Plane [1,.,.]
7.0000 8.0000
9.0000 10.000
11.000 12.000
15.4.4 getDims
with Arrays
Working
This function returns the number of dimensions of an array:
a = arrayinit(4|4|5|2,0);
print getdims(a);
4.00
15.4.5 getOrders
This function returns the sizes of each dimension of an array. The length of the vector returned by
getOrders is the dimension of the array:
a = arrayinit(4|4|5|2,0);
print getOrders(a);
4.00
4.00
5.00
2.00
15-27
GAUSS User Guide
15.4.6 arraytomat
a = arrayinit(2|2,0);
b = arraytomat(a);
type(a);
21.000
type(b);
6.0000
15.4.7 mattoarray
b = rndn(2,2);
a = mattoarray(b);
type(b);
6.0000
type(a);
21.000
Many of the GAUSS functions have been re-designed to work with arrays. There are two general
approaches to this implementation. There are exceptions, however, and you are urged to refer to
the documention if you are not sure how a particular GAUSS function handles array input.
15-28
Working with Arrays
In the first approach, the function returns an element-by-element result that is strictly conformable
to the input. For example, cdfnc returns an array of identical size and shape to the input array:
a = areshape(seqa(-2,.5,12),2|3|2);
b = cdfnc(a);
print b;
Plane [1,.,.]
0.9772 0.9331
with Arrays
Working
0.8413 0.6914
0.5000 0.3085
Plane [2,.,.]
0.1586 0.0668
0.0227 0.0062
0.0013 0.0002
In the second approach, which applies generally to GAUSS matrix functions, the function operates
on the matrix defined by the last two trailing dimensions of the array. Thus, given a 5×10×3 array,
moment returns a 5×3×3 array of five moment matrices computed from the five 10×3 matrices in
the input array.
Only the last two trailing dimensions matter; i.e., given a 2×3×4×5×10×6 array, moment returns a
2×3×4×5×6×6 array of moment matrices.
For example, in the following the result is a 2×3 array of 3×1 vectors of singular values of a 2×3
array of 6×3 matrices:
a = areshape(seqa(1,1,108),2|3|6|3);
b=svds(a);
print b;
15-29
GAUSS User Guide
Plane [1,1,.,.]
45.894532
1.6407053
1.2063156e-015
Plane [1,2,.,.]
118.72909
0.63421188
5.8652600e-015
Plane [1,3,.,.]
194.29063
0.38756064
1.7162751e-014
Plane [2,1,.,.]
270.30524
0.27857175
1.9012118e-014
Plane [2,2,.,.]
346.47504
0.21732995
1.4501098e-014
Plane [2,3,.,.]
422.71618
0.17813229
1.6612287e-014
It might be tempting to conclude from this example that, in general, a GAUSS function’s behavior
15-30
Working with Arrays
on the last two trailing dimensions of an array is strictly analogous to the GAUSS function’s
behavior on a matrix. This may be true with some of the functions, but not all. For example, the
GAUSS meanc function returns a column result for matrix input. However, the behavior for the
GAUSS amean function is not analogous. This function takes a second argument that specifies on
which dimension the mean is to be taken. That dimension is then collapsed to a size of 1. Thus:
a = areshape(seqa(1,1,24),2|3|4);
print a;
Plane [1,.,.]
with Arrays
Working
1.000 2.000 3.000 4.000
5.000 6.000 7.000 8.000
9.000 10.000 11.000 12.000
Plane [2,.,.]
/*
** means computed across rows
*/
b = amean(a,1);
print b;
Plane [1,.,.]
2.500
6.500
10.500
Plane [2,.,.]
14.500
15-31
GAUSS User Guide
18.500
22.500
/*
** means computed down columns
*/
c = amean(a,2);
print c;
Plane [1,.,.]
Plane [2,.,.]
/*
** means computed along 3rd dimension
*/
d = amean(a,3);
print d;
Plane [1,.,.]
Suppose we have N cases observed at T times. Let yit be an observation on a dependent variable
for the ith case at time t, Xit an observation of k independent variables for the ith case at time t, B, a
15-32
Working with Arrays
is a variance components model where µi is a random error term uncorrelated with it , but which is
correlated within cases. This implies an NT×NT residual moment matrix that is block diagonal
with N T×T moment matrices with the following form:
with Arrays
Working
σµ 2
σµ + σ . . .
2 2
σ2µ
.. .. .. ..
. . . .
σ2µ σ2µ . . . σ2µ + σ2
Using GAUSS arrays, we can compute the log-likelihood of this model without resorting to do
loops. Let Y be a 100×3×1 array of observations on the dependent variable, and X a 100×3×5
array of observations on the independent variables. Further let B be a 5×1 vector of coefficients,
and sigu and sige be the residual variances of µ and respectively. Then, in explicit steps we
compute
N = 100;
T = 3;
K = 5;
15-33
GAUSS User Guide
All of this can be made more efficient by nesting statements, which eliminates copying of
temporary intervening arrays to local arrays. It is also useful to add a check for the positive
definiteness of sigma:
N = 100;
T = 3;
K = 5;
const = -0.5*N*T*ln(2*pi);
oldt = trapchk(1);
trap 1,1;
sigmai = invpd(sigu*ones(T,T)+sige*eye(T));
trap oldt,1;
if not scalmiss(sigmai);
E = Y - amult(X,areshape(B,N|K|1));
lnl = const + 0.5*N*ln(detl)-
0.5*asum(amult(amult(atranspose(E,1|3|2),
areshape(sigmai,N|T|T)),E),3);
else;
lnl = error(0);
endif;
15-34
Working with Arrays
15.7 Appendix
This is an incomplete list of special functions for working with arrays. Many GAUSS functions
have been modified to handle arrays and are not listed here. For example, cdfnc computes the
complement of the Normal cdf for each element of an array just as it would for a matrix. See the
documentation for these GAUSS functions for information about their behavior with arrays.
with Arrays
Working
aeye Create an array of identity matrices.
15-35
GAUSS User Guide
loopnextindex Increment an index vector to the next logical index and jump to the
specified label if the index did not wrap to the beginning.
15-36
Structures 16
Structures
16.1 Basic Structures
The list of members can include scalars, arrays, matrices, strings, and string arrays, as well as
other structures. As a type, scalars are unique to structures and don’t otherwise exist.
For example, the following defines a structure containing the possible contents:
struct generic_example {
scalar x;
matrix y;
16-1
GAUSS User Guide
string s1;
string array s2
struct other_example t;
};
A useful convention is to put the structure definition into a file with a .sdf extension. Then, for
any command file or source code file that requires this definition, put
#include filename.sdf
For example:
#include example.sdf
These statements create structure definitions that persist until the workspace is cleared. They do
not create structures, only structure-type definitions. The next section describes how to create an
instance of a structure.
For example:
#include example.sdf
struct generic_example p0;
16-2
Structures
p0.x = rndn(20,3);
mn = meanc(p0.x);
Global structures are initialized at compile time. Each member of the structure is initialized
according to the following schedule:
Structures
scalar 0, a scalar zero
matrix , an empty matrix with zero rows and zero columns
array 0, a 1-dimensional array set to zero
string ””, a null string
string array ””, a 1×1 string array set to null
If a global already exists in memory, it will not be reinitialized. It may be the case in your program
that when it is rerun, the global variables may need to be reset to default values. That is, your
program may depend on certain members of a structure being set to default values that are set to
some other value later in the program. When you rerun this program, you will want to reinitialize
the global structure. To do this, make an assignment to at least one of the members. This can be
made convenient by writing a procedure that declares a structure and initializes one of its members
to a default value, and then returns it. For example:
/* ds.src */
#include ds.sdf
proc dsCreate;
16-3
GAUSS User Guide
struct DS d0;
d0.dataMatrix = 0;
retp(d0);
endp;
Calling this function after declaring an instance of the structure will ensure initialization to default
values each time your program is run:
struct DS d0;
d0 = dsCreate;
Local structures, which are structures defined inside procedures, are initialized at the first
assignment. The procedure may have been written in such a way that a subset of structures are
used an any one call, and in that case time is saved by not initializing the unused structures. They
will be initialized to default values only when the first assignment is made to one of its members.
#include ds.sdf
struct DS p0;
p0 = reshape(dsCreate,5,1);
This creates a 5×1 vector of instances of DS structures, with all of the members initialized to
default values.
When the instance members have been set to some other values, reshape will produce multiple
copies of that instance set to those values.
16-4
Structures
#include trade.sdf
struct option p0,p1,p2;
p0 = optionCreate;
p1 = optionCreate;
p2 = p1 | p0;
Structure indexing may be used to reference a particular element in a structure array. The syntax
follows that of matrix indexing. For example, given the following structure definition:
struct example1 {
matrix x;
matrix y;
string str;
Structures
};
In this example, e1a and e1b are concatenated to create a 2×1 array of example1 structures that is
assigned back to e1a. Then the y member of the [2,1] element of e1a is set to a random matrix.
Indexing of structure arrays can occur on multiple levels. For example, let’s define the following
structures:
struct example3 {
16-5
GAUSS User Guide
matrix w;
string array sa;
};
struct example2 {
matrix z;
struct example3 e3;
};
struct example1 {
matrix x;
matrix y;
string str;
struct example2 e2;
};
Let’s assume that we have an example1 structure e1 like the one displayed in Figure 16.1. We
could then index the structure as follows:
r = e1.e2[3,1].e3[2,1].w
You can also use indexing to reference the structure itself, rather than a member of that structure:
e3tmp = e1.e2[3,1].e3[.,1];
16-6
Structures
e1
matrix x
matrix y
string str
struct example2 e2
e2 is a 3x1 array
Structures
e3 is a 2x1 array e3 is a 1x1 array e3 is a 3x1 array
matrix w matrix w
string aray sa string aray sa
16-7
GAUSS User Guide
In this case, e3tmp would be an array of 3×1 example3 structures, since the [3,1] member of
e1.e2 contains a 3×1 array of example3 structures.
It is important to remember, however, that when indexing a structure array on multiple levels, only
the final index may resolve to an array of structures. For example:
e3tmp = e1.e2[.,1].e3[2,1];
Instances and vectors or matrices of instances of structures can be saved in a file on the disk, and
later loaded from the file onto the disk. The syntax for saving an instance to the disk is
ret = savestruct(instance,filename);
For example:
#include ds.sdf
struct DS p0;
p0 = reshape(dsCreate,2,1);
retc = saveStruct(p2,"p2");
This saves the vector of instances in a file called p2.fsr. retc will be zero if the save was
successful; otherwise, nonzero.
16-8
Structures
For example:
#include trade.sdf;
struct DS p3;
{ p3, retc } = loadstruct("p2","ds");
Structures
Structures or members of structures can be passed to procedures. When a structure is passed as an
argument to a procedure, it is passed by value. The structure becomes a local copy of the structure
that was passed. The data in the structure is not duplicated unless the local copy of the structure
has a new value assigned to one of its members. Structure arguments must be declared in the
procedure definition:
struct rectangle {
matrix ulx;
matrix uly;
matrix lrx;
matrix lry;
};
Local structures are defined using a struct statement inside the procedure definition:
16-9
GAUSS User Guide
A structure pointer is a separate data type that contains the address of a structure and is used to
reference that structure.
struct example_struct {
matrix x;
matrix y;
};
However, at this point, esp is not yet pointing at anything. It has only been defined to be the kind
of pointer that points at example_struct structures. To set it to point at a particular structure
instance, we must first create the structure instance:
16-10
Structures
and then we can set esp to point at es by setting esp to the address of es:
esp = &es;
copies the contents of the structure that esp is pointing at (i.e., the contents of es) to es2. It is the
same as
Structures
struct example_struct es2;
es2 = es;
To reference a member of a structure, we use a “dot” syntax. For example, we might use the
following code to set the x member of es.
es.x = rndn(3,3);
To reference a member of a structure using a pointer to that structure, we use an “arrow” syntax.
For example, we might use the following code to set the x member of es using the pointer esp:
esp->x = rndn(10,5);
16-11
GAUSS User Guide
This code will modify es, since esp is merely a pointer to es.
struct example_struct_2 {
matrix z;
struct example_struct *ep;
};
sp1->sp2->x;
nor
s.sp1->x;
will ever be valid (sp1 and sp2 are assumed to be structure pointers, s a structure instance, and x a
matrix). The “arrow” (->) will only be valid if it is used for the first (or furthest left) dereference,
as in:
sp1->st.x;
At this point we do not support indexing of structure pointers. Thus, a structure pointer should
point at a scalar structure instance, not a matrix of structures. However, you may index members
of that scalar structure instance. So, for example, let us suppose that you defined the following
structure types:
struct sb {
matrix y;
matrix z;
16-12
Structures
};
struct sa {
matrix x;
struct structb s;
};
and then created an instance of an sa structure, a0, setting a0.s to a 3×2 matrix of sb structures.
The following would be legal:
struct sa *sap
sap = &a0;
sap->s[3,1].y = rndn(3,3);
Structures
Structure pointers are especially useful in cases where structures are passed into and out of
procedures. If a procedure takes a structure as an argument and modifies any members of that
structure, then it makes a local copy of the entire structure before modifying it. Thus if you want
to have the modified copy of the structure after running the procedure, you need to pass the
structure out of the procedure as one of its return arguments. For example:
struct example_struct {
matrix x;
matrix y;
matrix z;
};
16-13
GAUSS User Guide
es1.x = rndn(1000,100);
es1.y = rndn(1000,1);
es1 = product(es1);
In this example, the structure es1 is passed into the procedure, copied and modified. The modified
structure is then passed out of the procedure and assigned back to es1.
Structure pointers allow you to avoid such excessive data copying and eliminate the need to pass a
structure back out of a procedure in cases like this. When you pass a structure pointer into a
procedure and then modify a member of the structure that it references, the actual structure is
modified rather than a local copy of it. Thus there is no need to pass the modifed structure back
out of the procedure. For example, the above example could be accomplished using structure
pointers as follows:
struct example_struct {
matrix x;
matrix y;
matrix z;
};
es1p = &es1;
es1.x = rndn(1000,100);
es1.y = rndn(1000,1);
product(es1p);
In this case, the procedure modifies the structure es1, which es1p is pointing at, instead of a local
copy of the structure.
16-14
Structures
There are three common types of structures that will be found in the GAUSS Run-Time Library
and applications.
The DS and PV structures are defined in the GAUSS Run-Time Library. Their definitions are
found in ds.sdf and pv.sdf, respectively, in the src source code subdirectory.
Before structures, many procedures in the Run-Time Library and all applications had global
variables serving a variety of purposes, such as setting and altering defaults. Currently, these
variables are being entered as members of “control” structures.
The DS structure, or “data” structure, is a very simple structure. It contains a member for each
GAUSS data type. The following is found in ds.sdf:
Structures
struct DS {
scalar type;
matrix dataMatrix;
array dataArray;
string dname;
string array vnames;
};
This structure was designed for use by the various optimization functions in GAUSS, in particular,
sqpSolvemt, as well as a set of gradient procedures, gradmt, hessmt, et al.
These procedures all require that the user provide a procedure computing a function (to be
optimized or take the derivative of, etc.), which takes the DS structure as an argument. The
Run-Time Library procedures such as sqpSolvemt take the DS structure as an argument and
pass it on to the user-provided procedure without modification. Thus, the user can put into that
structure whatever might be needed as data in the procedure.
16-15
GAUSS User Guide
#include ds.sdf
struct DS d0;
d0 = dsCreate;
The PV structure, or parameter vector structure, is used by various optimization, modelling, and
gradient procedures, in particular sqpSolvemt, for handling the parameter vector. The GAUSS
Run-Time Library contains special functions that work with this structure. They are prefixed by
“pv” and defined in pv.src. These functions store matrices and arrays with parameters in the
structure, and retrieve various kinds of information about the parameters and parameter vector
from it.
The various procedures in the Run-Time Library and applications for optimization, modelling,
derivatives, etc., all require a parameter vector. Parameters in complex models, however, often
come in matrices of various types, and it has been the responsibility of the programmer to generate
the parameter vector from the matrices and vice versa. The PV procedures make this problem
much more convenient to solve.
The typical situation involves two parts: first, “packing” the parameters into the PV structure,
which is then passed to the Run-Time Library procedure or application; and second,
“unpacking” the PV structure in the user-provided procedure for use in computing the objective
function. For example, to pack parameters into a PV structure:
#include sqpsolvemt.sdf
/* starting values */
16-16
Structures
struct PV p0;
p0 = pvPack(pvCreate,b0,"b0");
p0 = pvPack(p0,garch,"garch");
p0 = pvPack(p0,arch,"arch");
p0 = pvPack(p0,omega,"omega");
/* data */
z = loadd("tseries");
struct DS d0;
d0.dataMatrix = z;
Next, in the user-provided procedure for computing the objective function, in this case minus the
log-likelihood, the parameter vector is unpacked:
Structures
proc ll(struct PV p0, struct DS d0);
local b0,garch,arch,omega,p,q,h,u,vc,w;
b0 = pvUnpack(p0,"b0");
garch = pvUnpack(p0,"garch");
arch = pvUnpack(p0,"arch");
omega = pvUnpack(p0,"omega");
p = rows(garch);
q = rows(arch);
u = d0.dataMatrix - b0;
vc = moment(u,0)/rows(u);
w = omega + (zeros(q,q) | shiftr((u.*ones(1,q))’,
seqa(q-1,-1,q))) * arch;
h = recserar(w,vc*ones(p,1),garch);
logl = -0.5 * ((u.*u)./h + ln(2*pi) + ln(h));
retp(logl);
16-17
GAUSS User Guide
endp;
Masked Matrices
The pvUnpack function unpacks parameters into matrices or arrays for use in computations. The
first argument is a PV structure containing the parameter vector. Sometimes the matrix or vector is
partly parameters to be estimated (that is, a parameter to be entered in the parameter vector) and
partly fixed parameters. To distinguish between estimated and fixed parameters, an additional
argument is used in the packing function called a “mask”, which is strictly conformable to the
input matrix. Its elements are set to 1 for an estimated parameter and 0 for a fixed parameter. For
example:
p0 = pvPackm(p0,.1*eye(3),"theta",eye(3));
Here just the diagonal of a 3×3 matrix is added to the parameter vector.
When this matrix is unpacked, the entire matrix is returned with current values of the parameters
on the diagonal:
print pvUnpack(p0,"theta");
Symmetric Matrices
Symmetric matrices are a special case because even if the entire matrix is to be estimated, only the
nonredundant portion is to be put into the parameter vector. Thus, for them there are special
procedures. For example:
16-18
Structures
vc = { 1 .6 .4, .6 1 .2, .4 .2 1 };
p0 = pvPacks(p0,vc,"vc");
There is also a procedure for masking in case only a subset of the nonredundant elements are to be
included in the parameter vector:
vc = { 1 .6 .4, .6 1 .2, .4 .2 1 };
mask = { 1 1 0, 1 1 0, 0 0 1 };
p0 = pvPacksm(p0,vc,"vc",mask);
Fast Unpacking
When unpacking matrices using a matrix name, pvUnpack has to make a search through a list of
Structures
names, which is relatively time-consuming. This can be alleviated by using an index rather than a
name in unpacking. To do this, though, requires using a special pack procedure that establishes the
index:
p0 = pvPacki(p0,b0,"b0",1);
p0 = pvPacki(p0,garch,"garch",2);
p0 = pvPacki(p0,arch,"arch",3);
p0 = pvPacki(p0,omega,"omega",4);
Now they may be unpacked using the index number:
b0 = pvUnpack(p0,1);
garch = pvUnpack(p0,2);
arch = pvUnpack(p0,3);
omega = pvUnpack(p0,4);
When packed with an index number, they may be unpacked either by index or by name, but
unpacking by index is faster.
16-19
GAUSS User Guide
pvList
This procedure generates a list of the matrices or arrays packed into the structure:
p0 = pvPack(p0,b0,"b0");
p0 = pvPack(p0,garch,"garch");
p0 = pvPack(p0,arch,"arch");
p0 = pvPack(p0,omega,"omega");
print pvList(p0);
b0
garch
arch
omega
pvLength
print pvLength(p0);
6.0000
pvGetParNames
print pvGetParNames(p0);
16-20
Structures
b0[1,1]
garch[1,1]
garch[2,1]
arch[1,1]
arch[2,1]
omega[1,1]
pvGetParVector
print pvGetParVector(p0);
1.0000
0.1000
0.1000
Structures
0.1000
0.1000
1.0000
pvPutParVector
This procedure replaces the parameter vector with the one in the argument:
1.5000
0.2000
0.2000
0.3000
0.3000
0.8000
16-21
GAUSS User Guide
pvGetIndex
This procedure returns the indices in the parameter vector of the parameters in a matrix. These
indices are useful when setting linear constraints or bounds in sqpSolvemt. Bounds, for example,
are set by specifying a K×2 matrix where K is the length of the parameter vector and the first
column are the lower bounds and the second the upper bounds. To set the bounds for a particular
parameter, then, requires knowing where that parameter is in the parameter vector. This
information can be found using pvGetIndex. For example:
Another important class of structures is the “control” structure. Applications developed before
structures were introduced into GAUSS typically handled some program specifications by the use
of global variables which had some disadvantages, in particular, preventing the nesting of calls to
procedures.
Currently, the purposes served by global variables are now served by the use of a control structure.
For example, for sqpSolvemt:
struct sqpSolvemtControl {
matrix A;
matrix B;
matrix C;
matrix D;
scalar eqProc;
scalar ineqProc;
16-22
Structures
matrix bounds;
scalar gradProc;
scalar hessProc;
scalar maxIters;
scalar dirTol;
scalar CovType;
scalar feasibleTest;
scalar maxTries;
scalar randRadius;
scalar trustRadius;
scalar seed;
scalar output;
scalar printIters;
matrix weights;
};
Structures
16.4 sqpSolvemt
sqpSolvemt is a procedure in the GAUSS Run-Time Library that solves the general nonlinear
programming problem using a Sequential Quadratic Programming descent method, that is, it
solves
min f (θ)
subject to
Aθ = B linear equality
Cθ>=D linear inequality
H(θ) = 0 nonlinear equality
G(θ)>=0 nonlinear inequality
θlb <=θ<=θub bounds
16-23
GAUSS User Guide
The linear and bounds constraints are redundant with respect to the nonlinear constraints, but are
treated separately for computational convenience.
The call to sqpSolvemt has four input arguments and one output argument:
out = SQPsolveMT(&fct,P,D,C);
The first input argument is a pointer to the objective function to be minimized. The procedure
computing this objective function has two arguments: a PV structure containing the start values,
and a DS structure containing data, if any. For example:
Note that this procedure returns a vector rather than a scalar. When the objective function is a
properly defined log-likelihood, returning a vector of minus log-probabilities permits the
calculation of a QML covariance matrix of the parameters.
16-24
Structures
C an sqpSolvemtControl structure
The DS structure is optional. sqpSolvemt passes this argument on to the user-provided procedure
that &fct is pointing to without modification. If there is no data, a default structure can be passed
to it.
sqpSolvemtControl Structure
Structures
vector of the parameters.
16-25
GAUSS User Guide
Bounds 1×2 or K×2 matrix, bounds on parameters. If 1×2 all parameters have same
bounds.
Default = -1e256 1e256 .
GradProc scalar, pointer to a procedure that computes the gradient of the function
with respect to the parameters. When such a procedure has been provided,
it has two input arguments, instances of PV and DS structures, and one
output argument, the derivatives. If the function procedure returns a scalar,
the gradient procedure returns a 1×K row vector of derivatives. If function
procedure turns an N×1 vector, the gradient procedure returns an N×K
matrix of derivatives.
This procedure may compute a subset of the derivatives. sqpSolvemt will
compute numerical derivatives for all those elements set to missing values
in the return vector or matrix.
Default = .; i.e., no gradient procedure has been provided.
HessProc scalar, pointer to a procedure that computes the Hessian; i.e., the matrix of
second order partial derivatives of the function with respect to the
parameters. When such a procedure has been provided, it has two input
arguments, instances of PV and DS structures, and one output argument, a
vector of computed inequality constraints.
Default = .; i.e., no Hessian procedure has been provided.
Whether the objective function procedure returns a scalar or vector, the
Hessian procedure must return a K×K matrix. Elements set to missing
values will be computed numerically by sqpSolvemt.
16-26
Structures
FeasibleTest scalar, if nonzero, parameters are tested for feasibility before computing
function in line search. If function is defined outside inequality boundaries,
then this test can be turned off. Default = 1.
seed scalar, if nonzero, seeds random number generator for random search,
otherwise time in seconds from midnight is used.
trustRadius scalar, radius of the trust region. If scalar missing, trust region not applied.
The trust sets a maximum amount of the direction at each iteration. Default
= .001.
Structures
16.4.2 Output Argument
struct SQPsolveMTOut {
struct PV par;
scalar fct;
struct SQPsolveMTLagrange lagr;
scalar retcode;
matrix moment;
matrix hessian;
matrix xproduct;
};
16-27
GAUSS User Guide
0 normal convergence
1 forced exit
2 maximum number of iterations exceeded
3 function calculation failed
4 gradient calculation failed
5 Hessian calculation failed
6 line search failed
7 error with constraints
8 function complex
9 feasible direction couldn’t be found
16-28
Structures
16.4.3 Example
Define
Y = Λη + θ
where Λ is a K×L matrix of loadings, η an L×1 vector of unobserved “latent” variables, and θ a
K×1 vector of unobserved errors. Then
Σ = ΛΦΛ0Ψ
where Φ is the L×L covariance matrix of the latent variables, and Ψ is the K×K covariance matrix
Structures
of the errors.
Not all elements of Λ, Φ, and Ψ can be estimated. At least one element of each column of Λ must
be fixed to 1, and Ψ is usually a diagonal matrix.
Constraints
16-29
GAUSS User Guide
This command file can be found in the file sqpfact.e in the examples subdirectory:
#include sqpsolvemt.sdf
lmask = { 0 0,
1 0,
0 0,
0 1 };
tmask = { 1 0 0 0,
0 1 0 0,
0 0 1 0,
0 0 0 1 };
struct PV par0;
par0 = pvCreate;
par0 = pvPackm(par0,lambda,"lambda",lmask);
par0 = pvPacks(par0,phi,"phi");
par0 = pvPacksm(par0,psi,"psi",tmask);
16-30
Structures
c0 = sqpSolveMTcontrolCreate;
c0.bounds = ones(pvLength(par0),1).*(-1e250˜1e250);
c0.bounds[lind,1] = zeros(rows(lind),1);
c0.bounds[lind,2] = 10*ones(rows(lind),1);
c0.bounds[tind,1] = .001*ones(rows(tind),1);
c0.bounds[tind,2] = 100*ones(rows(tind),1);
c0.output = 1;
c0.printIters = 1;
c0.trustRadius = 1;
Structures
c0.ineqProc = &ineq;
c0.covType = 1;
struct DS d0;
d0 = dsCreate;
d0.dataMatrix = loadd("maxfact");
lambdahat = pvUnpack(out0.par,"lambda");
phihat = pvUnpack(out0.par,"phi");
psihat = pvUnpack(out0.par,"psi");
print "estimates";
print;
print "lambda" lambdahat;
16-31
GAUSS User Guide
print;
print "phi" phihat;
print;
print "psi" psihat;
struct PV stderr;
stderr = out0.par;
if not scalmiss(out0.moment);
stderr = pvPutParVector(stderr,sqrt(diag(out0.moment)));
lambdase = pvUnpack(stderr,"lambda");
phise = pvUnpack(stderr,"phi");
psise = pvUnpack(stderr,"psi");
print "standard errors";
print;
print "lambda" lambdase;
print;
print "phi" phise;
print;
print "psi" psise;
endif;
output off;
lambda = pvUnpack(par1,"lambda");
phi = pvUnpack(par1,"phi");
psi = pvUnpack(par1,"psi");
sigma = lambda*phi*lambda’ + psi;
logl = -lnpdfmvn(data1.dataMatrix,sigma);
retp(logl);
endp;
16-32
Structures
local lambda,phi,psi,sigma,e;
lambda = pvUnpack(par1,"lambda");
phi = pvUnpack(par1,"phi");
psi = pvUnpack(par1,"psi");
sigma = lambda*phi*lambda’ + psi;
e = eigh(sigma) - .001; /* eigenvalues of sigma */
e = e | eigh(phi) - .001; /* eigenvalues of phi */
retp(e);
endp;
Structures
16-33
Run-Time Library Structures 17
Two structures are used by several GAUSS Run-Time Library functions for handling parameter
vectors and data: the PV parameter structure and the DS data structure.
Structures
RTL
The members of an instance of structure of type PV are all “private,” that is, not accessible directly
to the user. It is designed to handle parameter vectors for threadsafe optimization functions.
Entering and receiving parameter vectors, and accessing properties of this vector, are
accomplished using special functions.
Suppose you are optimizing a function containing a K×L matrix of coefficients. The optimization
function requires a parameter vector but your function uses a K×L matrix. Your needs and the
needs of the optimization function can be both satisfied by an instance of the structure of type PV.
For example:
struct PV p1;
p1 = pvCreate;
17-1
GAUSS User Guide
The pvCreate function initializes p1 to default values. pvPack enters the 4×3 matrix stored
row-wise as a 12×1 parameter vector for the optimization function. The optimization program will
pass the instance of the structure of type PV to your objective function.
By calling pvUnpack your 4×3 coefficient matrix is retrieved from the parameter vector. For
example, in your procedure you have
x = pvUnpack(p1,"coefficients");
and now x is a 4×3 matrix of coefficients for your use in calculating the object function.
Suppose that your objective function has parameters to be estimated in a covariance matrix. The
covariance matrix is a symmetric matrix where only the lower left portion contains unique values
for estimation. To handle this, use pvPacks. For example:
struct PV p1;
p1 = pvCreate;
cov = { 1 .1 .1,
.1 1 .1,
.1 .1 1 };
p1 = pvPacks(p1,cov,"covariance");
Only the lower left portion of cov will be stored in the parameter vector. When the covariance
matrix is unpacked, the parameters in the parameter vector will be entered into both the lower and
upper portions of the matrix.
There may be cases where only a portion of a matrix being used to compute the objective function
are parameters to be estimated. In this case use pvPackm with a “mask” matrix that contains ones
where parameters are to be estimated and zeros otherwise. For example,
17-2
Run-Time Library Structures
struct PV p1;
p1 = pvCreate;
cov = { 1 .5,
.5 1 };
mask = { 0 1,
1 0 };
p1 = pvPacksm(p1,cov,"correlation",mask);
Here only the one element in the lower left of cov is stored in the parameter vector. Suppose the
optimization program sends a trial value for that parameter of, say, .45. When the matrix is
unpacked in your procedure it will contain the fixed values associated with the zeros in the mask
as well as the trial value in that part of the matrix associated with the ones. Thus,
print unpack(p1,"correlation");
1.0000 .4500
.4500 1.0000
A mask may also be used with general matrices to store a portion of a matrix in the parameter
Structures
RTL
vector.
struct PV p1;
p1 = pvCreate;
m = { 0 .5 1,
.5 0 .3 };
mask = { 0 1 1,
1 0 0};
p1 = pvPackm(p1,m,"coefficients",mask);
17-3
GAUSS User Guide
A PV instance can, of course, hold parameters from all these types of matrices: symmetric, masked
symmetric, rectangular, and masked rectangular. For example:
lmask = { 0 0,
1 0,
0 0,
0 1 };
tmask = { 1 0 0 0,
0 1 0 0,
0 0 1 0,
0 0 0 1 };
struct PV par0;
par0 = pvCreate;
par0 = pvPackm(par0,lambda,"lambda",lmask);
par0 = pvPacks(par0,phi,"phi");
par0 = pvPacksm(par0,theta,"theta",tmask);
It isn’t necessary to know where in the parameter vector the parameters are located in order to use
them in your procedure calculating the objective function. Thus:
lambda = pvUnpack(par1,"lambda");
17-4
Run-Time Library Structures
phi = pvUnpack(par1,"phi");
theta = pvUnpack(par1,"theta");
sigma = lambda*phi*lambda’ + theta;
Additional functions are available to retrieve information on the properties of the parameter vector.
pvGetParVector and pvPutParVector get and put parameter vector from and into the PV
instance, pvGetParNames retrieves names for the elements of the parameter vector, pvList
returns the list of matrix names in the PV instance, pvLength the length of the parameter vector.
struct PV p1;
p1 = pvCreate;
cov = { 1 .5,
.5 1 };
mask = { 0 1,
1 0 };
p1 = pvPacksm(p1,cov,"correlation",mask);
print pvGetParVector(p1);
.5000
Structures
RTL
p1 = pvPutParVector(p1,.8);
print pvGetParVector(p1);
.8000
print pvUnpack(p1,"correlation");
1.0000 .8000
.8000 1.0000
print pvGetParNames(p1);
correlation[2,1]
17-5
GAUSS User Guide
print pvLength(p1);
1.0000
Also, pvTest tests an instance to make sure it is properly constructed. pvCreate generates an
initialized instance, and pvGetIndex returns the indices of the parameters of an input matrix in
the parameter vector. This last function is most useful when constructing linear constraint indices
for the optimization programs.
Unpacking matrices using matrix names is slow because it requires a string search through a string
array of names. A set of special packing functions are provided that avoid the search altogether.
These functions use a “table” of indices that you specify to find the matrix in the PV instance. For
example:
struct PV p1;
p1 = pvCreate;
y = rndn(4,1);
x = rndn(4,4);
p1 = pvPacki(p1,y,"Y",1);
p1 = pvPacki(p1,x,"X",2);
print pvUnpack(p1,1);
.3422
.0407
.5611
.0953
print pvUnpack(p1,"Y");
17-6
Run-Time Library Structures
.3422
.0407
.5611
.0953
The call to pvPacki puts an entry in the table associating the matrix in its second argument with
the index 1. As indicated above the matrix can be unpacked either by index or by name.
Unpacking by index, however, is much faster than by name.
Note that the matrix can be unpacked using either the index or the matrix name.
There are index versions of all four of the packing functions, pvPacki, pvPackmi, pvPacksi, and
pvPacksmi.
struct DS d0;
Structures
RTL
d0.dataMatrix M×K matrix, data
d0.dataArray N-dimensional array, data
d0.type scalar
d0.dname string
d0.vnames string array
The definition and use of the elements of d0 are determined by the particular application and are
mostly up to the user. A typical use might use a vector of structures. For example, suppose the
objective function requires a vector of observations on a dependent variable as well as on K
independent variables. Then:
17-7
GAUSS User Guide
struct DS d0;
d0 = dsCreate;
y = rndn(20,1);
x = rndn(20,5);
d0 = reshape(d0,2,1);
d0[1].dataMatrix = y;
d0[2].dataMatrix = X;
The d0 instance would be passed to the optimization program which then passes it to your
procedure computing the objective function. For example:
A particular application may require setting other members of the DS instance for particular
purposes, but in general you may use them for your own purposes. For example, d0.dname could
be set to a GAUSS dataset name from which you read the data in the objective function procedure,
or d0.vnames could be set to the variable names of the columns of the data stored in
d0.dataMatrix, or d0.type could be an indicator variable for the elements of a vector of DS
instances.
The following are complete examples of the use of the PV and DS structures. The first example fits
a set of data to the Micherlitz model. It illustrates packing and unpacking by index.
#include sqpsolvemt.sdf
struct DS Y;
Y = dsCreate;
Y.dataMatrix = 3.183|
17-8
Run-Time Library Structures
3.059|
2.871|
2.622|
2.541|
2.184|
2.110|
2.075|
2.018|
1.903|
1.770|
1.762|
1.550;
struct DS X;
X = dsCreate;
X.dataMatrix = seqa(1,1,13);
struct DS Z;
Z = reshape(Z,2,1);
Z[1] = Y;
Z[2] = X;
Structures
RTL
/* default values */
struct PV par1;
par1 = pvCreate;
17-9
GAUSS User Guide
start = { 2, 4, 2 };
par1 = pvPacki(par1,start,"Parameters",1);
estimates = pvGetParVector(out1.par);
print " parameter estimates ";
print estimates;
print;
print " standard errors ";
print sqrt(diag(out1.moment));
p0 = pvUnpack(par1,1);
w = exp(-p0[3]*Z[2].dataMatrix);
e = z[1].dataMatrix - p0[1] - p0[2] * w;
s2 = moment(e,0) / rows(e);
e1 = - ones(rows(e),1);
e2 = -w;
e3 = p0[2]*Z[2].dataMatrix.*w;
w = (1 - e.*e / s2) / rows(e);
g = e.*e1 + w*(e’e1);
g = g ˜ (e.*e2 + w*(e’e2));
g = g ˜ (e.*e3 + w*(e’e3));
17-10
Run-Time Library Structures
retp(4*g/(rows(e)*s2));
endp;
\#include sqpsolvemt.sdf
lmask = { 0 0,
1 0,
0 0,
0 1 };
Structures
RTL
0.0 0.0 0.0 0.6 };
tmask = { 1 0 0 0,
0 1 0 0,
0 0 1 0,
0 0 0 1 };
struct PV par0;
par0 = pvCreate;
par0 = pvPackm(par0,lambda,"lambda",lmask);
par0 = pvPacks(par0,phi,"phi");
par0 = pvPacksm(par0,theta,"theta",tmask);
17-11
GAUSS User Guide
c0.bounds = ones(pvLength(par0),1).*(-1e250˜1e250);
c0.bounds[lind,1] = zeros(rows(lind),1);
c0.bounds[lind,2] = 10*ones(rows(lind),1);
c0.bounds[tind,1] = .001*ones(rows(tind),1);
c0.bounds[tind,2] = 100*ones(rows(tind),1);
c0.ineqProc = &ineq;
c0.covType = 1;
struct DS d0;
d0 = dsCreate;
d0.dataMatrix = loadd("maxfact");
lambdahat = pvUnpack(out0.par,"lambda");
phihat = pvUnpack(out0.par,"phi");
thetahat = pvUnpack(out0.par,"theta");
print "estimates";
print;
print "lambda" lambdahat;
print;
17-12
Run-Time Library Structures
struct PV stderr;
stderr = out0.par;
if not scalmiss(out0.moment);
stderr =
pvPutParVector(stderr,sqrt(diag(out0.moment)));
lambdase = pvUnpack(stderr,"lambda");
phise = pvUnpack(stderr,"phi");
thetase = pvUnpack(stderr,"theta");
Structures
RTL
proc lpr(struct PV par1, struct DS data1);
local lambda,phi,theta,sigma,logl;
lambda = pvUnpack(par1,"lambda");
phi = pvUnpack(par1,"phi");
theta = pvUnpack(par1,"theta");
sigma = lambda*phi*lambda’ + theta;
logl = -lnpdfmvn(data1.dataMatrix,sigma);
retp(logl);
endp;
17-13
GAUSS User Guide
lambda = pvUnpack(par1,"lambda");
phi = pvUnpack(par1,"phi");
theta = pvUnpack(par1,"theta");
sigma = lambda*phi*lambda’ + theta;
retp(e);
endp;
17-14
Multi-Threaded
Programming in GAUSS 18
The term thread comes from the phrase “thread of execution”—simply, it denotes a section of code
that you want to execute. A single-threaded program has only one thread of execution, i.e., the
program itself. A multi-threaded program is one that can have multiple threads—sections of
code—executing simultaneously. Since these threads are part of the same program, they share the
same workspace, and see and operate on the same symbols. Threads allow you to take full
advantage of the hardware processing resources available on hyper-threaded, multi-core, and
multi-processor systems, executing independent calculations simultaneously, combining and using
the results of their work when done.
18-1
GAUSS User Guide
ThreadJoin Completes the definition of a set of threads, waits until they are done.
ThreadStat n = m’m;
ThreadBegin;
y = x’x;
z = y’y;
ThreadEnd;
Finally, ThreadJoin completes the definition of a set of threads. It waits for the threads in a set to
finish and rejoin the creating (the parent) thread, which can then continue, making use of their
individual calculations:
18-2
Multi-Threaded Programming in GAUSS
ThreadBegin; // Thread 1
y = x’x;
z = y’y;
ThreadEnd;
ThreadBegin; // Thread 2
q = r’r;
r = q’q;
ThreadEnd;
ThreadStat n = m’m; // Thread 3
ThreadStat p = o’o; // Thread 4
ThreadJoin; // waits for Threads 1-4 to finish
b = z + r + n’p; // Using the results
This is really the one and only thing you need to know about threads: threads are separate sections
of the same program, executing simultaneously, operating on the same data. In fact, it’s so
fundamental it’s worth saying again: threads are separate sections of code in a program, running at
the same time, using the same workspace, referencing and operating on the same symbols.
This raises basic issues of workflow and data integrity. How do you manage the creation and
execution of threads, and make use of the work they do? And how do you maintain data integrity?
(You do not want two threads assigning to the same symbol at the same time.)
To handle thread workflow, GAUSS employs a split-and-join approach. At various points in your
program (as many as you like), you define a set of threads that will be created and run as a group.
When created, the threads in the set execute simultaneously, each doing useful work. The parent
thread waits for the created threads to complete, then continues, the results of their work now
available for further use.
Threads
To maintain data integrity, we introduce the writer-must-isolate (informally, the
any-thread-can-read-unless-some-thread-writes) programming rule. That is to say, symbols that
are read from but not assigned to can be referenced by as many threads in a set as you like.
Symbols that are assigned to, however, must be wholly owned by a single thread. No other thread
in the set can reference that symbol. They cannot assign to it, nor can they read from it. They
cannot refer to it at all.
18-3
GAUSS User Guide
Note: the writer-must-isolate rule only applies to the threads within a given set (including any
child thread sets they may create). It does not apply between thread sets that have no chance of
running simultaneously.
For threads defined in the main code, the writer-must-isolate rule applies to the global symbols.
For threads defined in procedures or keywords, it applies to the global symbols, local symbols, and
the procedure/keyword arguments.
One—you can define threads anywhere. You can define them in the main code, you can define
them in proc’s and keyword’s, and yes, you can define them inside other threads.
Two—you can call proc’s and keyword’s from threads. This is what really ties everything
together. You can call a proc from a thread, and that proc can create threads, and any of those
threads can call proc’s, and any of those proc’s can create threads, and ... you get the picture.
q = chol(b);
ThreadBegin;
x = q + m;
ThreadBegin;
y = x’x;
z = q’m;
ThreadEnd;
ThreadBegin;
a = b + x;
c = a + m;
ThreadEnd;
ThreadJoin;
q = m’c;
ThreadEnd;
18-4
Multi-Threaded Programming in GAUSS
ThreadBegin;
ThreadStat r = m’m;
ThreadStat s = m + inv(b);
ThreadJoin;
t = r’s;
ThreadEnd;
ThreadJoin;
x = r+s+q+z-t;
proc bef(x);
local y,t;
ThreadStat y = nof(x);
ThreadStat t = dof(x’x);
ThreadJoin;
t = t+y;
retp(t);
endp;
proc abr(m);
local x,y,z,a,b;
a = m’m;
ThreadStat x = inv(m);
ThreadStat y = bef(m);
ThreadStat z = dne(a);
ThreadJoin;
Threads
b = chut(x,y,z,a);
retp(inv(b));
endp;
s = rndn(500,500);
18-5
GAUSS User Guide
ThreadStat t = abr(s);
ThreadStat q = abr(sˆ2);
ThreadStat r = che(s);
ThreadJoin;
w = del(t,q,r);
print w[1:10,1:10];
This means you can multi-thread anything you want, and call it from anywhere. You can
multi-thread all the proc’s and keyword’s in your libraries, and call them freely anywhere in your
multi-threaded programs.
A few points on coding restrictions. First, you can’t interlace thread definition statements and
regular statements. You can’t do this:
ThreadStat a = b’b;
n = q;
ThreadStat c = d’d;
ThreadJoin;
Or this:
if k == 1;
ThreadStat a = b’b;
elseif k == 2;
ThreadStat a = c’c;
endif;
if j == 1;
ThreadStat d = e’e;
elseif j == 2;
ThreadStat d = f’f;
endif;
ThreadJoin;
18-6
Multi-Threaded Programming in GAUSS
Each set of threads is defined as a group, and always completed by a ThreadJoin, like this:
n = q;
ThreadStat a = b’b;
ThreadStat c = d’d;
ThreadJoin;
And this:
ThreadBegin;
if k == 1;
a = b’b;
elseif k == 2;
a = c’c;
endif;
ThreadEnd;
ThreadBegin;
if j == 1;
d = e’e;
elseif j == 2;
d = f’f;
endif;
ThreadEnd;
ThreadJoin;
Second—as stated above, you can reference read-only symbols in as many threads within a set as
you like, but any symbols that are assigned to must be wholly owned by a single thread. A symbol
that is assigned to by a thread cannot be written or read by any other thread in that set. This is the
Threads
writer-must-isolate rule.
ThreadStat x = y’y;
ThreadStat z = y+y;
18-7
GAUSS User Guide
ThreadStat a = b-y;
ThreadJoin;
ThreadStat x = y’y;
Threadstat z = x’x;
ThreadStat a = b-y;
ThreadJoin;
This is because the threads within a set run simultaneously. Thus, there is no way of knowing
when an assignment to a symbol has taken place, no way of knowing in one thread the “state” of a
symbol in another.
Let’s revisit the nested thread example for a minute and see how the writer-must-isolate rule
applies to it:
18-8
Multi-Threaded Programming in GAUSS
The main point here is that any symbols a thread or its children assign to must be isolated from all
the other threads (and their children) of the same nesting level in that set. On the other hand, the
children of a thread can freely read/write symbols that are read/written by their parent, because
there is no risk of simultaneity; they must only isolate written symbols from their siblings and
siblings’ offspring.
If you break the writer-must-isolate rule, your program (and probably GAUSS) will crash.
Worse, until it crashes, it will be happily producing indeterminate results.
Finally—the ThreadEnd command is what tells a thread to terminate, so you mustn’t write code
that keeps a thread from reaching it. For example, don’t retp from the middle of a thread:
ThreadStat m = imt( 9 );
ThreadBegin;
x = q[1];
if x = 1;
retp(z);
else;
r = z + 2;
endif;
ThreadEnd;
ThreadJoin;
And don’t use goto to jump into or out of the middle of a thread:
Threads
retry:
ThreadBegin;
{ err, x } = fna(q);
if err;
18-9
GAUSS User Guide
goto badidea;
endif;
x = fnb(x);
ThreadEnd;
ThreadStat y = fnb(y);
ThreadJoin;
z = fnc(x,y);
save z;
end;
badidea:
errorlog "Error computing fna(q)";
q = fnd(q);
goto retry;
Basically, don’t do anything that will keep a thread from reaching its ThreadEnd command.
That’s the only way it knows its work is done. end and stop are okay to call, though—they will
bring the program to an end as usual, and terminate all running threads in the process.
(You can use goto and labels to jump around within a thread—that is, within the confines of a
ThreadBegin/ThreadEnd pair.)
18-10
Libraries
Libraries 19
The GAUSS library system allows for the creation and maintenance of modular programs. The
user can create “libraries” of frequently used functions that the GAUSS system will automatically
find and compile whenever they are referenced in a program.
19.1 Autoloader
The autoloader resolves references to procedures, keywords, matrices, and strings that are not
defined in the program from which they are referenced. The autoloader automatically locates and
compiles the files containing the symbol definitions that are not resolved during the compilation of
the main file. The search path used by the autoloader is first the current directory, and then the
paths listed in the src_path configuration variable in the order they appear. src_path can be
defined in the GAUSS configuration file.
19-1
GAUSS User Guide
When the compiler encounters a symbol that has not previously been defined, that is called a
“forward reference”. GAUSS handles forward references in two ways, depending on whether they
are “left-hand side” or “right-hand side” references.
Left-Hand Side
A left-hand side reference is usually a reference to a symbol on the left-hand side of the equal sign
in an expression.
x = 5;
Left-hand side references, since they are assignments, are assumed to be matrices. In the statement
above, x is assumed to be a matrix and the code is compiled accordingly. If, at execution time, the
expression actually returns a string, the assignment is made and the type of the symbol x is forced
to string.
Some commands are implicit left-hand side assignments. There is an implicit left-hand side
reference to x in each statement below:
clear x;
load x;
open x = myfile;
Right-Hand Side
A right-hand side reference is usually a reference to a symbol on the right-hand side of the equal
sign in an expression such as:
19-2
Libraries
Libraries
z = 6;
y = z + dog;
print y;
In the program above, since dog is not previously known to the compiler, the autoloader will
search for it in the active libraries. If it is found, the file containing it will be compiled. If it is not
found in a library, the autoload/autodelete state will determine how it is handled.
If the autoloader is OFF, no forward references are allowed. Every procedure, matrix, and string
referenced by your program must be defined before it is referenced. An external statement can
be used above the first reference to a symbol, but the definition of the symbol must be in the main
file or in one of the files that are #include’d. No global symbols are deleted automatically.
If the autoloader is ON, GAUSS searches for unresolved symbol references during compilation
using a specific search path as outlined below. If the autoloader is OFF, an Undefined symbol
error message will result for right-hand side references to unknown symbols.
When autoload is ON, the autodelete state controls the handling of references to unknown
symbols.
The following search path will be followed to locate any symbols not previously defined:
Autodelete ON
1. user library
2. user-specified libraries.
3. gauss library
4. current directory, then src_path for files with a .g extension.
Forward references are allowed and .g files need not be in a library. If there are symbols that
cannot be found in any of the places listed above, an Undefined symbol error message will be
19-3
GAUSS User Guide
generated and all uninitialized variables and all procedures with global references will be deleted
from the global symbol table. This autodeletion process is transparent to the user, since the
symbols are automatically located by the autoloader the next time the program is run. This process
results in more compile time, which may or may not be significant, depending on the speed of the
computer and the size of the program.
Autodelete OFF
1. user library
2. user-specified libraries.
3. gauss library
All .g files must be listed in a library. Forward references to symbols that are not listed in an
active library are not allowed. For example:
x = rndn(10,10);
y = sym(x); /* Forward reference to sym */
proc sym(x);
retp(x+x’);
endp;
Use an external statement for anything referenced above its definition if autodelete is OFF:
x = rndn(10,10);
y = sym(x);
proc sym(x);
retp(x+x’);
endp;
19-4
Libraries
Libraries
When autodelete is OFF, symbols not found in an active library will not be added to the symbol
table. This prevents the creation of uninitialized procedures in the global symbol table. No
deletion of symbols from the global symbol table will take place.
Libraries
The first place GAUSS looks for a symbol definition is in the “active” libraries. A GAUSS library
is a text file that serves as a dictionary to the source files that contain the symbol definitions. When
a library is active, GAUSS will look in it whenever it is looking for a symbol it is trying to resolve.
The library statement is used to make a library active. Library files should be located in the
subdirectory listed in the lib_path configuration variable. Library files have an .lcg extension.
Suppose you have several procedures that are all related and you want them all defined in the same
file. You can create such a file, and, with the help of a library, the autoloader will be able to find
the procedures defined in that file whenever they are called.
First, create the file that is to contain your desired procedure definitions. By convention, this file is
usually named with a .src extension, but you may use any name and any file extension. In this
file, put all the definitions of related procedures you wish to use. Here is an example of such a file.
It is called norm.src:
/*
** norm.src
**
** This is a file containing the definitions of three
** procedures which return the norm of a matrix x.
** The three norms calculated are the 1-norm, the
** inf-norm and the E-norm.
*/
proc onenorm(x);
retp(maxc(sumc(abs(x))));
endp;
proc infnorm(x);
19-5
GAUSS User Guide
retp(maxc(sumc(abs(x’))));
endp;
proc Enorm(x);
retp(sumc(sumc(x.*x)));
endp;
Next, create a library file that contains the name of the file you want access to, and the list of
symbols defined in it. This can be done with the lib command. (For details, see lib in the
GAUSS L R.)
A library file entry has a filename that is flush left. The drive and path can be included to speed up
the autoloader. Indented below the filename are the symbols included in the file. There can be
multiple symbols listed on a line, with spaces between. The symbol type follows the symbol
name, with a colon delimiting it from the symbol name. The valid symbol types are:
keyword keyword.
proc procedure.
string string.
struct structure.
If the symbol type is missing, the colon must not be present and the symbol type is assumed to be
proc. Both library files below are valid:
Example 1
19-6
Libraries
Libraries
/*
** math
**
** This library lists files and procedures for mathematical routines.
*/
norm.src
onenorm:proc infnorm:proc Enorm:proc
complex.src
cmmult:proc cmdiv:proc cmadd:proc cmsoln:proc
poly.src
polychar:proc polyroot:proc polymult:proc
Example 2
/*
** math
**
** This library lists files and procedures for mathematical routines.
*/
c:\gauss\src\norm.src
onenorm : proc
infnorm : proc
Enorm : proc
c:\gauss\src\complex.src
cmmult : proc
cmdiv : proc
cmadd : proc
cmsoln : proc
c:\gauss\src\fcomp.src
feq : proc
fne : proc
flt : proc
fgt : proc
fle : proc
19-7
GAUSS User Guide
fge : proc
c:\gauss\src\fcomp.dec
_fcmptol : matrix
Once the autoloader finds, via the library, the file containing your procedure definition, everything
in that file will be compiled. For this reason, you should combine related procedures in the same
file in order to minimize the compiling of procedures not needed by your program. In other words,
you should not combine unrelated functions in one .src file because if one function in a .src file
is needed, the whole file will be compiled.
user Library
This is a library for user-created procedures. If the autoloader is ON, the user library is the first
place GAUSS looks when trying to resolve symbol references.
You can update the user library with the lib command as follows:
No user library is shipped with GAUSS. It will be created the first time you use the lib
command to update it.
For details on the parameters available with the lib command, see the GAUSS L
R.
.g Files
If autoload and autodelete are ON and a symbol is not found in a library, the autoloader will
assume it is a procedure and look for a file that has the same name as the symbol and a .g
extension. For example, if you have defined a procedure called square, you could put the
definition in a file called square.g in one of the subdirectories listed in your src_path. If
autodelete is OFF, the .g file must be listed in an active library; for example, in the user library.
19-8
Libraries
Libraries
19.2 Global Declaration Files
If your application makes use of several global variables, create a file containing declare
statements. Use files with the extension .dec to assign default values to global matrices and
strings with declare statements and to declare global N-dimensional arrays, sparse matrices,
and structures, which will be initialized as follows:
In order to declare structures in a .dec file, you must #include the file(s) containing the
definitions of the types of structures that you wish to declare at the top of your .dec file. For
example, if you have the following structure type definition in a file called mystruct.sdf:
struct mystruct {
matrix m;
array a;
scalar scal;
string array sa;
};
You could declare an instance of that structure type, called ms, in a .dec file as follows:
#include mystruct.sdf
See declare in the C R, Chapter 33, for more information.
A file with a .ext extension containing the same symbols in external statements can also be
created and #include’d at the top of any file that references these global variables. An
19-9
GAUSS User Guide
appropriate library file should contain the name of the .dec files and the names of the globals they
declare. This allows you to reference global variables across source files in an application.
Here is an example that illustrates the way in which .dec, .ext, .lcg and .src files work
together. Always begin the names of global matrices or strings with ‘_’ to distinguish them from
procedures.
.src File:
/*
** fcomp.src
**
** These functions use _fcmptol to fuzz the comparison operations
** to allow for roundoff error.
**
** The statement: y = feq(a,b);
**
** is equivalent to: y = a eq b;
**
** Returns a scalar result, 1 (true) or 0 (false)
**
** y = feq(a,b);
** y = fne(a,b);
*/
#include fcomp.ext
proc feq(a,b);
retp(abs(a-b) <= _fcmptol);
endp;
proc fne(a,b);
retp(abs(a-b) > _fcmptol);
endp;
.dec File:
19-10
Libraries
Libraries
/*
** fcomp.dec - global declaration file for fuzzy comparisons.
*/
.ext File:
/*
** fcomp.ext - external declaration file for fuzzy comparisons.
*/
.lcg File:
/*
** fcomp.lcg - fuzzy compare library
*/
fcomp.dec
_fcmptol:matrix
fcomp.src
feq:proc
fne:proc
With the exception of the library (.lcg) files, these files must be located along your src_path.
The library files must be on your lib_path. With these files in place, the autoloader will be able
to find everything needed to run the following programs:
library fcomp;
x = rndn(3,3);
xi = inv(x);
19-11
GAUSS User Guide
xix = xi*x;
if feq(xix,eye(3));
print "Inverse within tolerance.";
else;
print "Inverse not within tolerance.";
endif;
If the default tolerance of 1e-14 is too tight, the tolerance can be relaxed:
library fcomp;
x = rndn(3,3);
xi = inv(x);
xix = xi*x;
_fcmptol = 1e-12; /* reset tolerance */
if feq(xix,eye(3));
print "Inverse within tolerance.";
else;
print "Inverse not within tolerance.";
endif;
19.3 Troubleshooting
Below is a partial list of errors you may encounter in using the library system, followed by the
most probable cause.
19-12
Libraries
Libraries
The autoloader cannot find a file called prt.dec. Check for this file. It should exist
somewhere along your src_path, if you have it listed in prt.lcg.
Undefined symbols:
PRTVEC /gauss/src/tstprt.g(2)
The symbol prtvec could not be found. Check if the file containing prtvec is in the
src_path. You may have not activated the library that contains your symbol definition.
Do so in a library statement.
Below is some advice you are encouraged to follow when constructing your own library system:
• Whenever possible, declare variables in a file that contains only declare statements. When
your program is run again without clearing the workspace, the file containing the variable
declarations will not be compiled and declare warnings will be prevented.
19-13
GAUSS User Guide
• To avoid meaningless redefinition errors and declare warnings, never declare a global
more than once in any one file. Redefinition error messages and declare warnings are
meant to help you prevent name conflicts, and will be useless to you if your code generates
them normally.
By following these guidelines, any declare warnings and redefinition errors you get will be
meaningful. By knowing that such warnings and errors are significant, you will be able to debug
your programs more efficiently.
19-14
Compiler
Compiler 20
GAUSS allows you to compile your large, frequently used programs to a file that can be run over
and over with no compile time. The compiled image is usually smaller than the uncompiled
source. GAUSS is not a native code compiler; rather, it compiles to a form of pseudocode. The file
will have a .gcg extension.
The compile command will compile an entire program to a compiled file. An attempt to edit a
compiled file will cause the source code to be loaded into the editor if it is available to the system.
The run command assumes a compiled file if no extension is given, and that a file with a .gcg
extension is in the src_path. A saveall command is available to save the current contents of
memory in a compiled file for instant recall later. The use command will instantly load a
compiled program or set of procedures at the beginning of an ASCII program before compiling the
rest of the ASCII program file.
Since the compiled files are encoded binary files, the compiler is useful for developers who do not
want to distribute their source code.
20-1
GAUSS User Guide
Source code program files that can be run with the run command can be compiled to .gcg files
with the compile command:
compile qxy.e;
All procedures, global matrices, arrays, strings and string arrays, and the main program segment
will be saved in the compiled file. The compiled file can be run later using the run command. Any
libraries used in the program must be present and active during the compile, but not when the
program is run. If the program uses the dlibrary command, the .dll files must be present when
the program is run and the dlibrary path must be set to the correct subdirectory. This will be
handled automatically in your configuration file. If the program is run on a different computer than
it was compiled on, the .dll files must be present in the correct location. sysstate (case 24) can
be used to set the dlibrary path at run-time.
The simplest way to create a compiled file containing a set of frequently used procedures is to use
saveall and an external statement:
library pgraph;
external proc xy,logx,logy,loglog,hist;
saveall pgraph;
Just list the procedures you will be using in an external statement and follow it with a saveall
statement. It is not necessary to list procedures that you do not explicitly call, but are called from
20-2
Compiler
another procedure, because the autoloader will automatically find them before the saveall
command is executed. Nor is it necessary to list every procedure you will be calling, unless the
source will not be available when the compiled file is use’d.
Remember, the list of active libraries is NOT saved in the compiled file, so you may still need a
Compiler
library statement in a program that is use’ing a compiled file.
20.3 Debugging
If you are using compiled code in a development situation in which debugging is important,
compile the file with line number records. After the development is over, you can recompile
without line number records if the maximum possible execution speed is important. If you want to
guarantee that all procedures contain line number records, put a new statement at the top of your
program and turn line number tracking on.
20-3
File I/O 21
File I/O
The following is a partial list of the I/O commands in the GAUSS programming language:
21-1
GAUSS User Guide
21-2
File I/O
File I/O
GAUSS has facilities for reading and writing ASCII files. Since most software can also read and
write ASCII files, this provides one method of sharing data between GAUSS and many other
kinds of programs.
Reading
Files containing numeric data that are delimited with spaces or commas and are small enough to fit
into a single matrix or string can be read with load. Larger ASCII data files can be converted to
GAUSS data sets with the ATOG utility program (see ATOG, Chapter 28). ATOG can convert
packed ASCII files as well as delimited files.
For small delimited data files, the load statement can be used to load the data directly into a
GAUSS matrix. The resulting GAUSS matrix must be no larger than the limit for a single matrix.
For example,
will load the data in the file dat1.asc into an N×1 matrix x. This method is preferred because
rows(x) can be used to determine how many elements were actually loaded, and the matrix can
be reshape’d to the desired form:
21-3
GAUSS User Guide
This will load the data into a 100×5 matrix. If there are more or fewer than 500 numbers in the
data set, the matrix will automatically be reshaped to 100×5.
Writing
To write data to an ASCII file the print or printfm command is used to print to the auxiliary
output. The resulting files are standard ASCII files and can be edited with GAUSS’s editor or
another text editor.
The output and outwidth commands are used to control the auxiliary output. The print or
printfm command is used to control what is sent to the output file.
The window can be turned on and off using screen. When printing a large amount of data to the
auxiliary output, the window can be turned off using the command
screen off;
This will make the process much faster, especially if the auxiliary output is a disk file.
It is easy to forget to turn the window on again. Use the end statement to terminate your
programs; end will automatically perform screen on and output off.
The following commands can be used to control printing to the auxiliary output:
21-4
File I/O
File I/O
This example illustrates printing a matrix to a file:
The numbers in the matrix x will be printed with a field width of 8 spaces per number, and with 2
places beyond the decimal point. The resulting file will be an ASCII data file. It will have 132
column lines maximum.
A more extended example follows. This program will write the contents of the GAUSS file
mydata.dat into an ASCII file called mydata.asc. If there is an existing file by the name of
mydata.asc, it will be overwritten:
21-5
GAUSS User Guide
endo;
fp = close(fp);
end;
The output ... reset command will create an auxiliary output file called mydata.asc to
receive the output. The window is turned off to speed up the process. The GAUSS data file
mydata.dat is opened for reading and 200 rows are read per iteration until the end of the file is
reached. The data read are printed to the auxiliary output mydata.asc only, because the window
is off.
getf will read a file and return it in a string variable. Any kind of file can be read in this way as
long as it will fit into a single string variable.
To read files sequentially, use fopen to open the file and use fgets, fputs, and associated
functions to read and write the file. The current position in a file can be determined with ftell.
The following example uses these functions to copy an ASCII text file:
do until eof(fin);
21-6
File I/O
call close(fin);
call close(fout);
retp(0);
File I/O
endp;
GAUSS data sets are the preferred method of storing data contained in a single matrix for use
within GAUSS. Use of these data sets allows extremely fast reading and writing of data. Many
library functions are designed to read data from these data sets.
If you want to store multiple variables of various types in a single file, see GAUSS D A,
Section 21.3.
21.2.1 Layout
GAUSS data sets are arranged as matrices; that is, they are organized in terms of rows and
columns. The columns in a data file are assigned names, and these names are stored in the header,
or, in the case of the v89 format, in a separate header file.
The limit on the number of rows in a GAUSS data set is determined by disk size. The limit on the
number of columns is limited by RAM. Data can be stored in 2, 4, or 8 bytes per number, rather
than just 8 bytes as in the case of GAUSS matrix files.
21-7
GAUSS User Guide
Data sets can be created with the create or datacreate command. The names of the columns,
the type of data, etc., can be specified. (For details, see create in the GAUSS L
R.)
Data sets, unlike matrices, cannot change from real to complex, or vice-versa. Data sets are always
stored a row at a time. The rows of a complex data set, then, have the real and imaginary parts
interleaved, element by element. For this reason, you cannot write rows from a complex matrix to
a real data set—there is no way to interleave the data without rewriting the entire data set. If you
must, explicitly convert the rows of data first, using the real and imag functions (see the GAUSS
L R), and then write them to the data set. Rows from a real matrix CAN be
written to a complex data set; GAUSS simply supplies 0’s for the imaginary part.
To create a complex data set, include the complex flag in your create command.
The basic functions in GAUSS for reading data files are open and readr:
open f1 = dat1;
x = readr(f1,100);
The call to readr in this example will read in 100 rows from dat1.dat. The data will be assigned
to a matrix x.
loadd and saved can be used for loading and saving small data sets.
21-8
File I/O
The following example illustrates the creation of a GAUSS data file by merging (horizontally
concatenating) two existing data sets:
file1 = "dat1";
file2 = "dat2";
outfile = "daty";
open fin1 = ˆfile1 for read;
open fin2 = ˆfile2 for read;
varnames = getname(file1)|getname(file2);
otyp = maxc(typef(fin1)|typef(fin2));
create fout = ˆoutfile with ˆvarnames,0,otyp;
File I/O
nr = 400;
do until eof(fin1) or eof(fin2);
y1 = readr(fin1,nr);
y2 = readr(fin2,nr);
r = maxc(rows(y1)|rows(y2));
y = y1[1:r,.] ˜ y2[1:r,.];
call writer(fout,y);
endo;
closeall fin1,fin2,fout;
In this example, data sets dat1.dat and dat2.dat are opened for reading. The variable names
from each data set are read using getname, and combined in a single vector called varnames. A
variable called otyp is created, which will be equal to the larger of the two data types of the input
files. This will insure that the output is not rounded to less precision than the input files. A new
data set daty.dat is created using the create ... with ... command. Then, on every
iteration of the loop, 400 rows are read in from each of the two input data sets, horizontally
concatenated, and written out to daty.dat. When the end of one of the input files is reached,
reading and writing will stop. The closeall command is used to close all files.
Although GAUSS itself does not distinguish between numeric and character columns in a matrix
or data set, some of the GAUSS Application programs do. When creating a data set, it is important
to indicate the type of data in the various columns. The following discusses two ways of doing this.
21-9
GAUSS User Guide
The v89 data set format distinguished between character and numeric data in data sets by the case
of the variable names associated with the columns. The v96 data set format, however, stores this
type information separately, resulting in a much cleaner and more robust method of tracking
variable types, and greater freedom in the naming of data set variables.
When you create a data set, you can supply a vector indicating the type of data in each column of
the data set. For example:
data = { M 32 21500,
F 27 36000,
F 28 19500,
M 25 32000 };
vnames = { "Sex" "Age" "Pay" };
vtypes = { 0 1 1 };
create f = mydata with ˆvnames, 3, 8, vtypes;
call writer(f,data);
f = close(f);
The call to getnamef in this example returns a string array rather than a character vector, so you
can print it without the ‘$’ prefix.
21-10
File I/O
However, this is now obsolete; use vartypef and v96 data sets to be compatible with future
versions.
File I/O
The GAUSS Data Archive (GDA) is extremely powerful and flexible, giving you much greater
control over how you store your data. There is no limitation on the number of variables that can be
stored in a GDA, and the only size limitation is the amount of available disk space. Moreover,
GDA’s are designed to hold whatever type of data you want to store in them. You may write
matrices, arrays, strings, string arrays, sparse matrices, and structures to a GDA, and the GDA will
keep track of the type, size and location of each of the variables contained in it. Since GAUSS
now supports reading and writing to GDA’s that were created on other platforms, GDA’s provide a
simple solution to the problem of sharing data across platforms.
To create a GAUSS Data Archive, call gdaCreate, which creates a GDA containing only header
information. It is recommended that file names passed into gdaCreate have a .gda extension;
however, gdaCreate will not force an extension.
To write variables to the GDA, you must call gdaWrite. A single call to gdaWrite writes only
one variable to the GDA. Writing multiple variables requires multiple calls to gdaWrite.
ret = gdaCreate("myfile.gda",1);
21-11
GAUSS User Guide
ret = gdaWrite("myfile.gda",rndn(100,50),"x1");
ret = gdaWrite("myfile.gda","This is a string","str1");
ret = gdaWrite("myfile.gda",394,"x2");
1 x1 matrix 100 × 50
2 str1 string 16 chars
3 x2 matrix 1×1
The following table details the commands that you may use to read various types of variables from
a GAUSS Data Archive:
gdaRead, gdaReadSparse, and gdaReadStruct take a variable name and return the variable
data. gdaReadByIndex returns the variable data for a specified variable index.
For example, to get the variable x1 out of myfile.gda, you could call:
y = gdaRead("myfile.gda","x1");
or
y = gdaReadByIndex("myfile.gda",1);
21-12
File I/O
If you want to read only a part of a matrix, array, string, or string array from a GDA, call
gdaReadSome. Sparse matrices and structures may not be read in parts.
To overwrite an entire variable in a GDA, you may call gdaUpdate or gdaUpdateAndPack. If the
new variable is not the same size as the variable that it is replacing, gdaUpdate will leave empty
bytes in the file, while gdaUpdateAndPack will pack the file (from the location of the variable
that is being replaced to the end of the file) to remove those empty bytes.
File I/O
gdaUpdate is usually faster, since it does not move data in the file unnecessarily. However, calling
gdaUpdate several times for one file may result in a file with a large number of empty bytes.
On the other hand, gdaUpdateAndPack uses disk space efficiently, but it may be slow for large
files (especially if the variable to be updated is one of the first variables in the file).
If speed and disk space are both concerns and you are going to update several variables, it will be
most efficient to use gdaUpdate to update the variables and then call gdaPack once at the end to
pack the file.
ret = gdaUpdate("myfile.gda",rndn(1000,100),"x1");
ret = gdaUpdateAndPack("myfile.gda",rndn(1000,100),"x1");
21-13
GAUSS User Guide
The save command takes a matrix in memory, adds a header that contains information on the
number of rows and columns in the matrix, and stores it on disk. Numbers are stored in double
precision just as they are in matrices in memory. These files have the extension .fmt.
Matrix files can be no larger than a single matrix. No variable names are associated with matrix
files.
GAUSS matrix files can be load’ed into memory using the load or loadm command or they can
be opened with the open command and read with the readr command. With the readr
command, a subset of the rows can be read. With the load command, the entire matrix is load’ed.
GAUSS matrix files can be open’ed for read, but not for append, or for update.
If a matrix file has been opened and assigned a file handle, rowsf and colsf can be used to
determine how many rows and columns it has without actually reading it into memory. seekr and
readr can be used to jump to particular rows and to read them into memory. This is useful when
only a subset of rows is needed at any time. This procedure will save memory and be much faster
than load’ing the entire matrix into memory.
21-14
File I/O
File I/O
.dht
Extended Data Set v89 .dat, Obsolete, use v96.
.dht
Data Set v92 .dat Obsolete, use v96.
Universal Data Set v96 .dat Supported for read/write.
Matrix files are binary files, and cannot be read with a text editor. They are created with save.
Matrix files with up to 8190 elements have a .fmt extension and a 16-byte header formatted as
follows:
Offset Description
The body of the file starts at offset 16 and consists of IEEE format double precision floating point
numbers or character elements of up to 8 characters. Character elements take up 8 bytes and are
padded on the right with zeros. The size of the body of the file is 8*rows*cols rounded up to the
next 16-byte paragraph boundary. Numbers are stored row by row. A 2×3 real matrix will be
21-15
GAUSS User Guide
stored on disk in the following way, from the lowest addressed element to the highest addressed
element:
For complex matrices, the size of the body of the file is 16*rows*cols. The entire real part of the
matrix is stored first, then the entire imaginary part. A 2×3 complex matrix will be stored on disk
in the following way, from the lowest addressed element to the highest addressed element:
Matrices with more than 8190 elements are saved in an extended format. These files have a
16-byte header formatted as follows:
Offset Description
The size of the body of an extended matrix file is 8*rows*cols (not rounded up to a paragraph
boundary). Aside from this, the body is the same as the small matrix v89 file.
String files are created with save. String files with up to 65519 characters have a 16-byte header
formatted as follows:
21-16
File I/O
Offset Description
The body of the file starts at offset 16. It consists of the string terminated with a null byte. The size
of the file is the 16-byte header plus the length of the string and null byte rounded up to the next
16-byte paragraph boundary.
File I/O
21.5.4 Extended String v89 (Obsolete)
Strings with more than 65519 characters are saved in an extended format. These files have a
16-byte header formatted as follows:
Offset Description
The body of the file starts at offset 16. It consists of the string terminated with a null byte. The size
of the file is the 16-byte header plus the length of the string and null byte rounded up to the next
8-byte boundary.
All data sets are created with create. v89 data sets consist of two files; one .dht contains the
header information; the second (.dat) contains the binary data. The data will be one of three types:
21-17
GAUSS User Guide
The .dht file is used in conjunction with the .dat file as a descriptor file and as a place to store
names for the columns in the .dat file. Data sets with up to 8175 columns have a .dht file
formatted as follows:
Offset Description
Column names begin at offset 128 and are stored 8 bytes each in ASCII format. Names with less
than 8 characters are padded on the right with bytes of 0.
The number of rows in the .dat file is calculated in GAUSS using the file size, columns, and data
type. This means that users can modify the .dat file by adding or deleting rows with other
software without updating the header information.
Names for the columns should be lowercase for character data, to be able to distinguish them from
numeric data with vartype.
GAUSS currently examines only the 4’s bit of the control flags. This bit is set to 0 for real data
sets, 1 for complex data sets. All other bits are 0.
Data sets are always stored a row at a time. A real data set with 2 rows and 3 columns will be
stored on disk in the following way, from the lowest addressed element to the highest addressed
21-18
File I/O
element:
The rows of a complex data set are stored with the real and imaginary parts interleaved, element
by element. A 2×3 complex data set, then, will be stored on disk in the following way, from the
lowest addressed element to the highest addressed element:
File I/O
[1, 1]r [1, 1]i [1, 2]r [1, 2]i [1, 3]r [1, 3]i
[2, 1]r [2, 1]i [2, 2]r [2, 2]i [2, 3]r [2, 3]i
Data sets with more than 8175 columns are saved in an extended format that cannot be read by the
16-bit version. These files have a .dht descriptor file formatted as follows:
Offset Description
Aside from the differences in the descriptor file and the number of columns allowed in the data
file, extended data sets conform to the v89 data set description specified above.
21-19
GAUSS User Guide
Offset Description
0-3 always 0
4-7 always 0xEECDCDCD
8-11 reserved
12-15 reserved
16-19 reserved
20-23 0 - real matrix, 1 - complex matrix
24-27 number of dimensions
0 - scalar
1 - row vector
2 - column vector, matrix
28-31 header size, 128 + number of dimensions * 4, padded to 8-byte boundary
32-127 reserved
If the data is a scalar, the data will directly follow the header.
If the data is a row vector, an unsigned integer equaling the number of columns in the vector will
precede the data, along with 4 padding bytes.
If the data is a column vector or a matrix, there will be two unsigned integers preceding the data.
The first will represent the number of rows in the matrix and the second will represent the number
of columns.
The data area always begins on an even 8-byte boundary. Numbers are stored in double precision
(8 bytes per element, 16 if complex). For complex matrices, all of the real parts are stored first,
followed by all the imaginary parts.
Offset Description
0-3 always 0
4-7 always 0xEECFCFCF
21-20
File I/O
Offset Description
8-11 reserved
12-15 reserved
16-19 reserved
20-23 size of string in units of 8 bytes
24-27 length of string plus null terminator in bytes
28-127 reserved
The size of the data area is always divisible by 8, and is padded with nulls if the length of the string
is not evenly divisible by 8. If the length of the string is evenly divisible by 8, the data area will be
the length of the string plus 8. The data area follows immediately after the 128-byte header.
File I/O
21.5.9 Data Set v92 (Obsolete)
Offset Description
0-3 always 0
4-7 always 0xEECACACA
8-11 reserved
12-15 reserved
16-19 reserved
20-23 rows in data set
24-27 columns in data set
28-31 0 - real data set, 1 - complex data set
32-35 type of data in data set, 2, 4, or 8
36-39 header size in bytes is 128 + columns * 9
40-127 reserved
The variable names begin at offset 128 and are stored 8 bytes each in ASCII format. Each name
corresponds to one column of data. Names less than 8 characters are padded on the right with
bytes of zero.
The variable type flags immediately follow the variable names. They are 1-byte binary integers,
one per column, padded to an even 8-byte boundary. A 1 indicates a numeric variable and a 0
indicates a character variable.
21-21
GAUSS User Guide
The contents of the data set follow the header and start on an 8-byte boundary. Data is either 2-byte
signed integer, 4-byte single precision floating point or 8-byte double precision floating point.
Offset Description
If the data is a scalar, the data will directly follow the header.
If the data is a row vector, an unsigned integer equaling the number of columns in the vector will
21-22
File I/O
If the data is a column vector or a matrix, there will be two unsigned integers preceding the data.
The first will represent the number of rows in the matrix and the second will represent the number
of columns.
The data area always begins on an even 8-byte boundary. Numbers are stored in double precision
(8 bytes per element, 16 if complex). For complex matrices, all of the real parts are stored first,
followed by all the imaginary parts.
File I/O
21.5.11 Data Set v96
Offset Description
21-23
GAUSS User Guide
Offset Description
68-71 always 0
72-75 header size, 128 + columns * 33, padded to 8-byte boundary
76-79 reserved
80-83 rows in data set
84-87 columns in data set
88-127 reserved
The variable names begin at offset 128 and are stored 32 bytes each in ASCII format. Each name
corresponds to one column of data. Names less than 32 characters are padded on the right with
bytes of zero.
The variable type flags immediately follow the variable names. They are 1-byte binary integers,
one per column, padded to an even 8-byte boundary. A 1 indicates a numeric variable and a 0
indicates a character variable.
Contents of the data set follow the header and start on an 8-byte boundary. Data is either 2-byte
signed integer, 4-byte single precision floating point or 8-byte double precision floating point.
A GAUSS Data Archive consists of a header, followed by the variable data and, finally, an array of
variable descriptors containing information about each variable.
Header
21-24
File I/O
File I/O
44-55 32-bit unsigned integers reserved
56-63 64-bit unsigned integer number of variables
64-67 32-bit unsigned integer header size, 128
68-95 32-bit unsigned integers reserved
96-103 64-bit unsigned integer offset of variable descriptor table from end of
header
104-127 64-bit unsigned integers reserved
Variable Data
After the header comes the variable data. Matrices are laid out in row-major order, and strings are
written with a null-terminating byte.
For string arrays, an array of rows×columns struct satable’s is written out first, followed by the
string array data in row-major order with each element null terminated. A struct satable consists
of two members:
off size_t offset of element data from beginning of string array data
len size_t length of element data, including null-terminating byte
Arrays are written with the orders (sizes) of each dimension followed by the array data. For
example, the following 2×3×4 array:
21-25
GAUSS User Guide
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16
17 18 19 20
21 22 23 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Variable Structures
The variable data is followed by an array of variable descriptors. For each variable in the GDA,
there is a corresponding variable descriptor in this array. A variable descriptor is laid out as
follows:
21-26
File I/O
20-23 32-bit unsigned integer huge flag, indicates whether the variable is larger than
INT MAX
24-31 64-bit unsigned integer rows for matrices and string arrays
32-39 64-bit unsigned integer columns for matrices and string arrays, length for strings,
including null-terminating byte
40-47 64-bit unsigned integer index of the variable in the GDA
48-55 64-bit unsigned integer offset of variable data from end of header
56-63 64-bit unsigned integer length of variable data in bytes
64-143 string name of variable, null-terminated
File I/O
The variable type (bytes 0-3) may be any of the following:
20 array
30 matrix
40 string
50 string array
The size of pointer element (bytes 16-19) is the size of a pointer on the machine on which the
variable was written to the GDA. It will be set to 4 on 32-bit machines and 8 on 64-bit machines.
This element is used only for string array variables. If a GDA containing string arrays is created
on a 32-bit machine and then read on a 64-bit machine, or vice versa, then the size of pointer
element indicates how the members of the struct satable’s must be converted in order to be read
on the current machine.
The huge flag (bytes 20-23) is set to 1 if the variable size is greater than INT MAX, which is
defined as 2147483647. A variable for which the huge flag is set to 1 may not be read into GAUSS
on a 32-bit machine.
The variable index element (bytes 40-47) contains the index of the variable in the GDA. Although
the variable data is not necessarily ordered by index (see gdaUpdate), the variable descriptors are.
Therefore, the indices are always in ascending order.
21-27
Foreign Language Interface 22
FLI
The Foreign Language Interface (FLI) allows users to create functions written in C, FORTRAN, or
other languages, and call them from a GAUSS program. The functions are placed in dynamic
libraries (DLLs, also known as shared libraries or shared objects) and linked in at run-time as
needed. The FLI functions are:
GAUSS recognizes a default dynamic library directory, a directory where it will look for your
dynamic-link libraries when you call dlibrary. You can specify the default directory in
gauss.cfg by setting dlib_path. As it is shipped, gauss.cfg specifies $(GAUSSDIR)/dlib as
the default directory.
22-1
GAUSS User Guide
3. Return an integer.
In C syntax, then, your functions would take one of the following forms:
1. int func(void);
Functions can be written to take a list of up to 100 arguments, or a vector (in C terms, a
1-dimensional array) of up to 1000 arguments. This does not affect how the function is called from
GAUSS; the dllcall statement will always appear to pass the arguments in a list. That is, the
dllcall statement will always look as follows:
dllcall func(a,b,c,d[[,e...]]);
For details on calling your function, passing arguments to it, getting data back, and what the return
value means, see dllcall in the GAUSS L R.
22-2
Foreign Language Interface
The following describes how to build a dynamic library called hyp.dll (on Windows) or
libhyp.so (on UNIX/Linux) from the source file hyp.c.
As mentioned in the previous section, your FLI functions may take only pointers to doubles as
arguments. Therefore, you should define your FLI functions to be merely wrapper functions that
cast their arguments as necessary and then call the functions that actually do the work. This is
demonstrated in the source file hyp.c:
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
FLI
/* this does all the work, not exported */
static int hypo(double *x, double *y, double *h, int r, int c)
{
double *wx;
double *wy;
double *dp;
double *sp1;
double *sp2;
int i, elems;
elems = r*c;
22-3
GAUSS User Guide
free(wx);
return 30; /* out of memory */
}
dp = wx;
sp1 = x;
dp = wy;
sp2 = y;
dp = h;
sp1 = wx;
sp2 = wy;
22-4
Foreign Language Interface
return 0;
}
FLI
The following Makefiles contain the compile and link commands you would use to build the
dynamic library on various platforms. For explanations of the various flags used, see the
documentation for your compiler and linker.
Windows
hyp.dll: hyp.obj
link /dll /out:hyp.dll hyp.obj
hyp.obj: hyp.c
cl -c -MD -GX hyp.c
Solaris
22-5
GAUSS User Guide
libhyp.so: hyp.c
$(CC) -G $(CCOPTIONS) -o $@ hyp.c -lm
Linux
libhyp.so: hyp.cpp
$(CC) $(CCOPTIONS) -o $@ hyp.c
For details on linking your dynamic library, see dlibrary in the GAUSS L R.
22-6
Data Transformations 23
GAUSS allows expressions that directly reference variables (columns) of a data set. This is done
within the context of a data loop:
Data Loop
csed = ln(sqrt(csed));
select csed > 0.35 and married $=\,= "y";
make chfac = hcfac + wcfac;
keep csed chfac stid recsum voom;
endata;
GAUSS translates the data loop into a procedure that performs the required operations, and then
calls the procedure automatically at the location (in your program) of the data loop. It does this by
translating your main program file into a temporary file and then executing the temporary file.
A data loop may be placed only in the main program file. Data loops in files that are #include’d
or autoloaded are not recognized.
23-1
GAUSS User Guide
A data loop begins with a dataloop statement and ends with an endata statement. Inside a data
loop, the following statements are supported:
In any expression inside a data loop, all text symbols not immediately followed by a left
parenthesis ‘(’ are assumed to be data set variable (column) names. Text symbols followed by a
left parenthesis are assumed to be procedure names. Any symbol listed in an extern statement is
assumed to be a matrix or string already in memory.
23-2
Data Transformations
All program statements in the main file and not inside a data loop are passed through to the
temporary file without modification. Program statements within a data loop that are preceded by a
‘#’ are passed through to the temporary file without modification. The user familiar with the code
generated in the temporary file can use this to do out-of-the-ordinary operations inside the data
loop.
The translator that processes data loops can be turned on and off. When the translator is on, there
are three distinct phases in running a program:
Data Loop
In the translation phase, the main program file is translated into a temporary file. Each data loop is
translated into a procedure and a call to this procedure is placed in the temporary file at the same
location as the original data loop. The data loop itself is commented out in the temporary file. All
the data loop procedures are placed at the end of the temporary file.
Depending upon the status of line number tracking, error messages encountered in this phase will
be printed with the file name and line numbers corresponding to the main file.
In the compilation phase, the temporary file is compiled. Depending upon the status of line
number tracking, error messages encountered in this phase will be printed with the file name and
23-3
GAUSS User Guide
line numbers corresponding to both the main file and the temporary file.
In the execution phase, the compiled program is executed. Depending on the status of line number
tracking, error messages will include line number references from both the main file and the
temporary file.
The following local variables are created by the translator and used in the produced code:
These variables are reserved, and should not be used within a dataloop... endata section.
23-4
The GAUSS Profiler 24
GAUSS now includes a profiler, which enables you to determine exactly how much time your
programs are spending on each line and in each called procedure, thereby providing you with the
information you need to increase the efficiency of your programs. The GAUSS Profiler and
tcollect are both run from a command prompt window, not at a GAUSS prompt.
There are two steps to using the GAUSS Profiler: collection and analysis.
To collect profiling information, you must run your GAUSS program in tcollect, an executable
shipped with GAUSS that is identical to tgauss except that it generates a file containing profiling
information each time it is run:
24-1
GAUSS User Guide
tcollect -b myfile.e
The output displayed by tcollect includes the name of the output file containing the profiling
information. tcollect output files have a gaussprof prefix and a .gco extension.
Note that running tcollect on long programs may generate a very large .gco output file. Thus
you may want to delete the .gco files on your machine regularly.
24.1.2 Analysis
To analyze the information stored in the tcollect output file, you must run the gaussprof
executable, which is also shipped with GAUSS, on that file. gaussprof produces an organized
report, displaying the time usage by procedure and by line.
gaussprof gaussprof_001.gco
-h suppress headers
-sp order procedure call sort order where order contains one or more of the folllowing:
24-2
The GAUSS Profiler
e exclusive time
t total time
c number of times called
p procedure name
a ascending order
d descending order (default)
-sl order line number sort order where order contains one or more of the folllowing:
The default, with no flags, is: -pl -sp dep -sl dtf.
Profiler
24-3
PQG
Publication Quality Graphics 25
GAUSS Publication Quality Graphics (PQG) is a set of routines built on the graphics functions
in GraphiC by Scientific Endeavors Corporation.
The main graphics routines include xy, xyz, surface, polar and log plots, as well as histograms,
bar, and box graphs. Users can enhance their graphs by adding legends, changing fonts, and
adding extra lines, arrows, symbols and messages.
The user can create a single full size graph, inset a smaller graph into a larger one, tile a window
with several equally sized graphs or place several overlapping graphs in the window. Graphic
panel size and location are all completely under the user’s control.
GAUSS PQG consists of a set of main graphing procedures and several additional procedures and
global variables for customizing the output.
All of the actual output to the window happens during the call to these main routines:
25-1
GAUSS User Guide
hist Histogram.
xy Cartesian graph.
There are four basic parts to a graphics program. These elements should be in any program that
uses graphics routines. The four parts are the header, data setup, graphics format setup, and
graphics call.
25-2
Publication Quality Graphics
PQG
Header
In order to use the graphics procedures, the pgraph library must be activated. This is done in the
library statement at the top of your program or command file. The next line in your program will
typically be a command to reset the graphics global variables to their default state. For example:
Data Setup
x = seqa(1,1,50);
y = sin(x);
Most of the graphics elements contain defaults that allow the user to generate a plot without
modification. These defaults, however, may be overridden by the user through the use of global
variables and graphics procedures. Some of the elements that may be configured by the user are
axes numbering, labeling, cropping, scaling, line and symbol sizes and types, legends, and colors.
The graphics routines take as input the user data and global variables that have previously been
set. It is in these routines where the graphics file is created and displayed.
Following are three PQG examples. The first two programs are different versions of the same
graph. The variables that begin with _p are the global control variables used by the graphics
25-3
GAUSS User Guide
routines. (For a detailed description of these variables, see G C V, Section
25.6.
Example 1 The routine being called here is a simple XY plot. The entire window will be used.
Four sets of data will be plotted with the line and symbol attributes automatically selected. This
graph will include a legend, title, and a time/date stamp (time stamp is on by default):
Example 2 Here is the same graph with more of the graphics format controlled by the user. The
first two data sets will be plotted using symbols at data points only (observed data); the data points
in the second two sets will be connected with lines (predicted results):
25-4
Publication Quality Graphics
PQG
xlabel("X Axis"); /* X axis label */
title("Example xy Graph"); /* main title */
xy(x,y); /* call to main routine */
Example 3 In this example, two graphics panels are drawn. The first is a full-sized surface
representation, and the second is a half-sized inset containing a contour of the same data located in
the lower left corner of the window:
25-5
GAUSS User Guide
While the structure has changed somewhat, the four basic elements of the graphics program are all
here. The additional routines begwind, endwind, makewind, nextwind, and setwind are all
used to control the graphic panels.
As Example 3 illustrates, the code between graphic panel functions (that is, setwind or
nextwind) may include assignments to global variables, a call to graphset, or may set up new
data to be passed to the main graphics routines.
You are encouraged to run the example programs supplied with GAUSS. Analyzing these
programs is perhaps the best way to learn how to use the PQG system. The example programs are
located on the examples subdirectory.
PQG uses a 4190×3120 pixel resolution grid on a 9.0×6.855-inch printable area. There are three
units of measure supported with most of the graphics global elements:
Inch Coordinates
Inch coordinates are based on the dimensions of the full-size 9.0×6.855-inch output page. The
origin is (0,0) at the lower left corner of the page. If the picture is rotated, the origin is at the upper
left. (For more information, see I U G P, Section 25.3.5.)
Plot Coordinates
Plot coordinates refer to the coordinate system of the graph in the units of the user’s X, Y and Z
axes.
Pixel Coordinates
Pixel coordinates refer to the 4096×3120 pixel coordinates of the full-size output page. The origin
is (0,0) at the lower left corner of the page. If the picture is rotated, the origin is at the upper left.
25-6
Publication Quality Graphics
PQG
25.3 Graphic Panels
Multiple graphic panels for graphics are supported. These graphic panels allow the user to display
multiple graphs on one window or page.
A graphic panel is any rectangular subsection of the window or page. Graphc panels may be any
size and position on the window and may be tiled or overlapping, transparent or nontransparent.
Tiled graphic panels do not overlap. The window can easily be divided into any number of tiled
graphic panels with the window command. window takes three parameters: number of rows,
number of columns, and graphic panel attribute (1=transparent, 0=nontransparent).
This example will divide the window into six equally sized graphic panels. There will be two rows
of three graphic panels–three graphic panels in the upper half of the window and three in the lower
half. The attribute value of 0 is arbitrary since there are no other graphic panels beneath them.
window(nrows,ncols,attr);
window(2,3,0);
Overlapping graphic panels are laid on top of one another as they are created, much as if you were
using the cut and paste method to place several graphs together on one page. An overlapping
graphic panel is created with the makewind command.
In this example, makewind will create an overlapping graphic panel that is 4 inches wide by 2.5
inches tall, positioned 1 inch from the left edge of the page and 1.5 inches from the bottom of the
page. It will be nontransparent:
makewind(hsize,vsize,hpos,vpos,attr);
25-7
GAUSS User Guide
window(2,3,0);
makewind(4,2.5,1,1.5,0);
A nontransparent graphic panel is one that is blanked before graphics information is written to it.
Therefore, information in any previously drawn graphic panels that lie under it will not be visible.
A transparent graphic panel is one that is not blanked, allowing the graphic panel beneath it to
“show through”. Lines, symbols, arrows, error bars, and other graphics objects may extend from
one graphic panel to the next by using transparent graphic panels. First, create the desired graphic
panel configuration. Then create a full-window, transparent graphic panel using the makewind or
window command. Set the appropriate global variables to position the desired object on the
transparent graphic panel. Use the draw procedure to draw it. This graphic panel will act as a
transparent “overlay” on top of the other graphic panels. Transparent graphic panels can be used to
add text or to superimpose one graphic panel on top of another.
25-8
Publication Quality Graphics
PQG
nextwind Set to next available graphic panel number.
This example creates four tiled graphic panels and one graphic panel that overlaps the other four:
library pgraph;
graphset;
begwind;
25-9
GAUSS User Guide
Some global variables allow coordinates to be input in inches. If a coordinate value is in inches
and is being used in a graphic panel, that value will be scaled to “graphic panel inches” and
positioned relative to the lower left corner of the graphic panel. A “graphic panel inch” is a true
inch in size only if the graphic panel is scaled to the full window, otherwise X coordinates will be
scaled relative to the horizontal graphic panel size and Y coordinates will be scaled relative to the
vertical graphic panel size.
The functions savewind and loadwind allow the user to save graphic panel configurations. Once
graphic panels are created (using makewind and window), savewind may be called. This will
save to disk the global variables containing information about the current graphic panel
configuration. To load this configuration again, call loadwind. (See loadwind in the GAUSS
L R.
Graphics text elements, such as titles, messages, axes labels, axes numbering, and legends, can be
modified and enhanced by changing fonts and by adding superscripting, subscripting, and special
mathematical symbols.
To make these modifications and enhancements, the user can embed “escape codes” in the text
strings that are passed to title, xlabel, ylabel and asclabel or assigned to _pmsgstr and
_plegstr.
25-10
Publication Quality Graphics
PQG
The escape code \L (or \l) can be embedded into title strings to create a multiple line title:
or
Four fonts are supplied with the Publication Quality Graphics system. They are Simplex,
Complex, Simgrma, and Microb. (For a list of the characters available in each font, see Appendix
A.)
Fonts are loaded by passing to the fonts procedure a string containing the names of all fonts to be
loaded. For example, this statement will load all four fonts:
25-11
GAUSS User Guide
The fonts command must be called before any of the fonts can be used in text strings. A font can
then be selected by embedding an escape code of the form “\20n” in the string that is to be written
in the new font. The n will be 1, 2, 3 or 4, depending on the order in which the fonts were loaded
in fonts. If the fonts were loaded as in the previous example, the escape characters for each
would be:
\201 Simplex
\202 Complex
\203 Microb
\204 Simgrma
The following example demonstrates how to select a font for use in a string:
Once a font is selected, all succeeding text will use that font until another font is selected. If no
fonts are selected by the user, a default font (Simplex) is loaded and selected automatically for all
text work.
The following examples illustrate the use of the Simgrma font; they assume that Simgrma was the
fourth font loaded. (For the available Simgrma characters and their numbers, see Appendix A.)
The Simgrma characters are specified by either:
25-12
Publication Quality Graphics
PQG
R
For example, to get an integral sign “ ” in Simgrma, embed either a “\044” or a “,” in a string
that has been set to use the Simgrma font.
To produce the title f (x) = sin2 (πx), the following title string should be used:
title("\201f(x) = sin[2](\204p\201x)");
To number the major X axis tick marks with multiples of π/4, the following could be passed to
asclabel:
xtics is used to make sure that major tick marks are placed in the appropriate places.
This example will number the X axis tick marks with the labels µ−2 , µ−1 , 1, µ, and µ2 :
This example illustrates the use of several of the special Simgrma symbols:
This produces:
Z
√
e−µ /2 dµ
2
1/2π
25-13
GAUSS User Guide
25.5 Colors
The following global variables are used to control various graphics elements. Default values are
provided. Any or all of these variables can be set before calling one of the main graphing routines.
The default values can be modified by changing the declarations in pgraph.dec and the
statements in the procedure graphset in pgraph.src. graphset can be called whenever the
user wants to reset these variables to their default values.
_pageshf 2×1 vector, the graph will be shifted to the right and up if this is not 0. If this is
0, the graph will be centered on the output page. Default is 0.
Note: Used internally. (For the same functionality, see makewind in the
GAUSS L R.) This is used by the graphic panel routines. The
user must not set this when using the graphic panel procedures.
_pagesiz 2×1 vector, size of the graph in inches on the printer output. Maximum size is
9.0×6.855 inches (unrotated) or 6.855×9.0 inches (rotated). If this is 0, the
maximum size will be used. Default is 0.
Note: Used internally. (For the same functionality, see makewind in the
GAUSS L R). This is used by the graphic panel routines. The
user must not set this when using the graphic panel procedures.
_parrow M×11 matrix, draws one arrow per row of the input matrix (for total of M
arrows). If scalar zero, no arrows will be drawn.
25-14
Publication Quality Graphics
PQG
[M,1] x starting point.
[M,2] y starting point.
[M,3] x ending point.
[M,4] y ending point.
[M,5] ratio of the length of the arrow head to half its width.
[M,6] size of arrow head in inches.
[M,7] type and location of arrow heads. This integer number will be
interpreted as a decimal expansion mn, for example: if 10, then m = 1, n =
0.
25-15
GAUSS User Guide
[M,11] controls thickness of lines used to draw arrow. This value may be zero
or greater. A value of zero is normal line width.
_parrow = { 1 1 2 2 3 0.2 11 10 2 6 0,
3 4 2 2 3 0.2 11 10 2 6 0 };
_parrow3 M×12 matrix, draws one 3-D arrow per row of the input matrix (for a total of M
arrows). If scalar zero, no arrows will be drawn.
25-16
Publication Quality Graphics
PQG
[M,11] line type:
1 dashed
2 dotted
3 short dashes
4 closely spaced dots
5 dots and dashes
6 solid
[M,12] controls thickness of lines used to draw arrow. This value may be zero
or greater. A value of zero is normal line width.
_parrow3 = { 1 1 1 2 2 2 3 0.2 11 10 6 0,
3 4 5 2 2 2 3 0.2 11 10 6 0 };
_paxes scalar, 2×1, or 3×1 vector for independent control for each axis. The first
element controls the X axis, the second controls the Y axis, and the third (if set)
controls the Z axis. If 0 the axis will not be drawn. Default is 1.
If this is a scalar, it will be expanded to that value.
For example:
_paxht scalar, size of axes labels in inches. If 0, a default size will be computed.
Default is 0.
_pbartyp 1×2 or K×2 matrix. Controls bar shading and colors in bar graphs and
histograms.
The first column controls the bar shading:
0 no shading
1 dots
25-17
GAUSS User Guide
2 vertical cross-hatch
3 diagonal lines with positive slope
4 diagonal lines with negative slope
5 diagonal cross-hatch
6 solid
The second column controls the bar color, see C, Section 25.5.
_pbarwid scalar, width of bars in bar graphs and histograms. The valid range is 0-1. If 0,
the bars will be a single pixel wide. If 1, the bars will touch each other. Default
is 0.5, so the bars take up about half the space open to them.
_pbox scalar, draws a box (border) around the entire graph. Set to desired color of box
to be drawn. Use 0 if no box is desired. Default is 0.
_pboxctl 5×1 vector, controls box plot style, width, and color. Used by procedure box
only.
[1] box width between 0 and 1. If 0, the box plot is drawn as two vertical
lines representing the quartile ranges with a filled circle representing the
50th percentile.
[2] box color, see C, Section 25.5. If 0, the colors may be individually
controlled using global variable _pcolor.
[3] min/max style for the box symbol. One of the following:
1 minimum and maximum taken from the actual limits of the data.
Elements 4 and 5 are ignored.
2 statistical standard with the minimum and maximum calculated
according to interquartile range as follows:
intqrange = 75th − 25th
min = 25th − 1.5 intqrange
max = 75th + 1.5 intqrange
25-18
Publication Quality Graphics
PQG
_pboxlim 5×M output matrix containing computed percentile results from procedure box.
M corresponds to each column of input y data.
_pcolor scalar or K×1 vector, colors for main curves in xy, xyz and log graphs. To use
a single color set for all curves set this to a scalar color value. If 0, use default
colors. Default is 0.
The default colors come from a global vector called _pcsel. This vector can be
changed by editing pgraph.dec to change the default colors, see C,
Section 25.5 (_pcsel is not documented elsewhere).
_pcrop scalar or 1×5 vector, allows plot cropping for different graphic elements to be
individually controlled. Valid values are 0 (disabled) or 1 (enabled). If cropping
is enabled, any graphical data sent outside the axes area will not be drawn. If
this is a scalar, it is expanded to a 1×5 vector using the given value for all
elements. All cropping is enabled by default.
This example will crop main curves, and lines and circles drawn by _pline.
_pcrop = { 1 1 0 1 0 };
_pcross scalar. If 1, the axes will intersect at the (0,0) X-Y location if it is visible.
Default is 0, meaning the axes will be at the lowest end of the X-Y coordinates.
25-19
GAUSS User Guide
_pdate date string. If this contains characters, the date will be appended and printed.
The default is set as follows (the first character is a font selection escape code):
If this is set to a null string, no date will be printed. (For more information on
using fonts within strings, see G T E, Section 25.4.
_perrbar M×9 matrix, draws one error bar per row of the input matrix. If scalar 0, no
error bars will be drawn. Location values are in plot coordinates.
[M,1] x location.
[M,2] left end of error bar.
[M,3] right end of error bar.
[M,4] y location.
[M,5] bottom of error bar.
[M,6] top of error bar.
[M,7] line type:
1 dashed
2 dotted
3 short dashes
4 closely spaced dots
5 dots and dashes
6 solid
[M,8] color, see C, Section 25.5.
[M,9] line thickness.. This value may be 0 or greater. A value of 0 is normal
line width.
_perrbar = { 1 0 2 2 1 3 6 2 0 };
_pframe 2×1 vector, controls frame around axes area. On 3-D plots this is a cube
surrounding the 3-D workspace.
25-20
Publication Quality Graphics
PQG
[1] 1 frame on
0 frame off
[2] 1 tick marks on frame
0 no tick marks
_plctrl scalar or K×1 vector to control whether lines and/or symbols will be displayed
for the main curves. This also controls the frequency of symbols on main
curves. The number of rows (K) is equal to the number of individual curves to
be plotted in the graph. Default is 0.
This example draws a line for the first curve, draws a line and plots a symbol
every 10 data points for the second curve, and plots symbols only every 5 data
points for the third curve:
_plctrl = { 0, 10, -5 };
25-21
GAUSS User Guide
_plegctl = 1;
This example creates a smaller legend and positions it 2.5 inches from the left
and 1 inch from the bottom.
_plegctl = { 2 3 2.5 1 };
_plegstr string, legend entry text. Text for multiple curves is separated by a null byte
(“\000”).
For example:
_plev M×1 vector, user-defined contour levels for contour. Default is 0. (See
contour in the GAUSS L R.)
_pline M×9 matrix, to draw lines, circles, or radii. Each row controls one item to be
drawn. If this is a scalar zero, nothing will be drawn. Default is 0.
25-22
Publication Quality Graphics
PQG
1 line in plot coordinates
2 line in inch coordinates
3 line in pixel coordinates
4 circle in plot coordinates
5 circle in inch coordinates
6 radius in plot coordinates
7 radius in inch coordinates
[M,2] line type:
1 dashed
2 dotted
3 short dashes
4 closely spaced dots
5 dots and dashes
6 solid
[M,3-7] coordinates and dimensions:
if item type is line (1<=_pline[M,1]<=3):
[M,3] x starting point.
[M,4] y starting point.
[M,5] x ending point.
[M,6] y ending point.
[M,7] 0 if this is a continuation of a curve, 1 if this begins a new curve.
25-23
GAUSS User Guide
_pline3d M×9 matrix. Allows extra lines to be added to an xyz or surface graph in 3-D
plot coordinates.
_plotshf 2×1 vector, distance of plot from lower left corner of output page in inches.
[1] x distance.
[2] y distance.
If scalar 0, there will be no shift. Default is 0.
Note: Used internally. (For the same functionality, see axmargin in the
GAUSS L R.) This is used by the graphic panel routines. The
user must not set this when using the graphic panel procedures.
25-24
Publication Quality Graphics
PQG
_plotsiz 2×1 vector, size of the axes area in inches. If scalar 0, the maximum size will
be used.
Note: Used internally. (For the same functionality, see axmargin in the
GAUSS L R.) This is used by the graphic panel routines. The
user must not set this when using the graphic panel procedures.
_pltype scalar or K×1 vector, line type for the main curves. If this is a nonzero scalar, all
lines will be this type. If scalar 0, line types will be default styles. Default is 0.
1 dashed
2 dotted
3 short dashes
4 closely spaced dots
5 dots and dashes
6 solid
The default line types come from a global vector called _plsel. This vector
can be changed by editing pgraph.dec to change the default line types
(_plsel is not documented elsewhere.)
_plwidth scalar or K×1 vector, line thickness for main curves. This value may be zero or
greater. A value of zero is normal (single pixel) line width. Default is 0.
_pmcolor 9×1 vector, color values to use for plot, see C, Section 25.5.
[1] axes.
[2] axes numbers.
[3] X axis label.
[4] Y axis label.
[5] Z axis label.
[6] title.
[7] box.
[8] date.
[9] background.
25-25
GAUSS User Guide
_pmsgctl L×7 matrix of control information for printing the strings contained in
_pmsgstr.
_pmsgstr string, contains a set of messages to be printed on the plot. Each message is
separated from the next by a null byte (\000). The number of messages must
correspond to the number of rows in the _pmsgctl control matrix. This can be
created as follows:
_pnotify scalar, controls window output during the creation of the graph. Default is 1.
_pnum scalar, 2×1 or 3×1 vector for independent control for axes numbering. The first
element controls the X axis numbers, the second controls the Y axis numbers,
and the third (if set) controls the Z axis numbers. Default is 1.
If this value is scalar, it will be expanded to a vector.
25-26
Publication Quality Graphics
PQG
2 axes numbers displayed, horizontally oriented on axis
For example:
_pnumht scalar, size of axes numbers in inches. If 0, a size of .13 will be used. Default is
0.
_psilent scalar. If 0, a beep will sound when the graph is finished drawing to the
window. Default is 1 (no beep).
_pstype scalar or K×1 vector, controls symbol used at data points. To use a single
symbol type for all points, set this to one of the following scalar values:
If this is a vector, each line will have a different symbol. Symbols will repeat if
there are more lines than symbol types. Default is 0 (no symbols are shown).
25-27
GAUSS User Guide
[M,1] x location.
[M,2] y location.
[M,3] symbol type, see _pstype earlier.
[M,4] symbol height. If this is 0, a default height of 5.0 will be used.
[M,5] symbol color, see C, Section 25.5.
[M,6] type of coordinates:
1 plot coordinates
2 inch coordinates
[M,7] line thickness. A value of zero is normal line width.
_psym3d M×7 matrix for plotting extra symbols on a 3-D (surface or xyz) graph.
_psymsiz scalar or K×1 vector, symbol size for the symbols on the main curves. This is
NOT related to _psym. If 0, a default size of 5.0 is used.
_ptek string, name of Tektronix format graphics file. This must have a .tkf
extension. If this is set to a null string, the graphics file will be suppressed. The
default is graphic.tkf.
_ptitlht scalar, the height of the title characters in inches. If this is 0, a default height of
approx. 0.13 inch will be used.
25-28
Publication Quality Graphics
PQG
_pversno string, the graphics version number.
_pxpmax scalar, the maximum number of places to the right of the decimal point for the
X axis numbers. Default is 12.
_pxsci scalar, the threshold in digits above which the data for the X axis will be scaled
and a power of 10 scaling factor displayed. Default is 4.
_pypmax scalar, the maximum number of places to the right of the decimal point for the
Y axis numbers. Default is 12.
_pysci scalar, the threshold in digits above which the data for the Y axis will be scaled
and a power of 10 scaling factor displayed. Default is 4.
_pzclr scalar, row vector, or K×2 matrix, Z level color control for procedures surface
and contour. (See surface in the GAUSS L R.)
_pzoom 1×3 row vector, magnifies the graphics display for zooming in on detailed areas
of the graph. If scalar 0, no magnification is performed. Default is 0.
To see the upper left quarter of the screen magnified 2 times use:
_pzoom = { 2 25 75 };
_pzpmax scalar, the maximum number of places to the right of the decimal point for the
Z axis numbers. Default is 3.
_pzsci scalar, the threshold in digits above which the data for the Z axis will be scaled
and a power of 10 scaling factor displayed. Default is 4.
25-29
Editor
Graphics
Graphics Editor 26
26.1 Introduction to the Graphics Editor
The GAUSS graphics editor is a utility for composing pages containing GAUSS graphics files. Its
primary purpose is to provide the user with a toolbox for creating and annotating graphs created
by GAUSS using all of the fonts available on your Windows system. It is not meant to be a
full-featured publishing tool but rather a supplemental utility for dynamically importing and easily
arranging multiple graphics files on a single page.
26.1.1 Overview
The graphics editor allows the user to interactively create any number of graphical objects for
composing documents. It is launched by selecting Tools from the GAUSS menu bar, then
Graphics Editor, or by clicking on the Graphics Editor icon on the GAUSS toolbar.
Once the document has been created, it may be saved for later modification. All of the objects and
their respective properties contained in the document are preserved. The document may also be
exported to other formats.
26-1
GAUSS User Guide
The graphics editor workspace is a window allowing access to a single page with tools for
composing the document. The page is defined by user-defined properties such as page orientation
and margin settings.
It provides a dialog bar for user-selection of the current pen and brush properties.
Zoom capability is provided for detailed accuracy and accomodating a wide-variety of computer
display resolutions.
26.2.1 Toolbar
The toolbar is displayed across the top of the application window, below the menu bar. The toolbar
provides quick mouse access to many tools used in the graphics editor.
26-2
Graphics Editor
To hide or display the toolbar, choose Toolbar from the View menu (ALT, V, T).
Editor
Graphics
Figure 26.2: Graphics Editor Toolbar
The status bar is displayed at the bottom of the graphics editor window. To display or hide the
status bar, use the Status Bar command in the View menu.
The left area of the status bar describes actions of menu items as you use the arrow keys to
navigate through menus. This area similarly shows messages that describe the actions of toolbar
buttons as you depress them, before releasing them. If after viewing the description of the toolbar
command you wish not to execute the command, then release the mouse button while the pointer
is off the toolbar button.
In addition, the status bar provides short hints while using the graphical interface such as defining,
sizing, and moving objects.
Indicator Description
26-3
GAUSS User Guide
The x/y mouse position on the page (in units specified in the View Properties menu).
Import Imports a file of another format. Currently only the GAUSS graphics format
.tkf is supported.
Save Saves the active document to its current name and directory. When you save
a document for the first time, the graphics editor displays the Save As dialog
box so you can name your document (CTRL+S). If you want to change the
name and directory of an existing document before saving it, choose the
Save As command.
Save As Saves and names the active document. The graphics editor displays the Save
As dialog box so you can name your document.
Print Prints a document. This command presents a Print dialog box where you
may specify the ranges of pages to be printed, the number of copies, the
destination printer, and other printer setup options (CTRL+P).
Print Preview Displays the active document as it would appear when printed. When you
choose this command, the main window will be replaced with a print
preview window in which one or two pages will be displayed in their printed
format. The print preview toolbar offers you the option to view either one or
two pages at a time, move back and forth through the document, zoom in
and out of pages, and initiate a print job.
26-4
Graphics Editor
Print Setup Allows you to select a printer and printer connection. This command
presents a Print Setup dialog box where you specify the printer and its
connection.
Exit Ends the graphics editor session. The graphics editor prompts you to save
Editor
Graphics
documents with unsaved changes (ALT+F4).
Cut Removes the currently selected data from the document and put it on the
clipboard (CTRL+X). This command is unavailable if there is no data
currently selected. Cutting data to the clipboard replaces the contents
previously stored there.
Copy Copies currently selected data onto the clipboard (CTRL+C). This command
is unavailable if there is no data currently selected. Copying data to the
clipboard replaces the contents previously stored there.
Paste Inserts a copy of the clipboard contents at the insertion point (CTRL+V).
This command is unavailable if the clipboard is empty.
Toolbar Displays and hides the Toolbar, which includes buttons for some of the most
common commands such as File Open. A check mark appears next to the
26-5
GAUSS User Guide
menu item when the Toolbar is displayed. See T, Section ??, for help
on using the toolbar.
Status Bar Displays and hides the Status Bar, which describes the action to be executed
by the selected menu or depressed toolbar button and keyboard latch state. A
check mark appears next to the menu item when the Status Bar is displayed.
See S B, Section ??, for help on using the status bar.
Properties Allows you to change user-defined page/view settings. See S
P/V P, Section 26.2.11, for more information.
Zoom Allows you to change user-defined zoom control. See U Z
F, Section 26.2.11, for more information.
The Draw menu offers the following commands; see G O, Section 26.2.13, for more
detailed information on each.
26-6
Graphics Editor
The Export menu offers the following commands; see F M, Section 26.3, for more
detailed information on each.
Editor
Graphics
Encapsulated Writes an Encapsulated Postscript file.
Postscript
The Help menu offers the following commands, which provide you assistance with this
application:
Help Topics Displays the opening screen of Help. From the opening screen, you can
jump to step-by-step instructions for using the graphics editor and various
types of reference information. Once you open Help, you can click the
Contents button whenever you want to return to the opening screen.
About Displays the copyright notice and version number of this application.
Once an object has been selected, its action context menu may be displayed by right-clicking
inside the object.
26-7
GAUSS User Guide
Cut Removes the currently selected data from the document and put it on the
clipboard. This command is unavailable if there is no data currently selected.
Cutting data to the clipboard replaces the contents previously stored there.
Copy Copies currently selected data onto the clipboard. This command is
unavailable if there is no data currently selected.
Delete Deletes currently selected data. This command is unavailable if there is no
data currently selected.
Z-Order Changes the objects position in the z-order of the document’s list. The
z-order allows the user to control in what order the object is drawn on the
page. To change, select the Z-Order menu item from the action context menu
and select one of the following:
Move to Top - Moves the object to the top of the list.
Move to Bottom - Moves the object to the bottom of the list.
Edit Allows you to modify the object. This menu item currently applies to text
objects only.
Deselect De-selects the object.
Properties Opens the object’s property dialog.
The Page Context menu is displayed by pressing the right mouse button when no object is selected.
Paste Copies an object from the clipboard to the page if one is available.
Retain Aspect Check or uncheck the aspect ratio state. When checked, this forces the
object to retain its aspect ratio while sizing it from the top or sides; sizing
Ratio from the corners overrides this setting.
See M G O, Section 26.2.14, for more information.
26-8
Graphics Editor
The following describes how to set various page and view properties, including how to set the
page orientation and margins, use the zoom feature, and set the color options.
Editor
Graphics
Setting the Page Orientation and Margins
Page Orientation The document page orientation of landscape or portrait is set from the
Properties dialog under the View menu.
Reset Colors Pushing this button will reset the available colors to the initial IBM 16-color
scheme. Each color may be set to a custom color; see P/F P,
Section 26.2.12 for more information.
To set the zoom, click the zoom drop-down control on the toolbar or select Zoom in the View
menu.
You may reset the colors to the original IBM 16-color scheme by pressing the Reset Colors to
Initial button.
26-9
GAUSS User Guide
All drawing is done with a currently selected pen and brush. The current pen has a width attribute
and color attribute.
The current brush (for painting object backgrounds) has a color attribute.
Left-click in one of the color boxes in the dialog bar to set the current color. The dialog bar is
located to the left of the drawing area. The color is immediately displayed in the sample box at the
top of the color box area.
Double-clicking in the color box will allow you you to customize that particular color.
Left-click in one of the width boxes in the dialog bar shown below the color boxes.
Right-click in one of the color boxes in the dialog bar to set the current fill color. The dialog bar is
located to the left of the drawing area. The fill color is immediately displayed in the sample box at
the top of the color box area.
Transparent Fill To set a transparent fill color, right-click in the sample box at the top of the
color boxes. This will cause the drawing object to have no fill associated
with it and allow objects beneath it to show through.
Customizing Double-clicking in the color box will allow you to customize that particular
the Color color.
26-10
Graphics Editor
The graphics editor allows the user to interactively create any number of the following graphical
objects for composing your document:
Editor
Graphics
Creating a TKF Graphics window
To create a TKF graphics window, select the Graphics Window menu item from the Draw menu or
press the Create TKF window icon on the toobar.
A properties dialog is presented which allows the selection of a GAUSS-generated TKF graphics
file. There are two ways to do this. If one or more graphics files are currently being displayed in
GAUSS, those filenames will appear in the Active Graphs drop-down control and may be selected.
Otherwise, pressing the Browse button will present a common open file dialog for selecting the
file.
The graphics window border and fill colors may be set using their respective Color buttons.
Once the file has been selected and the user presses OK, the graphics window is created in a
default size located at the top-left corner of the page.
To create a text object, select the Text menu item from the Draw menu or press the Draw text
button on the toolbar.
Next, position the mouse where you want the top-left corner of your text then press and hold the
left mouse button. (You may also move the object by pressing and holding the right mouse button
while keeping the left button depressed). Drag the mouse to the bottom-left corner and let up on
the mouse button.
26-11
GAUSS User Guide
The text region will be redrawn in the current background fill color and contain a text cursor
inside. At this point you may enter your text at the cursor using the last selected text font.
If the text requires more lines than the current bounding box allows, the box will be resized as
needed.
To save your text when you have finished typing, press the OK icon on the text toolbar indicated
by a green check mark or press the SHIFT+ENTER key. Clicking the mouse button outside the
text window will also save the text and complete the operation.
The box will be redrawn with the proper font background, border color and margin settings.
To cancel out of the text and lose changes, press the Cancel icon on the text toolbar indicated by a
red X or press the ESCAPE key.
Note: The text object may be rotated at any angle from the Object Properties menu.
Creating a Line
To create a line, select the Line menu item from the Draw menu or press the Draw line button on
the toolbar. This puts the editor into the draw line state indicated by a crosshair cursor.
Next, position the mouse where you want the first end point of the line then press and hold the left
mouse button. (You may also move the line by pressing and holding the right mouse button while
keeping the left button depressed). Drag the mouse to the location for the second end point and
release the mouse button.
Note: Pressing the CTRL key while defining a line or arrow forces the line to be vertical or
horizontal.
At any time after the above process you may modify the line object.
26-12
Graphics Editor
Creating an Arrow
To create an arrow, select the Arrow menu item from the Draw menu or press the Draw arrow
button on the toolbar. This puts the editor into the draw arrow state indicated by a crosshair cursor.
Editor
Graphics
Arrow style controls for defining the arrow appear on the dialog bar to the left of the drawing area.
Defining the endpoints are identical to the steps for defining a line above. However, the arrow head
size and shape may be set using the additional arrow style controls in the dialog bar.
Once the second endpoint is defined, the arrow will be redrawn in the current pen color.
At any time after the above process you may modify the arrow object.
Arrow Styles
There are a combination of two styles of arrows: Open/Closed, and Single/Double-headed arrows.
A closed arrow is one whose arrowhead is filled in with the current pen color. An open one has no
fill.
A single-headed arrow is a line with an arrowhead on one end. A double-headed arrow has an
arrowhead at both ends.
Creating a Rectangle
To create a rectangle object, select the Rectangle menu item from the Draw menu or press the
Draw rectangle button on the toolbar. This puts the editor into the draw rectangle state indicated
by a crosshair cursor.
Next, position the mouse where you want the top-left corner then press and hold the left mouse
button. (You may also move the object by pressing and holding the right mouse button while
keeping the left button depressed). Drag the mouse to the bottom-left corner and release the mouse
button.
The rectangle will be redrawn in the current background fill and border color.
26-13
GAUSS User Guide
At any time after the above process you may modify the rectangle object.
Creating an Ellipse
Note: a circle is first created when defining an ellipse. After the circle has been defined it may be
dynamically reshaped into an ellipse of the desired size using the mouse.
To create the circle, select the Ellipse menu item from the Draw menu or press the Draw ellipse
button on the toolbar. This puts the editor into the draw ellipse state indicated by a crosshair cursor.
Next, position the mouse where you want the center of the circle then press and hold the left
mouse button. (You may also move the object by pressing and holding the right mouse button
while keeping the left button depressed). Drag the mouse to the desired radius and release the
mouse button.
The circle will be redrawn in the current background fill and border color.
At any time after the above process you may modify the circle to any other elliptical size and
shape.
First, ensure you are in selection mode by pressing the Select toolbar button or choosing the Select
item in the Draw menu. Selection mode is indicated with an arrow cursor.
Next, select the object you want to modify by left-clicking anywhere inside or on the object. It
will then become highlighted.
Once selected, it may be sized, moved, or modified with one of the actions listed in the object’s
action context menu.
26-14
Graphics Editor
Aspect Ratio
You can force an object’s aspect ratio to be retained while sizing it by checking this menu item
from the Page Context menu.
Editor
Graphics
When this is checked, all the sides of the object are sized by the same amount as the side being
moved, eliminating the need to resize all sides independently.
If the object is being sized by a corner point, this feature is ignored, eliminating the need for the
user to check/uncheck the aspect ratio menu item needlessly.
The aspect ratio feature has no effect for lines and arrows.
Object Properties
The properties dialog box allows you to modify various attributes of the object depending on its
type.
See G O, Section 26.2.13, for more information about object properties.
Sizing an Object
Next, click and hold the left mouse button in one of the object’s highlight points. Depending on
the type of object, the new size is defined by how you move the mouse. Lifting the mouse button
sets the new size and causes it to be redrawn.
If the object is a TKF graphics window, text object, rectangle, or ellipse, then the aspect ratio may
be retained depending on the selection state of the aspect ratio menu item state when grabbing one
of the four sides.
Grabbing the corner of an object allows you to size it in any direction ignoring the state of the
aspect ratio menu.
26-15
GAUSS User Guide
If the object is a line or arrow, then it may be moved during the sizing operation by pressing and
holding the right mouse button while still depressing the left mouse button.
Note: A rotated text object may not be sized. It may only be sized in a non-rotated state.
Moving an Object
Next, click and hold the left mouse button somewhere inside the object. Drag the object to the new
location and lift the mouse button. The object is then redrawn in the new location.
If the object is a line or arrow, then it may also be moved during the sizing operation.
The graphics editor stores the document as a list of vector-based graphical objects. These are
binary files and cannot be edited by hand. It uses a default extension of .pge.
You may export your document to other formats for inclusion in web pages, word-processors, and
publishing applications.
The Export menu enables you to easily export graphic files to some of the most frequently used
graphic formats.
26-16
Graphics Editor
To write an Encapsulated Postscript file, select the Encapsulated Postscript menu item from the
Export menu. This displays the entire file dialog.
Editor
Graphics
Filename Enter or browse to the desired output filename. The default extension is
.eps.
Convert lines Check this item if you want to convert all colors in the image to black.
to Black
Scale Factor Enter a scale factor if necessary. By default, the graphics editor uses a very
high internal resolution for the best possible quality. However, some
applications are unable to correctly scale the data when importing. This may
be worked around by scaling the data during the export stage.
Minimum line Enter the minimum line width value if you want to darken the lines in
width the exported file.
Although the fonts you select for your text box may appear fine in the graphics editor, it is possible
the target application importing it may not interpret them correctly. Every application has its own
EPS interpreter, and the availability of your font depends on it. If you are having problems of this
type, try using the Enhanced Metafile format conversion. This format has no such font limitations.
A JPEG image file is a widely used bitmap format for inclusion in web pages due to its
compression characteristics.
To write a JPEG file, select the JPEG Image File menu item from the Export menu.
Note: Because JPEG is a bitmap format, the image written to the file is exactly as seen on the
display. Thus you may be requried to zoom out on some lower-resolution displays to obtain an
image of the entire document.
26-17
GAUSS User Guide
A TIFF (Tag Image File Format) image file is an older but widely used bitmap format.
To write a TIFF file, select the TIFF Image File menu item from the Export menu.
Note: Because TIFF is a bitmap format, the image written to the file is exactly as seen on the
display. Thus you may be required to zoom out on some lower-resolution displays to obtain an
image of the entire document.
Windows Metafile
An enhanced metafile is a vector-based file and is considered the best method for export/import on
the Microsoft Windows platform.
To create a Windows Enhanced Metafile, select the Windows Metafile menu item from the Export
menu.
Autoscale
Checking this option forces the translated to automatically scale the Enhanced Metafile. This is
the best setting for most applications for importing. However, some applications require a more
precise format. If the importing application has trouble, uncheck this option.
26-18
Time and Date 27
In the year 1 AD the calendar in general use was the Julian calendar. The Gregorian calendar that
we use today was not invented until the late 1500’s. This new calendar changed the method of
calculating leap years on century marks. With the Julian system simply every fourth year was a
leap year. The Gregorian system made every fourth year a leap year with the exception of century
marks which are only leap years if divisible by 400. The British adoption of this calendar, which
the GAUSS date functions are based on, did not happen until the year 1752. In that year eleven
days were removed; September 2, 1752 was followed by September 14, 1752.
dtvnormal and utctodtv are accurate back to 1 AD. The rest of the GAUSS date functions
assume a normal Gregorian system regardless of year. Thus, they will not account for the days
taken out in September of 1752, nor will they account for all century marks being leap years
before the adoption of the Gregorian system in 1752.
The time is given by your operating system, daylight savings time is not automatically accounted
for by GAUSS in calculations.
27-1
GAUSS User Guide
The Time and Date formats in GAUSS fall into one of two major categories, matrix/vector and
string. The matrix/vector formats can be used for either calculations or if desired for output. The
string formats are, however, mostly for use as ouput. Some manipulation of strings is possible
with the use of the stof function.
d = date;
d;
1997.00 /* Year */
5.00000 /* Month */
29.0000 /* Day */
56.4700 /* Hundredths of a second since midnight */
t = time;
t;
These vectors can be written to a string of the desired form by passing them through the
corresponding function.
5/29/97
datestrymd(d);
27-2
Time and Date
19970529
10:17:33
A list and brief description of these, and other related functions is provided in the table in section
27.2.
Where:
dtvnormal normalizes a date. The last two elements are ignored for input, as shown in the
following example. They are set to the correct values on output. The input can be 1×8 or N×8.
dtv = { 1954 3 17 4 16 0 0 0 };
dtv = dtvnormal(dtv);
27-3
GAUSS User Guide
1954 3 17 4 16 0 3 75
1954 3 417 4 16 0 3 75
dtv = dtvnormal(dtv);
print dtv;
1955 4 21 4 16 0 4 110
Following is a partial listing of the time and date functions available in GAUSS.
27-4
Time and Date
di f = 1.1274488e + 09
di f = 130.00000
_isleap returns a matrix of ones and zeros, ones when the corresponding year is a leap year.
27-5
GAUSS User Guide
x = seqa(1970,1,20);
y = _isleap(x);
delif(x,abs(y-1));
x = seqa(1983,1,3);
y = _daypryr(x);
sumc(y);
1096.0000
g = { 1986, 2, 23, 0 };
dy = dayinyr(g);
sumc(y)+dy;
1150.0000
For more information on any of these functions see their respective pages in the command
reference.
Iterations of a program can be timed with the use of the hsec function in the following manner.
27-6
Time and Date
In the case of a program running from one day into the next you would need to replace the hsec
function with the date function. The ethsec function should be used to compute the time
difference; a straight subtraction as in the previous example will not give the desired result.
27-7
ATOG 28
ATOG
ATOG is a stand-alone conversion utility that converts ASCII files into GAUSS data sets. ATOG
can convert delimited and packed ASCII files into GAUSS data sets. ATOG can be run from a
batch file or the command line; it is not run from a GAUSS prompt but rather from a command
prompt window.
atog cmdfile
where cmdfile is the name of the command file. If no extension is given, .cmd will be assumed. If
no command file is specified, a command summary will be displayed.
28-1
GAUSS User Guide
The principle commands for converting an ASCII file that is delimited with spaces or commas are
given in the following example:
input agex.asc;
output agex;
invar $ race # age pay $ sex region;
outvar region age sex pay;
outtyp d;
In this example, a delimited ASCII file agex.asc is converted to a double precision GAUSS data
file agex.dat. The input file has five variables. The file will be interpreted as having five columns:
28-2
ATOG
The output file will have four columns since the first column of the input file (race) is not included
in the output variables. The columns of the output file are:
The variable names are saved in the file header. Unless preservecase has been specified, the
names of character variables will be saved in lowercase, and the names of numeric variables will
be saved in uppercase. The $ in the invar statement specifies that the variables that follow are
character type. The # specifies numeric. If $ and # are not used in an invar statement, the default
is numeric.
ATOG
28.2 Commands
append
append;
No assumptions are made regarding the format of the existing file. Make certain that the number,
order, and type of data converted match the existing file. ATOG creates v96 format data files, so
will only append to v96 format data files.
28-3
GAUSS User Guide
complex
Instructs ATOG to convert the ASCII file into a complex GAUSS data set:
complex;
Complex GAUSS data sets are stored by rows, with the real and imaginary parts interleaved,
element by element. ATOG assumes the same structure for the ASCII input file, and will thus read
TWO numbers out for EACH variable specified.
input
Specifies the file name of the ASCII file to be converted. The full path name can be used in the file
specification.
input data.raw;
The command:
input /research/data/myfile.asc;
28-4
ATOG
invar
Soft Delimited ASCII Files Soft delimited files may have spaces, commas, or cr/lf as delimiters
between elements. Two or more consecutive delimiters with no data between them are treated as
one delimiter. For example:
ATOG
6 VAR02 numeric
7 VAR03 numeric
8 VAR04 numeric
9 VAR05 numeric
10 VAR06 numeric
11 VAR07 numeric
12 VAR08 numeric
13 VAR09 numeric
14 VAR10 numeric
15 X001 numeric
16 X002 numeric
17 X003 numeric
18 X004 numeric
19 X005 numeric
As the input file is translated, the first 19 elements will be interpreted as the first row (observation),
the next 19 will be interpreted as the second row, and so on. If the number of elements in the file is
not evenly divisible by 19, the final incomplete row will be dropped and a warning message will
be given.
28-5
GAUSS User Guide
Hard Delimited ASCII Files Hard delimited files have a printable character as a delimiter
between elements. Two delimiters without intervening data between them will be interpreted as a
missing. If \n is specified as a delimiter, the file should have one element per line and blank lines
will be considered missings. Otherwise, delimiters must be printable characters. The dot ‘.’ is
illegal and will always be interpreted as a missing value. To specify the backslash as a delimiter,
use \\. If \r is specified as a delimiter, the file will be assumed to contain one case or record per
line with commas between elements and no comma at the end of the line.
For hard delimited files the delimit subcommand is used with the invar command. The
delimit subcommand has two optional parameters. The first parameter is the delimiter. The
default is a comma. The second parameter is an ‘N’. If the second parameter is present, ATOG
will expect N delimiters. If it is not present, ATOG will expect N-1 delimiters.
This example:
while
or
28-6
ATOG
If the invar statement specified three variables and N-1 delimiters, this file would be interpreted
as having three rows containing a missing in the [2,1] element and the [3,3] element like this:
ATOG
If N delimiters had been specified, this file would be interpreted as having two rows, and a final
incomplete row that is dropped:
The spaces were shown only for clarity and are not significant in delimited files so:
BILL,222.3,123.2,456.4,345.2,533.2,
STEVE,624.3,340.3,,624.3,639.5,
TOM,244.2,834.3,602.3,333.4,822.5
Linefeeds are significant only if \n is specified as the delimiter, or when using \r. This example:
28-7
GAUSS User Guide
will expect a file with no comma after the final element in each row:
Packed ASCII Files Packed ASCII files must have fixed length records. The record
subcommand is used to specify the record length, and variables are specified by giving their type,
starting position, length, and the position of an implicit decimal point if necessary.
outvar is not used with packed ASCII files. Instead, invar is used to specify only those variables
to be included in the output file.
For packed ASCII files the syntax of the invar command is as follows:
where,
reclen the total record length in bytes, including the final carriage return/line feed if
applicable. Records must be fixed length.
start starting position of the field in the record, 1 is the first position. The default
is 1.
length the length of the field in bytes. The default is 8.
prec optional; a decimal point will be inserted automatically prec places in from
the RIGHT edge of the field.
28-8
ATOG
If several variables are listed after a format definition, each succeeding field will be assumed to
start immediately after the preceding field. If an asterisk is used to specify the starting position, the
current logical default will be assumed. An asterisk in the length position will select the current
default for both length and prec. This is illegal: (3,8.*).
The type change characters $ and # are used to toggle between character and numeric data type.
The examples below assume a 32-byte record with a carriage return/line feed occupying the last 2
bytes of each record. The data below can be interpreted in different ways using different invar
statements:
ABCDEFGHIJ12345678901234567890<CR><LF>
| | | | | |
position 1 10 20 30 31 32
ATOG
This example:
This example:
28-9
GAUSS User Guide
msym
Specifies the character in the input file that is to be interpreted as a missing value. This example:
msym &;
defines the character ‘&’ as the missing value character. The default ‘.’ (dot) will always be
interpreted as a missing value unless it is part of a numeric value.
nocheck
Optional; suppresses automatic checking of packed ASCII record length and output data type. The
default is to increase the record length by 2 bytes if the second record in a packed file starts with
cr/lf, and any files that have explicitly defined character data will be output in double precision
regardless of the type specified.
output
The name of the GAUSS data set. A file will be created with the extension .dat. For example:
output /gauss/dat/test;
28-10
ATOG
outtyp
Selects the numerical accuracy of the output file. Use of this command should be dictated by the
accuracy of the input data and storage space limitations. The format is:
outtyp fmt;
D or 8 double precision
F or 4 single precision (default)
I or 2 integer
ATOG
2 integer 4 −32768<=X<=32767
4 single precision 6–7 8.43x10−37 <=|X|<=3.37x10+38
8 double precision 15–16 4.19x10−307 <=|X|<=1.67x10+308
If the output type is integer, the input numbers will be truncated to integers. If your data has more
than 6 or 7 significant digits, specify outtyp as double.
Character data require outtyp d. ATOG automatically selects double precision when character
data is specified in the invar statement, unless you have specified nocheck.
The precision of the storage selected does not affect the accuracy of GAUSS calculations using the
data. GAUSS converts all data to double precision when the file is read.
outvar
Selects the variables to be placed in the GAUSS data set. The outvar command needs only the
list of variables to be included in the output data set. They can be in any order. In this example:
28-11
GAUSS User Guide
preservecase
Optional; preserves the case of variable names. The default is nopreservcase, which will force
variable names for numeric variables to upper case and character variables to lower case.
28.3 Examples
Example 1 The first example is a soft delimited ASCII file called agex1.asc. The file contains
seven columns of ASCII data:
28-12
ATOG
input /gauss/agex1.asc;
output agex1;
invar $month #temp pres vol $true var[02];
outvar month true temp pres vol;
ATOG
case 1 Jan yes 167.3 822.4 6.34e+6
case 2 Feb no 165.8 987.3 5.63e+6
case 3 Mar yes 165.3 842.3 7.34e+6
type char char numeric numeric numeric
The data set is double precision since character data is explicitly specified.
Example 2 The second example is a packed ASCII file xlod.asc The file contains 32-character
records:
AEGDRFCSTy02345678960631567890<CR><LF>
EDJTAJPSTn12395863998064839561<CR><LF>
GWDNADMSTy19827845659725234451<CR><LF>
| | | | | |
position 1 10 20 30 31 32
28-13
GAUSS User Guide
input /gauss/dat/xlod.asc;
output xlod2;
invar record=32 $(1,3) client[2] zone (*,1) reg #(20,5) zip;
The data set is double precision since character data is explicitly specified.
Example 3 The third example is a hard delimited ASCII file called cplx.asc. The file contains
six columns of ASCII data:
input /gauss/cplx.asc;
output cplx;
invar delimit #cvar[3];
complex;
28-14
ATOG
The data set defaults to single precision, since no character data is present, and no outtyp
command is specified.
ATOG
atog - Can’t read temporary file
Notify Aptech Systems.
28-15
GAUSS User Guide
This is a warning that character data was found in a variable that was specified as
numeric. The conversion will continue.
28-16
ATOG
@ comment @
ATOG
atog - Statement too long
Command file statements must be less than 16384 bytes.
28-17
GAUSS User Guide
28-18
Error Messages 29
The following is a list of error messages intrinsic to the GAUSS programming language. Error
messages generated by library functions are not included here.
Messages
Error
load Input file too large.
getf Input file too large.
29-1
GAUSS User Guide
29-2
Error Messages
Messages
Error
A symbol has been referenced that has not been given a definition.
29-3
GAUSS User Guide
29-4
Error Messages
Messages
Error
G0049 Target matrix not complex
29-5
GAUSS User Guide
29-6
Error Messages
Messages
Error
You are using an argument of the wrong data type (e.g., inputting a matrix when a string
is called for).
G0073 Redefinition of
29-7
GAUSS User Guide
29-8
Error Messages
Messages
Error
A comment was never closed.
29-9
GAUSS User Guide
29-10
Error Messages
G0113 Unclosed (
G0114 Unclosed [
Messages
Error
G0118 Argument must be integer
File handles must be integral.
29-11
GAUSS User Guide
29-12
Error Messages
Messages
Error
G0147 Bad argument to trig function
29-13
GAUSS User Guide
29-14
Error Messages
Messages
Error
G0171 Called symbol is wrong type
You are attempting to call a local procedure using a pointer to something else.
29-15
GAUSS User Guide
29-16
Error Messages
The seed argument to rndns and rndus must be a simple local or global variable
reference. It cannot be an expression or constant. These functions are obsolete, please
use rndlcn and rndlcu
Messages
Error
G0294 Error opening temporary file
One of the temporary files could not be opened. The directory may be full.
29-17
GAUSS User Guide
29-18
Error Messages
Messages
Error
G0378 Cannot make External Variable
29-19
GAUSS User Guide
29-20
Error Messages
Messages
Error
G0409 outtyp already defaulted to 8
Character data has been found in the output data set before an outtyp 2 or outtyp 4
statement. This is a warning only.
29-21
GAUSS User Guide
29-22
Error Messages
Messages
Error
G0454 Failure opening printer for output
29-23
GAUSS User Guide
29-24
Error Messages
Messages
Error
G0510 Duplicate name for structure member
29-25
GAUSS User Guide
29-26
Error Messages
Messages
Error
29-27
Maximizing Performance 30
These hints will help you maximize the performance of your new GAUSS System.
Some temporary files are created during the autoloading process. If you have a tmp_path
configuration variable or a tmp environment string that defines a path on a RAM disk, the
temporary files will be placed on the RAM disk.
For example:
Performance
set tmp=f:\tmp
A disk cache will also help, as well as having your frequently used files in the first path in the
src_path.
30-1
GAUSS User Guide
You can optimize your library .lcg files by putting the correct drive and path on each file name
listed in the library. The lib command will do this for you.
Use the compile command to precompile your large frequently used programs. This will
completely eliminate compile time when the programs are rerun.
30.2 Loops
The use of the built-in matrix operators and functions rather than do loops will ensure that you are
utilizing the potential of GAUSS.
Here is an example:
x = rndn(8000,1);
you could get a count of the elements with an absolute value greater than 1 with a do loop, like
this:
c = 0;
i = 1;
do while i <= rows(x);
if abs(x[i]) > 1;
c = c+1;
endif;
i = i+1;
endo;
print c;
30-2
Maximizing Performance
Computers today can have large amounts of RAM. This doesn’t mean that large data sets should
be read entirely into memory. Many GAUSS procedures and applications are written to allow for
data sets to be read in sections rather than all at once. Even if you have enough RAM to store the
data set completely, you should consider taking advantage of this feature. The speed-ups using this
feature can be significant. For example, ols is called using a data set stored in a matrix versus
stored on the disk in a GAUSS data set. The computer is a 2.8 Megahertz computer with Windows
XP.
y = rndn(250000,1);
x = rndn(250000,100);
xlbl = 0$+"X"+ftocv(seqa(1,1,100),1,0);
lbl = "Y" | xlbl;
call saved(y˜x,"test",lbl);
__output = 0;
t0 = date;
call ols("",y,x);
t1 = date;
t2 = date;
call ols("test","Y",xlbl);
t3 = date;
Performance
print ethsec(t2,t3)/100 " seconds;
print;
print ethsec(t0,t1)/100 " seconds";
25.750000 seconds
9.6720000 seconds
30-3
GAUSS User Guide
This represents more than a 50% speedup by leaving the data on the disk.
maxvec,maxbytes
maxvec is a GAUSS procedure that returns the value of the global variable __maxvec that
determines the amount of data to be read in at a time from a GAUSS data set. This value can be
modified for a particular run by setting __maxvec in your command file to some other value. The
value returned by a call to maxvec can be permanently modified by editing system.dec and
changing the value of __maxvec. The value returned when running GAUSS Light is always 8192.
maxbytes is a GAUSS procedure that returns the value of a scalar global __maxbytes that sets
the amount of available RAM. This value can be modified for a particular run by setting
__maxbytes in your command file to some other value. The value returned by a call to maxbytes
can be permanently modified by editing system.dec and changing the value of __maxbytes.
If you wish to force GAUSS procedures and applications to read a GAUSS data set in its entirety,
set __maxvec and __maxbytes to very large values.
The hard disk used for the swap file should be optimized occasionally with a disk optimizer. Use a
disk maintenance program to ensure that the disk media is in good shape.
There is a line for cache size in the gauss.cfg file. Set it to the size of the CPU data cache for
your computer.
This affects the choice of algorithms used for matrix multiply functions.
This will not change the results you get, but it can radically affect performance for large matrices.
30-4
Fonts
Fonts A
There are four fonts available in the Publication Quality Graphics System:
The following tables show the characters available in each font and their ASCII values. (For
details on selecting fonts for your graph, see S F, Section 25.4.1.
A-1
GAUSS User Guide
A.1 Simplex
A-2
Fonts
Fonts
A.2 Simgrma
A-3
GAUSS User Guide
A.3 Microb
A-4
Fonts
Fonts
A.4 Complex
A-5
Words
Reserved
Reserved Words Appendix B
The following words are used for GAUSS functions. You cannot use these names for variables or
procedures in your programs:
abs AmericanBSCall
acf AmericanBSCall_Greeks
aconcat AmericanBSCall_ImpVol
acos AmericanBSPut
aeye AmericanBSPut_Greeks
amax AmericanBSPut_ImpVol
amean amin
AmericanBinomCall amult
AmericanBinomCall_Greeks and
AmericanBinomCall_ImpVol annualTradingDays
AmericanBinomPut arccos
AmericanBinomPut_Greeks arcsin
AmericanBinomPut_ImpVol arctan
B-1
GAUSS User Guide
arctan2 asin
areshape asum
arrayalloc atan
arrayindex atan2
arrayinit atranspose
arraytomat axmargin
asclabel
balance base10
band begwind
bandchol besselj
bandcholsol bessely
bandltsol box
bandrv boxcox
bandsolpd break
bar
calcbox cdfn
call cdfn2
callexe cdfnc
cdfbeta cdfni
cdfbvn cdftc
cdfbvn2 cdftci
cdfbvn2e cdftnc
cdfchic cdftvn
cdfchii cdir
cdfchinc ceil
cdffc cfft
cdffnc cffti
cdfgam changedir
cdfmvn chdir
B-2
Reserved Words Appendix
checkinterrupt cons
chol continue
choldn contour
cholsol conv
cholup convertsatostr
Words
Reserved
chrs convertstrtosa
cint coreleft
clear corrm
clearg corrms
close corrvc
closeall corrx
cls corrxs
cmsplit cos
cmsplit2 cosh
code counts
color countwts
cols create
colsf crossprd
combinate crout
combinated croutp
comlog csrcol
commandeerm csrlin
commandeersa csrtype
compile cumprodc
complex cumsumc
con curve
cond cvtos
conformed cvtosa
conj
datacreate dataopen
datacreatecomplex datasave
datalist date
dataload datestr
B-3
GAUSS User Guide
datestring dotfge
datestrymd dotfgemt
dayinyr dotfgt
dayOfWeek dotfgtmt
debug dotfle
declare dotflemt
delete dotflt
deletefile dotfltmt
delif dotfne
denseSubmat dotfnemt
design draw
det dsCreate
detl dstat
dfft dstatmt
dffti dstatmtControlCreate
dfree dtdate
diag dtday
diagrv dttime
digamma dttodtv
disable dttostr
dlibrary dttoutc
dllcall dtvnormal
do dtvtodt
dos dtvtoutc
doswincloseall dummy
doswinopen dummybr
dotfeq dummydn
dotfeqmt
ed eigcg2
edit eigch
editm eigch2
eig eigh
eigcg eighv
B-4
Reserved Words Appendix
eigrg erfc
eigrg2 error
eigrs errorlog
eigrs2 etdays
eigv ethsec
Words
Reserved
elapsedTradingDays etstr
else EuropeanBinomCall
elseif EuropeanBinomCall_Greeks
enable EuropeanBinomCall_ImpVol
end EuropeanBinomPut
endfor EuropeanBinomPut_Greeks
endif EuropeanBinomPut_ImpVol
endo EuropeanBSCall
endp EuropeanBSCall_Greeks
endwind EuropeanBSCall_ImpVol
envget EuropeanBSPut
eof EuropeanBSPut_Greeks
eq EuropeanBSPut_ImpVol
eqSolve exctsmpl
eqSolvemt exec
eqSolvemtControlCreate execbg
eqSolvemtOutCreate exp
eqSolveSet expr
eqv external
erf eye
fcheckerr fftmi
fclearerr fftn
feq fge
feqmt fgemt
fflush fgets
fft fgetsa
ffti fgetsat
fftm fgetst
B-5
GAUSS User Guide
fgt fonts
fgtmt fontunload
fileinfo fontunloadall
files fopen
filesa for
fix format
fle formatcv
flemt formatnv
floor fputs
flt fputst
fltmt fseek
fmod fstrerror
fn ftell
fne ftocv
fnemt ftos
font ftostrc
fontload
gamma gdaread
gammaii gdareadbyindex
gausset gdareadsome
gdaappend gdareportvarinfo
gdacreate gdaupdate
gdadstat gdaupdateandpack
gdadstatmat gdawrite
gdagetindex gdawritesome
gdagetname gdtfastcat
gdagetnames ge
gdagetorders getarray
gdagettype getdims
gdagettypes getf
gdagetvarinfo getmatrix
gdaiscplx getmatrix4d
gdapack getname
B-6
Reserved Words Appendix
getnamef gosub
getNextTradingDay goto
getNextWeekDay gradMT
getnr gradMTm
getnrmt gradp
Words
Reserved
getorders graph
getpath graphgpg
getPreviousTradingDay graphinit
getPreviousWeekDay graphprt
getscalar3d graphset
getscalar4d graphsev3
getwind gt
hardcopy hessMTm
hasimag hessMTmw
header hessMTw
headermt hessp
hess hist
hessMT histf
hessMTg histp
hessMTgw hsec
if indnv
imag indsav
indcv int
indexcat intgrat2
indices intgrat3
indices2 inthp
indicesf intHP1
indicesfn intHP2
B-7
GAUSS User Guide
intHP3 intsimp
intHP4 inv
inthpControlCreate invpd
intquad1 invswp
intquad2 iscplx
intquad3 iscplxf
intrleav isinfnanmiss
intrleavsa ismiss
intrsect isSparse
intrsectsa
key keyw
keyav keyword
keymatchmc
lag lapsvds
lag1 lapsvdusv
lagn le
lapeighb let
lapeighi lib
lapeighvb library
lapeighvi license_id
lapgeig line
lapgeigh linsolve
lapgeighv ln
lapgeigv lncdfbvn
lapgschur lncdfbvn2
lapgsvdcst lncdfmvn
lapgsvds lncdfn
lapgsvdst lncdfn2
lapsvdcusv lncdfnc
B-8
Reserved Words Appendix
lnfact loess
lngamma loessmt
lnpdfmvn loessmtControlCreate
lnpdfmvt log
lnpdfn loglog
Words
Reserved
lnpdft logx
load logy
loadarray loopnextindex
loadd lower
loadexe lowmat
loadf lowmat1
loadk lpos
loadm lprint
loadp lpwidth
loads lshow
loadstruct lt
loadwind ltrisol
local lu
locate lusol
machEpsilon mbesselei1
makevars mbesseli
makewind mbesseli0
margin mbesseli1
matalloc meanc
matinit median
matrix mergeby
mattoarray mergebysa
maxbytes mergevar
maxc minc
maxindc minindc
maxvec miss
mbesselei missex
mbesselei0 missrv
B-9
GAUSS User Guide
moment movingaveExpwgt
momentd movingaveWgt
movingave msym
nametype nextn
ndpchk nextnevn
ndpclex nextwind
ndpcntrl not
ne null
new null1
nextindex numCombinations
oldfft open
oldffti openpqg
ols optn
olsmt optnevn
olsmtControlCreate or
olsqr orth
olsqr2 output
olsqrmt outwidth
ones
pacf pi
packr pinv
parse pinvmt
pause plot
pdfn plotsym
B-10
Reserved Words Appendix
polar putarray
polychar putf
polyeval pvCreate
polyint pvgetindex
polymake pvgetparnames
Words
Reserved
polymat pvgetparvector
polymroot pvLength
polymult pvList
polyroot pvnumoffsets
pop pvoffsets
pqgwin pvPack
prcsn pvPacki
previousindex pvPackm
princomp pvPackmi
print pvPacks
printdos pvPacksi
printfm pvPacksm
printfmt pvPacksmi
proc pvputparvector
prodc pvtest
push pvunpack
QNewton qre
QNewtonmt qrep
QNewtonmtControlCreate qrsol
QNewtonmtOutCreate qrtsol
qnewtonset qtyr
QProg qtyre
QProgmt qtyrep
qprogMTInCreate quantile
qqr quantiled
qqre quantilem
qqrep quantilemd
qr qyr
B-11
GAUSS User Guide
qyre qyrep
rank rndkmi
rankindx rndkmn
readr rndKMnb
real rndKMp
recode rndkmu
recserar rndKMvm
recsercp rndLCbeta
recserrc rndLCgam
register_off rndlci
register_on rndlcn
register_reset rndLCnb
register_show rndLCp
renamefile rndlcu
replay rndLCvm
rerun rndmod
reshape rndmult
retp rndn
return rndnb
rev rndns
rfft rndp
rffti rndseed
rfftip rndu
rfftn rndus
rfftnp rndvm
rfftp rotater
rndbeta round
rndcon rows
rndgam rowsf
rndi rref
rndKMbeta run
rndKMgam
B-12
Reserved Words Appendix
satocv sinh
satostrC sleep
save solpd
Words
Reserved
saveall sortc
saved sortcc
savestruct sortd
savewind sorthc
scale sorthcc
scale3d sortind
scalerr sortindc
scalinfnanmiss sortindmc
scalmiss sortmc
schtoc sortr
schur sortrc
screen sparseCols
scroll sparseEye
searchsourcepath sparseFD
seekr sparseFP
selif sparseHConcat
seqa sparseNZE
seqm sparseOnes
setarray sparseRows
setcnvrt sparseScale
setdif sparseSet
setdifsa sparseSolve
setvars sparseSubmat
setvmode sparseTD
setvwrmode sparseTranspose
setwind sparseTrTD
shell sparseTscalar
shiftr sparseVConcat
show spline
showpqg spline1D
sin spline2D
singleindex sqpmt_feasible
B-13
GAUSS User Guide
sqpmt_meritFunct strtofcplx
sqpSolve strtriml
SQPsolveMT strtrimr
sqpSolveMTcontrolCreate strtrunc
sqpSolveMTlagrangeCreate strtruncl
sqpSolveMToutCreate strtruncpad
sqpSolveset strtruncr
sqrt struct
stdc submat
stocv subscat
stof substute
stop subvec
strcombine sumc
strindx sumr
string surface
strlen svd
strput svd1
strrindx svd2
strsect svdcusv
strsplit svds
strsplitpad svdusv
strtodt sysstate
strtodtd system
strtof
tab timedt
tan timestr
tanh timeutc
tempname title
ThreadBegin tkf2eps
ThreadEnd tkf2ps
ThreadJoin tkf2ps_margin
ThreadStat tocart
time todaydt
B-14
Reserved Words Appendix
toeplitz trim
token trimr
topolar trunc
trace type
trap typecv
Words
Reserved
trapchk typef
trigamma
union until
unionsa upmat
uniqindmc upmat1
uniqindx upper
uniqindxsa use
unique utctodt
uniquemc utctodtv
uniquesa utrisol
vals vec
varget vech
vargetl vecr
varmall vfor
varmares vget
varput view
varputl viewxyz
vartype vlist
vartypef vnamecv
vcm volume
vcms vput
vcx vread
vcxs vtypecv
B-15
GAUSS User Guide
wait winpan
waitc winprint
walkindex winprintpqg
while winrefresh
winclear winrefresharea
wincleararea winresize
winclearttylog winsetactive
winclose winsetbackground
wincloseall winsetcolor
winconvertpqg winsetcolorcells
window winsetcolormap
wingetactive winsetcursor
wingetattributes winsetforeground
wingetcolorcells winsetrefresh
wingetcursor winsettextwrap
winmove winwrite
winopenpqg winzoompqg
winopentext writer
winopentty
x_indcv xtics
xlabel xy
xor xyz
xpnd
ylabel ytics
B-16
Reserved Words Appendix
zeros ztics
zlabel
Words
Reserved
B-17
Singularity Tolerance Appendix C
Singularity
The tolerance used to determine whether or not a matrix is singular can be changed. The default
value is 1.0e-14 for both the LU and the Cholesky decompositions. The tolerance for each
decomposition can be changed separately. The following operators are affected by a change in the
tolerance:
Crout LU Decomposition
crout(x)
croutp(x)
inv(x)
det(x)
y/x when neither x nor y is scalar and x is square.
Cholesky Decomposition
chol(x)
invpd(x)
solpd(y,x)
y/x when neither x nor y is scalar and x is not square.
C-1
GAUSS User Guide
The tolerance value may be read or set using the sysstate function, cases 13 and 14.
There is no perfect tolerance for determining singularity. The default is 1.0e-14. You can adjust
this as necessary.
A numerically better method of determining singularity is to use cond to determine the condition
number of the matrix. If the equation
1 / cond(x) + 1 eq 1
is true, then the matrix is usually considered singular to machine precision. (See LAPACK for a
detailed discussion on the relationship between the matrix condition and the number of significant
figures of accuracy to be expected in the result.)
C-2
Index
Index
! , 11-6 / = , 11-11
∗ , 11-5 ./= , 11-12
*∼ , 11-7 = = , 10-40, 11-10
.* , 11-6 .= = , 11-12
> , 11-11
Index
.*. , 11-6
+ , 11-4 . > , 11-12
− , 11-4 >= , 11-11
/ , 11-5 . >= , 11-12
./ , 11-6 < , 11-10
% , 11-5 . < , 11-11
ˆ , 11-6, 11-19 <= , 11-10
.ˆ , 11-6 . <= , 11-12
Index-1
Index
Index-2
Index
Index
bandsolpd, 33-58 cdfCauchyinv, 33-76
bar shading, 25-17 cdfchic, 33-76
bar width, 25-18 cdfchii, 33-78
bar, 33-58 cdfchinc, 33-79
base10, 33-60 cdfexp, 33-80
batch mode, 3-1 cdfexpinv, 33-81
begwind, 33-61 cdffc, 33-81
besselj, 33-61 cdffnc, 33-83
bessely, 33-62 cdfgam, 33-84
beta function, 33-68 cdfGenPareto, 33-86
beta, 33-63 cdfLaplace, 33-87
binary file, loading, 33-371 cdfLaplaceinv, 33-88
binary files, 21-15 cdflogistic, 33-88
bivariate Normal, 33-70 cdflogisticinv, 33-89
blank lines, 10-38 cdfm.src, 33-91, 33-93, 33-94, 33-96,
bookmarks, 5-2 33-97, 33-99
Boolean operators, 11-13 cdfmvn, 33-90
box, 25-18 cdfmvn2e, 33-93
box, 33-64 cdfmvnce, 33-90
Index-3
Index
Index-4
Index
Index
control flow, 10-31 cursor, 33-165, 33-524
control structures, 16-22 curve, 33-168
conv, 33-147 cvtos, 33-169
conversion, character to ASCII value, 33-933
D
conversion, float to ASCII, 33-326, 33-327
conversion, string to floating point, 33-850 data coding, 32-38
convertsatostr, 33-148 data handling, 32-33
convertstrtosa, 33-148 data loop, 23-1
convolution, 33-147 data sets, 21-7, 32-36
coordinates, 25-6 data transformations, 23-1, 33-126, 33-192
correlation matrix, 33-149, 33-150, 33-573, data, writing, 33-953
33-579 datacreate, 33-170
corrm, 33-149 datacreatecomplex, 33-172
corrms, 33-150 datalist, 33-174
corrvc, 33-149 dataload, 33-175
corrx, 33-149 dataloop translator, 3-6
corrxs, 33-150 dataloop, 33-176
cos, 33-150 dataopen, 33-177
cosh, 33-151 datasave, 33-178
Index-5
Index
Index-6
Index
Index
edit windows, 5-1 erf, 33-259
edit, 33-239 erfc, 33-259
editing matrices, 6-1 erfccplx, 33-260
editor, 33-239 erfcplx, 33-260
editor properties, 5-2 error bar, 25-20
editor, alternate, 33-238 error code, 33-261, 33-772
Editor, Matrix, 6-1 error function, 33-259
eig, 33-239 error handling, 32-50
eigenvalues, 32-9, 33-239 error messages, 29-1, 33-262, 33-501
eigenvalues and eigenvectors, 33-243 Error Output Window, 5-7
eigh, 33-241 error trapping, 33-913
eighv, 33-242 error, 33-261
eigv, 33-243 errorlog, 33-262
elapsedTradingDays, 33-244 errorlogat, 33-262
element-by-element conformability, 11-1, escape character, 10-22
14-5 etdays, 33-263
element-by-element operators, 11-1 ethsec, 33-264
else, 33-420 etstr, 27-5, 33-265
elseif, 33-420 EuropeanBinomCall, 33-265
Index-7
Index
Index-8
Index
Index
ftos, 33-327 gdaSave, 33-359
ftostrC, 33-330 gdaUpdate, 33-361
function, 10-37, 33-467, 33-634 gdaUpdateAndPack, 33-363
functions, 32-46 gdaVars, 33-364
fuzzy conditional functions, 32-12 gdaWrite, 33-365
gdaWrite32, 33-366
G
gdaWriteSome, 33-367
gamma function, 33-332 ge, 11-11
gamma, 33-332 .ge, 11-12
gamma, incomplete, 33-84 generalized inverse, 33-459, 33-609, 33-610
gamma, log, 33-509 Generalized Pareto, 33-86, 33-604
gammacplx, 33-332 getarray, 33-369
gammaii, 33-333 getArray, 15-12
GAUSS Data Archives, 21-11, 21-24, 32-34 getdims, 33-370
GAUSS Source Browser, 8-1 getDims, 15-27
Gauss-Legendre quadrature, 33-448 getf, 33-371
gausset, 31-6, 33-334 getmatrix, 33-372
gdaAppend, 33-334 getMatrix, 15-13
gdaCreate, 33-336 getmatrix4D, 33-373
Index-9
Index
Index-10
Index
Index
indexing matrices, 10-40, 11-16 intrleav, 33-453
indexing procedures, 11-17 intrleavsa, 33-454
indexing, array, 14-3 intrsect, 33-455
indexing, structure, 16-5 intrsectsa, 33-456
indices, 33-426 intsimp, 33-456
indices2, 33-427 inv, 33-457
indicesf, 33-428 invar, ATOG command, 28-5
indicesfn, 33-429 inverse cosine, 33-29
indnv, 33-430 inverse sine, 33-30
indsav, 33-432 inverse, generalized, 33-459, 33-609, 33-610
infinity, 10-28 inverse, matrix, 33-457
initialize, 12-4 inverse, sweep, 33-459
initializing arrays, 15-1 invpd, 33-457
inner product, 11-5 invswp, 33-459
input, ATOG command, 28-4 iscplx, 33-460
input, console, 33-137 iscplxf, 33-461
input, keyboard, 33-137 isden, 33-461
installation, 2-1 isinfnanmiss, 33-462
installation, UNIX/Linux, 2-1 ismiss, 33-462
Index-11
Index
J lapgeigv, 33-478
lapgschur, 33-487
Jacobian, 33-394 lapgsvdcst, 33-479
K lapgsvds, 33-482
lapgsvdst, 33-484
keep (dataloop), 33-463 Laplace, 33-87, 33-88, 33-604
key, 33-464 lapsvdcusv, 33-488
keyav, 33-466 lapsvds, 33-490
keyboard input, 33-141 lapsvdusv, 33-491
keyboard, reading, 33-464 le, 11-10
keys, command, 5-17 .le, 11-12
keys, edit, 5-16 least squares, 11-5
keys, function, 5-18 least squares regression, 33-571, 33-575
keys, menu, 5-19 left-hand side, 19-2
keys, movement, 5-15 legend, 25-22
keys, text selection, 5-17 let, 33-492
keystroke macros, 5-2 lib, 33-497
keyw, 33-466 libraries, 19-1, 32-47
keyword, 12-1, 12-7 libraries, active, 33-499
keyword procedure, 33-467 Library Tool, 7-1
keyword, 33-467 library, 33-498
keywords, 32-46 line numbers, 33-501
Kronecker, 11-6 line thickness, 25-16, 25-20, 25-25
line type, 25-25
L
linear algebra, 32-5
label, 10-35, 10-39, 12-1, 33-385, 33-388 linear equation, 33-798
lag (dataloop), 33-468 linear equation solution, 11-5
lag1, 33-469 lines, 25-21, 25-22, 25-24
lagn, 33-469 #linesoff, 33-501
lambda, 33-79 #lineson, 33-501
lapeighb, 33-470 linsolve, 33-502
lapeighi, 33-471 listwise (dataloop), 33-503
lapeigvb, 33-472 listwise deletion, 33-220, 33-222, 33-559,
lapeigvi, 33-474 33-600
lapgeig, 33-475 literal, 10-23, 11-19
lapgeigh, 33-476 ln, 33-503
lapgeighv, 33-477 lncdfbvn, 33-504
Index-12
Index
Index
loadp, 33-514 M
loads, 33-514
loadstruct, 33-522 machEpsilon, 33-535
loadwind, 33-522 machine epsilon, 33-100, 33-891, 33-896
local variable declaration, 12-3 machine requirements, 2-2
local variables, 10-8, 12-3, 33-523 macros, 5-2
local, 12-2, 33-523 magnification, 25-29
locate, 33-524 make (dataloop), 33-536
loess, 33-524 makevars, 33-536
loessmt, 33-525 makewind, 33-538
loessmtControlCreate, 33-526 margin, 33-539
log coordinates, 33-528 matalloc, 33-540
log factorial, 33-509 matinit, 33-541
log gamma, 33-509 matrices, indexing, 10-40
log, 33-527 matrix conformability, 11-1
log, base 10, 33-527 Matrix Editor, 6-1
log, natural, 33-503 matrix files, 21-13
logging commands, 33-134 matrix manipulation, 32-24
logical operators, 11-13 matrix, creation, 33-492
Index-13
Index
Index-14
Index
O _pagesiz, 25-14
pairwise deletion, 11-5, 33-221, 33-222,
obsolete commands, 0-1 33-559
ols, 33-571 panel data, 15-32
olsmt, 33-575 _parrow, 25-14
olsmtControlCreate, 33-581 _parrow3, 25-16
olsqr, 33-581 parse, 33-601
olsqr2, 33-582 pause, 33-602
olsqrmt, 33-583 _paxes, 25-17
ones, 33-584 _paxht, 25-17
open, 33-585 _pbartyp, 25-17
operators, 10-1, 11-4 _pbarwid, 25-18
operators, element-by-element, 11-1 _pbox, 25-18
optimization, 32-17 _pboxlim, 25-19
optn, 33-590 _pcolor, 25-19
optnevn, 33-590 _pcrop, 25-19
or, 11-13, 11-14 _pcross, 25-19
.or, 11-15 _pdate, 25-20
Index
orth, 33-592 pdfCauchy, 33-602
orthogonal complement, 33-569 pdfexp, 33-603
orthonormal, 33-568, 33-592
pdfGenPareto, 33-604
outer product, 11-6
pdfLaplace, 33-604
output, 21-4
pdflogistic, 33-605
output functions, 32-56
pdfn, 33-606
output, 33-593
pdfRayleigh, 33-607
output, ATOG command, 28-10
pdfWeibull, 33-607
outtyp (dataloop), 33-596
_perrbar, 25-20
outtyp, ATOG command, 28-11
_pframe, 25-20
outvar, ATOG command, 28-11
_pgrid, 25-21
outwidth, 33-596
pi, 33-608
P pinv, 33-609
pinvmt, 33-610
pacf, 33-597 pixel coordinates, 25-6
packed ASCII, 28-1, 28-8 _plctrl, 25-21
packedToSp, 33-598 _plegctl, 25-22
packr, 33-600 _plegstr, 25-22
_pageshf, 25-14 _plev, 25-22
Index-15
Index
Index-16
Index
Index
pvPacks, 33-650 quadrature, 33-448
pvPacksi, 33-651 quantile, 33-694
pvPacksm, 33-652 quantiled, 33-695
pvPacksmi, 33-654 qyr, 33-697
pvPutParVector, 33-656 qyre, 33-698
pvTest, 33-658 qyrep, 33-700
pvUnpack, 33-658
_pxpmax, 25-29 R
_pxsci, 25-29
radii, 25-22
_pypmax, 25-29
random numbers, 32-11
_pysci, 25-29
rank of a matrix, 33-702
_pzclr, 25-29
rank, 33-702
_pzoom, 25-29
rankindx, 33-703
_pzpmax, 25-29
Rayleigh, 33-104, 33-105, 33-607
_pzsci, 25-29
readr, 33-704
Q real, 33-705
recode (dataloop), 33-708
QNewton, 33-659 recode, 33-706
QNewtonmt, 33-662 recserar, 33-709
Index-17
Index
Index-18
Index
Index
sequence functions, 32-21 sorting, 32-40
series functions, 32-21 sortmc, 33-804
set difference function, 33-785 sortr, sortrc, 33-805
setArray, 15-16 Source Browser, 8-1
setarray, 33-784 spaces, 11-16
setdif, 33-785 spaces, extraneous, 10-38, 11-16, 11-17
setdifsa, 33-786 sparse matrices, 32-28
setvars, 33-787 spBiconjGradSol, 33-806
setvwrmode, 33-788 spChol, 33-808
setwind, 33-788 spConjGradSol, 33-809
shell, 33-789 spCreate, 33-811
shiftr, 33-790 spDenseSubmat, 33-812
show, 33-791 spDiagRvMat, 33-813
Simpson’s method, 33-456 spEigv, 33-815
sin, 33-794 spEye, 33-817
sine, inverse, 33-30 spGetNZE, 33-818
singleindex, 33-795 spLDL, 33-820
singular value decomposition, 33-876, spline, 33-819
33-878, 33-880 spLU, 33-821
Index-19
Index
Index-20
Index
Index
svd, 33-875 timedt, 33-903
svd1, 33-876 timestr, 33-904
svd2, 33-877 timeutc, 33-905
svdcusv, 33-878 timing functions, 33-420
svds, 33-879 title, 33-905
svdusv, 33-880 tkf2eps, 33-906
sweep inverse, 33-459 tkf2ps, 33-907
symbol names, 10-39 tocart, 33-907
symbol table, 33-791 todaydt, 33-908
symbol table type, 33-918 Toeplitz matrix, 33-908
symbols, allocate maximum number, 33-565 toeplitz, 33-908
syntax, 10-37 token, 33-909
sysstate, 33-881 toolbar, graphics editor, 26-2
system, 33-896 toolbars, 4-10
T topolar, 33-910
trace program execution, 33-911
t distribution, Student’s, 33-105 trace, 33-911
tab, 33-897 translation phase, 23-3
table, 11-6 transpose, 11-8
Index-21
Index
Index-22
Index
waitc, 33-951 Z
walkindex, 33-951
watch variables, 6-3 zeros, 33-970
watch window, 5-10 zeta, 33-971
Weibull, 33-110, 33-111, 33-607 zlabel, 33-971
weighted count, 33-154 zooming graphs, 25-29
while, 33-206 ztics, 33-972
window, 21-4
window, 33-953
window, clear, 33-126
workbox, 33-947, 33-949
workspace, 33-190, 33-792
writer, 33-953
X
xlabel, 33-955
xlsGetSheetCount, 33-956
Index
xlsGetSheetSize, 33-956
xlsGetSheetTypes, 33-957
xlsMakeRange, 33-958
xlsReadM, 33-959
xlsReadSA, 33-960
xlsWrite, 33-962
xlsWriteM, 33-963
xlsWriteSA, 33-964
xor, 11-14
.xor, 11-15
xpnd, 33-966
xtics, 33-967
xy, 33-968
xyz, 33-968
Y
ylabel, 33-969
ytics, 33-969
Index-23