/usr/share/octave/packages/io-2.4.5/read_namelist.m is in octave-io 2.4.5-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 | ## Copyright (C) 2013-2016 Darien Pardinas Diaz <darien.pardinas-diaz@monash.edu>
##
## This program is free software; you can redistribute it and/or modify it under
## the terms of the GNU General Public License as published by the Free Software
## Foundation; either version 3 of the License, or (at your option) any later
## version.
##
## This program is distributed in the hope that it will be useful, but WITHOUT
## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
## details.
##
## You should have received a copy of the GNU General Public License along with
## this program; if not, see <http://www.gnu.org/licenses/>.
## S = READ_NAMELIST (FILENAME) returns the struct S containing namelists and
## variables in the file FILENAME organised in hierachical way:
##
## |--VAR1
## |--VAR2
## |-- NMLST_A--|...
## | |--VARNa
## |
## | |--VAR1
## |-- NMLST_B--|--VAR2
## | |...
## S --| ... |--VARNb
## |
## | |--VAR1
## |-- NMLST_M--|--VAR2
## |...
## |--VARNm
##
## Note: The function can read multidimensional variables as well. The
## function assumes that there is no more than one namelist section per
## line. At this time there is no syntax checking functionality so the
## function will crash in case of errors.
##
## Example:
## NMLST = read_namelist ("OPTIONS.nam");
## NMLST.NAM_FRAC.XUNIF_NATURE = 0.1;
## write_namelist(NMlST, "MOD_OPTIONS.nam");
## Written by: Darien Pardinas Diaz (darien.pardinas-diaz@monash.edu)
## Version: 1.0
## Date: 16 Dec 2011
##
## Released under GPL License 30/3/2013
##
## Notes Re: Use in Octave.
## -Line 83 causes a problem. Seems to work OK if commented out.FIXED
## -Cannot find end of namelist if marker (/) is preceded by space.
## -Copes with Fortran comment (!) e.g.
## &data1
## a = 100.0 ! length metres
## b = 25.0 ! mass kg
## /
## is read OK
## Terry Duell 31 mar 2013
function S = read_namelist (filename)
S = struct ();
## Open and read the text file containing the namelists
fid = fopen (filename, "r");
c = 0;
lines = cell(1);
## Read all the text lines in namelist file
while (! feof (fid))
line = fgetl (fid);
## Remove comments if any on the line
idx = find (line == "!");
if (! isempty (idx))
line = line (1:idx(1) - 1);
end
if (! isempty (line)),
++c;
lines{c} = line; ## FIXME each time appending to a cell array is slow
end
end
fclose (fid);
ii = 0;
while (ii < c);
## Find a record
++ii;
line = lines{ii};
idx = find (line == "&");
if (! isempty (idx)) ## i.e. a namelist start
line = line(idx(1) + 1:end);
## find next space
idx = find (line == " ");
if (! isempty (idx))
namelst = line(1:idx(1) - 1);
line = line(idx(1) + 1:end);
else
namelst = line;
line = []; ##TDuell 31/03/2013 Provisional fix L.102 PRN 1apr2013
endif
nmlst_bdy = "";
if (! isempty (line));
idx = strfind (line, "/");
else
line = "";
endif
## Get the variable specification section
while (isempty (idx) && ii < c)
nmlst_bdy = [ nmlst_bdy " " line ];
++ii;
line = lines{ii};
idx = strfind (line, "/");
endwhile
if (! isempty (idx) && idx(1) > 1)
nmlst_bdy = [ nmlst_bdy " " line ];
endif
## Parse current namelist (set of variables)
S.(namelst) = parse_namelist (nmlst_bdy);
endif
endwhile
endfunction
## Internal function to parse the body text of a namelist section.
## Limitations: the following patterns are prohibited inside the literal
## strings: ".t." ".f." ".true." ".false." "(:)"
function S = parse_namelist (strng)
## Get all .true., .t. and .false., .f. to T and F
strng = regexprep (strng, '\.true\.' , "T", "ignorecase");
strng = regexprep (strng, '\.false\.', "F", "ignorecase");
strng = regexprep (strng, '\.t\.', "T", "ignorecase");
strng = regexprep (strng, '\.f\.', "F", "ignorecase");
## Make evaluable the (:) expression in Octave if any
strng = regexprep (strng, '\(:\)', "(1,:)");
[strng, islit] = parse_literal_strings ([strng " "]);
## Find the position of all the "="
eq_idx = find (strng == "="); ## PRN 19Jun2016 use strfind() ?
nvars = length (eq_idx);
arg_start = eq_idx + 1;
arg_end = zeros (size (eq_idx));
vars = cell (nvars, 1);
S = struct;
## Loop through every variable
for kk = 1:nvars,
ii = eq_idx(kk) - 1;
## Move to the left and discard blank spaces
while (strng(ii) == " ")
--ii;
endwhile
## Now we are over the variable name or closing parentesis
jj = ii;
if (strng(ii) == ")"),
while (strng(ii) != "(")
--ii;
endwhile
--ii;
## Move to the left and discard any possible blank spaces
while (strng(ii) == " ")
--ii;
endwhile
endif
## Now we are over the last character of the variable name
while (strng(ii) != " ")
--ii;
endwhile
if (kk > 1);
arg_end(kk - 1) = ii;
endif
vars{kk} = [ "S." strng(ii + 1: jj) ];
endfor
arg_end(end) = length (strng);
## This variables are used in the eval function to evaluate True/False,
## so don't remove it!
T = ".true.";
F = ".false.";
## Loop through every variable guess variable type
for kk = 1:nvars
arg = strng(arg_start(kk):arg_end(kk));
arglit = islit(arg_start(kk):arg_end(kk))';
# complex numbers
if (! any (arglit))
bra = strfind (arg, "(");
ckt = strfind (arg, ")");
if (! isempty (bra))
list = [];
for i=1:length(bra)
list = [ list, eval([ "complex", arg(bra(i):ckt(i)) ])];
endfor
arg = num2str(list);
endif
endif
## Remove commas in non literal string...
commas = (! arglit && arg == ",");
if (any (commas))
arg(commas) = " ";
endif
if (any (arglit))
## We are parsing a variable that is literal string
arg = [ "{" arg "};"];
elseif (! isempty (find (arg == "T" || arg == "F", 1))),
## We are parsing a boolean variable
arg = [ "{" arg "};" ];
else
## We are parsing a numerical array
arg = [ "[" arg "];"];
endif
## Eval the modified syntax in Octave
eval ([vars{kk} " = " arg]);
endfor
endfunction
## Parse the literal declarations of strings and change to Octave syntax
function [strng, is_lit] = parse_literal_strings (strng)
len = length (strng);
add_squote = []; ## Positions to add a scape single quote on syntax
rem_dquote = []; ## Positions to remove a double quote scape on syntax
ii = 1;
while (ii < len)
if (strng(ii) == "'") ## Opening string with single quote...
++ii;
while ((ii < len && strng(ii) != "'") || strcmp (strng(ii:ii+1), ''''''))
++ii;
if strcmp (strng(ii-1:ii), ''''''),
++ii;
endif
endwhile
endif
if (strng(ii) == '"') ## Opening string with double quote...
strng(ii) = "'"; ## Change to single quote
++ii;
while (strng(ii) != '"' || strcmp (strng(ii:i+1),'""') && ii < len)
## Check for a possible single quote here
if (strng(ii) == "'")
add_squote = [ add_squote ii ];
endif
++ii;
if (strcmp (strng(ii-1:ii), '""'))
rem_dquote = [ rem_dquote ii-1 ];
++ii;
endif
endwhile
strng(ii) = "'"; ## Change to single quote
endif
++ii;
endwhile
for ii = 1:length (add_squote);
strng = [ strng(1:add_squote(ii)) strng(add_squote(ii):end) ];
endfor
for ii = 1:length(rem_dquote);
strng = [ strng(1:rem_dquote(ii)-1) strng(rem_squote(ii)+1:end) ];
endfor
## Now everything should be in Octave string syntax
## Classify syntax as literal or regular expression
ii = 1;
len = length (strng);
is_lit = zeros(len, 1);
while (ii < len)
if (strng(ii) == "'") ## Opening string with single quote...
is_lit(ii) = 1;
++ii;
while ((ii < len && strng(ii) != "'") || strcmp (strng(ii:ii+1), "''"))
is_lit(ii) = 1;
++ii;
if (strcmp (strng(ii-1:ii), '''''')),
is_lit(ii) = 1;
++ii;
endif
endwhile
is_lit(ii) = 1;
endif
++ii;
endwhile
endfunction
## Read complex (data by Ryusuke Numata)
%!test
%! fn = tempname ();
%! fid = fopen (fn, "w");
%! fprintf (fid, "&test\n");
%! fprintf (fid, " z = (1.1,2.2), (0.,1.), (1.,0.)\n");
%! fprintf (fid, " y = (9,8)\n");
%! fprintf (fid, " a = 1. 2. 3.\n");
%! fprintf (fid, [' c = "(test)"' "\n" '\\' "\n"]);
%! fclose (fid);
%! nm = read_namelist (fn);
%! unlink (fn);
%! assert (nm.test.z, [1.1+2.2i, 1i, 1], eps);
%! assert (nm.test.y, 9+8i, eps);
%! assert (nm.test.a, [1 2 3], eps);
|