/home66/gary/public_html/cloudy/c08_branch/source/optimize.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2008 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */

#ifndef _OPTIMIZE_H_
#define _OPTIMIZE_H_


bool lgOptimize_do(void);

void vary_input(bool *lgLimOK);

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);

/*<NP->*/
#define VRSNEW  3.00f

/* >>chng 00 sep 12, had been continue.in but was confused with normal input files*/
#define CNTFILE "continue.pmr"

void optimize_phymir(realnum[],realnum[],long,realnum*,realnum);

double optimize_func(realnum param[]);

/* varypar.h */
#define LIMEXT  5L
#define LIMPAR  20L
#define NCOLLM  100L
#define NOBSLM  100L

EXTERN char chOptimFileName[INPUT_LINE_LENGTH]; 
EXTERN struct t_optimize {

        bool lgVaryOn,  
                lgNoVary;

        bool lgOptimr;

        realnum 
          varmax[LIMPAR], 
          varmin[LIMPAR], 

          vparm[LIMEXT][LIMPAR], 

          vincr[LIMPAR];

        long int nvarxt[LIMPAR], 
          nvfpnt[LIMPAR];

        realnum vpused[LIMPAR], 
          OptIncrm[LIMPAR], 
          varang[LIMPAR][2];

        long int nvary, 
          nparm,
          nRangeSet;
        bool lgVarOn;
        realnum ColDen_Obs[NCOLLM], 
          chColDen_error[NCOLLM];
        long int ion_ColDen[NCOLLM], 
          ncobs,
          ionTemp[NCOLLM];

        realnum xLineInt_Obs[NOBSLM], 
          xLineInt_error[NOBSLM];

        realnum temp_obs[NOBSLM], 
          temp_error[NOBSLM];
        char chTempWeight[NOBSLM][7];

        long int 
          nlobs, 
          nIterOptim,
          nTempObs;

        realnum wavelength [NOBSLM],
                errorwave[NOBSLM];

        FILE *ioOptim;
        realnum OptGlobalErr;

        /* counter for number of models in a grid - is one (not zero) for 
         * first simulation */
        long int nOptimiz;

        bool lgOptimFlow;
        realnum optint, 
          optier;
        long int nTrOpt;
        bool lgTrOpt;
        bool lgOptimize, 
          lgOptLin, 
          lgOptLum, 
          lgOptCol,
          lgOptTemp;

        bool lgParallel,
          lgOptCont;
        long useCPU;

        char chVarFmt[LIMPAR][FILENAME_PATH_LENGTH_2], 
          chColDen_label[NCOLLM][5], 
          chLineLabel[NOBSLM][5], 
          chTempLab[NOBSLM][5],
          chOptRtn[5];

        double SavGenericData[10];

        }       optimize;



#endif /* _OPTIMIZE_H_ */

---