5#define NPLXMX (LIMPAR*(LIMPAR+6)+1)
40 if( strcmp(
optimize.chOptRtn,
"PHYM") == 0 )
59 else if( strcmp(
optimize.chOptRtn,
"SUBP") == 0 )
61 fprintf(
ioQQQ,
" Begin optimization with SUBPLEX\n" );
65 fprintf(
ioQQQ,
" Increase size of NPLXMX in parameter statements to handle this many variables.\n" );
66 fprintf(
ioQQQ,
" I need at least %5ld\n", need );
107 fprintf(
ioQQQ,
" SUBPLEX exceeding maximum iterations.\n"
108 " This can be reset with the OPTIMZE ITERATIONS command.\n" );
116 fprintf(
ioQQQ,
" trace return optimize_subplex:\n" );
119 fprintf(
ioQQQ,
" Values:" );
120 for( ii=1; ii <=
optimize.nvarxt[j]; ii++ )
122 fprintf(
ioQQQ,
"\n" );
131 called.lgTalkIsOK =
cpu.i().lgMPI_talk();
132 prt.lgFaintOn =
true;
136 fprintf(
ioQQQ,
" **************************************************\n" );
137 fprintf(
ioQQQ,
" **************************************************\n" );
138 fprintf(
ioQQQ,
" **************************************************\n" );
139 fprintf(
ioQQQ,
"\n Cloudy was called %4ld times.\n\n",
optimize.nOptimiz );
179 fprintf(
ioQQQ,
"The number of variable options on this line makes no sense to me.\n");
183 fprintf(
ioQQQ,
" Optimal command: %s\n",
input.chCardSav[np]);
184 fprintf(
ioQQQ,
" Smallest value:%10.2e Largest value:%10.2e Allowed range %10.2e to %10.2e\n",
190 if(
cpu.i().lgMaster() )
194 for( i=0; i <=
input.nSave; i++ )
195 fprintf( ioOptim,
"%s\n",
input.chCardSav[i]);
200 fprintf(
ioQQQ,
"\f" );
NORETURN void TotalInsanity(void)
#define DEBUG_ENTRY(funcname)
FILE * open_data(const char *fname, const char *mode, access_scheme scheme)
char chOptimFileName[INPUT_LINE_LENGTH]
void optimize_phymir(realnum[], const realnum[], long, chi2_type *, realnum)
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)
chi2_type optimize_func(const realnum param[], int grid_index=-1)