orthogonal.OLHD_Cioppa07

orthogonal.OLHD_Cioppa07(m)

Orthogonal Latin Hyercube Design. Based on the construction method of Cioppa and Lucas (2007)

Parameters

Name	Type	Description	Default
`m`	int	A positive integer, and it must be greater than or equal to 2	required

Raises

Type	Description
ValueError	If m is not greater than or equal to 2

Returns

Type	Description
numpy.numpy.ndarray	An orthogonal LHD with the following run size: (n=2^m + 1) and factor size: (k= m+ (m-1 choose 2))

Examples: Create an orthogonal LHD with m=4. So n=2^m+1=17 runs and k=4+3=7 factors

import pyLHD
pyLHD.OLHD_Cioppa07(m=4)

array([[ 1., -2., -4., -8.,  3.,  7.,  5.],
       [ 2.,  1., -3., -7., -4., -8.,  6.],
       [ 3., -4.,  2., -6., -1.,  5., -7.],
       [ 4.,  3.,  1., -5.,  2., -6., -8.],
       [ 5., -6., -8.,  4.,  7., -3., -1.],
       [ 6.,  5., -7.,  3., -8.,  4., -2.],
       [ 7., -8.,  6.,  2., -5., -1.,  3.],
       [ 8.,  7.,  5.,  1.,  6.,  2.,  4.],
       [ 0.,  0.,  0.,  0.,  0.,  0.,  0.],
       [-1.,  2.,  4.,  8., -3., -7., -5.],
       [-2., -1.,  3.,  7.,  4.,  8., -6.],
       [-3.,  4., -2.,  6.,  1., -5.,  7.],
       [-4., -3., -1.,  5., -2.,  6.,  8.],
       [-5.,  6.,  8., -4., -7.,  3.,  1.],
       [-6., -5.,  7., -3.,  8., -4.,  2.],
       [-7.,  8., -6., -2.,  5.,  1., -3.],
       [-8., -7., -5., -1., -6., -2., -4.]])

Create an orthogonal LHD with m=5. So n=2^m+1=33 runs and k=5+7=11 factors

import pyLHD
pyLHD.OLHD_Cioppa07(m=5)

array([[  1.,  -2.,  -4.,  -8., -16.,   3.,   7.,  15.,   5.,  13.,   9.],
       [  2.,   1.,  -3.,  -7., -15.,  -4.,  -8., -16.,   6.,  14.,  10.],
       [  3.,  -4.,   2.,  -6., -14.,  -1.,   5.,  13.,  -7., -15.,  11.],
       [  4.,   3.,   1.,  -5., -13.,   2.,  -6., -14.,  -8., -16.,  12.],
       [  5.,  -6.,  -8.,   4., -12.,   7.,  -3.,  11.,  -1.,   9., -13.],
       [  6.,   5.,  -7.,   3., -11.,  -8.,   4., -12.,  -2.,  10., -14.],
       [  7.,  -8.,   6.,   2., -10.,  -5.,  -1.,   9.,   3., -11., -15.],
       [  8.,   7.,   5.,   1.,  -9.,   6.,   2., -10.,   4., -12., -16.],
       [  9., -10., -12., -16.,   8.,  11.,  15.,  -7.,  13.,  -5.,  -1.],
       [ 10.,   9., -11., -15.,   7., -12., -16.,   8.,  14.,  -6.,  -2.],
       [ 11., -12.,  10., -14.,   6.,  -9.,  13.,  -5., -15.,   7.,  -3.],
       [ 12.,  11.,   9., -13.,   5.,  10., -14.,   6., -16.,   8.,  -4.],
       [ 13., -14., -16.,  12.,   4.,  15., -11.,  -3.,  -9.,  -1.,   5.],
       [ 14.,  13., -15.,  11.,   3., -16.,  12.,   4., -10.,  -2.,   6.],
       [ 15., -16.,  14.,  10.,   2., -13.,  -9.,  -1.,  11.,   3.,   7.],
       [ 16.,  15.,  13.,   9.,   1.,  14.,  10.,   2.,  12.,   4.,   8.],
       [  0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.],
       [ -1.,   2.,   4.,   8.,  16.,  -3.,  -7., -15.,  -5., -13.,  -9.],
       [ -2.,  -1.,   3.,   7.,  15.,   4.,   8.,  16.,  -6., -14., -10.],
       [ -3.,   4.,  -2.,   6.,  14.,   1.,  -5., -13.,   7.,  15., -11.],
       [ -4.,  -3.,  -1.,   5.,  13.,  -2.,   6.,  14.,   8.,  16., -12.],
       [ -5.,   6.,   8.,  -4.,  12.,  -7.,   3., -11.,   1.,  -9.,  13.],
       [ -6.,  -5.,   7.,  -3.,  11.,   8.,  -4.,  12.,   2., -10.,  14.],
       [ -7.,   8.,  -6.,  -2.,  10.,   5.,   1.,  -9.,  -3.,  11.,  15.],
       [ -8.,  -7.,  -5.,  -1.,   9.,  -6.,  -2.,  10.,  -4.,  12.,  16.],
       [ -9.,  10.,  12.,  16.,  -8., -11., -15.,   7., -13.,   5.,   1.],
       [-10.,  -9.,  11.,  15.,  -7.,  12.,  16.,  -8., -14.,   6.,   2.],
       [-11.,  12., -10.,  14.,  -6.,   9., -13.,   5.,  15.,  -7.,   3.],
       [-12., -11.,  -9.,  13.,  -5., -10.,  14.,  -6.,  16.,  -8.,   4.],
       [-13.,  14.,  16., -12.,  -4., -15.,  11.,   3.,   9.,   1.,  -5.],
       [-14., -13.,  15., -11.,  -3.,  16., -12.,  -4.,  10.,   2.,  -6.],
       [-15.,  16., -14., -10.,  -2.,  13.,   9.,   1., -11.,  -3.,  -7.],
       [-16., -15., -13.,  -9.,  -1., -14., -10.,  -2., -12.,  -4.,  -8.]])