1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
|
#! /usr/bin/env python
import openturns as ot
# Defining parameters
dimension = 5
size = 100
# Build OT LHS algorithm
lhs = ot.LHSExperiment(ot.JointDistribution([ot.Uniform()] * dimension), size)
# Generate design without optimization
design = lhs.generate()
# Defining space fillings
spaceFillingC2 = ot.SpaceFillingC2()
spaceFillingMinDist = ot.SpaceFillingMinDist()
spaceFillingPhiP = ot.SpaceFillingPhiP()
spaceFillingPhiP50 = ot.SpaceFillingPhiP(50)
# print the criteria on this design
print(
"C2=%f MinDist=%f PhiP=%f, PhiP(50)=%f"
% tuple(
[
sf.evaluate(design)
for sf in [
spaceFillingC2,
spaceFillingMinDist,
spaceFillingPhiP,
spaceFillingPhiP50,
]
]
)
)
|