--- /dev/null
+/* Copyright (c) 2007-2008 Massachusetts Institute of Technology
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining
+ * a copy of this software and associated documentation files (the
+ * "Software"), to deal in the Software without restriction, including
+ * without limitation the rights to use, copy, modify, merge, publish,
+ * distribute, sublicense, and/or sell copies of the Software, and to
+ * permit persons to whom the Software is furnished to do so, subject to
+ * the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be
+ * included in all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+ * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+ * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+ * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+ * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+ * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+ * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+ */
+
+#include <stdlib.h>
+
+#include "nlopt.h"
+#include "config.h"
+
+/*-----------------------------------------------------------------------*/
+/* wrappers around f77 procedures */
+
+typedef void (*nlopt_f77_func)(double *val, const int *n, const double *x,
+ double *gradient, const int *need_gradient,
+ void *func_data);
+
+typedef struct {
+ nlopt_f77_func f;
+ void *f_data;
+} f77_func_data;
+
+static double f77_func_wrap(int n, const double *x, double *grad, void *data)
+{
+ f77_func_data *d = (f77_func_data *) data;
+ double val;
+ int need_gradient = grad != 0;
+ d->f(&val, &n, x, grad, &need_gradient, d->f_data);
+ return val;
+}
+
+/*-----------------------------------------------------------------------*/
+/* rather than trying to detect the Fortran name-mangling scheme with
+ autoconf, we just include wrappers with all common name-mangling
+ schemes ... this avoids problems and also allows us to work with
+ multiple Fortran compilers on the same machine .
+
+ Note that our Fortran function names do not contain underscores;
+ otherwise, we would need to deal with the additional headache that
+ g77 appends two underscores in that case. */
+
+#ifndef WINDOWS_F77_MANGLING
+
+/* name + underscore is by far the most common (gfortran, g77, Intel, ...) */
+# define F77(a, A) a ## _
+# include "f77funcs.h"
+
+/* AIX and HPUX use just the lower-case name */
+# undef F77
+# define F77(a, A) a
+# include "f77funcs.h"
+
+/* old Cray UNICOS used just the upper-case name */
+# undef F77
+# define F77(a, A) A
+# include "f77funcs.h"
+
+#else /* WINDOWS_F77_MANGLING */
+
+/* Various mangling conventions common (?) under Windows. */
+
+/* name + underscore for gfortran, g77, ...? */
+# define F77(a, A) a ## _
+# include "f77funcs.h"
+
+/* Digital/Compaq/HP Visual Fortran, Intel Fortran. stdcall attribute
+ is apparently required to adjust for calling conventions (callee
+ pops stack in stdcall). See also:
+ http://msdn.microsoft.com/library/en-us/vccore98/html/_core_mixed.2d.language_programming.3a_.overview.asp
+*/
+# undef F77
+# if defined(__GNUC__)
+# define F77(a, A) __attribute__((stdcall)) A
+# elif defined(_MSC_VER) || defined(_ICC) || defined(_STDCALL_SUPPORTED)
+# define F77(a, A) __stdcall A
+# else
+# define F77(a, A) A /* oh well */
+# endif
+# include "f77funcs.h"
+
+#endif /* WINDOWS_F77_MANGLING */
--- /dev/null
+/* Copyright (c) 2007-2008 Massachusetts Institute of Technology
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining
+ * a copy of this software and associated documentation files (the
+ * "Software"), to deal in the Software without restriction, including
+ * without limitation the rights to use, copy, modify, merge, publish,
+ * distribute, sublicense, and/or sell copies of the Software, and to
+ * permit persons to whom the Software is furnished to do so, subject to
+ * the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be
+ * included in all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+ * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+ * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+ * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+ * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+ * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+ * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+ */
+
+/* Fortran API wrappers, using the F77 macro defined in f77api.c.
+ This header file is #included one or more times from f77api.c
+ in order to define verions of the Fortran API for various compilers.
+
+ In homage to Fortran 77, we stick with 6-character subroutine names.
+ The return value of the function is converted into the first argument
+ of the subroutine. */
+
+/* nlopt_minimize_constrained */
+void F77(nloptc,NLOPTC)(int *info,
+ const int *algorithm,
+ const int *n, nlopt_f77_func f, void *f_data,
+ const int *m, nlopt_f77_func fc,
+ char *fc_data, char *fc_second_datum,
+ const double *lb, const double *ub,
+ double *x, double *minf,
+ const double *minf_max,
+ const double *ftol_rel, const double *ftol_abs,
+ const double *xtol_rel, const double *xtol_abs,
+ const int *have_xtol_abs,
+ const int *maxeval, const double *maxtime)
+{
+ f77_func_data d, *dc;
+ int i;
+
+ d.f = f; d.f_data = f_data;
+ if (*m < 0) { *info = NLOPT_INVALID_ARGS; return; }
+ dc = (f77_func_data *) malloc(sizeof(f77_func_data) * *m);
+ if (*m > 0 && !dc) { *info = NLOPT_OUT_OF_MEMORY; return; }
+ for (i = 0; i < *m; ++i) {
+ dc[i].f = fc;
+ dc[i].f_data = fc_data + i * (fc_second_datum - fc_data);
+ }
+
+ *info = nlopt_minimize_constrained((nlopt_algorithm) *algorithm,
+ *n, f77_func_wrap, &d,
+ *m, f77_func_wrap,
+ dc, sizeof(f77_func_data),
+ lb, ub, x, minf,
+ *minf_max, *ftol_rel, *ftol_abs,
+ *xtol_rel,
+ *have_xtol_abs ? xtol_abs : 0,
+ *maxeval, *maxtime);
+
+ if (dc) free(dc);
+}
+
+/* nlopt_minimize */
+void F77(nloptm,NLOPTM)(int *info,
+ const int *algorithm,
+ const int *n, nlopt_f77_func f, void *f_data,
+ const double *lb, const double *ub,
+ double *x, double *minf,
+ const double *minf_max,
+ const double *ftol_rel, const double *ftol_abs,
+ const double *xtol_rel, const double *xtol_abs,
+ const int *have_xtol_abs,
+ const int *maxeval, const double *maxtime)
+{
+ int m = 0;
+ F77(nloptc,NLOPTC)(info, algorithm, n, f, f_data, &m, 0, 0, 0,
+ lb, ub, x, minf, minf_max, ftol_rel, ftol_abs,
+ xtol_rel, xtol_abs, have_xtol_abs, maxeval, maxtime);
+}
+
+void F77(nlosr,NLOSR)(const int *seed) { nlopt_srand((unsigned long) *seed); }
+void F77(nlosrt,NLOSRT)(void) { nlopt_srand_time(); }
+void F77(nloptv,NLOPTV)(int *major, int *minor, int *bugfix) {
+ nlopt_version(major, minor, bugfix);
+}
+void F77(nlogls,NLOGLS)(int *ideriv, int *inonderiv, int *maxeval)
+{
+ nlopt_algorithm deriv, nonderiv;
+ nlopt_get_local_search_algorithm(&deriv, &nonderiv, maxeval);
+ *ideriv = deriv;
+ *inonderiv = nonderiv;
+}
+void F77(nlosls,NLOSLS)(int *ideriv, int *inonderiv, int *maxeval)
+{
+ nlopt_algorithm deriv = (nlopt_algorithm) *ideriv;
+ nlopt_algorithm nonderiv = (nlopt_algorithm) *inonderiv;
+ nlopt_set_local_search_algorithm(deriv, nonderiv, *maxeval);
+}
~ianmdlvl / nlopt.git / commitdiff