BOX_BEHNKEN
Box-Behnken Sampling

BOX_BEHNKEN is a C library which returns a set of sample argument values to be used to analyze the behavior of a function.

Given DIM_NUM parameters, each of which is prescribed to lie between a minimum and maximum range, and a function taking these parameters as arguments, this program returns a table of sets of argument values that may be used to sample the function's behavior.

One sample point is the average. All other sample points are generated by setting a single spatial coordinate to its average value, and all other spatial coordinates to either the minimum or maximum. In this way, 1 + DIM_NUM*2^(DIM_NUM-1) points are generated.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Languages:

BOX_BEHNKEN is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

HALTON, a C library which computes elements of a Halton sequence.

HAMMERSLEY, a C library which computes elements of a Hammersley Quasi Monte Carlo (QMC) sequence, using a simple interface.

LATIN_RANDOM, a C library which computes elements of a Latin Hypercube dataset, choosing points at random.

UNIFORM, a C library which computes elements of a uniform pseudorandom sequence.

VAN_DER_CORPUT, a C library which computes elements of a (1 dimensional) van der Corput sequence.

Reference:

1. George Box, Donald Behnken,
   Some new three level designs for the study of quantitative variables,
   Technometrics,
   Volume 2, pages 455-475, 1960.

Source Code:

• box_behnken.c, the source code.
• box_behnken.h, the include file.

Examples and Tests:

• box_behnken_test.c, a sample calling program.
• box_behnken_test.txt, the output file.
• box_behnken_04_33.txt, contains a Box-Behnken dataset of 33 values for 4 parameters.

List of Routines:

• BOX_BEHNKEN returns a Box-Behnken design for the given number of factors.
• BOX_BEHNKEN_SIZE returns the size of a Box-Behnken design.
• I4_MAX returns the maximum of two I4's.
• I4_MIN returns the smaller of two I4's.
• I4_POWER returns the value of I^J.
• R8MAT_TRANSPOSE_PRINT prints an R8MAT, transposed.
• R8MAT_TRANSPOSE_PRINT_SOME prints some of an R8MAT, transposed.
• R8MAT_WRITE writes an R8MAT file.
• TIMESTAMP prints the current YMDHMS date as a time stamp.

You can go up one level to the C source codes.