luksan/luksan.h

   1 #ifndef LUKSAN_H
   2 #define LUKSAN_H
   3
   4 #include <nlopt.h>
   5 #include <nlopt-util.h>
   6
   7 #ifdef __cplusplus
   8 extern "C"
   9 {
  10 #endif /* __cplusplus */
  11
  12 nlopt_result luksan_plis(int n, nlopt_func f, void *f_data,
  13                   const double *lb, const double *ub, /* bounds */
  14                   double *x, /* in: initial guess, out: minimizer */
  15                   double *minf,
  16                   nlopt_stopping *stop);
  17
  18 nlopt_result luksan_plip(int n, nlopt_func f, void *f_data,
  19                          const double *lb, const double *ub, /* bounds */
  20                          double *x, /* in: initial guess, out: minimizer */
  21                          double *minf,
  22                          nlopt_stopping *stop,
  23                          int method);
  24
  25 /*****************************  internal routines *************************/
  26
  27 /* mssubs.c: */
  28 void luksan_mxdcmd__(int *n, int *m, double *a,
  29                      double *x, double *alf, double *y, double *z__);
  30 void luksan_mxdrcb__(int *n, int *m, double *a,
  31                      double *b, double *u, double *v, double *x, int *
  32                      ix, int *job);
  33 void luksan_mxdrcf__(int *n, int *m, double *a,
  34                      double *b, double *u, double *v, double *x, int *
  35                      ix, int *job);
  36 void luksan_mxdrmm__(int *n, int *m, double *a,
  37                      double *x, double *y);
  38 void luksan_mxdrsu__(int *n, int *m, double *a,
  39                      double *b, double *u);
  40 void luksan_mxucop__(int *n, double *x, double *y,
  41                      int *ix, int *job);
  42 void luksan_mxudir__(int *n, double *a, double *x,
  43                      double *y, double *z__, int *ix, int *job);
  44 void luksan_mxuneg__(int *n, double *x, double *y,
  45                      int *ix, int *job);
  46 void luksan_mxuzer__(int *n, double *x, int *ix,
  47                      int *job);
  48 void luksan_mxvcop__(int *n, double *x, double *y);
  49 void luksan_mxvdif__(int *n, double *x, double *y,
  50                      double *z__);
  51 void luksan_mxvneg__(int *n, double *x, double *y);
  52 void luksan_mxvscl__(int *n, double *a, double *x,
  53                      double *y);
  54 void luksan_mxvset__(int *n, double *a, double *x);
  55 double luksan_mxudot__(int *n, double *x, double *y, int *ix,
  56                        int *job);
  57 double luksan_mxvdot__(int *n, double *x, double *y);
  58 void luksan_mxvdir__(int *n, double *a, double *x,
  59                      double *y, double *z__);
  60 void luksan_mxdcmu__(int *n, int *m, double *a,
  61                      double *alf, double *x, double *y);
  62 void luksan_mxvlin__(int *n, double *a, double *x,
  63                      double *b, double *y, double *z__);
  64 void luksan_mxdcmv__(int *n, int *m, double *a,
  65                      double *alf, double *x, double *u, double *bet,
  66                      double *y, double *v);
  67 void luksan_mxvsav__(int *n, double *x, double *y);
  68 void luksan_mxvine__(int *n, int *ix);
  69 double luksan_mxvmax__(int *n, double *x);
  70
  71 /* pssubs.c: */
  72 void luksan_pcbs04__(int *nf, double *x, int *ix,
  73                      double *xl, double *xu, double *eps9, int *kbf);
  74 void luksan_ps1l01__(double *r__, double *rp,
  75                      double *f, double *fo, double *fp, double *p,
  76                      double *po, double *pp, double *minf, double *fmax,
  77                      double *rmin, double *rmax, double *tols, double *
  78                      tolp, double *par1, double *par2, int *kd, int *ld,
  79                      int *nit, int *kit, int *nred, int *mred, int *
  80                      maxst, int *iest, int *inits, int *iters, int *kters,
  81                      int *mes, int *isys);
  82 void luksan_pulsp3__(int *n, int *m, int *mf,
  83                      double *xm, double *gr, double *xo, double *go,
  84                      double *r__, double *po, double *sig, int *iterh,
  85                      int *met3);
  86 void luksan_pulvp3__(int *n, int *m, double *xm,
  87                      double *xr, double *gr, double *s, double *so,
  88                      double *xo, double *go, double *r__, double *po,
  89                      double *sig, int *iterh, int *met2, int *met3,
  90                      int *met5);
  91 void luksan_pyadc0__(int *nf, int *n, double *x,
  92                      int *ix, double *xl, double *xu, int *inew);
  93 void luksan_pyfut1__(int *n, double *f, double *
  94                      fo, double *umax, double *gmax, double *dmax__,
  95                      double *tolx, double *tolf, double *tolb, double *
  96                      tolg, int *kd, int *nit, int *kit, int *mit, int *
  97                      nfv, int *mfv, int *nfg, int *mfg, int *ntesx,
  98                      int *mtesx, int *ntesf, int *mtesf, int *ites,
  99                      int *ires1, int *ires2, int *irest, int *iters,
 100                      int *iterm);
 101 void luksan_pyrmc0__(int *nf, int *n, int *ix,
 102                      double *g, double *eps8, double *umax, double *gmax,
 103                      double *rmax, int *iold, int *irest);
 104 void luksan_pytrcd__(int *nf, double *x, int *ix,
 105                      double *xo, double *g, double *go, double *r__,
 106                      double *f, double *fo, double *p, double *po,
 107                      double *dmax__, int *kbf, int *kd, int *ld, int *
 108                      iters);
 109 void luksan_pytrcg__(int *nf, int *n, int *ix,
 110                      double *g, double *umax, double *gmax, int *kbf,
 111                      int *iold);
 112 void luksan_pytrcs__(int *nf, double *x, int *ix,
 113                      double *xo, double *xl, double *xu, double *g,
 114                      double *go, double *s, double *ro, double *fp,
 115                      double *fo, double *f, double *po, double *p,
 116                      double *rmax, double *eta9, int *kbf);
 117 void luksan_pnint1__(double *rl, double *ru, double *fl,
 118                      double *fu, double *pl, double *pu, double *r__,
 119                      int *mode, int *mtyp, int *merr);
 120
 121 /* Common Block Declarations */
 122 typedef struct {
 123      int nres, ndec, nin, nit;
 124      /* int nfv;   -- now stored in stop->nevals */
 125      int nfg, nfh;
 126 } stat_common;
 127
 128 /* number of double variables that can be stored in scratch memory
 129    ... it's >= 2007, and this is in the context of scientific computation,
 130    so assume that at least 10M are available, and that sizeof(double)==8 */
 131 #define MEMAVAIL 1310720
 132
 133 #ifdef __cplusplus
 134 }  /* extern "C" */
 135 #endif /* __cplusplus */
 136
 137 #endif /* LUKSAN_H */