mathlib-odrpack.c
1 //! @file mathlib-odrpack.c
2 //
3 //! @section copyright
4 //
5 // This file is part of VIF - vintage fortran compiler.
6 // Copyright 2020-2025 J. Marcel van der Veer <algol68g@xs4all.nl>.
7 //
8 //! @section license
9 //
10 // This program is free software; you can redistribute it and/or modify it
11 // under the terms of the gnu general public license as published by the
12 // free software foundation; either version 3 of the license, or
13 // (at your option) any later version.
14 //
15 // This program is distributed in the hope that it will be useful, but
16 // without any warranty; without even the implied warranty of merchantability
17 // or fitness for a particular purpose. See the GNU general public license for
18 // more details. you should have received a copy of the GNU general public
19 // license along with this program. If not, see <http://www.gnu.org/licenses/>.
20 //
21 //! @Synopsis
22 //
23 //! Subprograms from ODRPACK
24
25 // Compiled from Fortran source code by VIF.
26
27 // ODRPACK 2.01 -- Software for Weighted Orthogonal Distance Regression
28 //
29 // ODRPACK originates from the National Institute of Standards and Technology.
30 // The source code for ODRPACK was retrieved from:
31 //
32 // netlib.org/odrpack
33 //
34 // ODRPACK currently lives on for instance in Python library scipy.odr.
35 //
36 // VIF only includes the double precision version of ODRPACK.
37 //
38 // The code needed trivial adaptation for use with VIF - the only modification
39 // needed was to replace the variable name 'UNIT' that is a reserved word in VIF,
40 // with 'FILENO'.
41 //
42
43 // The license for ODRPACK Fortran source code is:
44 //
45 // NIST-developed software is provided by NIST as a public service. You may use,
46 // copy and distribute copies of the software in any medium, provided that you
47 // keep intact this entire notice. You may improve, modify and create derivative
48 // works of the software or any portion of the software, and you may copy and
49 // distribute such modifications or works. Modified works should carry a notice
50 // stating that you changed the software and should note the date and nature of
51 // any such change. Please explicitly acknowledge the National Institute of
52 // Standards and Technology as the source of the software.
53 //
54 // NIST-developed software is expressly provided "AS IS." NIST MAKES NO WARRANTY
55 // OF ANY KIND, EXPRESS, IMPLIED, IN FACT OR ARISING BY OPERATION OF LAW,
56 // INCLUDING, WITHOUT LIMITATION, THE IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
57 // FOR A PARTICULAR PURPOSE, NON-INFRINGEMENT AND DATA ACCURACY. NIST NEITHER
58 // REPRESENTS NOR WARRANTS THAT THE OPERATION OF THE SOFTWARE WILL BE
59 // UNINTERRUPTED OR ERROR-FREE, OR THAT ANY DEFECTS WILL BE CORRECTED. NIST DOES
60 // NOT WARRANT OR MAKE ANY REPRESENTATIONS REGARDING THE USE OF THE SOFTWARE OR
61 // THE RESULTS THEREOF, INCLUDING BUT NOT LIMITED TO THE CORRECTNESS, ACCURACY,
62 // RELIABILITY, OR USEFULNESS OF THE SOFTWARE.
63 //
64 // You are solely responsible for determining the appropriateness of using and
65 // distributing the software and you assume all risks associated with its use,
66 // including but not limited to the risks and costs of program errors, compliance
67 // with applicable laws, damage to or loss of data, programs or equipment, and the
68 // unavailability or interruption of operation. This software is not intended to
69 // be used in any situation where a failure could cause risk of injury or damage
70 // to property. The software developed by NIST employees is not subject to
71 // copyright protection within the United States.
72
73 // VIF THU 20 NOV 2025 22:21:00 ** MARCEL ** ODRPACK PAGE 00001
74
75
76
77 // M M A RRRRRR CCCCC EEEEEEE L
78 // MM MM A A R R C C E L
79 // M M M M A A R R C E L
80 // M M M A A RRRRRR C EEEEE L
81 // M M AAAAAAA R R C E L
82 // M M A A R R C C E L
83 // M M A A R R CCCCC EEEEEEE LLLLLLL
84
85 // OOOOO DDDDDD RRRRRR PPPPPP A CCCCC K K
86 // O O D D R R P P A A C C K K
87 // O O D D R R P P A A C K K
88 // O O D D RRRRRR PPPPPP A A C KKK
89 // O O D D R R P AAAAAAA C K K
90 // O O D D R R P A A C C K K
91 // OOOOO DDDDDD R R P A A CCCCC K K
92
93
94
95
96 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DIAGNOSTICS PAGE 00002
97
98 // ** linker ** object size 1243120 bytes
99 // ** statistics ** 56 subprograms, no errors, no warnings
100
101
102 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** JOB CONTROL PAGE 00003
103
104
105 // Line JCL *...*....1....|....2....|....3....|....4....|....5....|....6....|....7..*.|....8
106
107
108
109 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00004
110
111 /*
112 Generated by VIF - experimental VIntage Fortran compiler.
113 VIF release 1.4.1
114 */
115
116 #if defined (__GNUC__)
117 #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
118 #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
119 #if (__GNUC__ >= 14)
120 #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
121 #pragma GCC diagnostic ignored "-Wimplicit-int"
122 #pragma GCC diagnostic ignored "-Wint-conversion"
123 #pragma GCC diagnostic ignored "-Wreturn-mismatch"
124 #endif
125 #else
126 #error VIF requires GCC
127 #endif
128
129 #include <vif.h>
130
131 static int_4 _km1 = -1, _k0 = 0, _k1 = 1;
132
133
134 #define _dc_0 " no"
135 #define _dc_1 " yes"
136 #define _dc_2 " unchecked"
137 #define _dc_3 " verified"
138 #define _dc_4 " questionable"
139 #define _dc_5 " "
140 #define _dc_6 "1,"
141 #define _dc_7 "n,"
142 #define _dc_8 " --- problem size:"
143 #define _dc_9 " n = "
144 #define _dc_10 " (number with nonzero weight = "
145 #define _dc_11 ")"
146 #define _dc_12 " nq = "
147 #define _dc_13 " m = "
148 #define _dc_14 " np = "
149 #define _dc_15 " (number unfixed = "
150 #define _dc_16 " --- control values:"
151 #define _dc_17 " job = "
152 #define _dc_18 " = abcde, where"
153 #define _dc_19 " a="
154 #define _dc_20 " ==> fit is a restart."
155 #define _dc_21 " ==> fit is not a restart."
156 #define _dc_22 " b="
157 #define _dc_23 " ==> deltas are initialized"
158 #define _dc_24 " to zero."
159 #define _dc_25 " by user."
160 #define _dc_26 " ==> deltas are fixed at"
161 #define _dc_27 " zero since e="
162 #define _dc_28 "."
163 #define _dc_29 " c="
164 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00005
165
166 #define _dc_30 " ==> covariance matrix will"
167 #define _dc_31 " be computed using"
168 #define _dc_32 " derivatives re-"
169 #define _dc_33 "evaluated at the solution."
170 #define _dc_34 " derivatives from the"
171 #define _dc_35 " last iteration."
172 #define _dc_36 " not be computed."
173 #define _dc_37 " d="
174 #define _dc_38 " ==> derivatives are"
175 #define _dc_39 " supplied by user."
176 #define _dc_40 " derivatives were checked."
177 #define _dc_41 " results appear questionable."
178 #define _dc_42 " results appear correct."
179 #define _dc_43 " derivatives were not"
180 #define _dc_44 " checked."
181 #define _dc_45 " estimated by central"
182 #define _dc_46 " differences."
183 #define _dc_47 " estimated by forward"
184 #define _dc_48 " e="
185 #define _dc_49 " ==> method is implicit odr."
186 #define _dc_50 " ==> method is explicit odr."
187 #define _dc_51 " ==> method is explicit ols."
188 #define _dc_52 " ndigit = "
189 #define _dc_53 " (estimated by odrpack)"
190 #define _dc_54 " (supplied by user)"
191 #define _dc_55 " taufac = "
192 #define _dc_56 " --- stopping criteria:"
193 #define _dc_57 " sstol = "
194 #define _dc_58 " (sum of squares stopping tolerance)"
195 #define _dc_59 " partol = "
196 #define _dc_60 " (parameter stopping tolerance)"
197 #define _dc_61 " maxit = "
198 #define _dc_62 " (maximum number of iterations)"
199 #define _dc_63 " --- initial sum of squared weighted deltas ="
200 #define _dc_64 " "
201 #define _dc_65 " initial penalty function value ="
202 #define _dc_66 " penalty term ="
203 #define _dc_67 " penalty parameter ="
204 #define _dc_68 " --- initial weighted sum of squares ="
205 #define _dc_69 " sum of squared weighted deltas ="
206 #define _dc_70 " sum of squared weighted epsilons ="
207 #define _dc_71 " --- explanatory variable and delta weight summary:"
208 #define _dc_72 " --- explanatory variable summary:"
209 #define _dc_73 " index x(i,j) delta(i,j) fixed"
210 #define _dc_74 " scale weight derivative"
211 #define _dc_75 " "
212 #define _dc_76 " assessment"
213 #define _dc_77 " (i,j) (ifixx)"
214 #define _dc_78 " (scld) (wd) "
215 #define _dc_79 " scale weight "
216 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00006
217
218 #define _dc_80 " "
219 #define _dc_81 " step size"
220 #define _dc_82 " (scld) (wd) (stpd)"
221 #define _dc_83 " index x(i,j)"
222 #define _dc_84 " (i,j) "
223 #define _dc_85 " --- response variable and epsilon error weight"
224 #define _dc_86 " summary:"
225 #define _dc_87 " index y(i,l) weight"
226 #define _dc_88 " (i,l) (we)"
227 #define _dc_89 " --- function parameter summary:"
228 #define _dc_90 " index beta(k) fixed scale"
229 #define _dc_91 " derivative"
230 #define _dc_92 " "
231 #define _dc_93 " assessment"
232 #define _dc_94 " (k) (ifixb) (sclb)"
233 #define _dc_95 " "
234 #define _dc_96 " step size"
235 #define _dc_97 " (stpb)"
236 #define _dc_98 " "
237 #define _dc_99 " "
238 #define _dc_100 " "
239 #define _dc_101 " "
240 #define _dc_102 " "
241 #define _dc_103 "yes"
242 #define _dc_104 " no"
243 #define _dc_105 " cum. penalty act. rel. pred. rel."
244 #define _dc_106 " it. no. fn function sum-of-sqs sum-of-sqs"
245 #define _dc_107 " g-n"
246 #define _dc_108 " num. evals value reduction reduction"
247 #define _dc_109 " tau/pnorm step"
248 #define _dc_110 " ---- ------ ----------- ----------- -----------"
249 #define _dc_111 " --------- ----"
250 #define _dc_112 " cum. act. rel. pred. rel."
251 #define _dc_113 " it. no. fn weighted sum-of-sqs sum-of-sqs"
252 #define _dc_114 " num. evals sum-of-sqs reduction reduction"
253 #define _dc_115 " g-n beta -------------->"
254 #define _dc_116 " tau/pnorm step index value"
255 #define _dc_117 " --------- ---- ----- -----"
256 #define _dc_118 " penalty parameter value = "
257 #define _dc_119 " "
258 #define _dc_120 " to"
259 #define _dc_121 " "
260 #define _dc_122 " "
261 #define _dc_123 " --- stopping conditions:"
262 #define _dc_124 " info = "
263 #define _dc_125 " ==> sum of squares convergence."
264 #define _dc_126 " ==> parameter convergence."
265 #define _dc_127 " ==> sum of squares convergence and"
266 #define _dc_128 " parameter convergence."
267 #define _dc_129 " ==> iteration limit reached."
268 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00007
269
270 #define _dc_130 " ==> unexpected value,"
271 #define _dc_131 " probably indicating"
272 #define _dc_132 " incorrectly specified"
273 #define _dc_133 " user input."
274 #define _dc_134 " = abcd, where a nonzero value for digit a,"
275 #define _dc_135 " b, or c indicates why"
276 #define _dc_136 " the results might be questionable,"
277 #define _dc_137 " and digit d indicates"
278 #define _dc_138 " the actual stopping condition."
279 #define _dc_139 " a=1 ==> derivatives are"
280 #define _dc_140 " questionable."
281 #define _dc_141 " b=1 ==> user set istop to"
282 #define _dc_142 " nonzero value during last"
283 #define _dc_143 " call to subroutine fcn."
284 #define _dc_144 " c=1 ==> derivatives are not"
285 #define _dc_145 " full rank at the solution."
286 #define _dc_146 " c=2 ==> derivatives are zero"
287 #define _dc_147 " rank at the solution."
288 #define _dc_148 " d=1 ==> sum of squares convergence."
289 #define _dc_149 " d=2 ==> parameter convergence."
290 #define _dc_150 " d=3 ==> sum of squares convergence"
291 #define _dc_151 " and parameter convergence."
292 #define _dc_152 " d=4 ==> iteration limit reached."
293 #define _dc_153 " incorrectly specified"
294 #define _dc_154 " = abcde, where a nonzero value for a given"
295 #define _dc_155 " digit indicates an"
296 #define _dc_156 " abnormal stopping condition."
297 #define _dc_157 " a=5 ==> user stopped computations"
298 #define _dc_158 " in subroutine fcn."
299 #define _dc_159 " ==> computations were"
300 #define _dc_160 " stopped during the"
301 #define _dc_161 " function evaluation."
302 #define _dc_162 " stopped because"
303 #define _dc_163 " derivatives with"
304 #define _dc_164 " respect to delta were"
305 #define _dc_165 " computed by"
306 #define _dc_166 " subroutine fcn when"
307 #define _dc_167 " fit is ols."
308 #define _dc_168 " jacobian evaluation."
309 #define _dc_169 " a=6 ==> numerical instabilities"
310 #define _dc_170 " have been detected,"
311 #define _dc_171 " possibly indicating"
312 #define _dc_172 " a discontinuity in the"
313 #define _dc_173 " derivatives or a poor"
314 #define _dc_174 " poor choice of problem"
315 #define _dc_175 " scale or weights."
316 #define _dc_176 " niter = "
317 #define _dc_177 " (number of iterations)"
318 #define _dc_178 " nfev = "
319 #define _dc_179 " (number of function evaluations)"
320 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00008
321
322 #define _dc_180 " njev = "
323 #define _dc_181 " (number of jacobian evaluations)"
324 #define _dc_182 " irank = "
325 #define _dc_183 " (rank deficiency)"
326 #define _dc_184 " rcond = "
327 #define _dc_185 " (inverse condition number)"
328 #define _dc_186 " istop = "
329 #define _dc_187 " (returned by user from"
330 #define _dc_188 " subroutine fcn)"
331 #define _dc_189 " --- final sum of squared weighted deltas = "
332 #define _dc_190 " final penalty function value = "
333 #define _dc_191 " penalty term = "
334 #define _dc_192 " penalty parameter = "
335 #define _dc_193 " --- final weighted sums of squares = "
336 #define _dc_194 " sum of squared weighted deltas = "
337 #define _dc_195 " sum of squared weighted epsilons = "
338 #define _dc_196 " --- residual standard deviation = "
339 #define _dc_197 " degrees of freedom ="
340 #define _dc_198 " --- estimated beta(j), j = 1, ..., np:"
341 #define _dc_199 " --- estimated delta(i,*), i = 1, ..., n:"
342 #define _dc_200 " --- estimated epsilon(i) and delta(i,*), i = 1, ..., n:"
343 #define _dc_201 " --- estimated epsilon(i), i = 1, ..., n:"
344 #define _dc_202 " "
345 #define _dc_203 " --- estimated epsilon(i,"
346 #define _dc_204 "), i = 1, ..., n:"
347 #define _dc_205 " --- estimated delta(i,"
348 #define _dc_206 " index value"
349 #define _dc_207 " index value -------------->"
350 #define _dc_208 " beta s.d. beta"
351 #define _dc_209 " ---- 95% confidence interval ----"
352 #define _dc_210 " n.b. standard errors and confidence intervals are"
353 #define _dc_211 " computed using"
354 #define _dc_212 " derivatives calculated at the beginning"
355 #define _dc_213 " of the last iteration,"
356 #define _dc_214 " and not using derivatives re-evaluated at the"
357 #define _dc_215 " final solution."
358 #define _dc_216 " n.b. the standard errors of the estimated betas were"
359 #define _dc_217 " not computed because"
360 #define _dc_218 " the derivatives were not available. either maxit"
361 #define _dc_219 " is 0 and the third"
362 #define _dc_220 " digit of job is greater than 1, or the most"
363 #define _dc_221 " recently tried values of"
364 #define _dc_222 " beta and/or x+delta were identified as"
365 #define _dc_223 " unacceptable by user supplied"
366 #define _dc_224 " subroutine fcn."
367 #define _dc_225 " not computed."
368 #define _dc_226 " (see info above.)"
369 #define _dc_227 " beta status"
370 #define _dc_228 " "
371 #define _dc_229 "to"
372 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00009
373
374 #define _dc_230 " "
375 #define _dc_231 "estimated"
376 #define _dc_232 " fixed"
377 #define _dc_233 " dropped"
378 #define _dc_234 " n.b. no parameters were fixed by the user or"
379 #define _dc_235 " dropped at the last"
380 #define _dc_236 " iteration because they caused the model to be"
381 #define _dc_237 " rank deficient."
382 #define _dc_238 " n.b. no change was made to the user supplied parameter"
383 #define _dc_239 " values because"
384 #define _dc_240 " maxit=0."
385 #define _dc_241 "(/' i',"
386 #define _dc_242 "(' delta(i,',i1,')')/)"
387 #define _dc_243 "(' epsilon(i,',i1,')'),"
388 #define _dc_244 "(' epsilon(i,',i1,')')/)"
389 #define _dc_245 "odr"
390 #define _dc_246 "ols"
391 #define _dc_247 " *** initial summary for fit by method of "
392 #define _dc_248 " ***"
393 #define _dc_249 " *** iteration reports for fit by method of "
394 #define _dc_250 " *** final summary for fit by method of "
395 #define _dc_251 " error : n is less than one."
396 #define _dc_252 " error : m is less than one."
397 #define _dc_253 " error : np is less than one"
398 #define _dc_254 " or np is greater than n."
399 #define _dc_255 " error : nq is less than one."
400 #define _dc_256 " error : ldx is less than n."
401 #define _dc_257 " error : ldy is less than n."
402 #define _dc_258 " error : ldifx is less than n"
403 #define _dc_259 " and ldifx is not equal to one."
404 #define _dc_260 " error : ldscld is less than n"
405 #define _dc_261 " and ldscld is not equal to one."
406 #define _dc_262 " error : ldstpd is less than n"
407 #define _dc_263 " and ldstpd is not equal to one."
408 #define _dc_264 " error : ldwe is less than n"
409 #define _dc_265 " and ldwe is not equal to one or"
410 #define _dc_266 " or"
411 #define _dc_267 " ld2we is less than nq"
412 #define _dc_268 " and ld2we is not equal to one."
413 #define _dc_269 " error : ldwd is less than n"
414 #define _dc_270 " and ldwd is not equal to one."
415 #define _dc_271 " error : lwork is less than "
416 #define _dc_272 ","
417 #define _dc_273 " the smallest acceptable dimension of array work."
418 #define _dc_274 " error : liwork is less than "
419 #define _dc_275 " the smallest acceptable dimension of array"
420 #define _dc_276 " iwork."
421 #define _dc_277 " error : scld(i,j) is less than or equal to zero"
422 #define _dc_278 " for some i = 1, ..., n and j = 1, ..., m."
423 #define _dc_279 " when scld(1,1) is greater than zero"
424 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00010
425
426 #define _dc_280 " and ldscld is greater than or equal to n then"
427 #define _dc_281 " each of the n by m elements of"
428 #define _dc_282 " scld must be greater than zero."
429 #define _dc_283 " error : scld(1,j) is less than or equal to zero"
430 #define _dc_284 " for some j = 1, ..., m."
431 #define _dc_285 " and ldscld is equal to one then"
432 #define _dc_286 " each of the 1 by m elements of"
433 #define _dc_287 " error : sclb(k) is less than or equal to zero"
434 #define _dc_288 " for some k = 1, ..., np."
435 #define _dc_289 " all np elements of"
436 #define _dc_290 " sclb must be greater than zero."
437 #define _dc_291 " error : stpd(i,j) is less than or equal to zero"
438 #define _dc_292 " when stpd(1,1) is greater than zero"
439 #define _dc_293 " and ldstpd is greater than or equal to n then"
440 #define _dc_294 " stpd must be greater than zero."
441 #define _dc_295 " error : stpd(1,j) is less than or equal to zero"
442 #define _dc_296 " and ldstpd is equal to one then"
443 #define _dc_297 " error : stpb(k) is less than or equal to zero"
444 #define _dc_298 " stpb must be greater than zero."
445 #define _dc_299 " error : at least one of the (nq by nq) arrays starting"
446 #define _dc_300 " in we(i,1,1), i = 1, ..., n, is not positive"
447 #define _dc_301 " semidefinite. when we(1,1,1) is greater than"
448 #define _dc_302 " or equal to zero, and ldwe is greater than or"
449 #define _dc_303 " equal to n, and ld2we is greater than or equal"
450 #define _dc_304 " to nq, then each of the (nq by nq) arrays in we"
451 #define _dc_305 " must be positive semidefinite."
452 #define _dc_306 " error : at least one of the (1 by nq) arrays starting"
453 #define _dc_307 " in we(i,1,1), i = 1, ..., n, has a negative"
454 #define _dc_308 " element. when we(1,1,1) is greater than or"
455 #define _dc_309 " equal to zero, and ldwe is greater than or equal"
456 #define _dc_310 " to n, and ld2we is equal to 1, then each of the"
457 #define _dc_311 " (1 by nq) arrays in we must have only non-"
458 #define _dc_312 " negative elements."
459 #define _dc_313 " error : the (nq by nq) array starting in we(1,1,1) is"
460 #define _dc_314 " not positive semidefinite. when we(1,1,1) is"
461 #define _dc_315 " greater than or equal to zero, and ldwe is equal"
462 #define _dc_316 " to 1, and ld2we is greater than or equal to nq,"
463 #define _dc_317 " then the (nq by nq) array in we must be positive"
464 #define _dc_318 " semidefinite."
465 #define _dc_319 " error : the (1 by nq) array starting in we(1,1,1) has"
466 #define _dc_320 " a negative element. when we(1,1,1) is greater"
467 #define _dc_321 " than or equal to zero, and ldwe is equal to 1,"
468 #define _dc_322 " and ld2we is equal to 1, then the (1 by nq)"
469 #define _dc_323 " array in we must have only nonnegative elements."
470 #define _dc_324 " error : the number of nonzero arrays in array we is"
471 #define _dc_325 " less than np."
472 #define _dc_326 " error : at least one of the (m by m) arrays starting"
473 #define _dc_327 " in wd(i,1,1), i = 1, ..., n, is not positive"
474 #define _dc_328 " definite. when wd(1,1,1) is greater than zero,"
475 #define _dc_329 " and ldwd is greater than or equal to n, and"
476 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00011
477
478 #define _dc_330 " ld2wd is greater than or equal to m, then each"
479 #define _dc_331 " of the (m by m) arrays in wd must be positive"
480 #define _dc_332 " definite."
481 #define _dc_333 " error : at least one of the (1 by m) arrays starting"
482 #define _dc_334 " in wd(i,1,1), i = 1, ..., n, has a nonpositive"
483 #define _dc_335 " element. when wd(1,1,1) is greater than zero,"
484 #define _dc_336 " ld2wd is equal to 1, then each of the (1 by m)"
485 #define _dc_337 " arrays in wd must have only positive elements."
486 #define _dc_338 " error : the (m by m) array starting in wd(1,1,1) is"
487 #define _dc_339 " not positive definite. when wd(1,1,1) is"
488 #define _dc_340 " greater than zero, and ldwd is equal to 1, and"
489 #define _dc_341 " ld2wd is greater than or equal to m, then the"
490 #define _dc_342 " (m by m) array in wd must be positive definite."
491 #define _dc_343 " error : the (1 by m) array starting in wd(1,1,1) has a"
492 #define _dc_344 " nonpositive element. when wd(1,1,1) is greater"
493 #define _dc_345 " than zero, and ldwd is equal to 1, and ld2wd is"
494 #define _dc_346 " equal to 1, then the (1 by m) array in wd must"
495 #define _dc_347 " have only positive elements."
496 #define _dc_348 "*"
497 #define _dc_349 " *** derivative checking report for fit by method of "
498 #define _dc_350 " for response "
499 #define _dc_351 " of observation "
500 #define _dc_352 " "
501 #define _dc_353 " user"
502 #define _dc_354 " "
503 #define _dc_355 " "
504 #define _dc_356 " supplied"
505 #define _dc_357 " relative"
506 #define _dc_358 " derivative "
507 #define _dc_359 " derivative wrt"
508 #define _dc_360 " value"
509 #define _dc_361 " difference"
510 #define _dc_362 " assessment "
511 #define _dc_363 " beta("
512 #define _dc_364 " --- "
513 #define _dc_365 "verified"
514 #define _dc_366 "questionable (see note "
515 #define _dc_367 " delta("
516 #define _dc_368 " notes:"
517 #define _dc_369 " (1) user supplied and finite difference derivatives"
518 #define _dc_370 " agree, but"
519 #define _dc_371 " results are questionable because both are zero."
520 #define _dc_372 " (2) user supplied and finite difference derivatives"
521 #define _dc_373 " results are questionable because one is"
522 #define _dc_374 " identically zero"
523 #define _dc_375 " and the other is only approximately zero."
524 #define _dc_376 " (3) user supplied and finite difference derivatives"
525 #define _dc_377 " disagree, but"
526 #define _dc_378 " and the other is not."
527 #define _dc_379 " (4) user supplied and finite difference derivatives"
528 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00012
529
530 #define _dc_380 " finite difference derivative is questionable"
531 #define _dc_381 " because either"
532 #define _dc_382 " the ratio of relative curvature to relative"
533 #define _dc_383 " slope is too high"
534 #define _dc_384 " or the scale is wrong."
535 #define _dc_385 " (5) user supplied and finite difference derivatives"
536 #define _dc_386 " because the"
537 #define _dc_387 " ratio of relative curvature to relative slope is"
538 #define _dc_388 " too high."
539 #define _dc_389 " (6) user supplied and finite difference derivatives"
540 #define _dc_390 " have at least 2 digits in common."
541 #define _dc_391 " (7) user supplied and finite difference derivatives"
542 #define _dc_392 " disagree, and"
543 #define _dc_393 " have fewer than 2 digits in common. derivative"
544 #define _dc_394 " checking must"
545 #define _dc_395 " be turned off in order to proceed."
546 #define _dc_396 " number of reliable digits in function results "
547 #define _dc_397 " (estimated by odrpack)"
548 #define _dc_398 " (supplied by user)"
549 #define _dc_399 " number of digits of agreement required between "
550 #define _dc_400 " user supplied and finite difference derivative for "
551 #define _dc_401 " user supplied derivative to be considered verified "
552 #define _dc_402 " row number at which derivatives were checked "
553 #define _dc_403 " -values of the explanatory variables at this row"
554 #define _dc_404 " "
555 #define _dc_405 "x("
556 #define _dc_406 " variable istop has been returned with a nonzero value "
557 #define _dc_407 " from user supplied subroutine fcn when invoked using the"
558 #define _dc_408 " initial estimates of beta and delta supplied by the "
559 #define _dc_409 " user. the initial estimates must be adjusted to allow "
560 #define _dc_410 " proper evaluation of subroutine fcn before the "
561 #define _dc_411 " regression procedure can continue."
562 #define _dc_412 " from user supplied subroutine fcn. this occurred during"
563 #define _dc_413 " the computation of the number of reliable digits in the "
564 #define _dc_414 " predicted values (f) returned from subroutine fcn, indi-"
565 #define _dc_415 " cating that changes in the initial estimates of beta(k),"
566 #define _dc_416 " k=1,np, as small as 2*beta(k)*sqrt(machine precision), "
567 #define _dc_417 " where machine precision is defined as the smallest value"
568 #define _dc_418 " e such that 1+e>1 on the computer being used, prevent "
569 #define _dc_419 " subroutine fcn from being properly evaluated. the "
570 #define _dc_420 " initial estimates must be adjusted to allow proper "
571 #define _dc_421 " evaluation of subroutine fcn during these computations "
572 #define _dc_422 " before the regression procedure can continue."
573 #define _dc_423 " the derivative checking procedure, indicating that "
574 #define _dc_424 " changes in the initial estimates of beta(k), k=1,np, as "
575 #define _dc_425 " small as max[beta(k),1/sclb(k)]*10**(-neta/2), and/or "
576 #define _dc_426 " of delta(i,j), i=1,n and j=1,m, as small as "
577 #define _dc_427 " max[delta(i,j),1/scld(i,j)]*10**(-neta/2), where neta "
578 #define _dc_428 " is defined to be the number of reliable digits in "
579 #define _dc_429 " predicted values (f) returned from subroutine fcn, "
580 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00013
581
582 #define _dc_430 " prevent subroutine fcn from being properly evaluated. "
583 #define _dc_431 " the initial estimates must be adjusted to allow proper "
584 #define _dc_432 " from user supplied subroutine fcn when invoked for "
585 #define _dc_433 " derivative evaluations using the initial estimates of "
586 #define _dc_434 " beta and delta supplied by the user. the initial "
587 #define _dc_435 " estimates must be adjusted to allow proper evaluation "
588 #define _dc_436 " of subroutine fcn before the regression procedure can "
589 #define _dc_437 " continue."
590 #define _dc_438 " the correct form of the call statement is "
591 #define _dc_439 " call dodr"
592 #define _dc_440 " + (fcn,"
593 #define _dc_441 " + n,m,np,nq,"
594 #define _dc_442 " + beta,"
595 #define _dc_443 " + y,ldy,x,ldx,"
596 #define _dc_444 " + we,ldwe,ld2we,wd,ldwd,ld2wd,"
597 #define _dc_445 " + job,"
598 #define _dc_446 " + iprint,lunerr,lunrpt,"
599 #define _dc_447 " + work,lwork,iwork,liwork,"
600 #define _dc_448 " + info)"
601 #define _dc_449 " call dodrc"
602 #define _dc_450 " + ifixb,ifixx,ldifx,"
603 #define _dc_451 " + job,ndigit,taufac,"
604 #define _dc_452 " + sstol,partol,maxit,"
605 #define _dc_453 " + stpb,stpd,ldstpd,"
606 #define _dc_454 " + sclb,scld,ldscld,"
607 #define _dc_455 " ******************************************************* "
608 #define _dc_456 " * odrpack version 2.01 of 06-19-92 (double precision) * "
609
610
611
612 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** TYPEDEFS PAGE 00014
613
614 static FORMAT *__fmt_a = NULL;
615 typedef char char_1[2];
616 typedef char char_3[4];
617 typedef char char_7[8];
618 typedef char char_15[16];
619 typedef char char_127[128];
620
621
622 // VIF THU 20 NOV 2025 22:21:00 ** GLOBAL SCOPE ** FREQUENCY TABLE PAGE 00015
623
624 #define __ncalls 57
625 static CALLS __calls[__ncalls] = {
626 { // ** body not listed **"dacces", 0}, // subroutine
628 { // ** body not listed **"detaf", 0}, // subroutine
630 { // ** body not listed **"dfctr", 0}, // subroutine
632 { // ** body not listed **"dflags", 0}, // subroutine
634 { // ** body not listed **"difix", 0}, // subroutine
636 { // ** body not listed **"diwinf", 0}, // subroutine
638 { // ** body not listed **"djacfd", 0}, // subroutine
640 { // ** body not listed **"djck", 0}, // subroutine
642 { // ** body not listed **"djckm", 0}, // subroutine
644 { // ** body not listed **"dluno", 0}, // subroutine
646 { // ** body not listed **"dmprec", 0}, // real*8 function
648 { // ** body not listed **"dodcnt", 0}, // subroutine
650 { // ** body not listed **"dodlm", 0}, // subroutine
652 { // ** body not listed **"dodpc1", 0}, // subroutine
654 { // ** body not listed **"dodpc3", 0}, // subroutine
656 { // ** body not listed **"dodpe1", 0}, // subroutine
658 { // ** body not listed **"dodpe3", 0}, // subroutine
660 { // ** body not listed **"dodphd", 0}, // subroutine
662 { // ** body not listed **"dodr", 0}, // subroutine
664 { // ** body not listed **"dodvcv", 0}, // subroutine
666 { // ** body not listed **"dppnml", 0}, // real*8 function
668 { // ** body not listed **"dpvb", 0}, // subroutine
670 { // ** body not listed **"dscale", 0}, // subroutine
672 { // ** body not listed **"dscld", 0}, // subroutine
676 { // ** body not listed **"dsolve", 0}, // subroutine
678 { // ** body not listed **"dvevtr", 0}, // subroutine
680 { // ** body not listed **"dwinf", 0}, // subroutine
682 { // ** body not listed **"dxpy", 0}, // subroutine
684 { // ** body not listed **NULL, 0}
685 };
967 { // ** body not listed **
1055 }
1056
1057 // VIF THU 20 NOV 2025 22:21:00 ** DESUBI ** DESUBI PAGE 00027
1058
1059 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1060 // S S U U B B R R O O U U T I NN N E
1061 // S U U B B R R O O U U T I N N N E
1062 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1063 // S U U B B R R O O U U T I N N N E
1064 // S S U U B B R R O O U U T I N NN E
1065 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1066
1067 // DDDDDD EEEEEEE SSSSS U U BBBBBB III
1068 // D D E S S U U B B I
1069 // D D E S U U B B I
1070 // D D EEEEE SSSSS U U BBBBBB I
1071 // D D E S U U B B I
1072 // D D E S S U U B B I
1073 // DDDDDD EEEEEEE SSSSS UUUUU BBBBBB III
1074
1075
1076 // VIF THU 20 NOV 2025 22:21:00 ** DESUBI ** DIAGNOSTICS PAGE 00028
1077
1078
1079 // VIF THU 20 NOV 2025 22:21:00 ** DESUBI ** GENERATED CODE PAGE 00029
1080
1081 int_4 _desubi (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, real_8 _p_ alpha_, real_8
1082 _p_ tt_, int_4 _p_ ldtt_, int_4 _p_ i_, real_8 _p_ e_)
1083 { // ** body not listed **
1204 }
1205
1206 // VIF THU 20 NOV 2025 22:21:00 ** DETAF ** DETAF PAGE 00032
1207
1208 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1209 // S S U U B B R R O O U U T I NN N E
1210 // S U U B B R R O O U U T I N N N E
1211 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1212 // S U U B B R R O O U U T I N N N E
1213 // S S U U B B R R O O U U T I N NN E
1214 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1215
1216 // DDDDDD EEEEEEE TTTTTTT A FFFFFFF
1217 // D D E T A A F
1218 // D D E T A A F
1219 // D D EEEEE T A A FFFFF
1220 // D D E T AAAAAAA F
1221 // D D E T A A F
1222 // DDDDDD EEEEEEE T A A F
1223
1224
1225 // VIF THU 20 NOV 2025 22:21:00 ** DETAF ** DIAGNOSTICS PAGE 00033
1226
1227
1228 // VIF THU 20 NOV 2025 22:21:00 ** DETAF ** GENERATED CODE PAGE 00034
1229
1230 int_4 _detaf (int_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ xplusd_, real_8 _p_
1231 beta_, real_8 _p_ epsmac_, int_4 _p_ nrow_, real_8 _p_ partmp_, real_8 _p_ pv0_, int_4 _p_ ifixb_, int_4 _p_ ifixx_,
1232 int_4 _p_ ldifx_, int_4 _p_ istop_, int_4 _p_ nfev_, real_8 _p_ eta_, int_4 _p_ neta_, real_8 _p_ wrk1_, real_8 _p_
1233 wrk2_, real_8 _p_ wrk6_, real_8 _p_ wrk7_)
1234 { // ** body not listed **
1354 }
1355
1356 // VIF THU 20 NOV 2025 22:21:00 ** DEVJAC ** DEVJAC PAGE 00037
1357
1358 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1359 // S S U U B B R R O O U U T I NN N E
1360 // S U U B B R R O O U U T I N N N E
1361 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1362 // S U U B B R R O O U U T I N N N E
1363 // S S U U B B R R O O U U T I N NN E
1364 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1365
1366 // DDDDDD EEEEEEE V V J A CCCCC
1367 // D D E V V J A A C C
1368 // D D E V V J A A C
1369 // D D EEEEE V V J A A C
1370 // D D E V V J J AAAAAAA C
1371 // D D E V V J J A A C C
1372 // DDDDDD EEEEEEE V JJJJJ A A CCCCC
1373
1374
1375 // VIF THU 20 NOV 2025 22:21:00 ** DEVJAC ** DIAGNOSTICS PAGE 00038
1376
1377
1378 // VIF THU 20 NOV 2025 22:21:00 ** DEVJAC ** GENERATED CODE PAGE 00039
1379
1380 int_4 _devjac (int_4 (*_fcn)(), logical_4 _p_ anajac_, logical_4 _p_ cdjac_, int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_,
1381 int_4 _p_ nq_, real_8 _p_ betac_, real_8 _p_ beta_, real_8 _p_ stpb_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_
1382 ldifx_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ delta_, real_8 _p_ xplusd_, real_8 _p_ stpd_, int_4 _p_ ldstpd_,
1383 real_8 _p_ ssf_, real_8 _p_ tt_, int_4 _p_ ldtt_, int_4 _p_ neta_, real_8 _p_ fn_, real_8 _p_ stp_, real_8 _p_ wrk1_,
1384 real_8 _p_ wrk2_, real_8 _p_ wrk3_, real_8 _p_ wrk6_, real_8 _p_ fjacb_, logical_4 _p_ isodr_, real_8 _p_ fjacd_,
1385 real_8 _p_ we1_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, int_4 _p_ njev_, int_4 _p_ nfev_, int_4 _p_ istop_, int_4 _p_ info_
1386 )
1387 { // ** body not listed **
1475 }
1476
1477 // VIF THU 20 NOV 2025 22:21:00 ** DFCTR ** DFCTR PAGE 00041
1478
1479 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1480 // S S U U B B R R O O U U T I NN N E
1481 // S U U B B R R O O U U T I N N N E
1482 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1483 // S U U B B R R O O U U T I N N N E
1484 // S S U U B B R R O O U U T I N NN E
1485 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1486
1487 // DDDDDD FFFFFFF CCCCC TTTTTTT RRRRRR
1488 // D D F C C T R R
1489 // D D F C T R R
1490 // D D FFFFF C T RRRRRR
1491 // D D F C T R R
1492 // D D F C C T R R
1493 // DDDDDD F CCCCC T R R
1494
1495
1496 // VIF THU 20 NOV 2025 22:21:00 ** DFCTR ** DIAGNOSTICS PAGE 00042
1497
1498
1499 // VIF THU 20 NOV 2025 22:21:00 ** DFCTR ** GENERATED CODE PAGE 00043
1500
1501 int_4 _dfctr (logical_4 _p_ oksemi_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, int_4 _p_ info_)
1502 { // ** body not listed **
1567 }
1568
1569 // VIF THU 20 NOV 2025 22:21:00 ** DFCTRW ** DFCTRW PAGE 00045
1570
1571 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1572 // S S U U B B R R O O U U T I NN N E
1573 // S U U B B R R O O U U T I N N N E
1574 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1575 // S U U B B R R O O U U T I N N N E
1576 // S S U U B B R R O O U U T I N NN E
1577 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1578
1579 // DDDDDD FFFFFFF CCCCC TTTTTTT RRRRRR W W
1580 // D D F C C T R R W W
1581 // D D F C T R R W W
1582 // D D FFFFF C T RRRRRR W W W
1583 // D D F C T R R W W W
1584 // D D F C C T R R W W W
1585 // DDDDDD F CCCCC T R R WW WW
1586
1587
1588 // VIF THU 20 NOV 2025 22:21:00 ** DFCTRW ** DIAGNOSTICS PAGE 00046
1589
1590
1591 // VIF THU 20 NOV 2025 22:21:00 ** DFCTRW ** GENERATED CODE PAGE 00047
1592
1593 int_4 _dfctrw (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ nq_, int_4 _p_ npp_, logical_4 _p_ isodr_, real_8 _p_ we_, int_4
1594 _p_ ldwe_, int_4 _p_ ld2we_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, real_8 _p_ wrk0_, real_8 _p_ wrk4_,
1595 real_8 _p_ we1_, int_4 _p_ nnzw_, int_4 _p_ info_)
1596 { // ** body not listed **
1773 }
1774
1775 // VIF THU 20 NOV 2025 22:21:00 ** DFLAGS ** DFLAGS PAGE 00051
1776
1777 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1778 // S S U U B B R R O O U U T I NN N E
1779 // S U U B B R R O O U U T I N N N E
1780 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1781 // S U U B B R R O O U U T I N N N E
1782 // S S U U B B R R O O U U T I N NN E
1783 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1784
1785 // DDDDDD FFFFFFF L A GGGGG SSSSS
1786 // D D F L A A G G S S
1787 // D D F L A A G S
1788 // D D FFFFF L A A G GGGG SSSSS
1789 // D D F L AAAAAAA G G S
1790 // D D F L A A G G S S
1791 // DDDDDD F LLLLLLL A A GGGGG SSSSS
1792
1793
1794 // VIF THU 20 NOV 2025 22:21:00 ** DFLAGS ** DIAGNOSTICS PAGE 00052
1795
1796
1797 // VIF THU 20 NOV 2025 22:21:00 ** DFLAGS ** GENERATED CODE PAGE 00053
1798
1799 int_4 _dflags (int_4 _p_ job_, logical_4 _p_ restrt_, logical_4 _p_ initd_, logical_4 _p_ dovcv_, logical_4 _p_ redoj_,
1800 logical_4 _p_ anajac_, logical_4 _p_ cdjac_, logical_4 _p_ chkjac_, logical_4 _p_ isodr_, logical_4 _p_ implct_)
1801 { // ** body not listed **
1870 }
1871
1872 // VIF THU 20 NOV 2025 22:21:00 ** DHSTEP ** DHSTEP PAGE 00055
1873
1874 // RRRRRR EEEEEEE A L 88888
1875 // R R E A A L * * 8 8
1876 // R R E A A L * * 8 8
1877 // RRRRRR EEEEE A A L ******* 88888
1878 // R R E AAAAAAA L * * 8 8
1879 // R R E A A L * * 8 8
1880 // R R EEEEEEE A A LLLLLLL 88888
1881
1882 // FFFFFFF U U N N CCCCC TTTTTTT III OOOOO N N
1883 // F U U NN N C C T I O O NN N
1884 // F U U N N N C T I O O N N N
1885 // FFFFF U U N N N C T I O O N N N
1886 // F U U N N N C T I O O N N N
1887 // F U U N NN C C T I O O N NN
1888 // F UUUUU N N CCCCC T III OOOOO N N
1889
1890 // DDDDDD H H SSSSS TTTTTTT EEEEEEE PPPPPP
1891 // D D H H S S T E P P
1892 // D D H H S T E P P
1893 // D D HHHHHHH SSSSS T EEEEE PPPPPP
1894 // D D H H S T E P
1895 // D D H H S S T E P
1896 // DDDDDD H H SSSSS T EEEEEEE P
1897
1898
1899 // VIF THU 20 NOV 2025 22:21:00 ** DHSTEP ** DIAGNOSTICS PAGE 00056
1900
1901
1902 // VIF THU 20 NOV 2025 22:21:00 ** DHSTEP ** GENERATED CODE PAGE 00057
1903
1904 real_8 _dhstep (int_4 _p_ itype_, int_4 _p_ neta_, int_4 _p_ i_, int_4 _p_ j_, real_8 _p_ stp_, int_4 _p_ ldstp_)
1905 { // ** body not listed **
1957 }
1958
1959 // VIF THU 20 NOV 2025 22:21:00 ** DIFIX ** DIFIX PAGE 00059
1960
1961 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1962 // S S U U B B R R O O U U T I NN N E
1963 // S U U B B R R O O U U T I N N N E
1964 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1965 // S U U B B R R O O U U T I N N N E
1966 // S S U U B B R R O O U U T I N NN E
1967 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1968
1969 // DDDDDD III FFFFFFF III X X
1970 // D D I F I X X
1971 // D D I F I X X
1972 // D D I FFFFF I X
1973 // D D I F I X X
1974 // D D I F I X X
1975 // DDDDDD III F III X X
1976
1977
1978 // VIF THU 20 NOV 2025 22:21:00 ** DIFIX ** DIAGNOSTICS PAGE 00060
1979
1980
1981 // VIF THU 20 NOV 2025 22:21:00 ** DIFIX ** GENERATED CODE PAGE 00061
1982
1983 int_4 _difix (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ ifix_, int_4 _p_ ldifix_, real_8 _p_ t_, int_4 _p_ ldt_, real_8 _p_
1984 tfix_, int_4 _p_ ldtfix_)
1985 { // ** body not listed **
2038 }
2039
2040 // VIF THU 20 NOV 2025 22:21:00 ** DINIWK ** DINIWK PAGE 00063
2041
2042 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
2043 // S S U U B B R R O O U U T I NN N E
2044 // S U U B B R R O O U U T I N N N E
2045 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
2046 // S U U B B R R O O U U T I N N N E
2047 // S S U U B B R R O O U U T I N NN E
2048 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
2049
2050 // DDDDDD III N N III W W K K
2051 // D D I NN N I W W K K
2052 // D D I N N N I W W K K
2053 // D D I N N N I W W W KKK
2054 // D D I N N N I W W W K K
2055 // D D I N NN I W W W K K
2056 // DDDDDD III N N III WW WW K K
2057
2058
2059 // VIF THU 20 NOV 2025 22:21:00 ** DINIWK ** DIAGNOSTICS PAGE 00064
2060
2061
2062 // VIF THU 20 NOV 2025 22:21:00 ** DINIWK ** GENERATED CODE PAGE 00065
2063
2064 int_4 _diniwk (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4
2065 _p_ liwork_, real_8 _p_ x_, int_4 _p_ ldx_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8 _p_ scld_, int_4 _p_ ldscld_,
2066 real_8 _p_ beta_, real_8 _p_ sclb_, real_8 _p_ sstol_, real_8 _p_ partol_, int_4 _p_ maxit_, real_8 _p_ taufac_, int_4
2067 _p_ job_, int_4 _p_ iprint_, int_4 _p_ lunerr_, int_4 _p_ lunrpt_, int_4 _p_ epsmai_, int_4 _p_ sstoli_, int_4 _p_
2068 partli_, int_4 _p_ maxiti_, int_4 _p_ taufci_, int_4 _p_ jobi_, int_4 _p_ iprini_, int_4 _p_ luneri_, int_4 _p_ lunrpi_
2069 , int_4 _p_ ssfi_, int_4 _p_ tti_, int_4 _p_ ldtti_, int_4 _p_ deltai_)
2070 { // ** body not listed **
2217 }
2218 // VIF THU 20 NOV 2025 22:21:00 ** DINIWK ** GENERATED CODE PAGE 00068
2219
2220
2221 // VIF THU 20 NOV 2025 22:21:00 ** DIWINF ** DIWINF PAGE 00069
2222
2223 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
2224 // S S U U B B R R O O U U T I NN N E
2225 // S U U B B R R O O U U T I N N N E
2226 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
2227 // S U U B B R R O O U U T I N N N E
2228 // S S U U B B R R O O U U T I N NN E
2229 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
2230
2231 // DDDDDD III W W III N N FFFFFFF
2232 // D D I W W I NN N F
2233 // D D I W W I N N N F
2234 // D D I W W W I N N N FFFFF
2235 // D D I W W W I N N N F
2236 // D D I W W W I N NN F
2237 // DDDDDD III WW WW III N N F
2238
2239
2240 // VIF THU 20 NOV 2025 22:21:00 ** DIWINF ** DIAGNOSTICS PAGE 00070
2241
2242
2243 // VIF THU 20 NOV 2025 22:21:00 ** DIWINF ** GENERATED CODE PAGE 00071
2244
2245 int_4 _diwinf (int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ msgbi_, int_4 _p_ msgdi_, int_4 _p_ ifix2i_, int_4
2246 _p_ istopi_, int_4 _p_ nnzwi_, int_4 _p_ nppi_, int_4 _p_ idfi_, int_4 _p_ jobi_, int_4 _p_ iprini_, int_4 _p_ luneri_,
2247 int_4 _p_ lunrpi_, int_4 _p_ nrowi_, int_4 _p_ ntoli_, int_4 _p_ netai_, int_4 _p_ maxiti_, int_4 _p_ niteri_, int_4
2248 _p_ nfevi_, int_4 _p_ njevi_, int_4 _p_ int2i_, int_4 _p_ iranki_, int_4 _p_ ldtti_, int_4 _p_ liwkmn_)
2249 { // ** body not listed **
2302 }
2303
2304 // VIF THU 20 NOV 2025 22:21:00 ** DJACCD ** DJACCD PAGE 00073
2305
2306 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
2307 // S S U U B B R R O O U U T I NN N E
2308 // S U U B B R R O O U U T I N N N E
2309 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
2310 // S U U B B R R O O U U T I N N N E
2311 // S S U U B B R R O O U U T I N NN E
2312 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
2313
2314 // DDDDDD J A CCCCC CCCCC DDDDDD
2315 // D D J A A C C C C D D
2316 // D D J A A C C D D
2317 // D D J A A C C D D
2318 // D D J J AAAAAAA C C D D
2319 // D D J J A A C C C C D D
2320 // DDDDDD JJJJJ A A CCCCC CCCCC DDDDDD
2321
2322
2323 // VIF THU 20 NOV 2025 22:21:00 ** DJACCD ** DIAGNOSTICS PAGE 00074
2324
2325
2326 // VIF THU 20 NOV 2025 22:21:00 ** DJACCD ** GENERATED CODE PAGE 00075
2327
2328 int_4 _djaccd (int_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
2329 x_, int_4 _p_ ldx_, real_8 _p_ delta_, real_8 _p_ xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8
2330 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_, real_8 _p_ ssf_, real_8 _p_ tt_, int_4 _p_ ldtt_, int_4 _p_ neta_,
2331 real_8 _p_ stp_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk3_, real_8 _p_ wrk6_, real_8 _p_ fjacb_, logical_4
2332 _p_ isodr_, real_8 _p_ fjacd_, int_4 _p_ nfev_, int_4 _p_ istop_)
2333 { // ** body not listed **
2546 }
2547
2548 // VIF THU 20 NOV 2025 22:21:00 ** DJACFD ** DJACFD PAGE 00080
2549
2550 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
2551 // S S U U B B R R O O U U T I NN N E
2552 // S U U B B R R O O U U T I N N N E
2553 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
2554 // S U U B B R R O O U U T I N N N E
2555 // S S U U B B R R O O U U T I N NN E
2556 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
2557
2558 // DDDDDD J A CCCCC FFFFFFF DDDDDD
2559 // D D J A A C C F D D
2560 // D D J A A C F D D
2561 // D D J A A C FFFFF D D
2562 // D D J J AAAAAAA C F D D
2563 // D D J J A A C C F D D
2564 // DDDDDD JJJJJ A A CCCCC F DDDDDD
2565
2566
2567 // VIF THU 20 NOV 2025 22:21:00 ** DJACFD ** DIAGNOSTICS PAGE 00081
2568
2569
2570 // VIF THU 20 NOV 2025 22:21:00 ** DJACFD ** GENERATED CODE PAGE 00082
2571
2572 int_4 _djacfd (int_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
2573 x_, int_4 _p_ ldx_, real_8 _p_ delta_, real_8 _p_ xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8
2574 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_, real_8 _p_ ssf_, real_8 _p_ tt_, int_4 _p_ ldtt_, int_4 _p_ neta_,
2575 real_8 _p_ fn_, real_8 _p_ stp_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk3_, real_8 _p_ wrk6_, real_8 _p_
2576 fjacb_, logical_4 _p_ isodr_, real_8 _p_ fjacd_, int_4 _p_ nfev_, int_4 _p_ istop_)
2577 { // ** body not listed **
2760 }
2761
2762 // VIF THU 20 NOV 2025 22:21:00 ** DJCKC ** DJCKC PAGE 00086
2763
2764 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
2765 // S S U U B B R R O O U U T I NN N E
2766 // S U U B B R R O O U U T I N N N E
2767 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
2768 // S U U B B R R O O U U T I N N N E
2769 // S S U U B B R R O O U U T I N NN E
2770 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
2771
2772 // DDDDDD J CCCCC K K CCCCC
2773 // D D J C C K K C C
2774 // D D J C K K C
2775 // D D J C KKK C
2776 // D D J J C K K C
2777 // D D J J C C K K C C
2778 // DDDDDD JJJJJ CCCCC K K CCCCC
2779
2780
2781 // VIF THU 20 NOV 2025 22:21:00 ** DJCKC ** DIAGNOSTICS PAGE 00087
2782
2783
2784 // VIF THU 20 NOV 2025 22:21:00 ** DJCKC ** GENERATED CODE PAGE 00088
2785
2786 int_4 _djckc (real_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
2787 xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8 _p_ eta_, real_8 _p_ tol_, int_4 _p_ nrow_,
2788 real_8 _p_ epsmac_, int_4 _p_ j_, int_4 _p_ lq_, real_8 _p_ hc_, logical_4 _p_ iswrtb_, real_8 _p_ fd_, real_8 _p_
2789 typj_, real_8 _p_ pvpstp_, real_8 _p_ stp0_, real_8 _p_ pv_, real_8 _p_ d_, real_8 _p_ diffj_, int_4 _p_ msg_, int_4
2790 _p_ istop_, int_4 _p_ nfev_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk6_)
2791 { // ** body not listed **
2901 }
2902
2903 // VIF THU 20 NOV 2025 22:21:00 ** DJCK ** DJCK PAGE 00091
2904
2905 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
2906 // S S U U B B R R O O U U T I NN N E
2907 // S U U B B R R O O U U T I N N N E
2908 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
2909 // S U U B B R R O O U U T I N N N E
2910 // S S U U B B R R O O U U T I N NN E
2911 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
2912
2913 // DDDDDD J CCCCC K K
2914 // D D J C C K K
2915 // D D J C K K
2916 // D D J C KKK
2917 // D D J J C K K
2918 // D D J J C C K K
2919 // DDDDDD JJJJJ CCCCC K K
2920
2921
2922 // VIF THU 20 NOV 2025 22:21:00 ** DJCK ** DIAGNOSTICS PAGE 00092
2923
2924
2925 // VIF THU 20 NOV 2025 22:21:00 ** DJCK ** GENERATED CODE PAGE 00093
2926
2927 int_4 _djck (int_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
2928 xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_,
2929 real_8 _p_ ssf_, real_8 _p_ tt_, int_4 _p_ ldtt_, real_8 _p_ eta_, int_4 _p_ neta_, int_4 _p_ ntol_, int_4 _p_ nrow_,
2930 logical_4 _p_ isodr_, real_8 _p_ epsmac_, real_8 _p_ pv0_, real_8 _p_ fjacb_, real_8 _p_ fjacd_, int_4 _p_ msgb_, int_4
2931 _p_ msgd_, real_8 _p_ diff_, int_4 _p_ istop_, int_4 _p_ nfev_, int_4 _p_ njev_, real_8 _p_ wrk1_, real_8 _p_ wrk2_,
2932 real_8 _p_ wrk6_)
2933 { // ** body not listed **
3088 }
3089
3090 // VIF THU 20 NOV 2025 22:21:00 ** DJCKF ** DJCKF PAGE 00097
3091
3092 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3093 // S S U U B B R R O O U U T I NN N E
3094 // S U U B B R R O O U U T I N N N E
3095 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3096 // S U U B B R R O O U U T I N N N E
3097 // S S U U B B R R O O U U T I N NN E
3098 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3099
3100 // DDDDDD J CCCCC K K FFFFFFF
3101 // D D J C C K K F
3102 // D D J C K K F
3103 // D D J C KKK FFFFF
3104 // D D J J C K K F
3105 // D D J J C C K K F
3106 // DDDDDD JJJJJ CCCCC K K F
3107
3108
3109 // VIF THU 20 NOV 2025 22:21:00 ** DJCKF ** DIAGNOSTICS PAGE 00098
3110
3111
3112 // VIF THU 20 NOV 2025 22:21:00 ** DJCKF ** GENERATED CODE PAGE 00099
3113
3114 int_4 _djckf (real_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
3115 xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8 _p_ eta_, real_8 _p_ tol_, int_4 _p_ nrow_, int_4
3116 _p_ j_, int_4 _p_ lq_, logical_4 _p_ iswrtb_, real_8 _p_ fd_, real_8 _p_ typj_, real_8 _p_ pvpstp_, real_8 _p_ stp0_,
3117 real_8 _p_ curve_, real_8 _p_ pv_, real_8 _p_ d_, real_8 _p_ diffj_, int_4 _p_ msg_, int_4 _p_ istop_, int_4 _p_ nfev_,
3118 real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk6_)
3119 { // ** body not listed **
3198 }
3199
3200 // VIF THU 20 NOV 2025 22:21:00 ** DJCKM ** DJCKM PAGE 00101
3201
3202 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3203 // S S U U B B R R O O U U T I NN N E
3204 // S U U B B R R O O U U T I N N N E
3205 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3206 // S U U B B R R O O U U T I N N N E
3207 // S S U U B B R R O O U U T I N NN E
3208 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3209
3210 // DDDDDD J CCCCC K K M M
3211 // D D J C C K K MM MM
3212 // D D J C K K M M M M
3213 // D D J C KKK M M M
3214 // D D J J C K K M M
3215 // D D J J C C K K M M
3216 // DDDDDD JJJJJ CCCCC K K M M
3217
3218
3219 // VIF THU 20 NOV 2025 22:21:00 ** DJCKM ** DIAGNOSTICS PAGE 00102
3220
3221
3222 // VIF THU 20 NOV 2025 22:21:00 ** DJCKM ** GENERATED CODE PAGE 00103
3223
3224 int_4 _djckm (real_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
3225 xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8 _p_ eta_, real_8 _p_ tol_, int_4 _p_ nrow_,
3226 real_8 _p_ epsmac_, int_4 _p_ j_, int_4 _p_ lq_, real_8 _p_ typj_, real_8 _p_ h0_, real_8 _p_ hc0_, logical_4 _p_
3227 iswrtb_, real_8 _p_ pv_, real_8 _p_ d_, real_8 _p_ diffj_, int_4 _p_ msg1_, int_4 _p_ msg_, int_4 _p_ istop_, int_4 _p_
3228 nfev_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk6_)
3229 { // ** body not listed **
3389 }
3390
3391 // VIF THU 20 NOV 2025 22:21:00 ** DJCKZ ** DJCKZ PAGE 00107
3392
3393 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3394 // S S U U B B R R O O U U T I NN N E
3395 // S U U B B R R O O U U T I N N N E
3396 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3397 // S U U B B R R O O U U T I N N N E
3398 // S S U U B B R R O O U U T I N NN E
3399 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3400
3401 // DDDDDD J CCCCC K K ZZZZZZZ
3402 // D D J C C K K Z
3403 // D D J C K K Z
3404 // D D J C KKK Z
3405 // D D J J C K K Z
3406 // D D J J C C K K *
3407 // DDDDDD JJJJJ CCCCC K K ZZZZZZZ
3408
3409
3410 // VIF THU 20 NOV 2025 22:21:00 ** DJCKZ ** DIAGNOSTICS PAGE 00108
3411
3412
3413 // VIF THU 20 NOV 2025 22:21:00 ** DJCKZ ** GENERATED CODE PAGE 00109
3414
3415 int_4 _djckz (real_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
3416 xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, int_4 _p_ nrow_, real_8 _p_ epsmac_, int_4 _p_ j_, int_4
3417 _p_ lq_, logical_4 _p_ iswrtb_, real_8 _p_ tol_, real_8 _p_ d_, real_8 _p_ fd_, real_8 _p_ typj_, real_8 _p_ pvpstp_,
3418 real_8 _p_ stp0_, real_8 _p_ pv_, real_8 _p_ diffj_, int_4 _p_ msg_, int_4 _p_ istop_, int_4 _p_ nfev_, real_8 _p_
3419 wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk6_)
3420 { // ** body not listed **
3488 }
3489
3490 // VIF THU 20 NOV 2025 22:21:00 ** DLUNO ** DLUNO PAGE 00111
3491
3492 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3493 // S S U U B B R R O O U U T I NN N E
3494 // S U U B B R R O O U U T I N N N E
3495 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3496 // S U U B B R R O O U U T I N N N E
3497 // S S U U B B R R O O U U T I N NN E
3498 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3499
3500 // DDDDDD L U U N N OOOOO
3501 // D D L U U NN N O O
3502 // D D L U U N N N O O
3503 // D D L U U N N N O O
3504 // D D L U U N N N O O
3505 // D D L U U N NN O O
3506 // DDDDDD LLLLLLL UUUUU N N OOOOO
3507
3508
3509 // VIF THU 20 NOV 2025 22:21:00 ** DLUNO ** DIAGNOSTICS PAGE 00112
3510
3511
3512 // VIF THU 20 NOV 2025 22:21:00 ** DLUNO ** GENERATED CODE PAGE 00113
3513
3514 int_4 _dluno (int_4 _p_ lun_, char _p_ fn_)
3515 { // ** body not listed **
3521 }
3522
3523 // VIF THU 20 NOV 2025 22:21:00 ** DLUNC ** DLUNC PAGE 00114
3524
3525 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3526 // S S U U B B R R O O U U T I NN N E
3527 // S U U B B R R O O U U T I N N N E
3528 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3529 // S U U B B R R O O U U T I N N N E
3530 // S S U U B B R R O O U U T I N NN E
3531 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3532
3533 // DDDDDD L U U N N CCCCC
3534 // D D L U U NN N C C
3535 // D D L U U N N N C
3536 // D D L U U N N N C
3537 // D D L U U N N N C
3538 // D D L U U N NN C C
3539 // DDDDDD LLLLLLL UUUUU N N CCCCC
3540
3541
3542 // VIF THU 20 NOV 2025 22:21:00 ** DLUNC ** DIAGNOSTICS PAGE 00115
3543
3544
3545 // VIF THU 20 NOV 2025 22:21:00 ** DLUNC ** GENERATED CODE PAGE 00116
3546
3547 int_4 _dlunc (int_4 _p_ lun_)
3548 { // ** body not listed **
3560 }
3561
3562 // VIF THU 20 NOV 2025 22:21:00 ** DMPREC ** DMPREC PAGE 00117
3563
3564 // RRRRRR EEEEEEE A L 88888
3565 // R R E A A L * * 8 8
3566 // R R E A A L * * 8 8
3567 // RRRRRR EEEEE A A L ******* 88888
3568 // R R E AAAAAAA L * * 8 8
3569 // R R E A A L * * 8 8
3570 // R R EEEEEEE A A LLLLLLL 88888
3571
3572 // FFFFFFF U U N N CCCCC TTTTTTT III OOOOO N N
3573 // F U U NN N C C T I O O NN N
3574 // F U U N N N C T I O O N N N
3575 // FFFFF U U N N N C T I O O N N N
3576 // F U U N N N C T I O O N N N
3577 // F U U N NN C C T I O O N NN
3578 // F UUUUU N N CCCCC T III OOOOO N N
3579
3580 // DDDDDD M M PPPPPP RRRRRR EEEEEEE CCCCC
3581 // D D MM MM P P R R E C C
3582 // D D M M M M P P R R E C
3583 // D D M M M PPPPPP RRRRRR EEEEE C
3584 // D D M M P R R E C
3585 // D D M M P R R E C C
3586 // DDDDDD M M P R R EEEEEEE CCCCC
3587
3588
3589 // VIF THU 20 NOV 2025 22:21:00 ** DMPREC ** DIAGNOSTICS PAGE 00118
3590
3591
3592 // VIF THU 20 NOV 2025 22:21:00 ** DMPREC ** GENERATED CODE PAGE 00119
3593
3594 real_8 _dmprec (void)
3595 { // ** body not listed **
3604 }
3605
3606 // VIF THU 20 NOV 2025 22:21:00 ** DODCHK ** DODCHK PAGE 00120
3607
3608 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3609 // S S U U B B R R O O U U T I NN N E
3610 // S U U B B R R O O U U T I N N N E
3611 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3612 // S U U B B R R O O U U T I N N N E
3613 // S S U U B B R R O O U U T I N NN E
3614 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3615
3616 // DDDDDD OOOOO DDDDDD CCCCC H H K K
3617 // D D O O D D C C H H K K
3618 // D D O O D D C H H K K
3619 // D D O O D D C HHHHHHH KKK
3620 // D D O O D D C H H K K
3621 // D D O O D D C C H H K K
3622 // DDDDDD OOOOO DDDDDD CCCCC H H K K
3623
3624
3625 // VIF THU 20 NOV 2025 22:21:00 ** DODCHK ** DIAGNOSTICS PAGE 00121
3626
3627
3628 // VIF THU 20 NOV 2025 22:21:00 ** DODCHK ** GENERATED CODE PAGE 00122
3629
3630 int_4 _dodchk (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, logical_4 _p_ isodr_, logical_4 _p_ anajac_,
3631 logical_4 _p_ implct_, int_4 _p_ ifixb_, int_4 _p_ ldx_, int_4 _p_ ldifx_, int_4 _p_ ldscld_, int_4 _p_ ldstpd_, int_4
3632 _p_ ldwe_, int_4 _p_ ld2we_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, int_4 _p_ ldy_, int_4 _p_ lwork_, int_4 _p_ lwkmn_,
3633 int_4 _p_ liwork_, int_4 _p_ liwkmn_, real_8 _p_ sclb_, real_8 _p_ scld_, real_8 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_
3634 info_)
3635 { // ** body not listed **
3773 }
3774
3775 // VIF THU 20 NOV 2025 22:21:00 ** DODCNT ** DODCNT PAGE 00125
3776
3777 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3778 // S S U U B B R R O O U U T I NN N E
3779 // S U U B B R R O O U U T I N N N E
3780 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3781 // S U U B B R R O O U U T I N N N E
3782 // S S U U B B R R O O U U T I N NN E
3783 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3784
3785 // DDDDDD OOOOO DDDDDD CCCCC N N TTTTTTT
3786 // D D O O D D C C NN N T
3787 // D D O O D D C N N N T
3788 // D D O O D D C N N N T
3789 // D D O O D D C N N N T
3790 // D D O O D D C C N NN T
3791 // DDDDDD OOOOO DDDDDD CCCCC N N T
3792
3793
3794 // VIF THU 20 NOV 2025 22:21:00 ** DODCNT ** DIAGNOSTICS PAGE 00126
3795
3796
3797 // VIF THU 20 NOV 2025 22:21:00 ** DODCNT ** GENERATED CODE PAGE 00127
3798
3799 int_4 _dodcnt (logical_4 _p_ short_, real_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8
3800 _p_ beta_, real_8 _p_ y_, int_4 _p_ ldy_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ we_, int_4 _p_ ldwe_, int_4 _p_
3801 ld2we_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, int_4
3802 _p_ job_, int_4 _p_ ndigit_, real_8 _p_ taufac_, real_8 _p_ sstol_, real_8 _p_ partol_, int_4 _p_ maxit_, int_4 _p_
3803 iprint_, int_4 _p_ lunerr_, int_4 _p_ lunrpt_, real_8 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_, real_8 _p_ sclb_,
3804 real_8 _p_ scld_, int_4 _p_ ldscld_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_
3805 info_)
3806 { // ** body not listed **
3940 }
3941
3942 // VIF THU 20 NOV 2025 22:21:00 ** DODDRV ** DODDRV PAGE 00130
3943
3944 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3945 // S S U U B B R R O O U U T I NN N E
3946 // S U U B B R R O O U U T I N N N E
3947 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3948 // S U U B B R R O O U U T I N N N E
3949 // S S U U B B R R O O U U T I N NN E
3950 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3951
3952 // DDDDDD OOOOO DDDDDD DDDDDD RRRRRR V V
3953 // D D O O D D D D R R V V
3954 // D D O O D D D D R R V V
3955 // D D O O D D D D RRRRRR V V
3956 // D D O O D D D D R R V V
3957 // D D O O D D D D R R V V
3958 // DDDDDD OOOOO DDDDDD DDDDDD R R V
3959
3960
3961 // VIF THU 20 NOV 2025 22:21:00 ** DODDRV ** DIAGNOSTICS PAGE 00131
3962
3963
3964 // VIF THU 20 NOV 2025 22:21:00 ** DODDRV ** GENERATED CODE PAGE 00132
3965
3966 int_4 _doddrv (logical_4 _p_ short_, logical_4 _p_ head_, logical_4 _p_ fstitr_, logical_4 _p_ prtpen_, int_4 (*_fcn)()
3967 , int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_ y_, int_4 _p_ ldy_, real_8 _p_
3968 x_, int_4 _p_ ldx_, real_8 _p_ we_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_
3969 ld2wd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, int_4 _p_ job_, int_4 _p_ ndigit_, real_8 _p_ taufac_,
3970 real_8 _p_ sstol_, real_8 _p_ partol_, int_4 _p_ maxit_, int_4 _p_ iprint_, int_4 _p_ lunerr_, int_4 _p_ lunrpt_,
3971 real_8 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_, real_8 _p_ sclb_, real_8 _p_ scld_, int_4 _p_ ldscld_, real_8
3972 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ maxit1_, real_8 _p_ tstimp_, int_4 _p_
3973 info_)
3974 { // ** body not listed **
4248 }
4249
4250 // VIF THU 20 NOV 2025 22:21:00 ** DODLM ** DODLM PAGE 00138
4251
4252 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
4253 // S S U U B B R R O O U U T I NN N E
4254 // S U U B B R R O O U U T I N N N E
4255 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
4256 // S U U B B R R O O U U T I N N N E
4257 // S S U U B B R R O O U U T I N NN E
4258 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
4259
4260 // DDDDDD OOOOO DDDDDD L M M
4261 // D D O O D D L MM MM
4262 // D D O O D D L M M M M
4263 // D D O O D D L M M M
4264 // D D O O D D L M M
4265 // D D O O D D L M M
4266 // DDDDDD OOOOO DDDDDD LLLLLLL M M
4267
4268
4269 // VIF THU 20 NOV 2025 22:21:00 ** DODLM ** DIAGNOSTICS PAGE 00139
4270
4271
4272 // VIF THU 20 NOV 2025 22:21:00 ** DODLM ** GENERATED CODE PAGE 00140
4273
4274 int_4 _dodlm (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ npp_, real_8 _p_ f_, real_8 _p_
4275 fjacb_, real_8 _p_ fjacd_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, real_8 _p_ ss_, real_8 _p_ tt_, int_4 _p_
4276 ldtt_, real_8 _p_ delta_, real_8 _p_ alpha2_, real_8 _p_ tau_, real_8 _p_ epsfcn_, logical_4 _p_ isodr_, real_8 _p_
4277 tfjacb_, real_8 _p_ omega_, real_8 _p_ u_, real_8 _p_ qraux_, int_4 _p_ jpvt_, real_8 _p_ s_, real_8 _p_ t_, int_4 _p_
4278 nlms_, real_8 _p_ rcond_, int_4 _p_ irank_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk3_, real_8 _p_ wrk4_,
4279 real_8 _p_ wrk5_, real_8 _p_ wrk_, int_4 _p_ lwrk_, int_4 _p_ istopc_)
4280 { // ** body not listed **
4409 }
4410
4411 // VIF THU 20 NOV 2025 22:21:00 ** DODMN ** DODMN PAGE 00143
4412
4413 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
4414 // S S U U B B R R O O U U T I NN N E
4415 // S U U B B R R O O U U T I N N N E
4416 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
4417 // S U U B B R R O O U U T I N N N E
4418 // S S U U B B R R O O U U T I N NN E
4419 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
4420
4421 // DDDDDD OOOOO DDDDDD M M N N
4422 // D D O O D D MM MM NN N
4423 // D D O O D D M M M M N N N
4424 // D D O O D D M M M N N N
4425 // D D O O D D M M N N N
4426 // D D O O D D M M N NN
4427 // DDDDDD OOOOO DDDDDD M M N N
4428
4429
4430 // VIF THU 20 NOV 2025 22:21:00 ** DODMN ** DIAGNOSTICS PAGE 00144
4431
4432
4433 // VIF THU 20 NOV 2025 22:21:00 ** DODMN ** GENERATED CODE PAGE 00145
4434
4435 int_4 _dodmn (logical_4 _p_ head_, logical_4 _p_ fstitr_, logical_4 _p_ prtpen_, int_4 (*_fcn)(), int_4 _p_ n_, int_4
4436 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ job_, real_8 _p_ beta_, real_8 _p_ y_, int_4 _p_ ldy_, real_8 _p_ x_,
4437 int_4 _p_ ldx_, real_8 _p_ we_, real_8 _p_ we1_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, real_8 _p_ wd_, int_4 _p_ ldwd_,
4438 int_4 _p_ ld2wd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8 _p_ betac_, real_8 _p_ betan_, real_8
4439 _p_ betas_, real_8 _p_ s_, real_8 _p_ delta_, real_8 _p_ deltan_, real_8 _p_ deltas_, real_8 _p_ t_, real_8 _p_ f_,
4440 real_8 _p_ fn_, real_8 _p_ fs_, real_8 _p_ fjacb_, int_4 _p_ msgb_, real_8 _p_ fjacd_, int_4 _p_ msgd_, real_8 _p_ ssf_
4441 , real_8 _p_ ss_, real_8 _p_ tt_, int_4 _p_ ldtt_, real_8 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_, real_8 _p_
4442 xplusd_, real_8 _p_ wrk_, int_4 _p_ lwrk_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_,
4443 int_4 _p_ info_)
4444 { // ** body not listed **
4984 }
4985
4986 // VIF THU 20 NOV 2025 22:21:00 ** DODPC1 ** DODPC1 PAGE 00156
4987
4988 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
4989 // S S U U B B R R O O U U T I NN N E
4990 // S U U B B R R O O U U T I N N N E
4991 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
4992 // S U U B B R R O O U U T I N N N E
4993 // S S U U B B R R O O U U T I N NN E
4994 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
4995
4996 // DDDDDD OOOOO DDDDDD PPPPPP CCCCC 1
4997 // D D O O D D P P C C 11
4998 // D D O O D D P P C 1 1
4999 // D D O O D D PPPPPP C 1
5000 // D D O O D D P C 1
5001 // D D O O D D P C C 1
5002 // DDDDDD OOOOO DDDDDD P CCCCC 11111
5003
5004
5005 // VIF THU 20 NOV 2025 22:21:00 ** DODPC1 ** DIAGNOSTICS PAGE 00157
5006
5007
5008 // VIF THU 20 NOV 2025 22:21:00 ** DODPC1 ** GENERATED CODE PAGE 00158
5009
5010 int_4 _dodpc1 (int_4 _p_ ipr_, int_4 _p_ lunrpt_, logical_4 _p_ anajac_, logical_4 _p_ cdjac_, logical_4 _p_ chkjac_,
5011 logical_4 _p_ initd_, logical_4 _p_ restrt_, logical_4 _p_ isodr_, logical_4 _p_ implct_, logical_4 _p_ dovcv_,
5012 logical_4 _p_ redoj_, int_4 _p_ msgb1_, int_4 _p_ msgb_, int_4 _p_ msgd1_, int_4 _p_ msgd_, int_4 _p_ n_, int_4 _p_ m_,
5013 int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ npp_, int_4 _p_ nnzw_, real_8 _p_ x_, int_4 _p_ ldx_, int_4 _p_ ifixx_, int_4
5014 _p_ ldifx_, real_8 _p_ delta_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, real_8 _p_ tt_, int_4 _p_ ldtt_,
5015 real_8 _p_ stpd_, int_4 _p_ ldstpd_, real_8 _p_ y_, int_4 _p_ ldy_, real_8 _p_ we_, int_4 _p_ ldwe_, int_4 _p_ ld2we_,
5016 real_8 _p_ pnlty_, real_8 _p_ beta_, int_4 _p_ ifixb_, real_8 _p_ ssf_, real_8 _p_ stpb_, int_4 _p_ job_, int_4 _p_
5017 neta_, real_8 _p_ taufac_, real_8 _p_ sstol_, real_8 _p_ partol_, int_4 _p_ maxit_, real_8 _p_ wss_, real_8 _p_ wssdel_
5018 , real_8 _p_ wsseps_)
5019 { // ** body not listed **
8249 }
8250
8251 // VIF THU 20 NOV 2025 22:21:00 ** DODPC2 ** DODPC2 PAGE 00221
8252
8253 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
8254 // S S U U B B R R O O U U T I NN N E
8255 // S U U B B R R O O U U T I N N N E
8256 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
8257 // S U U B B R R O O U U T I N N N E
8258 // S S U U B B R R O O U U T I N NN E
8259 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
8260
8261 // DDDDDD OOOOO DDDDDD PPPPPP CCCCC 22222
8262 // D D O O D D P P C C 2 2
8263 // D D O O D D P P C 2
8264 // D D O O D D PPPPPP C 22222
8265 // D D O O D D P C 2
8266 // D D O O D D P C C 2
8267 // DDDDDD OOOOO DDDDDD P CCCCC 2222222
8268
8269
8270 // VIF THU 20 NOV 2025 22:21:00 ** DODPC2 ** DIAGNOSTICS PAGE 00222
8271
8272
8273 // VIF THU 20 NOV 2025 22:21:00 ** DODPC2 ** GENERATED CODE PAGE 00223
8274
8275 int_4 _dodpc2 (int_4 _p_ ipr_, int_4 _p_ lunrpt_, logical_4 _p_ fstitr_, logical_4 _p_ implct_, logical_4 _p_ prtpen_,
8276 real_8 _p_ pnlty_, int_4 _p_ niter_, int_4 _p_ nfev_, real_8 _p_ wss_, real_8 _p_ actred_, real_8 _p_ prered_, real_8
8277 _p_ alpha_, real_8 _p_ tau_, real_8 _p_ pnorm_, int_4 _p_ np_, real_8 _p_ beta_)
8278 { // ** body not listed **
9264 }
9265
9266 // VIF THU 20 NOV 2025 22:21:00 ** DODPC3 ** DODPC3 PAGE 00243
9267
9268 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
9269 // S S U U B B R R O O U U T I NN N E
9270 // S U U B B R R O O U U T I N N N E
9271 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
9272 // S U U B B R R O O U U T I N N N E
9273 // S S U U B B R R O O U U T I N NN E
9274 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
9275
9276 // DDDDDD OOOOO DDDDDD PPPPPP CCCCC 33333
9277 // D D O O D D P P C C 3 3
9278 // D D O O D D P P C 3
9279 // D D O O D D PPPPPP C 33333
9280 // D D O O D D P C 3
9281 // D D O O D D P C C 3 3
9282 // DDDDDD OOOOO DDDDDD P CCCCC 33333
9283
9284
9285 // VIF THU 20 NOV 2025 22:21:00 ** DODPC3 ** DIAGNOSTICS PAGE 00244
9286
9287
9288 // VIF THU 20 NOV 2025 22:21:00 ** DODPC3 ** GENERATED CODE PAGE 00245
9289
9290 int_4 _dodpc3 (int_4 _p_ ipr_, int_4 _p_ lunrpt_, logical_4 _p_ isodr_, logical_4 _p_ implct_, logical_4 _p_ didvcv_,
9291 logical_4 _p_ dovcv_, logical_4 _p_ redoj_, logical_4 _p_ anajac_, int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_
9292 nq_, int_4 _p_ npp_, int_4 _p_ info_, int_4 _p_ niter_, int_4 _p_ nfev_, int_4 _p_ njev_, int_4 _p_ irank_, real_8 _p_
9293 rcond_, int_4 _p_ istop_, real_8 _p_ wss_, real_8 _p_ wssdel_, real_8 _p_ wsseps_, real_8 _p_ pnlty_, real_8 _p_ rvar_,
9294 int_4 _p_ idf_, real_8 _p_ beta_, real_8 _p_ sdbeta_, int_4 _p_ ifixb2_, real_8 _p_ f_, real_8 _p_ delta_)
9295 { // ** body not listed **
12364 }
12365
12366 // VIF THU 20 NOV 2025 22:21:00 ** DODPCR ** DODPCR PAGE 00305
12367
12368 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
12369 // S S U U B B R R O O U U T I NN N E
12370 // S U U B B R R O O U U T I N N N E
12371 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
12372 // S U U B B R R O O U U T I N N N E
12373 // S S U U B B R R O O U U T I N NN E
12374 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
12375
12376 // DDDDDD OOOOO DDDDDD PPPPPP CCCCC RRRRRR
12377 // D D O O D D P P C C R R
12378 // D D O O D D P P C R R
12379 // D D O O D D PPPPPP C RRRRRR
12380 // D D O O D D P C R R
12381 // D D O O D D P C C R R
12382 // DDDDDD OOOOO DDDDDD P CCCCC R R
12383
12384
12385 // VIF THU 20 NOV 2025 22:21:00 ** DODPCR ** DIAGNOSTICS PAGE 00306
12386
12387
12388 // VIF THU 20 NOV 2025 22:21:00 ** DODPCR ** GENERATED CODE PAGE 00307
12389
12390 int_4 _dodpcr (int_4 _p_ ipr_, int_4 _p_ lunrpt_, logical_4 _p_ head_, logical_4 _p_ prtpen_, logical_4 _p_ fstitr_,
12391 logical_4 _p_ didvcv_, int_4 _p_ iflag_, int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ npp_,
12392 int_4 _p_ nnzw_, int_4 _p_ msgb_, int_4 _p_ msgd_, real_8 _p_ beta_, real_8 _p_ y_, int_4 _p_ ldy_, real_8 _p_ x_,
12393 int_4 _p_ ldx_, real_8 _p_ delta_, real_8 _p_ we_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, real_8 _p_ wd_, int_4 _p_ ldwd_,
12394 int_4 _p_ ld2wd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, real_8 _p_ ssf_, real_8 _p_ tt_, int_4 _p_
12395 ldtt_, real_8 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_, int_4 _p_ job_, int_4 _p_ neta_, real_8 _p_ taufac_,
12396 real_8 _p_ sstol_, real_8 _p_ partol_, int_4 _p_ maxit_, real_8 _p_ wss_, real_8 _p_ rvar_, int_4 _p_ idf_, real_8 _p_
12397 sdbeta_, int_4 _p_ niter_, int_4 _p_ nfev_, int_4 _p_ njev_, real_8 _p_ actred_, real_8 _p_ prered_, real_8 _p_ tau_,
12398 real_8 _p_ pnorm_, real_8 _p_ alpha_, real_8 _p_ f_, real_8 _p_ rcond_, int_4 _p_ irank_, int_4 _p_ info_, int_4 _p_
12399 istop_)
12400 { // ** body not listed **
12557 }
12558
12559 // VIF THU 20 NOV 2025 22:21:00 ** DODPE1 ** DODPE1 PAGE 00311
12560
12561 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
12562 // S S U U B B R R O O U U T I NN N E
12563 // S U U B B R R O O U U T I N N N E
12564 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
12565 // S U U B B R R O O U U T I N N N E
12566 // S S U U B B R R O O U U T I N NN E
12567 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
12568
12569 // DDDDDD OOOOO DDDDDD PPPPPP EEEEEEE 1
12570 // D D O O D D P P E 11
12571 // D D O O D D P P E 1 1
12572 // D D O O D D PPPPPP EEEEE 1
12573 // D D O O D D P E 1
12574 // D D O O D D P E 1
12575 // DDDDDD OOOOO DDDDDD P EEEEEEE 11111
12576
12577
12578 // VIF THU 20 NOV 2025 22:21:00 ** DODPE1 ** DIAGNOSTICS PAGE 00312
12579
12580
12581 // VIF THU 20 NOV 2025 22:21:00 ** DODPE1 ** GENERATED CODE PAGE 00313
12582
12583 int_4 _dodpe1 (int_4 _p_ fileno_, int_4 _p_ d1_, int_4 _p_ d2_, int_4 _p_ d3_, int_4 _p_ d4_, int_4 _p_ d5_, int_4 _p_
12584 n_, int_4 _p_ m_, int_4 _p_ nq_, int_4 _p_ ldscld_, int_4 _p_ ldstpd_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, int_4 _p_
12585 ldwd_, int_4 _p_ ld2wd_, int_4 _p_ lwkmn_, int_4 _p_ liwkmn_)
12586 { // ** body not listed **
13367 }
13368
13369 // VIF THU 20 NOV 2025 22:21:00 ** DODPE2 ** DODPE2 PAGE 00329
13370
13371 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
13372 // S S U U B B R R O O U U T I NN N E
13373 // S U U B B R R O O U U T I N N N E
13374 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
13375 // S U U B B R R O O U U T I N N N E
13376 // S S U U B B R R O O U U T I N NN E
13377 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
13378
13379 // DDDDDD OOOOO DDDDDD PPPPPP EEEEEEE 22222
13380 // D D O O D D P P E 2 2
13381 // D D O O D D P P E 2
13382 // D D O O D D PPPPPP EEEEE 22222
13383 // D D O O D D P E 2
13384 // D D O O D D P E 2
13385 // DDDDDD OOOOO DDDDDD P EEEEEEE 2222222
13386
13387
13388 // VIF THU 20 NOV 2025 22:21:00 ** DODPE2 ** DIAGNOSTICS PAGE 00330
13389
13390
13391 // VIF THU 20 NOV 2025 22:21:00 ** DODPE2 ** GENERATED CODE PAGE 00331
13392
13393 int_4 _dodpe2 (int_4 _p_ fileno_, int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ fjacb_, real_8
13394 _p_ fjacd_, real_8 _p_ diff_, int_4 _p_ msgb1_, int_4 _p_ msgb_, logical_4 _p_ isodr_, int_4 _p_ msgd1_, int_4 _p_
13395 msgd_, real_8 _p_ xplusd_, int_4 _p_ nrow_, int_4 _p_ neta_, int_4 _p_ ntol_)
13396 { // ** body not listed **
14687 }
14688
14689 // VIF THU 20 NOV 2025 22:21:00 ** DODPE3 ** DODPE3 PAGE 00356
14690
14691 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
14692 // S S U U B B R R O O U U T I NN N E
14693 // S U U B B R R O O U U T I N N N E
14694 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
14695 // S U U B B R R O O U U T I N N N E
14696 // S S U U B B R R O O U U T I N NN E
14697 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
14698
14699 // DDDDDD OOOOO DDDDDD PPPPPP EEEEEEE 33333
14700 // D D O O D D P P E 3 3
14701 // D D O O D D P P E 3
14702 // D D O O D D PPPPPP EEEEE 33333
14703 // D D O O D D P E 3
14704 // D D O O D D P E 3 3
14705 // DDDDDD OOOOO DDDDDD P EEEEEEE 33333
14706
14707
14708 // VIF THU 20 NOV 2025 22:21:00 ** DODPE3 ** DIAGNOSTICS PAGE 00357
14709
14710
14711 // VIF THU 20 NOV 2025 22:21:00 ** DODPE3 ** GENERATED CODE PAGE 00358
14712
14713 int_4 _dodpe3 (int_4 _p_ fileno_, int_4 _p_ d2_, int_4 _p_ d3_)
14714 { // ** body not listed **
14866 }
14867 // VIF THU 20 NOV 2025 22:21:00 ** DODPE3 ** GENERATED CODE PAGE 00361
14868
14869
14870 // VIF THU 20 NOV 2025 22:21:00 ** DODPER ** DODPER PAGE 00362
14871
14872 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
14873 // S S U U B B R R O O U U T I NN N E
14874 // S U U B B R R O O U U T I N N N E
14875 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
14876 // S U U B B R R O O U U T I N N N E
14877 // S S U U B B R R O O U U T I N NN E
14878 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
14879
14880 // DDDDDD OOOOO DDDDDD PPPPPP EEEEEEE RRRRRR
14881 // D D O O D D P P E R R
14882 // D D O O D D P P E R R
14883 // D D O O D D PPPPPP EEEEE RRRRRR
14884 // D D O O D D P E R R
14885 // D D O O D D P E R R
14886 // DDDDDD OOOOO DDDDDD P EEEEEEE R R
14887
14888
14889 // VIF THU 20 NOV 2025 22:21:00 ** DODPER ** DIAGNOSTICS PAGE 00363
14890
14891
14892 // VIF THU 20 NOV 2025 22:21:00 ** DODPER ** GENERATED CODE PAGE 00364
14893
14894 int_4 _dodper (int_4 _p_ info_, int_4 _p_ lunerr_, logical_4 _p_ short_, int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_,
14895 int_4 _p_ nq_, int_4 _p_ ldscld_, int_4 _p_ ldstpd_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, int_4 _p_ ldwd_, int_4 _p_
14896 ld2wd_, int_4 _p_ lwkmn_, int_4 _p_ liwkmn_, real_8 _p_ fjacb_, real_8 _p_ fjacd_, real_8 _p_ diff_, int_4 _p_ msgb_,
14897 logical_4 _p_ isodr_, int_4 _p_ msgd_, real_8 _p_ xplusd_, int_4 _p_ nrow_, int_4 _p_ neta_, int_4 _p_ ntol_)
14898 { // ** body not listed **
15028 }
15029
15030 // VIF THU 20 NOV 2025 22:21:00 ** DODPHD ** DODPHD PAGE 00367
15031
15032 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
15033 // S S U U B B R R O O U U T I NN N E
15034 // S U U B B R R O O U U T I N N N E
15035 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
15036 // S U U B B R R O O U U T I N N N E
15037 // S S U U B B R R O O U U T I N NN E
15038 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
15039
15040 // DDDDDD OOOOO DDDDDD PPPPPP H H DDDDDD
15041 // D D O O D D P P H H D D
15042 // D D O O D D P P H H D D
15043 // D D O O D D PPPPPP HHHHHHH D D
15044 // D D O O D D P H H D D
15045 // D D O O D D P H H D D
15046 // DDDDDD OOOOO DDDDDD P H H DDDDDD
15047
15048
15049 // VIF THU 20 NOV 2025 22:21:00 ** DODPHD ** DIAGNOSTICS PAGE 00368
15050
15051
15052 // VIF THU 20 NOV 2025 22:21:00 ** DODPHD ** GENERATED CODE PAGE 00369
15053
15054 int_4 _dodphd (logical_4 _p_ head_, int_4 _p_ fileno_)
15055 { // ** body not listed **
15084 }
15085
15086 // VIF THU 20 NOV 2025 22:21:00 ** DODRC ** DODRC PAGE 00370
15087
15088 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
15089 // S S U U B B R R O O U U T I NN N E
15090 // S U U B B R R O O U U T I N N N E
15091 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
15092 // S U U B B R R O O U U T I N N N E
15093 // S S U U B B R R O O U U T I N NN E
15094 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
15095
15096 // DDDDDD OOOOO DDDDDD RRRRRR CCCCC
15097 // D D O O D D R R C C
15098 // D D O O D D R R C
15099 // D D O O D D RRRRRR C
15100 // D D O O D D R R C
15101 // D D O O D D R R C C
15102 // DDDDDD OOOOO DDDDDD R R CCCCC
15103
15104
15105 // VIF THU 20 NOV 2025 22:21:00 ** DODRC ** DIAGNOSTICS PAGE 00371
15106
15107
15108 // VIF THU 20 NOV 2025 22:21:00 ** DODRC ** GENERATED CODE PAGE 00372
15109
15110 int_4 _dodrc (real_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
15111 y_, int_4 _p_ ldy_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ we_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, real_8 _p_ wd_,
15112 int_4 _p_ ldwd_, int_4 _p_ ld2wd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, int_4 _p_ job_, int_4 _p_
15113 ndigit_, real_8 _p_ taufac_, real_8 _p_ sstol_, real_8 _p_ partol_, int_4 _p_ maxit_, int_4 _p_ iprint_, int_4 _p_
15114 lunerr_, int_4 _p_ lunrpt_, real_8 _p_ stpb_, real_8 _p_ stpd_, int_4 _p_ ldstpd_, real_8 _p_ sclb_, real_8 _p_ scld_,
15115 int_4 _p_ ldscld_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
15116 { // ** body not listed **
15152 }
15153
15154 // VIF THU 20 NOV 2025 22:21:00 ** DODR ** DODR PAGE 00373
15155
15156 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
15157 // S S U U B B R R O O U U T I NN N E
15158 // S U U B B R R O O U U T I N N N E
15159 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
15160 // S U U B B R R O O U U T I N N N E
15161 // S S U U B B R R O O U U T I N NN E
15162 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
15163
15164 // DDDDDD OOOOO DDDDDD RRRRRR
15165 // D D O O D D R R
15166 // D D O O D D R R
15167 // D D O O D D RRRRRR
15168 // D D O O D D R R
15169 // D D O O D D R R
15170 // DDDDDD OOOOO DDDDDD R R
15171
15172
15173 // VIF THU 20 NOV 2025 22:21:00 ** DODR ** DIAGNOSTICS PAGE 00374
15174
15175
15176 // VIF THU 20 NOV 2025 22:21:00 ** DODR ** GENERATED CODE PAGE 00375
15177
15178 int_4 _dodr (real_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
15179 y_, int_4 _p_ ldy_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ we_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, real_8 _p_ wd_,
15180 int_4 _p_ ldwd_, int_4 _p_ ld2wd_, int_4 _p_ job_, int_4 _p_ iprint_, int_4 _p_ lunerr_, int_4 _p_ lunrpt_, real_8 _p_
15181 work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
15182 { // ** body not listed **
15235 }
15236
15237 // VIF THU 20 NOV 2025 22:21:00 ** DODSTP ** DODSTP PAGE 00377
15238
15239 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
15240 // S S U U B B R R O O U U T I NN N E
15241 // S U U B B R R O O U U T I N N N E
15242 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
15243 // S U U B B R R O O U U T I N N N E
15244 // S S U U B B R R O O U U T I N NN E
15245 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
15246
15247 // DDDDDD OOOOO DDDDDD SSSSS TTTTTTT PPPPPP
15248 // D D O O D D S S T P P
15249 // D D O O D D S T P P
15250 // D D O O D D SSSSS T PPPPPP
15251 // D D O O D D S T P
15252 // D D O O D D S S T P
15253 // DDDDDD OOOOO DDDDDD SSSSS T P
15254
15255
15256 // VIF THU 20 NOV 2025 22:21:00 ** DODSTP ** DIAGNOSTICS PAGE 00378
15257
15258
15259 // VIF THU 20 NOV 2025 22:21:00 ** DODSTP ** GENERATED CODE PAGE 00379
15260
15261 int_4 _dodstp (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ npp_, real_8 _p_ f_, real_8 _p_
15262 fjacb_, real_8 _p_ fjacd_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, real_8 _p_ ss_, real_8 _p_ tt_, int_4 _p_
15263 ldtt_, real_8 _p_ delta_, real_8 _p_ alpha_, real_8 _p_ epsfcn_, logical_4 _p_ isodr_, real_8 _p_ tfjacb_, real_8 _p_
15264 omega_, real_8 _p_ u_, real_8 _p_ qraux_, int_4 _p_ kpvt_, real_8 _p_ s_, real_8 _p_ t_, real_8 _p_ phi_, int_4 _p_
15265 irank_, real_8 _p_ rcond_, logical_4 _p_ forvcv_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk3_, real_8 _p_
15266 wrk4_, real_8 _p_ wrk5_, real_8 _p_ wrk_, int_4 _p_ lwrk_, int_4 _p_ istopc_)
15267 { // ** body not listed **
15539 }
15540
15541 // VIF THU 20 NOV 2025 22:21:00 ** DODVCV ** DODVCV PAGE 00385
15542
15543 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
15544 // S S U U B B R R O O U U T I NN N E
15545 // S U U B B R R O O U U T I N N N E
15546 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
15547 // S U U B B R R O O U U T I N N N E
15548 // S S U U B B R R O O U U T I N NN E
15549 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
15550
15551 // DDDDDD OOOOO DDDDDD V V CCCCC V V
15552 // D D O O D D V V C C V V
15553 // D D O O D D V V C V V
15554 // D D O O D D V V C V V
15555 // D D O O D D V V C V V
15556 // D D O O D D V V C C V V
15557 // DDDDDD OOOOO DDDDDD V CCCCC V
15558
15559
15560 // VIF THU 20 NOV 2025 22:21:00 ** DODVCV ** DIAGNOSTICS PAGE 00386
15561
15562
15563 // VIF THU 20 NOV 2025 22:21:00 ** DODVCV ** GENERATED CODE PAGE 00387
15564
15565 int_4 _dodvcv (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ npp_, real_8 _p_ f_, real_8 _p_
15566 fjacb_, real_8 _p_ fjacd_, real_8 _p_ wd_, int_4 _p_ ldwd_, int_4 _p_ ld2wd_, real_8 _p_ ssf_, real_8 _p_ ss_, real_8
15567 _p_ tt_, int_4 _p_ ldtt_, real_8 _p_ delta_, real_8 _p_ epsfcn_, logical_4 _p_ isodr_, real_8 _p_ vcv_, real_8 _p_ sd_,
15568 real_8 _p_ wrk6_, real_8 _p_ omega_, real_8 _p_ u_, real_8 _p_ qraux_, int_4 _p_ jpvt_, real_8 _p_ s_, real_8 _p_ t_,
15569 int_4 _p_ irank_, real_8 _p_ rcond_, real_8 _p_ rss_, int_4 _p_ idf_, real_8 _p_ rvar_, int_4 _p_ ifixb_, real_8 _p_
15570 wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk3_, real_8 _p_ wrk4_, real_8 _p_ wrk5_, real_8 _p_ wrk_, int_4 _p_ lwrk_, int_4
15571 _p_ istopc_)
15572 { // ** body not listed **
15718 }
15719 // VIF THU 20 NOV 2025 22:21:00 ** DODVCV ** GENERATED CODE PAGE 00390
15720
15721
15722 // VIF THU 20 NOV 2025 22:21:00 ** DPACK ** DPACK PAGE 00391
15723
15724 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
15725 // S S U U B B R R O O U U T I NN N E
15726 // S U U B B R R O O U U T I N N N E
15727 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
15728 // S U U B B R R O O U U T I N N N E
15729 // S S U U B B R R O O U U T I N NN E
15730 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
15731
15732 // DDDDDD PPPPPP A CCCCC K K
15733 // D D P P A A C C K K
15734 // D D P P A A C K K
15735 // D D PPPPPP A A C KKK
15736 // D D P AAAAAAA C K K
15737 // D D P A A C C K K
15738 // DDDDDD P A A CCCCC K K
15739
15740
15741 // VIF THU 20 NOV 2025 22:21:00 ** DPACK ** DIAGNOSTICS PAGE 00392
15742
15743
15744 // VIF THU 20 NOV 2025 22:21:00 ** DPACK ** GENERATED CODE PAGE 00393
15745
15746 int_4 _dpack (int_4 _p_ n2_, int_4 _p_ n1_, real_8 _p_ v1_, real_8 _p_ v2_, int_4 _p_ ifix_)
15747 { // ** body not listed **
15765 }
15766
15767 // VIF THU 20 NOV 2025 22:21:00 ** DPPNML ** DPPNML PAGE 00394
15768
15769 // RRRRRR EEEEEEE A L 88888
15770 // R R E A A L * * 8 8
15771 // R R E A A L * * 8 8
15772 // RRRRRR EEEEE A A L ******* 88888
15773 // R R E AAAAAAA L * * 8 8
15774 // R R E A A L * * 8 8
15775 // R R EEEEEEE A A LLLLLLL 88888
15776
15777 // FFFFFFF U U N N CCCCC TTTTTTT III OOOOO N N
15778 // F U U NN N C C T I O O NN N
15779 // F U U N N N C T I O O N N N
15780 // FFFFF U U N N N C T I O O N N N
15781 // F U U N N N C T I O O N N N
15782 // F U U N NN C C T I O O N NN
15783 // F UUUUU N N CCCCC T III OOOOO N N
15784
15785 // DDDDDD PPPPPP PPPPPP N N M M L
15786 // D D P P P P NN N MM MM L
15787 // D D P P P P N N N M M M M L
15788 // D D PPPPPP PPPPPP N N N M M M L
15789 // D D P P N N N M M L
15790 // D D P P N NN M M L
15791 // DDDDDD P P N N M M LLLLLLL
15792
15793
15794 // VIF THU 20 NOV 2025 22:21:00 ** DPPNML ** DIAGNOSTICS PAGE 00395
15795
15796
15797 // VIF THU 20 NOV 2025 22:21:00 ** DPPNML ** GENERATED CODE PAGE 00396
15798
15799 real_8 _dppnml (real_8 _p_ p_)
15800 { // ** body not listed **
15930 }
15931
15932 // VIF THU 20 NOV 2025 22:21:00 ** DPPT ** DPPT PAGE 00399
15933
15934 // RRRRRR EEEEEEE A L 88888
15935 // R R E A A L * * 8 8
15936 // R R E A A L * * 8 8
15937 // RRRRRR EEEEE A A L ******* 88888
15938 // R R E AAAAAAA L * * 8 8
15939 // R R E A A L * * 8 8
15940 // R R EEEEEEE A A LLLLLLL 88888
15941
15942 // FFFFFFF U U N N CCCCC TTTTTTT III OOOOO N N
15943 // F U U NN N C C T I O O NN N
15944 // F U U N N N C T I O O N N N
15945 // FFFFF U U N N N C T I O O N N N
15946 // F U U N N N C T I O O N N N
15947 // F U U N NN C C T I O O N NN
15948 // F UUUUU N N CCCCC T III OOOOO N N
15949
15950 // DDDDDD PPPPPP PPPPPP TTTTTTT
15951 // D D P P P P T
15952 // D D P P P P T
15953 // D D PPPPPP PPPPPP T
15954 // D D P P T
15955 // D D P P T
15956 // DDDDDD P P T
15957
15958
15959 // VIF THU 20 NOV 2025 22:21:00 ** DPPT ** DIAGNOSTICS PAGE 00400
15960
15961
15962 // VIF THU 20 NOV 2025 22:21:00 ** DPPT ** GENERATED CODE PAGE 00401
15963
15964 real_8 _dppt (real_8 _p_ p_, int_4 _p_ idf_)
15965 { // ** body not listed **
16234 }
16235
16236 // VIF THU 20 NOV 2025 22:21:00 ** DPVB ** DPVB PAGE 00407
16237
16238 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16239 // S S U U B B R R O O U U T I NN N E
16240 // S U U B B R R O O U U T I N N N E
16241 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16242 // S U U B B R R O O U U T I N N N E
16243 // S S U U B B R R O O U U T I N NN E
16244 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16245
16246 // DDDDDD PPPPPP V V BBBBBB
16247 // D D P P V V B B
16248 // D D P P V V B B
16249 // D D PPPPPP V V BBBBBB
16250 // D D P V V B B
16251 // D D P V V B B
16252 // DDDDDD P V BBBBBB
16253
16254
16255 // VIF THU 20 NOV 2025 22:21:00 ** DPVB ** DIAGNOSTICS PAGE 00408
16256
16257
16258 // VIF THU 20 NOV 2025 22:21:00 ** DPVB ** GENERATED CODE PAGE 00409
16259
16260 int_4 _dpvb (int_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
16261 xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, int_4 _p_ nrow_, int_4 _p_ j_, int_4 _p_ lq_, real_8 _p_
16262 stp_, int_4 _p_ istop_, int_4 _p_ nfev_, real_8 _p_ pvb_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk6_)
16263 { // ** body not listed **
16282 }
16283
16284 // VIF THU 20 NOV 2025 22:21:00 ** DPVD ** DPVD PAGE 00410
16285
16286 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16287 // S S U U B B R R O O U U T I NN N E
16288 // S U U B B R R O O U U T I N N N E
16289 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16290 // S U U B B R R O O U U T I N N N E
16291 // S S U U B B R R O O U U T I N NN E
16292 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16293
16294 // DDDDDD PPPPPP V V DDDDDD
16295 // D D P P V V D D
16296 // D D P P V V D D
16297 // D D PPPPPP V V D D
16298 // D D P V V D D
16299 // D D P V V D D
16300 // DDDDDD P V DDDDDD
16301
16302
16303 // VIF THU 20 NOV 2025 22:21:00 ** DPVD ** DIAGNOSTICS PAGE 00411
16304
16305
16306 // VIF THU 20 NOV 2025 22:21:00 ** DPVD ** GENERATED CODE PAGE 00412
16307
16308 int_4 _dpvd (int_4 (*_fcn)(), int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, real_8 _p_ beta_, real_8 _p_
16309 xplusd_, int_4 _p_ ifixb_, int_4 _p_ ifixx_, int_4 _p_ ldifx_, int_4 _p_ nrow_, int_4 _p_ j_, int_4 _p_ lq_, real_8 _p_
16310 stp_, int_4 _p_ istop_, int_4 _p_ nfev_, real_8 _p_ pvd_, real_8 _p_ wrk1_, real_8 _p_ wrk2_, real_8 _p_ wrk6_)
16311 { // ** body not listed **
16330 }
16331
16332 // VIF THU 20 NOV 2025 22:21:00 ** DSCALE ** DSCALE PAGE 00413
16333
16334 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16335 // S S U U B B R R O O U U T I NN N E
16336 // S U U B B R R O O U U T I N N N E
16337 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16338 // S U U B B R R O O U U T I N N N E
16339 // S S U U B B R R O O U U T I N NN E
16340 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16341
16342 // DDDDDD SSSSS CCCCC A L EEEEEEE
16343 // D D S S C C A A L E
16344 // D D S C A A L E
16345 // D D SSSSS C A A L EEEEE
16346 // D D S C AAAAAAA L E
16347 // D D S S C C A A L E
16348 // DDDDDD SSSSS CCCCC A A LLLLLLL EEEEEEE
16349
16350
16351 // VIF THU 20 NOV 2025 22:21:00 ** DSCALE ** DIAGNOSTICS PAGE 00414
16352
16353
16354 // VIF THU 20 NOV 2025 22:21:00 ** DSCALE ** GENERATED CODE PAGE 00415
16355
16356 int_4 _dscale (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ scl_, int_4 _p_ ldscl_, real_8 _p_ t_, int_4 _p_ ldt_, real_8 _p_
16357 sclt_, int_4 _p_ ldsclt_)
16358 { // ** body not listed **
16415 }
16416
16417 // VIF THU 20 NOV 2025 22:21:00 ** DSCLB ** DSCLB PAGE 00417
16418
16419 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16420 // S S U U B B R R O O U U T I NN N E
16421 // S U U B B R R O O U U T I N N N E
16422 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16423 // S U U B B R R O O U U T I N N N E
16424 // S S U U B B R R O O U U T I N NN E
16425 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16426
16427 // DDDDDD SSSSS CCCCC L BBBBBB
16428 // D D S S C C L B B
16429 // D D S C L B B
16430 // D D SSSSS C L BBBBBB
16431 // D D S C L B B
16432 // D D S S C C L B B
16433 // DDDDDD SSSSS CCCCC LLLLLLL BBBBBB
16434
16435
16436 // VIF THU 20 NOV 2025 22:21:00 ** DSCLB ** DIAGNOSTICS PAGE 00418
16437
16438
16439 // VIF THU 20 NOV 2025 22:21:00 ** DSCLB ** GENERATED CODE PAGE 00419
16440
16441 int_4 _dsclb (int_4 _p_ np_, real_8 _p_ beta_, real_8 _p_ ssf_)
16442 { // ** body not listed **
16511 }
16512
16513 // VIF THU 20 NOV 2025 22:21:00 ** DSCLD ** DSCLD PAGE 00421
16514
16515 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16516 // S S U U B B R R O O U U T I NN N E
16517 // S U U B B R R O O U U T I N N N E
16518 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16519 // S U U B B R R O O U U T I N N N E
16520 // S S U U B B R R O O U U T I N NN E
16521 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16522
16523 // DDDDDD SSSSS CCCCC L DDDDDD
16524 // D D S S C C L D D
16525 // D D S C L D D
16526 // D D SSSSS C L D D
16527 // D D S C L D D
16528 // D D S S C C L D D
16529 // DDDDDD SSSSS CCCCC LLLLLLL DDDDDD
16530
16531
16532 // VIF THU 20 NOV 2025 22:21:00 ** DSCLD ** DIAGNOSTICS PAGE 00422
16533
16534
16535 // VIF THU 20 NOV 2025 22:21:00 ** DSCLD ** GENERATED CODE PAGE 00423
16536
16537 int_4 _dscld (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ tt_, int_4 _p_ ldtt_)
16538 { // ** body not listed **
16608 }
16609
16610 // VIF THU 20 NOV 2025 22:21:00 ** DSETN ** DSETN PAGE 00425
16611
16612 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16613 // S S U U B B R R O O U U T I NN N E
16614 // S U U B B R R O O U U T I N N N E
16615 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16616 // S U U B B R R O O U U T I N N N E
16617 // S S U U B B R R O O U U T I N NN E
16618 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16619
16620 // DDDDDD SSSSS EEEEEEE TTTTTTT N N
16621 // D D S S E T NN N
16622 // D D S E T N N N
16623 // D D SSSSS EEEEE T N N N
16624 // D D S E T N N N
16625 // D D S S E T N NN
16626 // DDDDDD SSSSS EEEEEEE T N N
16627
16628
16629 // VIF THU 20 NOV 2025 22:21:00 ** DSETN ** DIAGNOSTICS PAGE 00426
16630
16631
16632 // VIF THU 20 NOV 2025 22:21:00 ** DSETN ** GENERATED CODE PAGE 00427
16633
16634 int_4 _dsetn (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ x_, int_4 _p_ ldx_, int_4 _p_ nrow_)
16635 { // ** body not listed **
16654 }
16655
16656 // VIF THU 20 NOV 2025 22:21:00 ** DSOLVE ** DSOLVE PAGE 00428
16657
16658 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16659 // S S U U B B R R O O U U T I NN N E
16660 // S U U B B R R O O U U T I N N N E
16661 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16662 // S U U B B R R O O U U T I N N N E
16663 // S S U U B B R R O O U U T I N NN E
16664 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16665
16666 // DDDDDD SSSSS OOOOO L V V EEEEEEE
16667 // D D S S O O L V V E
16668 // D D S O O L V V E
16669 // D D SSSSS O O L V V EEEEE
16670 // D D S O O L V V E
16671 // D D S S O O L V V E
16672 // DDDDDD SSSSS OOOOO LLLLLLL V EEEEEEE
16673
16674
16675 // VIF THU 20 NOV 2025 22:21:00 ** DSOLVE ** DIAGNOSTICS PAGE 00429
16676
16677
16678 // VIF THU 20 NOV 2025 22:21:00 ** DSOLVE ** GENERATED CODE PAGE 00430
16679
16680 int_4 _dsolve (int_4 _p_ n_, real_8 _p_ t_, int_4 _p_ ldt_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ job_)
16681 { // ** body not listed **
16778 }
16779
16780 // VIF THU 20 NOV 2025 22:21:00 ** DUNPAC ** DUNPAC PAGE 00432
16781
16782 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16783 // S S U U B B R R O O U U T I NN N E
16784 // S U U B B R R O O U U T I N N N E
16785 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16786 // S U U B B R R O O U U T I N N N E
16787 // S S U U B B R R O O U U T I N NN E
16788 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16789
16790 // DDDDDD U U N N PPPPPP A CCCCC
16791 // D D U U NN N P P A A C C
16792 // D D U U N N N P P A A C
16793 // D D U U N N N PPPPPP A A C
16794 // D D U U N N N P AAAAAAA C
16795 // D D U U N NN P A A C C
16796 // DDDDDD UUUUU N N P A A CCCCC
16797
16798
16799 // VIF THU 20 NOV 2025 22:21:00 ** DUNPAC ** DIAGNOSTICS PAGE 00433
16800
16801
16802 // VIF THU 20 NOV 2025 22:21:00 ** DUNPAC ** GENERATED CODE PAGE 00434
16803
16804 int_4 _dunpac (int_4 _p_ n2_, real_8 _p_ v1_, real_8 _p_ v2_, int_4 _p_ ifix_)
16805 { // ** body not listed **
16823 }
16824
16825 // VIF THU 20 NOV 2025 22:21:00 ** DVEVTR ** DVEVTR PAGE 00435
16826
16827 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16828 // S S U U B B R R O O U U T I NN N E
16829 // S U U B B R R O O U U T I N N N E
16830 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16831 // S U U B B R R O O U U T I N N N E
16832 // S S U U B B R R O O U U T I N NN E
16833 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16834
16835 // DDDDDD V V EEEEEEE V V TTTTTTT RRRRRR
16836 // D D V V E V V T R R
16837 // D D V V E V V T R R
16838 // D D V V EEEEE V V T RRRRRR
16839 // D D V V E V V T R R
16840 // D D V V E V V T R R
16841 // DDDDDD V EEEEEEE V T R R
16842
16843
16844 // VIF THU 20 NOV 2025 22:21:00 ** DVEVTR ** DIAGNOSTICS PAGE 00436
16845
16846
16847 // VIF THU 20 NOV 2025 22:21:00 ** DVEVTR ** GENERATED CODE PAGE 00437
16848
16849 int_4 _dvevtr (int_4 _p_ m_, int_4 _p_ nq_, int_4 _p_ indx_, real_8 _p_ v_, int_4 _p_ ldv_, int_4 _p_ ld2v_, real_8 _p_
16850 e_, int_4 _p_ lde_, real_8 _p_ ve_, int_4 _p_ ldve_, int_4 _p_ ld2ve_, real_8 _p_ vev_, int_4 _p_ ldvev_, real_8 _p_
16851 wrk5_)
16852 { // ** body not listed **
16893 }
16894
16895 // VIF THU 20 NOV 2025 22:21:00 ** DWGHT ** DWGHT PAGE 00438
16896
16897 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
16898 // S S U U B B R R O O U U T I NN N E
16899 // S U U B B R R O O U U T I N N N E
16900 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
16901 // S U U B B R R O O U U T I N N N E
16902 // S S U U B B R R O O U U T I N NN E
16903 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
16904
16905 // DDDDDD W W GGGGG H H TTTTTTT
16906 // D D W W G G H H T
16907 // D D W W G H H T
16908 // D D W W W G GGGG HHHHHHH T
16909 // D D W W W G G H H T
16910 // D D W W W G G H H T
16911 // DDDDDD WW WW GGGGG H H T
16912
16913
16914 // VIF THU 20 NOV 2025 22:21:00 ** DWGHT ** DIAGNOSTICS PAGE 00439
16915
16916
16917 // VIF THU 20 NOV 2025 22:21:00 ** DWGHT ** GENERATED CODE PAGE 00440
16918
16919 int_4 _dwght (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ wt_, int_4 _p_ ldwt_, int_4 _p_ ld2wt_, real_8 _p_ t_, int_4 _p_
16920 ldt_, real_8 _p_ wtt_, int_4 _p_ ldwtt_)
16921 { // ** body not listed **
16999 }
17000
17001 // VIF THU 20 NOV 2025 22:21:00 ** DWINF ** DWINF PAGE 00442
17002
17003 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
17004 // S S U U B B R R O O U U T I NN N E
17005 // S U U B B R R O O U U T I N N N E
17006 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
17007 // S U U B B R R O O U U T I N N N E
17008 // S S U U B B R R O O U U T I N NN E
17009 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
17010
17011 // DDDDDD W W III N N FFFFFFF
17012 // D D W W I NN N F
17013 // D D W W I N N N F
17014 // D D W W W I N N N FFFFF
17015 // D D W W W I N N N F
17016 // D D W W W I N NN F
17017 // DDDDDD WW WW III N N F
17018
17019
17020 // VIF THU 20 NOV 2025 22:21:00 ** DWINF ** DIAGNOSTICS PAGE 00443
17021
17022
17023 // VIF THU 20 NOV 2025 22:21:00 ** DWINF ** GENERATED CODE PAGE 00444
17024
17025 int_4 _dwinf (int_4 _p_ n_, int_4 _p_ m_, int_4 _p_ np_, int_4 _p_ nq_, int_4 _p_ ldwe_, int_4 _p_ ld2we_, logical_4
17026 _p_ isodr_, int_4 _p_ deltai_, int_4 _p_ epsi_, int_4 _p_ xplusi_, int_4 _p_ fni_, int_4 _p_ sdi_, int_4 _p_ vcvi_,
17027 int_4 _p_ rvari_, int_4 _p_ wssi_, int_4 _p_ wssdei_, int_4 _p_ wssepi_, int_4 _p_ rcondi_, int_4 _p_ etai_, int_4 _p_
17028 olmavi_, int_4 _p_ taui_, int_4 _p_ alphai_, int_4 _p_ actrsi_, int_4 _p_ pnormi_, int_4 _p_ rnorsi_, int_4 _p_ prersi_
17029 , int_4 _p_ partli_, int_4 _p_ sstoli_, int_4 _p_ taufci_, int_4 _p_ epsmai_, int_4 _p_ beta0i_, int_4 _p_ betaci_,
17030 int_4 _p_ betasi_, int_4 _p_ betani_, int_4 _p_ si_, int_4 _p_ ssi_, int_4 _p_ ssfi_, int_4 _p_ qrauxi_, int_4 _p_ ui_,
17031 int_4 _p_ fsi_, int_4 _p_ fjacbi_, int_4 _p_ we1i_, int_4 _p_ diffi_, int_4 _p_ deltsi_, int_4 _p_ deltni_, int_4 _p_
17032 ti_, int_4 _p_ tti_, int_4 _p_ omegai_, int_4 _p_ fjacdi_, int_4 _p_ wrk1i_, int_4 _p_ wrk2i_, int_4 _p_ wrk3i_, int_4
17033 _p_ wrk4i_, int_4 _p_ wrk5i_, int_4 _p_ wrk6i_, int_4 _p_ wrk7i_, int_4 _p_ lwkmn_)
17034 { // ** body not listed **
17160 }
17161
17162 // VIF THU 20 NOV 2025 22:21:00 ** DXMY ** DXMY PAGE 00447
17163
17164 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
17165 // S S U U B B R R O O U U T I NN N E
17166 // S U U B B R R O O U U T I N N N E
17167 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
17168 // S U U B B R R O O U U T I N N N E
17169 // S S U U B B R R O O U U T I N NN E
17170 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
17171
17172 // DDDDDD X X M M Y Y
17173 // D D X X MM MM Y Y
17174 // D D X X M M M M Y Y
17175 // D D X M M M Y
17176 // D D X X M M Y
17177 // D D X X M M Y
17178 // DDDDDD X X M M Y
17179
17180
17181 // VIF THU 20 NOV 2025 22:21:00 ** DXMY ** DIAGNOSTICS PAGE 00448
17182
17183
17184 // VIF THU 20 NOV 2025 22:21:00 ** DXMY ** GENERATED CODE PAGE 00449
17185
17186 int_4 _dxmy (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ y_, int_4 _p_ ldy_, real_8 _p_ xmy_,
17187 int_4 _p_ ldxmy_)
17188 { // ** body not listed **
17197 }
17198
17199 // VIF THU 20 NOV 2025 22:21:00 ** DXPY ** DXPY PAGE 00450
17200
17201 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
17202 // S S U U B B R R O O U U T I NN N E
17203 // S U U B B R R O O U U T I N N N E
17204 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
17205 // S U U B B R R O O U U T I N N N E
17206 // S S U U B B R R O O U U T I N NN E
17207 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
17208
17209 // DDDDDD X X PPPPPP Y Y
17210 // D D X X P P Y Y
17211 // D D X X P P Y Y
17212 // D D X PPPPPP Y
17213 // D D X X P Y
17214 // D D X X P Y
17215 // DDDDDD X X P Y
17216
17217
17218 // VIF THU 20 NOV 2025 22:21:00 ** DXPY ** DIAGNOSTICS PAGE 00451
17219
17220
17221 // VIF THU 20 NOV 2025 22:21:00 ** DXPY ** GENERATED CODE PAGE 00452
17222
17223 int_4 _dxpy (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ y_, int_4 _p_ ldy_, real_8 _p_ xpy_,
17224 int_4 _p_ ldxpy_)
17225 { // ** body not listed **
17234 }
17235
17236 // VIF THU 20 NOV 2025 22:21:00 ** DZERO ** DZERO PAGE 00453
17237
17238 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
17239 // S S U U B B R R O O U U T I NN N E
17240 // S U U B B R R O O U U T I N N N E
17241 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
17242 // S U U B B R R O O U U T I N N N E
17243 // S S U U B B R R O O U U T I N NN E
17244 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
17245
17246 // DDDDDD ZZZZZZZ EEEEEEE RRRRRR OOOOO
17247 // D D Z E R R O O
17248 // D D Z E R R O O
17249 // D D Z EEEEE RRRRRR O O
17250 // D D Z E R R O O
17251 // D D * E R R O O
17252 // DDDDDD ZZZZZZZ EEEEEEE R R OOOOO
17253
17254
17255 // VIF THU 20 NOV 2025 22:21:00 ** DZERO ** DIAGNOSTICS PAGE 00454
17256
17257
17258 // VIF THU 20 NOV 2025 22:21:00 ** DZERO ** GENERATED CODE PAGE 00455
17259
17260 int_4 _dzero (int_4 _p_ n_, int_4 _p_ m_, real_8 _p_ a_, int_4 _p_ lda_)
17261 { // ** body not listed **
17284 }
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
