all functions - m
mac_primitives
mac_primitives, file sets FILE primitive data types to be native to MacIntosh, 8 byte double. Interpreted function, defined at i0/std.i line 2913
macl_primitives
macl_primitives, file sets FILE primitive data types to be native to MacIntosh, long double. Interpreted function, defined at i0/std.i line 2921
make
make or make, path or make, path, template=1 or make, path, template="Makefile_name" creates or updates the Makefile for a compiled package (plugin). If PATH is given, it is the directory containing the source code for the package; by default, PATH is ".", the current working directory. With the template= keyword, a template Makefile is written, which you can use as a starting point, filling in all the necessary information by hand. In this case, make makes no further checks on the contents of the directory. In the first form, the name of the Makefile template is "Makefile", but if the template= keyword is a string, that name is used instead. The template Makefile will not function at all until you edit it. Alternatively, without the template= keyword, make looks at the contents of the directory, and fills in as much of the Makefile as it can using what it sees. There are two cases: 1. A Makefile (or makefile) already exists. Provided the file contains a Y_MAKEDIR= line, make merely updates Y_MAKEDIR, Y_EXE, and Y_EXE_PKGS to reflect this running yorick, and attempts no other modifiactions of the existing Makefile. 2. No such Makefile exists. Then make writes a template Makefile, and reads all the .i interpreted source files, .c and .h C source files, and .f or .F or .m Fortran source files in the directory, filling in as much of the Makefile as possible. In simple cases, the resulting Makefile will be perfectly adequate. (If there is any Fortran source, however, it will never be adequate.) In case 2, the following steps are taken: A. Any .i files containing lines of the form plug_in, "pkname" or if (!is_void(plug_in)) plug_in, "pkgname" are identified as the "package include file(s)" which define the interpreted API for the package. The make function takes the name of the package from the "pkgname" in these lines. At least one .i file must contain at least one plug_in line; at most one "pkgname" must occur in all the plug_in lines in all the .i files in the directory. This information is used to fill in the PKG_NAME macro and the PKG_I macro in the Makefile (the latter becomes a list of every .i file containing a plug_in statement). B. All .c, .f, .F, or .m files are assumed to be required compiled source, and the OBJS macro in the Makefile is set to a list of all corresponding .o object files. C. If there are any .h C header files in the directory, every .c C source file is scanned for an include directive which reads the header. An approriate dependency line for each .o file which is discovered to depend on local .h files in this manner. These steps will be sufficient to produce a finished Makefile in many simple cases; an example is given in the extend/ directory of the yorick distribution (a compiled complementary error function). However, it is impossible for yorick to automatically discover complicated Makefile rules your package may require. The most common example is dependency on a third party library, which requires a -l loader option, and very possibly a -L option as well. You need to add these flags to PKG_DEPLIBS in the Makefile by hand. If you want such a package to be buildable by non-experts, you will need to provide a configure script in order to fill in the required information across a wide variety of platforms. This can be exceptionally challenging, which is why yorick itself has no dependencies beyond what is guaranteed to be present on every system where it can run (libc, libm, and, on Unix systems, libX11). Once the Makefile has been built it becomes part of the source code of your package, and should not be removed. (You may want it to have an include directive to get the results of any configure script which must be run, however.) Thereafter, use yorick -batch make.i (or simply make in a running yorick) to set Y_MAKEDIR, Y_EXE, and Y_EXE_PKGS appropriately for the particular platform where you are building (case 1 above). Read Y_HOME/Makepkg for a description of the various make targets available. In a nutshell: make debug builds a debugging version (never a plugin) make builds a release version (a plugin if possible) make install installs the release version (may require root) make clean removes all results, leaving only original source Interpreted function, defined at i/make.i line 6SEE ALSO: plug_in, autoload, cd
make_despace
make_despace Interpreted function, defined at i/make.i line 339
make_sphere
make_sphere(radius, [imax,jmax,kmax], [phi1, phi2], [theta1, theta2]) return a mesh (see hex_mesh) representing the given section of the sphere of given RADIUS. IMAX, JMAX, and KMAX are the number of nodes (cells+1) in the radial, longitude (phi), and colatitude (theta) directions, respectively. Note that for a right handed coordinate system, phi1SEE ALSO: hex_meshtheta2. Interpreted function, defined at i0/hex.i line 539
make_subst
make_subst Interpreted function, defined at i/make.i line 289
matchbc
matchbc Interpreted function, defined at i/demo1.i line 110
max
max(x) or max(x, y, z, ...) returns the scalar maximum value of its array argument, or, if more than one argument is supplied, returns an array of the maximum value for each array element among the several arguments. In the multi-argument case, the arguments must be conformable. Builtin function, documented at i0/std.i line 811SEE ALSO: min, sum, avg
md5
state = md5(data) state = md5(data, state) digest = md5(state) digest = md5(state, hex=1) compute MD5 digest of data, an array of 16 char with the hex=1 keyword, returned digest is a string (in hex) Interpreted function, defined at i/md5.i line 36SEE ALSO: md5sum
md5sum
digest = md5sum(filename) digest = md5sum(filename, hex=1) compute MD5 digest of a file, an array of 16 char with the hex=1 keyword, returned digest is a string (in hex) Interpreted function, defined at i/md5.i line 10SEE ALSO: md5
median
median(x) or median(x, which) returns the median of the array X. The search for the median takes place along the dimension of X specified by WHICH. WHICH defaults to 1, meaning the first index of X. The median function returns an array with one fewer dimension than its argument X (the WHICH dimension of X is missing in the result), in exact analogy with rank reducing index range functions. If dimsof(X)(WHICH) is odd, the result will have the same data type as X; if even, the result will be a float or a double, since the median is defined as the arithmetic mean between the two central values in that case. Interpreted function, defined at i0/std.i line 1231SEE ALSO: sort
merge
merge(true_expr, false_expr, condition) returns the values TRUE_EXPR or FALSE_EXPR where CONDITION is non-zero or zero, respectively. The result has the data type of TRUE_EXPR or FALSE_EXPR, promoted to the higher arithmetic type if necessary. The result has the dimensions of CONDITION. The number of elements in TRUE_EXPR must match the number of non-zero elements of CONDITION, and the number of elements in FALSE_EXPR must match the number of zero elements of CONDITION. (TRUE_EXPR or FALSE_EXPR should be nil if there are no such elements of CONDITION. Normally, TRUE_EXPR and FALSE_EXPR should be 1-D arrays if they are not nil.) This function is intended for vectorizing a function whose domain is divided into two or more parts, as in: func f(x) { big= (x>=threshhold); wb= where(big); ws= where(!big); if (is_array(wb)) { xx= x(wb); fb=SEE ALSO: mergef, merge2, where} if (is_array(ws)) { xx= x(ws); fs= } return merge(fb, fs, big); } Builtin function, documented at i0/std.i line 894
merge2
merge2(true_expr, false_expr, condition) returns the values TRUE_EXPR or FALSE_EXPR where CONDITION is non-zero or zero, respectively. The result has the data type of TRUE_EXPR or FALSE_EXPR, promoted to the higher arithmetic type if necessary. Unlike the merge function, TRUE_EXPR and FALSE_EXPR must be conformable with each other, and with the CONDITION. Interpreted function, defined at i0/std.i line 925SEE ALSO: merge, where, mergef
merge_n
merge_n(val1,mask1, val2,mask2, ...) return array with shape of MASKi (which must all have same shape) and values VALi where MASKi is true. Unspecified values will be zero; the data type of the result is the data type of the first non-nil VALi. Each VALi must be a 1D array of length sum(MASKi!=0). Interpreted function, defined at i/dawson.i line 161SEE ALSO: any_in, merge
mergef
y = mergef(x, f1, cond1, f2, cond2, ... felse) Evaluate F1(X(where(COND1))), F2(X(where(COND2))), and so on, until FELSE(X(where(!(COND1 | COND2 | ...)))) and merge all the results back into an array Y with the same dimensions as X. Each of the CONDi must have the same dimensions as X, and they must be mutally exclusive. During the evaluation of Fi, note that all of the local variables of the caller of mergef are available. The Fi are invoked as Fi(X(mergel)) and the variable mergel = where(CONDi) is available to the Fi, in case they need to extract any additional parameters. If noneof(CONDi) then Fi will not be called at all, otherwise, the Fi are invoked in order. The return value of Fi must have the same shape as its argument (which will be a 1D array or scalar). Use mergeg to construct secondary results the same shape as X and Y. Interpreted function, defined at i0/std.i line 947SEE ALSO: mergeg, merge
mergeg
z = mergeg(z, value) or z = mergeg(z) If secondary results are to be returned from a mergef, besides its return value, the Fi may construct them using the second form of mergef: z = mergeg(z, value) where z is a variable in the original caller of mergef, and value is its value where(CONDi). Note that the variable name of the first parameter must be the same as the variable name it is assigned to in this construction -- that variable is being used to hold the state of z as it is built. After the outer mergef returns, the caller needs to invoke z = mergeg(z) one final time to complete each secondary return value. z = []; y = mergef(x, f1, cond, f2); z = mergeg(z); ... func f1(x) {SEE ALSO: mergef, mergez = mergeg(z, exprz(x)); return expry(x); } func f2(x) { z = mergeg(z, exprz(x)); return expry(x); } Interpreted function, defined at i0/std.i line 993
mesh3
mesh3(x,y,z) or mesh3(x,y,z, f1,f2,...) or mesh3(xyz, f1,f2,...) or mesh3(nxnynz, dxdydz, x0y0z0, f1,f2,...) make mesh3 argument for slice3, xyz3, getv3, etc., functions. X, Y, and Z are each 3D coordinate arrays. The optional F1, F2, etc. are 3D arrays of function values (e.g. density, temperature) which have one less value along each dimension than the coordinate arrays. The "index" of each zone in the returned mesh3 is the index in these cell-centered Fi arrays, so every index from one through the total number of cells indicates one real cell. The Fi arrays can also have the same dimensions as X, Y, or Z in order to represent point-centered quantities. If X has four dimensions and the length of the first is 3, then it is interpreted as XYZ (which is the quantity actually stored in the returned cell list). If X is a vector of 3 integers, it is interpreted as [nx,ny,nz] of a uniform 3D mesh, and the second and third arguments are [dx,dy,dz] and [x0,y0,z0] respectively. (DXDYDZ represent the size of the entire mesh, not the size of one cell, and NXNYNZ are the number of cells, not the number of points.) Interpreted function, defined at i/slice3.i line 66SEE ALSO: slice3, xyz3, getv3, getc3
mesh_loc
mesh_loc(y0, x0) or mesh_loc(y0, x0, y, x) or mesh_loc(y0, x0, y, x, ireg) returns the zone index (=i+imax*(j-1)) of the zone of the mesh (X,Y) (with optional region number array IREG) containing the point (X0,Y0). If (X0,Y0) lies outside the mesh, returns 0. Thus, eg- ireg(mesh_loc(x0, y0, y, x, ireg)) is the region number of the region containing (x0,y0). If no mesh specified, uses default. X0 and Y0 may be arrays as long as they are conformable. For mesh_loc wrappers to duplicate the functionality of the digitize and interp functions in 2D, see the library file digit2.i. After #include "digit2.i", type: help,digit2 Builtin function, documented at i0/graph.i line 1352SEE ALSO: plmesh, moush, mouse
min
min(x) or min(x, y, z, ...) returns the scalar minimum value of its array argument, or, if more than one argument is supplied, returns an array of the minimum value for each array element among the several arguments. In the multi-argument case, the arguments must be conformable. Builtin function, documented at i0/std.i line 801SEE ALSO: max, sum, avg
mkdir
mkdir, directory_name rmdir, directory_name Create DIRECTORY_NAME with mkdir, or remove it with rmdir. The rmdir function only works if the directory is empty. Builtin function, documented at i0/std.i line 2478SEE ALSO: cd, lsdir, get_cwd, get_home
mkdoc
mkdoc, filename or mkdoc, filename, outname, lpp alphabetizes and indexes the DOCUMENT comments in FILENAME, and formats into "dictionary-like" pages for printing. If OUTNAME is not given or nil, the output file will be FILENAME with ".doc" replacing the ".i". If LPP is not given, it defaults to 58 -- the maximum number of available lines per page of output. (Use 55 to be able to print with "lpr -p" style page headings.) FILENAME can be an array of strings to combine several include files into a single document. Interpreted function, defined at i/mkdoc.i line 10SEE ALSO: help
mkdoc_cross
mkdoc_cross Interpreted function, defined at i/mkdoc.i line 206
mkdoc_page
mkdoc_page Interpreted function, defined at i/mkdoc.i line 258
mkdoc_scan
mkdoc_scan Interpreted function, defined at i/mkdoc.i line 144
mkdoc_title
mkdoc_title Interpreted function, defined at i/mkdoc.i line 218
mnbrent
fmin= mnbrent(f, x0, x1, x2) or fmin= mnbrent(f, x0, x1, x2, xmin) or fmin= mnbrent(f, x0, x1, x2, xmin, xerr) returns the minimum of the function F (of a single argument x), given three points X0, X1, and X2 such that F(X1) is less than either F(X0) or F(X2), and X1 is between X0 and X2. If the XMIN argument is provided, it is set to the x value which produced FMIN. If XERR is supplied, the search stops when a fractional error of XERR in x is reached; note that XERR smaller than the square root of the machine precision (or omitted) will cause convergence to machine precision in FMIN. The algorithm is Brent's method - a combination of inverse parabolic interpolation and golden section search - as adapted from Numerical Recipes Ch. 10 (Press, et. al.). Interpreted function, defined at i/roots.i line 235SEE ALSO: mxbrent, nraphson, f_inverse
mod_midpt
mod_midpt Interpreted function, defined at i/rkutta.i line 320
moon
xyz = moon(time) return position XYZ of the moon relative to center of earth at TIME; the XYZ has leading dimension 3; x is along the vernal equinox, z is ecliptic north. The corrections to the lunar orbit are from Schlyter (see sch_moon). Claimed accurate to 2 arc minutes over some reasonable time. TIME is in days since 0/Jan/00 (that is, 0000 UT 31/Dec/99). This is 1.5 days earlier than the J2000 epoch. Interpreted function, defined at i/kepler.i line 261SEE ALSO: solar_system, sch_moon, kepler
mouse
result= mouse(system, style, prompt) displays a PROMPT, then waits for a mouse button to be pressed, then released. Returns array of eleven doubles: result= [x_pressed, y_pressed, x_released, y_released, xndc_pressed, yndc_pressed, xndc_released, yndc_released, system, button, modifiers] If SYSTEM>=0, the first four coordinate values will be relative to that coordinate system. For SYSTEM<0, the first four coordinate values will be relative to the coordinate system under the mouse when the button was pressed. The second four coordinates are always normalized device coordinates, which start at (0,0) in the lower left corner of the 8.5x11 sheet of paper the picture will be printed on, with 0.0013 NDC unit being 1/72.27 inch (1.0 point). Look in the style sheet for the location of the viewport in NDC coordinates (see the style keyword). If STYLE is 0, there will be no visual cues that the mouse command has been called; this is intended for a simple click. If STYLE is 1, a rubber band box will be drawn; if STYLE is 2, a rubber band line will be drawn. These disappear when the button is released. Clicking a second button before releasing the first cancels the mouse function, which will then return nil. Ordinary text input also cancels the mouse function, which again returns nil. The left button reverses forground for background (by XOR) in order to draw the rubber band (if any). The middle and right buttons use other masks, in case the rubber band is not visible with the left button. long(result(9)) is the coordinate system in which the first four coordinates are to be interpreted. long(result(10)) is the button which was pressed, 1 for left, 2 for middle, and 3 for right (4 and 5 are also possible). long(result(11)) is a mask representing the modifier keys which were pressed during the operation: 1 for shift, 2 for shift lock, 4 for control, 8 for mod1 (alt or meta), 16 for mod2, 32 for mod3, 64 for mod4, and 128 for mod5. Holding the shift key and pressing the left mouse button is equivalent to pressing the middle mouse button. Similarly, pressing meta-left is equivalent to the right button. This permits access to the middle and right button functions on machines (e.g.- most laptops) with two button or one button mice. The long(result(10)) value returned by mouse() reflects this convention, returning 2 or 3 for those cases, even though it is button 1 that is actually being pressed. Therefore, there is no way to distinguish shift-left from shift-middle, because the long(result(11)) mask indicates tht the shift button is pressed in either case. (And on a machine without a middle button, there would be no way to emulate shift-middle anyway.) Builtin function, documented at i0/graph.i line 1370SEE ALSO: moush
moush
moush() or moush(y, x, ireg) returns the 1-origin zone index for the point clicked in for the default mesh, or for the mesh (X,Y) (region array IREG). Interpreted function, defined at i0/graph.i line 1429
mov3
mov3, xa,ya,za move the current 3D plot by XA along viewer's x-axis, YA along viewer's y-axis, and ZA along viewer's z-axis. Interpreted function, defined at i/pl3d.i line 85SEE ALSO: rot3, orient3, setz3, undo3, save3, restore3, light3
movie
movie, draw_frame or movie, draw_frame, time_limit or movie, draw_frame, time_limit, min_interframe runs a movie based on the given DRAW_FRAME function. The movie stops after a total elapsed time of TIME_LIMIT seconds, which defaults to 60 (one minute), or when the DRAW_FRAME function returns zero. func draw_frame(i) { // Input argument i is the frame number. // draw_frame should return non-zero if there are more // frames in this movie. A zero return will stop the // movie. // draw_frame must NOT include any fma command if the // making_movie variable is set (movie sets this variable // before calling draw_frame) } If MIN_INTERFRAME is specified, a pauses will be added as necessary to slow down the movie. MIN_INTERFRAME is a time in seconds (default 0). The keyword bracket_time= (again a time in seconds) can be used to adjust the duration of the pauses after the first and last frames. It may also be a two element array [beg, end]. If the pause at the end is greater than five seconds, you will be prompted to explain that hittingSEE ALSO: movie_statswill abort the final pause. If every frame of your movie has the same limits, use the limits command to fix the limits before you call movie. BUG: If you hit to start a movie early, it will not pause at the end of the movie at all. You probably should not use long initial pauses. Interpreted function, defined at i/movie.i line 10
movie_stats
movie_stats or movie_stats, timing prints statistics from the last movie command, or from the command which produced TIMING. TIMING is the contents of the movie_timing external variable after the movie command completes. Interpreted function, defined at i/movie.i line 116SEE ALSO: movie
mpeg_close
mpeg_close, mpeg Close the mpeg file corresponding to the MPEG encoder. Actually, this merely destroys the reference to the encoder; the file will remain open until the final reference is destroyed. Interpreted function, defined at i0/mpeg.i line 45SEE ALSO: mpeg_create, mpeg_write, mpeg_movie
mpeg_create
mpeg = mpeg_create(filename) or mpeg = mpeg_create(filename, params) Create an mpeg-1 movie file FILENAME. Write frames with mpeg_write, close with mpeg_close. The return value is an mpeg encoder object. If given, PARAMS is [bit_rate, frame_rate, gop_size, max_b_frames] which default to [ 400000, 25, 10, 1 ] The rates are per second, the gop_size is the number of frames before an I-frame is emitted, and max_b_frames is the largest number of consecutive B-frames. (The third kind of frame is the P-frame; generally, the encoder emits B-frames until it is forced to emit a P-frame by max_b_frames, or an I-frame by gop_size. The smaller these numbers, the higher quality the movie, but the lower the compression.) Any of the four PARAMS values may be zero to get the default value, except for max_b_frames, which should be <0 to get the default value. Builtin function, documented at i0/mpeg.i line 9SEE ALSO: mpeg_write, mpeg_close, mpeg_movie
mpeg_movie
mpeg_movie, filename, draw_frame or mpeg_movie, filename, draw_frame, time_limit or mpeg_movie, filename, draw_frame, time_limit, min_interframe An extension of the movie function (#include "movie.i") that generates an mpeg file FILENAME. The other arguments are the same as the movie function (which see). The draw_frame function is: func draw_frame(i) { // Input argument i is the frame number. // draw_frame should return non-zero if there are more // frames in this movie. A zero return will stop the // movie. // draw_frame must NOT include any fma command if the // making_movie variable is set (movie sets this variable // before calling draw_frame) } Interpreted function, defined at i0/mpeg.i line 58SEE ALSO: movie, mpeg_create, mpeg_write, mpeg_close
mpeg_write
mpeg_write, mpeg, rgb Write a frame RGB into the mpeg file corresponding to the MPEG encoder returned by mpeg_create. RGB is a 3-by-width-by-height array of char. Every frame must have the same width and height. To finish the movie and close the file, call mpeg_close. Note that you may have trouble rendering the resulting mpeg file if the image width and height are note multiples of 8. Builtin function, documented at i0/mpeg.i line 31SEE ALSO: mpeg_create, mpeg_close, mpeg_movie
msort
msort(x1, x2, x3, ...) returns an index list which sorts the array X1 into increasing order. Where X1 values are equal, the list will sort X2 into increasing order. Where both X1 and X2 are equal, X3 will be in increasing order, and so on. Finally, where all of the keys are equal, the returned list will leave the order unchanged from the input keys. The Xi may be numbers or strings (e.g.- X1 could be an integer while X2 was a string, and X3 was a real). The Xi must all be conformable, and each dimension of X1 must be as large as the corresponding dimension of any otehr Xi. Hence, msort(x) will return the same list as sort(x), except where the values of x are equal, in which case msort leaves the order unchanged, while sort non-deterministically permutes equal elements. This feature may cost a factor of two in speed, so don't use it unless you really need it. In general, msort will call sort up to twice per input argument. Interpreted function, defined at i/msort.i line 10SEE ALSO: sort, msort_rank
msort_rank
msort_rank(x) msort_rank(x, list) returns a list of longs the same size and shape as X, whose values are the "rank" of the corresponding element of X among all the elements of X -- the smallest element has rank 0 and the largest has the largest rank, which is equal to one less than the number of distinct values in the array X. If LIST is present, it is set to the order list returned by sort(x(*)). Interpreted function, defined at i/msort.i line 55SEE ALSO: msort, sort
multi_bav
multi_bav(gav) returns bin boundaries for the bin centers gav. The bin boundaries are taken at the geometric means between consecutive gav(i), with the endpoints extended slightly beyond the endpoints of gav. Interpreted function, defined at i/multi.i line 1141
multi_bins
multi_bins(mf) The MF parameter is an array of MultiFiles, each created by multif. Automatically generates the bin structure which will be used by multi_streak (if the GB keyword is not specified). Interpreted function, defined at i/multi.i line 432
multi_line
gb= multi_line(nbins, hnu0, dhnu, dhnu_min) returns 2*NBINS+1 bin boundary energies for 2*NBINS bins cenetered around a spectral line at HNU0 of width DHNU. The result begins at HNU0-DHNU and ends at HNU0+DHNU. The finest two bins (nearest HNU0) has width DHNU_MIN, and the remaining bins have equal ratio widths as you move away from HNU0. Interpreted function, defined at i/multi.i line 1229
multi_no_dups
xnd= multi_no_dups(x) returns its input vector X with any duplicate values removed. X must be non-decreasing and of length at least two. Interpreted function, defined at i/multi.i line 1253
multi_streak
result= multi_streak(mf, rays, slimits, gb=common_bins) like the streak function, but allows opacity to be built up from "slave files", in addition to the "master file" MF(1). The MF parameter is an array of MultiFiles, each created by multif. The master file MF(1) contains the mesh, and the master list of dump times. Only dump times which are present in this master list, and in every slave file, will be processed. The master file MF(1) need not contain any opacity or emissivity data at all; each of the slave files MF(2:0) must contain data for at least one zone. The emissivities and opacities from each file are interpolated onto a common group structure. This common group structure can be provided via the GB keyword to multi_streak. If it is not provided, GB is computed by examining the group boundary (or center) arrays from the master and every slave file, and building a group structure which is at least as fine as every component group structure, at every point in the spectrum. Example: File family "prob_p00" contains the mesh and opacities and emissivities for all zones. Family "pp_h00" contains post processed opacities and emissivities on a much finer spectral mesh, but only for zones in regions 1 and 2 of the original problem. You want to transport the emission from the inner regions 1 and 2 through the overlying material: restore, openb("prob_p00"), ireg; master= multif("prob_p00", zoneuse=where(ireg>2)); slave= multif("pp_h00", zonelist=where(ireg==1|ireg==2)); rays= ... slimits= ... drat_start= ... drat_stop= ... result= multi_streak([master,slave], rays, slimits); multic, master; multic, slave; Interpreted function, defined at i/multi.i line 251SEE ALSO: multio, multic, multif, MultiFile, multi_opac, multi_emiss,
multi_srcf, multi_gb, multi_gav, multi_zonelist, multi_times,
multi_bins
multi_streak_save
multi_streak_save, outname, mf, rays, slimits, gb=common_bins or multi_streak_save, outfile, mf, rays, slimits, gb=common_bins like the streak function, but allows opacity to be built up from "slave files", in addition to the "master file" MF(1) and saves the streak in a PDB history file. The MF parameter is an array of MultiFiles, each created by multif. The master file MF(1) contains the mesh, and the master list of dump times. Only dump times which are present in this master list, and in every slave file, will be processed. The master file MF(1) need not contain any opacity or emissivity data at all; each of the slave files MF(2:0) must contain data for at least one zone. If the first argument is OUTFILE, a file variable instead of a file name, then that file is used for output. You can create OUTFILE and add static variables to it with save (but do NOT call add_record) which streak_save otherwise wouldn't know about. The output file has history records at the same times as the input file. Each record contains "time" (a double scalar), and the two arrays "transp", the transparency (between 0 and 1), and "selfem", the self emission (which has the same units as ekap in the file F). The dimensions of transp and selfem are ngroup-by-2-by-nrays (where nrays represents zero or more dimensions, copied from the RAYS input array). The RAYS and SLIMITS inputs are placed into the output file as non-record variables, and any variables in the drat_static option are copied form F to the output file. The gb and gav variables are copied from F into the output file as well. If the drat_glist option is present, that is stored in the output file also. The emissivities and opacities from each file are interpolated onto a common group structure. This common group structure can be provided via the GB keyword to multi_streak. If it is not provided, GB is computed by examining the group boundary (or center) arrays from the master and every slave file, and building a group structure which is at least as fine as every component group structure, at every point in the spectrum. Example: File family "prob_p00" contains the mesh and opacities and emissivities for all zones. Family "pp_h00" contains post processed opacities and emissivities on a much finer spectral mesh, but only for zones in regions 1 and 2 of the original problem. File "prob_strk" contains the streak history. You want to transport the emission from the inner regions 1 and 2 through the overlying material: restore, openb("prob_p00"), ireg; master= multif("prob_p00", zoneuse=where(ireg>2)); slave= multif("pp_h00", zonelist=where(ireg==1|ireg==2)); fout= openb("prob_strk"); save, fout, kmax, lmax; rays= ... slimits= ... drat_start= ... drat_stop= ... result= multi_streak_save(fout, [master,slave], rays, slimits); multic, master; multic, slave; Interpreted function, defined at i/multi.i line 312SEE ALSO: multio, multic, multif, MultiFile, multi_streak, multi_opac,
multi_emiss, multi_srcf, multi_gb, multi_gav, multi_zonelist,
multi_times, multi_bins
multi_times
times= multi_times(mf) returns the list of times which will be used by multi_streak. This is the subset of streak_times(mf(1)) which occur in all of the slave files. The drat_start and drat_stop times work as usual. Interpreted function, defined at i/multi.i line 471
multic
multic, mf or multic, [mf1, mf2, mf3, ...] closes a MultiFile created with multif. Presented with an array of multifiles, closes them all. Interpreted function, defined at i/multi.i line 217SEE ALSO: multio, multif
multif
multif(mf) returns an ordinary file pointer for the MultiFile MF. Do not use close to close this pointer; just set it to [] when you are done. Use multic to properly close the MF. Interpreted function, defined at i/multi.i line 240SEE ALSO: multio, multic
multio
mf= multio(filename) or mf= multio(file) opens file FILENAME for use with the multi_streak function. The file MUST be subsequently closed using multic, since this function produces a hidden reference to the file. The function multif can be used to return an ordinary file pointer, given the returned MF structure. If the argument is already a stream FILE, that file will be used. The call still produces a hidden copy of FILE, so you may set your copy of the FILE variable to [], but do not close the file. The following keywords can be used to allow for variations in the variable names or units, and to specify the correspondence between the zones in this file, and the zones in the master file: zonelist=index_list -or- zonelist=zonelist_name is an index list into the (rt,zt) mesh arrays of the master file. If ireg is the region number array (having the same dimensions as rt or zt, and with its first row and column all 0), and if FILENAME contains opacity data only for zones with region numbers 1 and 2, you could open the file using: mf= multio(filename, zonelist=where(ireg==1 | ireg==2)) The zonelist should be nil only if the spatial dimensions of the opacity and emissivity in this file exactly match those of rt or zt in the master file. If zonelist is a string, it replaces the default name for the zonelist variable stored in the file (see multi_zonelist). zoneuse=index_list The zonelist specifies how the zones in this file correspond with those in the master file. The zoneuse list allows you to specify that only some of the zones actually present in the opacity and emissivity arrays of this file are to contribute to the total. This might be necessary to avoid double counting in a region covered by more than one file. Hence zoneuse is a list of indices into the spatial dimension(s) of the opacity and emissivity arrays in this file. If nil, all zones in this file will contribute. If present, and if zonelist is supplied as an array (rather than out of the file), zonelist should have the same length as zoneuse. As a special case, if zoneuse is a scalar 0, no opacity or emissivity will come from this file; this makes sense only if this is the master file. opac=oname, emiss=ename, srcf=sname specify non-default names for the opacity, emissivity, and osource function arrays. The defaults are given by the global variables mutli_opac, multi_emiss, and multi_srcf (see help). If the emissivity array is present in the file, it is preferred to the source function array, which will then be ignored. oscale=opacity_unit, escale=emissivity_unit are optional conversion factors to bring the units of the opac and emiss (or srcf) arrays into agreement among the various files which are to be used in a single run. The default value is 1.0 (i.e.- all files are expected to have the same units). gb=gbname, gav=gavname, gexist=gexistname -or- gexist=group_existence_map specify non-default names for the group boundary, group energy, and group existence arrays. The defaults are given by the global variables mutli_gb, multi_gav, and multi_gexist (see help). If the group boundary array is present in the file, it is preferred to the group energy array, which will then be ignored. The file should specify group boundaries if its opacity and emissivity are averaged over finite width bins; group energies if its opacity and emissivity are computed at points. The group existence map, if present, allows several disjoint spectral regions to exist in a single file. If the data type of gexist is not "string", it should be an array of length one less than gb, if gb is present, or gav, otherwise. By this means you can ignore spectral regions which are present in the file. gscale=photon_energy_unit is an optional conversion factor to bring the units of the gb (or gav) arrays into agreement among the various files which are to be used in a single run. The default value is 1.0 (i.e.- all files are expected to have the same units). tscale=time_unit is an optional conversion factor to bring the units of the time into agreement among the various files which are to be used in a single run. The default value is 1.0 (i.e.- all files are expected to have the same units). noextrap=1 if present and non-zero prevents the opacity and emissivity data from this file from being extrapolated as 1/hnu^3 in master bins at energies above the highest energy bin in this file. freqfirst=0 if present and non-zero means the frequency index is first for the opacity and emissivity arrays, instead of the default of frequency index last. Interpreted function, defined at i/multi.i line 62SEE ALSO: multic, multif, multi_streak, MultiFile, multi_opac, multi_emiss,
multi_srcf, multi_gb, multi_gav, multi_zonelist
mxbrent
fmax= mxbrent(f, x0, x1, x2) or fmax= mxbrent(f, x0, x1, x2, xmax) or fmax= mxbrent(f, x0, x1, x2, xmax, xerr) returns the maximum of the function F (of a single argument x), given three points X0, X1, and X2 such that F(X1) is greater than either F(X0) or F(X2), and X1 is between X0 and X2. If the XMAX argument is provided, it is set to the x value which produced FMAX. If XERR is supplied, the search stops when a fractional error of XERR in x is reached; note that XERR smaller than the square root of the machine precision (or omitted) will cause convergence to machine precision in FMAX. The algorithm is Brent's method - a combination of inverse parabolic interpolation and golden section search - as adapted from Numerical Recipes Ch. 10 (Press, et. al.). Interpreted function, defined at i/roots.i line 210SEE ALSO: mxbrent, nraphson, f_inverse