<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.5"/>
<title>cloudy: /home/iurt/rpmbuild/BUILD/c13.01/source/optimize_do.cpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
$(document).ready(function() { searchBox.OnSelectItem(0); });
</script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td style="padding-left: 0.5em;">
<div id="projectname">cloudy
 <span id="projectnumber">trunk</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.5 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="namespaces.html"><span>Namespaces</span></a></li>
<li><a href="annotated.html"><span>Data Structures</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
<li>
<div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.png"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
</span>
</div>
</li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>Globals</span></a></li>
</ul>
</div>
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
<a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark"> </span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark"> </span>Data Structures</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark"> </span>Namespaces</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark"> </span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark"> </span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark"> </span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(6)"><span class="SelectionMark"> </span>Typedefs</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(7)"><span class="SelectionMark"> </span>Enumerations</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(8)"><span class="SelectionMark"> </span>Enumerator</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(9)"><span class="SelectionMark"> </span>Friends</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(10)"><span class="SelectionMark"> </span>Macros</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(11)"><span class="SelectionMark"> </span>Pages</a></div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_b2f33c71d4aa5e7af42a1ca61ff5af1b.html">source</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">optimize_do.cpp</div> </div>
</div><!--header-->
<div class="contents">
<a href="optimize__do_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"> * others. For conditions of distribution and use see copyright notice in license.txt */</span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">/*lgOptimize_do main driver for optimization runs*/</span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="preprocessor">#include "<a class="code" href="cddefines_8h.html">cddefines.h</a>"</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"><a class="line" href="optimize__do_8cpp.html#a5646b62b1e61c0258f70ffa48bc03e02"> 5</a></span> <span class="preprocessor">#define NPLXMX (LIMPAR*(LIMPAR+6)+1)</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="preprocessor"></span><span class="preprocessor">#include "<a class="code" href="input_8h.html">input.h</a>"</span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="preprocessor">#include "<a class="code" href="called_8h.html">called.h</a>"</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="preprocessor">#include "<a class="code" href="prt_8h.html">prt.h</a>"</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="preprocessor">#include "<a class="code" href="save_8h.html">save.h</a>"</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="preprocessor">#include "<a class="code" href="optimize_8h.html">optimize.h</a>"</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> </div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment">/* called by cdDrive, this returns false if things went ok, true for disaster */</span></div>
<div class="line"><a name="l00013"></a><span class="lineno"><a class="line" href="optimize__do_8cpp.html#a854163e18fe32cf189c36238a3565bdf"> 13</a></span> <span class="keywordtype">bool</span> <a class="code" href="optimize_8h.html#a854163e18fe32cf189c36238a3565bdf">lgOptimize_do</a>(<span class="keywordtype">void</span>)</div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> {</div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span>  <span class="keywordtype">long</span> <span class="keywordtype">int</span> i, </div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span>  iflag, </div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span>  ii, </div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span>  iworke[<a class="code" href="optimize__do_8cpp.html#a5646b62b1e61c0258f70ffa48bc03e02">NPLXMX</a>], </div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span>  j, </div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span>  mode, </div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span>  need, </div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span>  nfe, </div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  np;</div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span>  <a class="code" href="optimize_8h.html#a5846ad020ddaf38b0a183086a05196bf">chi2_type</a> fret;</div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span>  <a class="code" href="cddefines_8h.html#a9aeda5ccf905f7d5241c0fe35434c820">realnum</a> fx,</div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span>  ptem[<a class="code" href="optimize_8h.html#ae79a5460d3745f29d3a230d4c19ae6c4">LIMPAR</a>], </div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span>  delta[<a class="code" href="optimize_8h.html#ae79a5460d3745f29d3a230d4c19ae6c4">LIMPAR</a>], </div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span>  toler,</div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span>  worke[<a class="code" href="optimize__do_8cpp.html#a5646b62b1e61c0258f70ffa48bc03e02">NPLXMX</a>];</div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> </div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span>  <a class="code" href="cddefines_8h.html#ada0cfa6927ba3683adc971d17643f2e5">DEBUG_ENTRY</a>( <span class="stringliteral">"lgOptimize_do()"</span> );</div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> </div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  <span class="comment">/* main driver for optimization runs</span></div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment"> * Drives cloudy to optimize variables;*/</span></div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> </div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span>  <span class="comment">/* code originally written by R.F. Carswell, IOA Cambridge */</span></div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span> </div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  toler = (<a class="code" href="cddefines_8h.html#a9aeda5ccf905f7d5241c0fe35434c820">realnum</a>)log10(1. + <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#afbec1aab225ace1f64e3d4d63ce8cae5">OptGlobalErr</a>);</div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> </div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span>  <span class="keywordflow">if</span>( strcmp(<a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aa195ca512625c1a5d4a244ba64ab3054">chOptRtn</a>,<span class="stringliteral">"PHYM"</span>) == 0 )</div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  {</div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span>  <span class="comment">/* Phymir method */</span></div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  <span class="keywordflow">for</span>( j=0; j < <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; j++ )</div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span>  {</div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  ptem[j] = <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][j];</div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  delta[j] = <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a18a9d1fdfbe0db4a629c22cb53c20ab3">vincr</a>[j];</div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span>  }</div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  <span class="comment">/* >>chng 06 jan 02, fix uninitialized var problem detected by valgrind/purify, PvH */</span></div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="keywordflow">for</span>( j=<a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; j < <a class="code" href="optimize_8h.html#ae79a5460d3745f29d3a230d4c19ae6c4">LIMPAR</a>; j++ )</div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  {</div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  ptem[j] = -FLT_MAX;</div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span>  delta[j] = -FLT_MAX;</div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span>  }</div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  <a class="code" href="optimize_8h.html#a44cc64a3fe98c5d8403cf75cd7a8f7c5">optimize_phymir</a>(ptem,delta,<a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>,&fret,toler);</div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  <span class="keywordflow">for</span>( j=0; j < <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; j++ )</div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][j] = ptem[j];</div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  }</div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span> </div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span>( strcmp(<a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aa195ca512625c1a5d4a244ba64ab3054">chOptRtn</a>,<span class="stringliteral">"SUBP"</span>) == 0 )</div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  {</div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" Begin optimization with SUBPLEX\n"</span> );</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  need = 2*<a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a> + <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>*(<a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a> + 4) + 1;</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  <span class="keywordflow">if</span>( need > <a class="code" href="optimize__do_8cpp.html#a5646b62b1e61c0258f70ffa48bc03e02">NPLXMX</a> )</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  {</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" Increase size of NPLXMX in parameter statements to handle this many variables.\n"</span> );</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" I need at least %5ld\n"</span>, need );</div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  <a class="code" href="cddefines_8h.html#a937636e18a3a517d1aeb2c111adb5160">cdEXIT</a>(<a class="code" href="cddefines_8h.html#a73efe787c131b385070f25d18b7c9aa4">EXIT_FAILURE</a>);</div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  }</div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  <span class="keywordflow">for</span>( j=0; j < <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; j++ )</div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  ptem[j] = <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][j];</div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span> </div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <span class="comment">/* The subroutine SUBPLX input into cloudy 8/4/94.</span></div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span> <span class="comment"> * The program itself is very well commented.</span></div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span> <span class="comment"> * The mode must set to 0 for the default values.</span></div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span> <span class="comment"> * The switch iflag tells if the program terminated normally. */</span></div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  mode = 0;</div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span> </div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  <span class="comment">/* >>chng 97 dec 08, remove first arg, optimize_func, since not used in routines */</span></div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  <a class="code" href="optimize_8h.html#a7ae56509b76849056e201c7c1fb3cc9e">optimize_subplex</a>(</div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  <span class="comment">/* the number of parameters to vary */</span></div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>,</div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  <span class="comment">/* the relative error, single number */</span></div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  toler,</div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  <span class="comment">/* maximum number of function evaluations before giving up */</span></div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a7ce2a373fbad4b9c2981728f14ddcb98">nIterOptim</a>,</div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <span class="comment">/* mode of operation, we simply set to zero */</span></div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  mode,</div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  <span class="comment">/* the initial changes in the guessed best coefficients, typically 0.2 to 1 */</span></div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a18a9d1fdfbe0db4a629c22cb53c20ab3">vincr</a>,</div>
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  <span class="comment">/* a vector of nvary initial parameters that are the starting guesses for the parameters */</span></div>
<div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  ptem,</div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  <span class="comment">/* a realnum, this is simply ignored */</span></div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  &fx,</div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <span class="comment">/* another parameter that is simply ignored, a long int */</span></div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  &nfe,</div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <span class="comment">/* a realnum that is NPLXMX long, used for working space by the routine */</span></div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  <span class="comment">/* an array that is 20*26 + 1 elements long, used for working space */</span></div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  worke,</div>
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  <span class="comment">/* a long int that is NPLXMX long, used for working space by the routine */</span></div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  <span class="comment">/* an array that is 20*26 + 1 elements long, used for working space */</span></div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  iworke,</div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  <span class="comment">/* a long int - says what happened, if -1 then exceeded nIterOptim iteration */</span></div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  &iflag);</div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span> </div>
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  <span class="keywordflow">if</span>( iflag == -1 )</div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  {</div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" SUBPLEX exceeding maximum iterations.\n"</span></div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  <span class="stringliteral">" This can be reset with the OPTIMZE ITERATIONS command.\n"</span> );</div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  }</div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span> </div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  <span class="keywordflow">for</span>( j=0; j < <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; j++ )</div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][j] = ptem[j];</div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span> </div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  <span class="keywordflow">if</span>( <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a95aa7de0bd5cd5d7a1b802831f77e08c">lgOptimFlow</a> )</div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  {</div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" trace return optimize_subplex:\n"</span> );</div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  <span class="keywordflow">for</span>( j=0; j < <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; j++ )</div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  {</div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" Values:"</span> );</div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  <span class="keywordflow">for</span>( ii=1; ii <= <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a967997b4c7075c5ae851e43cd22be93e">nvarxt</a>[j]; ii++ )</div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" %.2e"</span>, <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[ii-1][j] );</div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">"\n"</span> );</div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  }</div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  }</div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  }</div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  <a class="code" href="cddefines_8h.html#ac55e0608660d9a2b8dc1cc14d7f59081">TotalInsanity</a>();</div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span> </div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  <span class="comment">// turn output back on for final model</span></div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  <a class="code" href="called_8cpp.html#aee8b9d066c56b63ee5195666abd874f3">called</a>.<a class="code" href="structt__called.html#a637dca49e1651bde0439e7da1e74d155">lgTalk</a> = <a class="code" href="cpu_8h.html#a0ef714d18b90c8132aa7a5dc5288e899">cpu</a>.<a class="code" href="classt__cpu.html#a6abb28da25284366b6d540d89a8d0a3e">i</a>().<a class="code" href="classt__cpu__i.html#a48ca43968ec45a03abd16dc42247e1d7">lgMPI_talk</a>();</div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  <a class="code" href="called_8cpp.html#aee8b9d066c56b63ee5195666abd874f3">called</a>.<a class="code" href="structt__called.html#ae59ebcae64a1a8001c877224f750061d">lgTalkIsOK</a> = <a class="code" href="cpu_8h.html#a0ef714d18b90c8132aa7a5dc5288e899">cpu</a>.<a class="code" href="classt__cpu.html#a6abb28da25284366b6d540d89a8d0a3e">i</a>().<a class="code" href="classt__cpu__i.html#a48ca43968ec45a03abd16dc42247e1d7">lgMPI_talk</a>();</div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <a class="code" href="prt_8cpp.html#ab8e4d65eaa169fb4e03c49c9a41982fd">prt</a>.<a class="code" href="structt__prt.html#affdbed5161b4fa78d6ca8d5692c61a92">lgFaintOn</a> = <span class="keyword">true</span>;</div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span> </div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  <span class="keywordflow">if</span>( <a class="code" href="called_8cpp.html#aee8b9d066c56b63ee5195666abd874f3">called</a>.<a class="code" href="structt__called.html#a637dca49e1651bde0439e7da1e74d155">lgTalk</a> )</div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  {</div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" **************************************************\n"</span> );</div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" **************************************************\n"</span> );</div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" **************************************************\n"</span> );</div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">"\n Cloudy was called %4ld times.\n\n"</span>, <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a392f241ceb4d737760f27c75d30de897">nOptimiz</a> );</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span> </div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  <span class="keywordflow">for</span>( i=0; i < <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; i++ )</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  {</div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  np = <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aeb3b28d6646039bf4960d15ce4774e24">nvfpnt</a>[i];</div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span> </div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  <span class="comment">/* now generate the actual command with parameter */</span></div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  <span class="keywordflow">if</span>( <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a967997b4c7075c5ae851e43cd22be93e">nvarxt</a>[i] == 1 )</div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  {</div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  sprintf( <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">chCardSav</a>[np], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aad77d2eede243752d675904e31677f2b">chVarFmt</a>[i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][i] );</div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  }</div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span> </div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span>( <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a967997b4c7075c5ae851e43cd22be93e">nvarxt</a>[i] == 2 )</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  {</div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  sprintf( <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">chCardSav</a>[np], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aad77d2eede243752d675904e31677f2b">chVarFmt</a>[i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][i],</div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[1][i]);</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  }</div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span> </div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span>( <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a967997b4c7075c5ae851e43cd22be93e">nvarxt</a>[i] == 3 )</div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  {</div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  sprintf( <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">chCardSav</a>[np], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aad77d2eede243752d675904e31677f2b">chVarFmt</a>[i], </div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[1][i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[2][i] );</div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  }</div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span> </div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span>( <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a967997b4c7075c5ae851e43cd22be93e">nvarxt</a>[i] == 4 )</div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  {</div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  sprintf( <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">chCardSav</a>[np], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aad77d2eede243752d675904e31677f2b">chVarFmt</a>[i], </div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[1][i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[2][i],</div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[3][i] );</div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  }</div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span> </div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span>( <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a967997b4c7075c5ae851e43cd22be93e">nvarxt</a>[i] == 5 )</div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span>  {</div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  sprintf( <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">chCardSav</a>[np], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aad77d2eede243752d675904e31677f2b">chVarFmt</a>[i], </div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[1][i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[2][i],</div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[3][i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[4][i]);</div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  }</div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span> </div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  {</div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  fprintf(<a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>,<span class="stringliteral">"The number of variable options on this line makes no sense to me.\n"</span>);</div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <a class="code" href="cddefines_8h.html#a937636e18a3a517d1aeb2c111adb5160">cdEXIT</a>(<a class="code" href="cddefines_8h.html#a73efe787c131b385070f25d18b7c9aa4">EXIT_FAILURE</a>);</div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  }</div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span> </div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" Optimal command: %s\n"</span>, <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">chCardSav</a>[np]);</div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" Smallest value:%10.2e Largest value:%10.2e Allowed range %10.2e to %10.2e\n"</span>, </div>
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a307cb548866e5040574e935d2d978708">varmin</a>[i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a1fbd0f95091419455715057b4d6f375a">varmax</a>[i], <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aea1f04b029ebc3bc1d4743dae1665378">varang</a>[i][0], </div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#aea1f04b029ebc3bc1d4743dae1665378">varang</a>[i][1] );</div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  }</div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  }</div>
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span> </div>
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  <span class="keywordflow">if</span>( <a class="code" href="cpu_8h.html#a0ef714d18b90c8132aa7a5dc5288e899">cpu</a>.<a class="code" href="classt__cpu.html#a6abb28da25284366b6d540d89a8d0a3e">i</a>().<a class="code" href="classt__cpu__i.html#a1afb0be4fe53ab062a4ff138c6659ca8">lgMaster</a>() )</div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  {</div>
<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <span class="comment">/* save optimal parameters */</span></div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  FILE *ioOptim = <a class="code" href="cpu_8cpp.html#a365356a809aa94ee905d59cb878bba55">open_data</a>( <a class="code" href="optimize_8cpp.html#aa92b1264d51fbc90d2fb54addea21c74">chOptimFileName</a>, <span class="stringliteral">"w"</span>, <a class="code" href="cpu_8h.html#a2fd14f46f7beb7281273fda3c33aaa0ea04cd741ad7d3f65d83e84b4671bf695d">AS_LOCAL_ONLY</a> );</div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  <span class="keywordflow">for</span>( i=0; i <= <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ab83580320beca85226e7c9e487f1e6d5">nSave</a>; i++ )</div>
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  fprintf( ioOptim, <span class="stringliteral">"%s\n"</span>, <a class="code" href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a>.<a class="code" href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">chCardSav</a>[i]);</div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  fclose( ioOptim );</div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span> </div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  <span class="comment">/* recalculate values in cloudy for the best fit, and print out</span></div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span> <span class="comment"> * all the information */</span></div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">"\f"</span> );</div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span> </div>
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <span class="keywordflow">for</span>( i=0; i < <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">nvary</a>; i++ )</div>
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  ptem[i] = <a class="code" href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a>.<a class="code" href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">vparm</a>[0][i];</div>
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span> </div>
<div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  (<span class="keywordtype">void</span>)<a class="code" href="optimize_8h.html#aca3e55fbbb82e671d1786aed89b6c782">optimize_func</a>(ptem);</div>
<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  }</div>
<div class="line"><a name="l00207"></a><span class="lineno"> 207</span> </div>
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  <span class="keywordflow">return</span> <a class="code" href="cddefines_8cpp.html#ad63889bfb7ec3be761d72998e3eb439e">lgAbort</a>;</div>
<div class="line"><a name="l00209"></a><span class="lineno"> 209</span> }</div>
<div class="ttc" id="structt__input_html_ab83580320beca85226e7c9e487f1e6d5"><div class="ttname"><a href="structt__input.html#ab83580320beca85226e7c9e487f1e6d5">t_input::nSave</a></div><div class="ttdeci">long int nSave</div><div class="ttdef"><b>Definition:</b> <a href="input_8h_source.html#l00046">input.h:46</a></div></div>
<div class="ttc" id="cpu_8cpp_html_a365356a809aa94ee905d59cb878bba55"><div class="ttname"><a href="cpu_8cpp.html#a365356a809aa94ee905d59cb878bba55">open_data</a></div><div class="ttdeci">FILE * open_data(const char *fname, const char *mode, access_scheme scheme)</div><div class="ttdef"><b>Definition:</b> <a href="cpu_8cpp_source.html#l00616">cpu.cpp:616</a></div></div>
<div class="ttc" id="structt__optimize_html_a1fbd0f95091419455715057b4d6f375a"><div class="ttname"><a href="structt__optimize.html#a1fbd0f95091419455715057b4d6f375a">t_optimize::varmax</a></div><div class="ttdeci">realnum varmax[LIMPAR]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00183">optimize.h:183</a></div></div>
<div class="ttc" id="input_8h_html"><div class="ttname"><a href="input_8h.html">input.h</a></div></div>
<div class="ttc" id="cddefines_8h_html_ac55e0608660d9a2b8dc1cc14d7f59081"><div class="ttname"><a href="cddefines_8h.html#ac55e0608660d9a2b8dc1cc14d7f59081">TotalInsanity</a></div><div class="ttdeci">NORETURN void TotalInsanity(void)</div><div class="ttdef"><b>Definition:</b> <a href="service_8cpp_source.html#l00886">service.cpp:886</a></div></div>
<div class="ttc" id="input_8cpp_html_a79ef9860d57798efe65aeb1e7ff8ea01"><div class="ttname"><a href="input_8cpp.html#a79ef9860d57798efe65aeb1e7ff8ea01">input</a></div><div class="ttdeci">t_input input</div><div class="ttdef"><b>Definition:</b> <a href="input_8cpp_source.html#l00012">input.cpp:12</a></div></div>
<div class="ttc" id="structt__optimize_html_aeb3b28d6646039bf4960d15ce4774e24"><div class="ttname"><a href="structt__optimize.html#aeb3b28d6646039bf4960d15ce4774e24">t_optimize::nvfpnt</a></div><div class="ttdeci">long int nvfpnt[LIMPAR]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00194">optimize.h:194</a></div></div>
<div class="ttc" id="classt__cpu_html_a6abb28da25284366b6d540d89a8d0a3e"><div class="ttname"><a href="classt__cpu.html#a6abb28da25284366b6d540d89a8d0a3e">t_cpu::i</a></div><div class="ttdeci">t_cpu_i & i()</div><div class="ttdef"><b>Definition:</b> <a href="cpu_8h_source.html#l00334">cpu.h:334</a></div></div>
<div class="ttc" id="optimize_8h_html_aca3e55fbbb82e671d1786aed89b6c782"><div class="ttname"><a href="optimize_8h.html#aca3e55fbbb82e671d1786aed89b6c782">optimize_func</a></div><div class="ttdeci">chi2_type optimize_func(const realnum param[], int grid_index=-1)</div><div class="ttdef"><b>Definition:</b> <a href="optimize__func_8cpp_source.html#l00020">optimize_func.cpp:20</a></div></div>
<div class="ttc" id="structt__optimize_html_a392f241ceb4d737760f27c75d30de897"><div class="ttname"><a href="structt__optimize.html#a392f241ceb4d737760f27c75d30de897">t_optimize::nOptimiz</a></div><div class="ttdeci">long int nOptimiz</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00246">optimize.h:246</a></div></div>
<div class="ttc" id="structt__optimize_html_aea1f04b029ebc3bc1d4743dae1665378"><div class="ttname"><a href="structt__optimize.html#aea1f04b029ebc3bc1d4743dae1665378">t_optimize::varang</a></div><div class="ttdeci">realnum varang[LIMPAR][2]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00197">optimize.h:197</a></div></div>
<div class="ttc" id="optimize_8h_html_a44cc64a3fe98c5d8403cf75cd7a8f7c5"><div class="ttname"><a href="optimize_8h.html#a44cc64a3fe98c5d8403cf75cd7a8f7c5">optimize_phymir</a></div><div class="ttdeci">void optimize_phymir(realnum[], const realnum[], long, chi2_type *, realnum)</div></div>
<div class="ttc" id="structt__optimize_html_aad77d2eede243752d675904e31677f2b"><div class="ttname"><a href="structt__optimize.html#aad77d2eede243752d675904e31677f2b">t_optimize::chVarFmt</a></div><div class="ttdeci">char chVarFmt[LIMPAR][FILENAME_PATH_LENGTH_2]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00263">optimize.h:263</a></div></div>
<div class="ttc" id="save_8h_html"><div class="ttname"><a href="save_8h.html">save.h</a></div></div>
<div class="ttc" id="cddefines_8cpp_html_a94bc49f6569c17375ac7fbf6a66f6dd0"><div class="ttname"><a href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a></div><div class="ttdeci">FILE * ioQQQ</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8cpp_source.html#l00007">cddefines.cpp:7</a></div></div>
<div class="ttc" id="structt__optimize_html_a61fab5b0f2116787dde71f83a323e98f"><div class="ttname"><a href="structt__optimize.html#a61fab5b0f2116787dde71f83a323e98f">t_optimize::vparm</a></div><div class="ttdeci">realnum vparm[LIMEXT][LIMPAR]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00188">optimize.h:188</a></div></div>
<div class="ttc" id="structt__called_html_a637dca49e1651bde0439e7da1e74d155"><div class="ttname"><a href="structt__called.html#a637dca49e1651bde0439e7da1e74d155">t_called::lgTalk</a></div><div class="ttdeci">bool lgTalk</div><div class="ttdef"><b>Definition:</b> <a href="called_8h_source.html#l00012">called.h:12</a></div></div>
<div class="ttc" id="optimize_8h_html_a7ae56509b76849056e201c7c1fb3cc9e"><div class="ttname"><a href="optimize_8h.html#a7ae56509b76849056e201c7c1fb3cc9e">optimize_subplex</a></div><div class="ttdeci">void optimize_subplex(long int n, double tol, long int maxnfe, long int mode, realnum scale[], realnum x[], realnum *fx, long int *nfe, realnum work[], long int iwork[], long int *iflag)</div><div class="ttdef"><b>Definition:</b> <a href="optimize__subplx_8cpp_source.html#l00064">optimize_subplx.cpp:64</a></div></div>
<div class="ttc" id="structt__optimize_html_a7ce2a373fbad4b9c2981728f14ddcb98"><div class="ttname"><a href="structt__optimize.html#a7ce2a373fbad4b9c2981728f14ddcb98">t_optimize::nIterOptim</a></div><div class="ttdeci">long int nIterOptim</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00205">optimize.h:205</a></div></div>
<div class="ttc" id="optimize__do_8cpp_html_a5646b62b1e61c0258f70ffa48bc03e02"><div class="ttname"><a href="optimize__do_8cpp.html#a5646b62b1e61c0258f70ffa48bc03e02">NPLXMX</a></div><div class="ttdeci">#define NPLXMX</div><div class="ttdef"><b>Definition:</b> <a href="optimize__do_8cpp_source.html#l00005">optimize_do.cpp:5</a></div></div>
<div class="ttc" id="called_8h_html"><div class="ttname"><a href="called_8h.html">called.h</a></div></div>
<div class="ttc" id="optimize_8h_html"><div class="ttname"><a href="optimize_8h.html">optimize.h</a></div></div>
<div class="ttc" id="cddefines_8h_html_a9aeda5ccf905f7d5241c0fe35434c820"><div class="ttname"><a href="cddefines_8h.html#a9aeda5ccf905f7d5241c0fe35434c820">realnum</a></div><div class="ttdeci">float realnum</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00107">cddefines.h:107</a></div></div>
<div class="ttc" id="cddefines_8h_html_a73efe787c131b385070f25d18b7c9aa4"><div class="ttname"><a href="cddefines_8h.html#a73efe787c131b385070f25d18b7c9aa4">EXIT_FAILURE</a></div><div class="ttdeci">#define EXIT_FAILURE</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00144">cddefines.h:144</a></div></div>
<div class="ttc" id="classt__cpu__i_html_a1afb0be4fe53ab062a4ff138c6659ca8"><div class="ttname"><a href="classt__cpu__i.html#a1afb0be4fe53ab062a4ff138c6659ca8">t_cpu_i::lgMaster</a></div><div class="ttdeci">bool lgMaster() const </div><div class="ttdef"><b>Definition:</b> <a href="cpu_8h_source.html#l00314">cpu.h:314</a></div></div>
<div class="ttc" id="structt__prt_html_affdbed5161b4fa78d6ca8d5692c61a92"><div class="ttname"><a href="structt__prt.html#affdbed5161b4fa78d6ca8d5692c61a92">t_prt::lgFaintOn</a></div><div class="ttdeci">bool lgFaintOn</div><div class="ttdef"><b>Definition:</b> <a href="prt_8h_source.html#l00202">prt.h:202</a></div></div>
<div class="ttc" id="cddefines_8h_html_a937636e18a3a517d1aeb2c111adb5160"><div class="ttname"><a href="cddefines_8h.html#a937636e18a3a517d1aeb2c111adb5160">cdEXIT</a></div><div class="ttdeci">#define cdEXIT(FAIL)</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00438">cddefines.h:438</a></div></div>
<div class="ttc" id="classt__cpu__i_html_a48ca43968ec45a03abd16dc42247e1d7"><div class="ttname"><a href="classt__cpu__i.html#a48ca43968ec45a03abd16dc42247e1d7">t_cpu_i::lgMPI_talk</a></div><div class="ttdeci">bool lgMPI_talk() const </div><div class="ttdef"><b>Definition:</b> <a href="cpu_8h_source.html#l00315">cpu.h:315</a></div></div>
<div class="ttc" id="optimize_8h_html_ae79a5460d3745f29d3a230d4c19ae6c4"><div class="ttname"><a href="optimize_8h.html#ae79a5460d3745f29d3a230d4c19ae6c4">LIMPAR</a></div><div class="ttdeci">const long LIMPAR</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00061">optimize.h:61</a></div></div>
<div class="ttc" id="optimize_8cpp_html_ad9c1c492528e865ade1aeae90a429054"><div class="ttname"><a href="optimize_8cpp.html#ad9c1c492528e865ade1aeae90a429054">optimize</a></div><div class="ttdeci">t_optimize optimize</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8cpp_source.html#l00005">optimize.cpp:5</a></div></div>
<div class="ttc" id="structt__optimize_html_a18a9d1fdfbe0db4a629c22cb53c20ab3"><div class="ttname"><a href="structt__optimize.html#a18a9d1fdfbe0db4a629c22cb53c20ab3">t_optimize::vincr</a></div><div class="ttdeci">realnum vincr[LIMPAR]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00191">optimize.h:191</a></div></div>
<div class="ttc" id="prt_8cpp_html_ab8e4d65eaa169fb4e03c49c9a41982fd"><div class="ttname"><a href="prt_8cpp.html#ab8e4d65eaa169fb4e03c49c9a41982fd">prt</a></div><div class="ttdeci">t_prt prt</div><div class="ttdef"><b>Definition:</b> <a href="prt_8cpp_source.html#l00010">prt.cpp:10</a></div></div>
<div class="ttc" id="structt__optimize_html_aa195ca512625c1a5d4a244ba64ab3054"><div class="ttname"><a href="structt__optimize.html#aa195ca512625c1a5d4a244ba64ab3054">t_optimize::chOptRtn</a></div><div class="ttdeci">char chOptRtn[5]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00264">optimize.h:264</a></div></div>
<div class="ttc" id="optimize_8h_html_a5846ad020ddaf38b0a183086a05196bf"><div class="ttname"><a href="optimize_8h.html#a5846ad020ddaf38b0a183086a05196bf">chi2_type</a></div><div class="ttdeci">double chi2_type</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00049">optimize.h:49</a></div></div>
<div class="ttc" id="cddefines_8h_html"><div class="ttname"><a href="cddefines_8h.html">cddefines.h</a></div></div>
<div class="ttc" id="cddefines_8h_html_ada0cfa6927ba3683adc971d17643f2e5"><div class="ttname"><a href="cddefines_8h.html#ada0cfa6927ba3683adc971d17643f2e5">DEBUG_ENTRY</a></div><div class="ttdeci">#define DEBUG_ENTRY(funcname)</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00688">cddefines.h:688</a></div></div>
<div class="ttc" id="cpu_8h_html_a2fd14f46f7beb7281273fda3c33aaa0ea04cd741ad7d3f65d83e84b4671bf695d"><div class="ttname"><a href="cpu_8h.html#a2fd14f46f7beb7281273fda3c33aaa0ea04cd741ad7d3f65d83e84b4671bf695d">AS_LOCAL_ONLY</a></div><div class="ttdef"><b>Definition:</b> <a href="cpu_8h_source.html#l00195">cpu.h:195</a></div></div>
<div class="ttc" id="structt__optimize_html_a307cb548866e5040574e935d2d978708"><div class="ttname"><a href="structt__optimize.html#a307cb548866e5040574e935d2d978708">t_optimize::varmin</a></div><div class="ttdeci">realnum varmin[LIMPAR]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00184">optimize.h:184</a></div></div>
<div class="ttc" id="optimize_8h_html_a854163e18fe32cf189c36238a3565bdf"><div class="ttname"><a href="optimize_8h.html#a854163e18fe32cf189c36238a3565bdf">lgOptimize_do</a></div><div class="ttdeci">bool lgOptimize_do(void)</div><div class="ttdef"><b>Definition:</b> <a href="optimize__do_8cpp_source.html#l00013">optimize_do.cpp:13</a></div></div>
<div class="ttc" id="structt__input_html_ae1282ce92a528597572b8ad0f8216382"><div class="ttname"><a href="structt__input.html#ae1282ce92a528597572b8ad0f8216382">t_input::chCardSav</a></div><div class="ttdeci">char chCardSav[NKRD][INPUT_LINE_LENGTH]</div><div class="ttdef"><b>Definition:</b> <a href="input_8h_source.html#l00032">input.h:32</a></div></div>
<div class="ttc" id="prt_8h_html"><div class="ttname"><a href="prt_8h.html">prt.h</a></div></div>
<div class="ttc" id="structt__optimize_html_afbec1aab225ace1f64e3d4d63ce8cae5"><div class="ttname"><a href="structt__optimize.html#afbec1aab225ace1f64e3d4d63ce8cae5">t_optimize::OptGlobalErr</a></div><div class="ttdeci">realnum OptGlobalErr</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00243">optimize.h:243</a></div></div>
<div class="ttc" id="cpu_8h_html_a0ef714d18b90c8132aa7a5dc5288e899"><div class="ttname"><a href="cpu_8h.html#a0ef714d18b90c8132aa7a5dc5288e899">cpu</a></div><div class="ttdeci">static t_cpu cpu</div><div class="ttdef"><b>Definition:</b> <a href="cpu_8h_source.html#l00342">cpu.h:342</a></div></div>
<div class="ttc" id="structt__optimize_html_a95aa7de0bd5cd5d7a1b802831f77e08c"><div class="ttname"><a href="structt__optimize.html#a95aa7de0bd5cd5d7a1b802831f77e08c">t_optimize::lgOptimFlow</a></div><div class="ttdeci">bool lgOptimFlow</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00248">optimize.h:248</a></div></div>
<div class="ttc" id="optimize_8cpp_html_aa92b1264d51fbc90d2fb54addea21c74"><div class="ttname"><a href="optimize_8cpp.html#aa92b1264d51fbc90d2fb54addea21c74">chOptimFileName</a></div><div class="ttdeci">char chOptimFileName[INPUT_LINE_LENGTH]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8cpp_source.html#l00006">optimize.cpp:6</a></div></div>
<div class="ttc" id="structt__optimize_html_a967997b4c7075c5ae851e43cd22be93e"><div class="ttname"><a href="structt__optimize.html#a967997b4c7075c5ae851e43cd22be93e">t_optimize::nvarxt</a></div><div class="ttdeci">long int nvarxt[LIMPAR]</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00194">optimize.h:194</a></div></div>
<div class="ttc" id="called_8cpp_html_aee8b9d066c56b63ee5195666abd874f3"><div class="ttname"><a href="called_8cpp.html#aee8b9d066c56b63ee5195666abd874f3">called</a></div><div class="ttdeci">t_called called</div><div class="ttdef"><b>Definition:</b> <a href="called_8cpp_source.html#l00005">called.cpp:5</a></div></div>
<div class="ttc" id="structt__optimize_html_a5be42c3cf9b64d79ddea853de6d5f10f"><div class="ttname"><a href="structt__optimize.html#a5be42c3cf9b64d79ddea853de6d5f10f">t_optimize::nvary</a></div><div class="ttdeci">long int nvary</div><div class="ttdef"><b>Definition:</b> <a href="optimize_8h_source.html#l00200">optimize.h:200</a></div></div>
<div class="ttc" id="structt__called_html_ae59ebcae64a1a8001c877224f750061d"><div class="ttname"><a href="structt__called.html#ae59ebcae64a1a8001c877224f750061d">t_called::lgTalkIsOK</a></div><div class="ttdeci">bool lgTalkIsOK</div><div class="ttdef"><b>Definition:</b> <a href="called_8h_source.html#l00023">called.h:23</a></div></div>
<div class="ttc" id="cddefines_8cpp_html_ad63889bfb7ec3be761d72998e3eb439e"><div class="ttname"><a href="cddefines_8cpp.html#ad63889bfb7ec3be761d72998e3eb439e">lgAbort</a></div><div class="ttdeci">bool lgAbort</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8cpp_source.html#l00010">cddefines.cpp:10</a></div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated on Sat Oct 19 2013 10:43:49 for cloudy by  <a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/>
</a> 1.8.5
</small></address>
</body>
</html>
