Pysmo Basics¶
PySMO Tutorial¶
Python-based Surrogate Modelling Objects (PySMO) provides tools for generating different types of reduced order models. PySMO currently provides tools for sampling and surrogate model generation.
Installation¶
PySMO is installed by default as part of IDAES. For instructions on installing IDAES, see the online documentation.
One-Shot Sampling with PySMO¶
The PySMO package offers five common sampling methods for one-shot design:
- Latin Hypercube Sampling (LHS)
- Full-Factorial Sampling
- Halton Sampling
- Hammersley Sampling
- Centroidal voronoi tessellation (CVT) sampling
PySMO provides two modes for data sampling: creation and selection.
- In creation mode, PySMO creates a specified number of sample points from the bounds provided by the user.
- In selection mode, PySMO selects a specified number of data points from a user-supplied dataset or file.
Generating samples:¶
For demonstration purposes, let us consider a problem for which we need twenty-five (25) samples of temperature and pressure from within the ranges T = 273K - 373K, P = 1 MPa - 50 MPa. Let us generate these samples in PySMO.
Step 1: Import PySMO's sampling tool¶
For this demonstration, we will attempt to generate the samples using the Hammersley sampling method.
In [1]:
from idaes.core.surrogate.pysmo.sampling import HammersleySampling
Step 2: Specify sampling information and initialize class¶
All the sampling tools (except full-factorial sampling) require the same keyword arguments:
- data_input : must be a list of lists containing the problem bounds (when creating points),
or an input dataset (when selecting points from a dataset)
- number_of_samples : number of samples to be created or selected.
- sampling_type : "creation" or "selection".
For full factorial sampling, the user needs to enter a list of points in each dimension in place of the number of samples. Full-factorial sampling requires other inputs - details may be found in the documentation.
For our example, we will create the bounds and then initialize the class with the number of samples.
In [2]:
bounds_info = [[273, 1],[373, 50]]
init_data = HammersleySampling(data_input=bounds_info, number_of_samples=25, sampling_type="creation")
Sampling type: creation
Step 3: Create the samples¶
The samples are created by calling the sample_points method on the initialized class.
In [3]:
samples = init_data.sample_points()
print(samples)
[[273. 1. ] [277. 25.5 ] [281. 13.25 ] [285. 37.75 ] [289. 7.125 ] [293. 31.625 ] [297. 19.375 ] [301. 43.875 ] [305. 4.0625 ] [309. 28.5625 ] [313. 16.3125 ] [317. 40.8125 ] [321. 10.1875 ] [325. 34.6875 ] [329. 22.4375 ] [333. 46.9375 ] [337. 2.53125] [341. 27.03125] [345. 14.78125] [349. 39.28125] [353. 8.65625] [357. 33.15625] [361. 20.90625] [365. 45.40625] [369. 5.59375]]
Simple as that, the samples have been created!
Now, let us visualize the samples in a 2-D plot.
Step 4: Visualize samples with matplotlib¶
In [4]:
from matplotlib import pyplot as plt
plt.plot(samples[:, 0], samples[:, 1], 'o')
plt.xlabel(r'Temperature', fontsize=11)
plt.xlabel(r'Pressure', fontsize=11)
plt.xlim(272, 374)
plt.ylim(0, 50)
plt.show()
Generating surrogates with PySMO¶
PySMO currently provides tools for generating three types of surrogates:
- Polynomial surrogates
- Radial basis function (RBF) surrogates, and
- Kriging surrogates
Details about thee various methods may be found in the documentation.
Generating polynomial models¶
The PolynomialRegression class trains polynomial models from data.
As an example, let us generate a surrogate for the Brainin function.
The true Brainin function is given by the expression:
\begin{gather} \hat{y}(x_{1},x_{2})=\left(x_{2}-\frac{5.1x_{1}^{2}}{4\pi^{2}}+\frac{5x_{1}}{\pi}-6\right)^{2}+10\left[\left(1-\frac{1}{8\pi}\right)\cos\left(x_{1}\right)+1\right]+5x_{1}\nonumber \\ x_{1}\in\left[-5,10\right];x_{2}\in\left[0,15\right] \end{gather}
We have generated 30 points from the function and saved the information in a text file called "brainin_30.txt". We will use this data to train a simple polynomial model. The data is in XY format, with the outputs $y$ in the third column.
Step 1: Import and visualize the data¶
In [5]:
import numpy as np
brainin_data = np.loadtxt('brainin_30.txt')
print(brainin_data,'\n\nDataset shape:', brainin_data.shape)
[[ 3.15107413e+00 4.17554078e+00 4.03815759e+00] [ 1.36776386e+00 1.26716420e+01 8.60127077e+01] [-4.92921716e+00 1.82353681e+00 2.41888089e+02] [ 5.06123627e+00 1.23877913e+01 1.37242426e+02] [-2.94940115e+00 8.62639994e+00 1.07587601e+01] [ 8.36982931e+00 3.13803183e+00 7.24964869e+00] [-2.22007671e+00 1.62565336e+00 7.71988248e+01] [-1.70453761e+00 1.46793568e+01 3.99820398e+01] [ 7.17524724e+00 2.57911519e+00 1.78453527e+01] [ 7.24337123e+00 4.11110621e+00 2.36937839e+01] [ 1.47556275e+00 1.41004473e+01 1.14292050e+02] [ 7.26474068e+00 5.04167925e+00 2.96703147e+01] [-2.36884319e+00 5.59248069e+00 2.71582680e+01] [-4.91467239e+00 3.78639530e+00 1.84993457e+02] [ 6.93493763e+00 2.28824569e-01 1.85302363e+01] [ 3.98265065e+00 9.05706809e+00 5.75698037e+01] [-3.42278472e+00 5.72915167e+00 5.30702398e+01] [-4.45285915e+00 1.33561735e+01 1.27882794e+01] [ 9.71381286e+00 8.99129832e-01 4.14567857e+00] [ 8.35818917e+00 8.65352249e+00 5.34055688e+01] [ 6.13719534e+00 9.45275905e+00 8.93001840e+01] [ 3.72763289e+00 3.06586980e-01 4.42072359e+00] [-1.84960133e+00 8.17027317e+00 8.83431235e+00] [ 6.53672757e+00 3.76042844e+00 2.62862798e+01] [-7.11564644e-01 1.27859263e+01 4.84986100e+01] [ 9.62509740e+00 1.32727994e+01 1.13457908e+02] [ 3.92617659e-01 8.98288419e+00 3.17439666e+01] [ 3.21934175e-01 5.10285323e+00 1.92673649e+01] [-2.32878516e+00 3.56541313e+00 5.02034858e+01] [-4.32706576e+00 7.58147144e+00 6.60457721e+01]] Dataset shape: (30, 3)
Let us visualize the data:
In [6]:
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import pyplot as plt
fig1 = plt.figure(figsize=(6, 4), tight_layout=True)
ax = fig1.add_subplot(111, projection='3d')
ax.scatter3D(brainin_data[:, 0], brainin_data[:, 1], brainin_data[:, 2], cmap=brainin_data[:, 2])
ax.set_xlabel(r'$x_{1}$', fontsize=11)
ax.set_ylabel(r'$x_{2}$', fontsize=11)
ax.set_zlabel(r'$y$', fontsize=11)
/tmp/ipykernel_115827/1699643538.py:5: UserWarning: No data for colormapping provided via 'c'. Parameters 'cmap' will be ignored ax.scatter3D(brainin_data[:, 0], brainin_data[:, 1], brainin_data[:, 2], cmap=brainin_data[:, 2])
Out[6]:
Text(0.5, 0, '$y$')
Step 2: Import the polynomial model tool¶
In [7]:
from idaes.core.surrogate.pysmo.polynomial_regression import PolynomialRegression
Step 3: Specify the regression settings and initialize the PolynomialRegression class¶
The PolynomialRegression class takes a keyword arguments:
- original_data_input : The dataset for regression training. training_data is expected to contain xy_data,
with the output values (y) in the last column.
- regression_data_input : same as above
- maximum_polynomial_order : maximum order of the polynomial to be generated
It also takes a number of optional arguments:
- multinomials : True/False option for specifying second-order bi-variate terms. default is False
- training_split : The training/cross-validation split of training data. Must be between 0 and 1.
Default is 0.75
- fname : Filename for saving results (.pickle extension).
- overwrite : Option determining whether any existing file with the same name supplied in 'fname'
should be overwritten.
For this example, let us consider a 4th order polynomial with interaction terms. We will split the data 80/20 betweeen training and cross-validation.
In [8]:
poly_class = PolynomialRegression(original_data_input=brainin_data, regression_data_input=brainin_data,
maximum_polynomial_order=4, multinomials=1, training_split=0.8, number_of_crossvalidations=10, overwrite=True)
===========================Polynomial Regression=============================================== No iterations will be run. Default parameter estimation method is used. Parameter estimation method: pyomo
Step 4: Extract variable names¶
Next, we extract Pyomo variable names from the dataset. This should be done always.
In [9]:
vars = poly_class.get_feature_vector()
We can view the variables using Pyomo's pprint function:
In [10]:
vars.pprint()
IndexedParam : Size=2, Index={0, 1}, Domain=Any, Default=None, Mutable=True
Key : Value
0 : 0
1 : 0
Step 5: Specify additional regression terms, if required.¶
This is one of the unique features of PySMO - it allows the user to specify additional regression features if they want.
The additional features must be specified in terms of the Pyomo variables created when calling the get_feature_vector()
For this example, let us create three additional features: $x_{1}^{2}x_{2}^{2}$, $exp(x_1)$ and $exp(x_2)$. We do this by calling the set_additional_terms function:
In [11]:
from pyomo.environ import exp
poly_class.set_additional_terms([vars[0] * vars[0] * vars[1] * vars[1], exp(vars[0]), exp(vars[1])])
That's it - those features will now exist in the model.
Note that set_additional_terms an optional call - the regression process works just fine without it.
Step 6: Train the surrogate and view results¶
Next, we train the polynomial surrogate by calling training:
In [12]:
poly_class.training()
No iterations will be run. Best surrogate model is of order 3 with a cross-val S.S. Error of 19.737617 ------------------------------------------------------------ The final coefficients of the regression terms are: k | 57.607003 (x_ 1 )^ 1 | -19.078561 (x_ 2 )^ 1 | -17.191563 (x_ 1 )^ 2 | 3.889924 (x_ 2 )^ 2 | 1.163328 (x_ 1 )^ 3 | -0.351928 (x_ 2 )^ 3 | 0.018398 x_ 1 .x_ 2 | 3.100135 The coefficients of the extra terms in additional_regression_features are: Coeff. additional_regression_features[ 1 ]: -0.017933 Coeff. additional_regression_features[ 2 ]: 0.004746 Coeff. additional_regression_features[ 3 ]: -1.3e-05 Regression model performance on training data: Order: 3 / MAE: 6.154513 / MSE: 65.995248 / R^2: 0.978531 Results saved in solution.pickle
Out[12]:
======================================================================================================================== Results of polynomial regression run: Polynomial order : 3 Number of terms in polynomial model: 11 Polynomial Expression: -------------------------- 57.60700349216609 -19.07856084658467*IndexedParam[0] -17.19156303535964*IndexedParam[1] + 3.88992417499501*IndexedParam[0]**2 + 1.1633280318838002*IndexedParam[1]**2 -0.35192841345409565*IndexedParam[0]**3 + 0.018398180908132653*IndexedParam[1]**3 + 3.1001348356693885*(IndexedParam[1]*IndexedParam[0]) -0.017933037213871822*(IndexedParam[0]*IndexedParam[0]*IndexedParam[1]*IndexedParam[1]) + 0.004745711043424125*exp(IndexedParam[0]) -1.2975320262079149e-05*exp(IndexedParam[1]) -------------------------- Model training errors: ----------------------- Mean Squared Error (MSE) : 65.99524781736321 Root Mean Squared Error (RMSE) : 8.123745922747904 Mean Absolute error (MSE) : 6.154513199272381 Goodness of fit (R2) : 0.9785305338453677 ========================================================================================================================
The polynomial model seems to fit well based on the $R^2$. It should be noted that the metrics are only an indication of how well of how well the model fit the training data - the user needs to verify the model's performance on a test data set if possible.
We can view the parity and residual plots for the fit:
In [13]:
poly_class.parity_residual_plots()
PySMO is also able to compute the confidence intervals on the regression coefficients obtained by calling confint_regression(). This is left as an exercise for the user.
Step 8 (Optional): Generate Pyomo expression¶
If the user wishes, they can generate the Pyomo expression for the polynomial fit using PySMO's generate_expression. To do this, the user must pass in a list of Pyomo variables corresponding to each variable in the input dataset.
As a demonstration, let us create the variables $x_1$ and $x_2$ and generate the pyomo expression based on them:
In [14]:
from pyomo.environ import Var, ConcreteModel
m = ConcreteModel()
m.x = Var([1, 2])
print(poly_class.generate_expression([m.x[1], m.x[2]]))
57.60700349216609 - 19.07856084658467*x[1] - 17.19156303535964*x[2] + 3.88992417499501*x[1]**2 + 1.1633280318838002*x[2]**2 - 0.35192841345409565*x[1]**3 + 0.018398180908132653*x[2]**3 + 3.1001348356693885*(x[2]*x[1]) - 0.017933037213871822*(x[1]*x[1]*x[2]*x[2]) + 0.004745711043424125*exp(x[1]) - 1.2975320262079149e-05*exp(x[2])
Step 9 (Optional): Predict output at any unsampled point¶
Based on the model we trained, we can predict the surrogate value at any previously unsampled point.
Let us evaluate the surrogate at three points:
- $x_{1}=5$, $x_{2}=8$ (true function value: 57.9908)
- $x_{1}=-3$, $x_{2}=10$ (true function value: 4.2461)
- $x_{1}=-2$, $x_{2}=3$. (true function value: 50.8899)
We will pass the points in as an array.
In [15]:
unsampled_points = np.array([ [5, 8], [-3, 10], [-2, 3]])
ys = poly_class.predict_output(unsampled_points)
print(ys)
[[57.7897913 ] [12.74008178] [54.28524655]]
The model performs fairly well in predicting the value at two of our sampled points but is off on the value at [-3, 2]. For better model performance, additional training data is needed in this region. We will leave this to the user to try.
Further information about using PySMO's polynomial regression tool can be found in the documentation.
Generating RBF models¶
The RadialBasisFunction class trains RBF models from data. For details about RBF models, the user should consult the documentation.
As an example, we will again consider the Brainin function. The same dataset loaded previously will be used.
Step 1: Import the data and the RBF tool¶
In [16]:
import numpy as np
brainin_data = np.loadtxt('brainin_30.txt')
from idaes.core.surrogate.pysmo.radial_basis_function import RadialBasisFunctions
Step 2: Specify the RBF settings and initialize the RadialBasisFunctions class¶
The RadialBasisFunctions class takes a number of keyword arguments:
- XY_data : The dataset forRBF training. training_data is expected to contain xy_data,
with the output values (y) in the last column.
It also takes a number of optional arguments:
- regularization : Boolean variable determining whether regularization is done. Default is True.
- basis_function : Basis function transformation to be applied to the training data. PySMO offers
six basis function types including the Gaussian and Spline transformations. User
should consult documentation for full list of options.
- fname : Filename for saving (.pickle extension)
- overwrite : Option determining whether any existing file with the same name supplied in 'fname'
should be overwritten.
For this demonstration, we will train an RBF model with a Gaussian basis function:
In [17]:
rbf_class = RadialBasisFunctions(XY_data=brainin_data, basis_function='gaussian', overwrite=True)
Warning: solution.pickle already exists; previous file will be overwritten. Default parameter estimation method is used. Parameter estimation method: algebraic Basis function: gaussian Regularization done: True
Step 3: Extract variable names¶
Next, we extract Pyomo variable names from the dataset.
In [18]:
vars = rbf_class.get_feature_vector()
Step 4: Train the RBF surrogate¶
Next, we train the RBF surrogate by calling training:
In [19]:
rbf_class.training()
=========================================================================================================== 0.001 | 1e-05 | 1.3301889998012952 | 1.5846892219716487e+19 | 3518.7169222165 | 2999999.79962304 | 6.66133770282471e-10 0.001 | 2e-05 | 1.3231321486466063 | 1.5846892219716487e+19 | 3518.7169222165 | 1500000.399777789 | 3.3306699615604816e-10 0.001 | 5e-05 | 1.3204882439267729 | 1.5846892219716487e+19 | 3518.7169222165 | 600000.7598786094 | 1.332269316819644e-10 0.001 | 7.5e-05 | 1.3203184940700354 | 1.5846892219716487e+19 | 3518.7169222165 | 400000.83991739847 | 8.881802846913944e-11 0.001 | 0.0001 | 1.3203209127604238 | 1.5846892219716487e+19 | 3518.7169222165 | 300000.8799381776 | 6.66135768630344e-11 0.001 | 0.0002 | 1.3204538291268846 | 1.5846892219716487e+19 | 3518.7169222165 | 150000.93996800907 | 3.330689945357991e-11 0.001 | 0.0005 | 1.3206180157500538 | 1.5846892219716487e+19 | 3518.7169222165 | 60000.97598693101 | 1.3322893008134386e-11 0.001 | 0.00075 | 1.3206637488446558 | 1.5846892219716487e+19 | 3518.7169222165 | 40000.98399128094 | 8.882002686956478e-12 0.001 | 0.001 | 1.3206877851837313 | 1.5846892219716487e+19 | 3518.7169222165 | 30000.987993460352 | 6.661557526368512e-12 0.001 | 0.002 | 1.3207246211571395 | 1.5846892219716487e+19 | 3518.7169222165 | 15000.993996716828 | 3.3308897854837546e-12 0.001 | 0.005 | 1.320744101068588 | 1.5846892219716487e+19 | 3518.7169222165 | 6000.997598685179 | 1.332489140956112e-12 0.001 | 0.0075 | 1.3207457685491835 | 1.5846892219716487e+19 | 3518.7169222165 | 4000.9983991224544 | 8.884001088388281e-13 0.001 | 0.01 | 1.3207447379585622 | 1.5846892219716487e+19 | 3518.7169222165 | 3000.9987993419977 | 6.663555927803872e-13 0.001 | 0.02 | 1.3207338807217226 | 1.5846892219716487e+19 | 3518.7169222165 | 1500.999399670947 | 3.332888186926446e-13 0.001 | 0.05 | 1.3206922424549317 | 1.5846892219716487e+19 | 3518.7169222165 | 600.9997598683274 | 1.3344875424000144e-13 0.001 | 0.075 | 1.3206569542287037 | 1.5846892219716487e+19 | 3518.7169222165 | 400.9998399122136 | 8.903985102830827e-14 0.001 | 0.1 | 1.3206222275255648 | 1.5846892219716487e+19 | 3518.7169222165 | 300.9998799341613 | 6.68353994224627e-14 0.001 | 0.2 | 1.3204909253945514 | 1.5846892219716487e+19 | 3518.7169222165 | 150.99993996707875 | 3.3528722013693445e-14 0.001 | 0.5 | 1.320173489756647 | 1.5846892219716487e+19 | 3518.7169222165 | 60.99997598683119 | 1.3544715568432328e-14 0.001 | 0.75 | 1.3199938819612358 | 1.5846892219716487e+19 | 3518.7169222165 | 40.99998399122082 | 9.103825247263236e-15 0.001 | 1 | 1.3198876280547238 | 1.5846892219716487e+19 | 3518.7169222165 | 30.999987993415616 | 6.8833800866786844e-15 0.002 | 1e-05 | 1.362387822291202 | 5.682709296695643e+18 | 1261.8149406885866 | 2999996.1977211484 | 6.661329704995885e-10 0.002 | 2e-05 | 1.3439946308757786 | 5.682709296695643e+18 | 1261.8149406885866 | 1499998.59873096 | 3.3306659624331657e-10 0.002 | 5e-05 | 1.32722464857336 | 5.682709296695643e+18 | 1261.8149406885866 | 600000.0394952487 | 1.3322677172472568e-10 0.002 | 7.5e-05 | 1.3235440099224118 | 5.682709296695643e+18 | 1261.8149406885866 | 400000.3596569235 | 8.881792182989201e-11 0.002 | 0.0001 | 1.3220037168869367 | 5.682709296695643e+18 | 1261.8149406885866 | 300000.5197424581 | 6.661349688351817e-11 0.002 | 0.0002 | 1.3204880868273354 | 5.682709296695643e+18 | 1261.8149406885866 | 150000.75986908065 | 3.3306859463584503e-11 0.002 | 0.0005 | 1.3203538501484549 | 5.682709296695643e+18 | 1261.8149406885866 | 60000.90394757265 | 1.3322877012183521e-11 0.002 | 0.00075 | 1.3204389361296947 | 5.682709296695643e+18 | 1261.8149406885866 | 40000.93596503106 | 8.881992022986797e-12 0.002 | 0.001 | 1.3205000397654638 | 5.682709296695643e+18 | 1261.8149406885866 | 30000.951973759784 | 6.6615495283883294e-12 0.002 | 0.002 | 1.3206158464558628 | 5.682709296695643e+18 | 1261.8149406885866 | 15000.97598686735 | 3.3308857864938424e-12 0.002 | 0.005 | 1.3206969267656399 | 5.682709296695643e+18 | 1261.8149406885866 | 6000.990394744821 | 1.3324875413600215e-12 0.002 | 0.0075 | 1.3207137929011332 | 5.682709296695643e+18 | 1261.8149406885866 | 4000.9935964960405 | 8.883990424415434e-13 0.002 | 0.01 | 1.3207205715037869 | 5.682709296695643e+18 | 1261.8149406885866 | 3000.995197371845 | 6.663547929823476e-13 0.002 | 0.02 | 1.3207216922930969 | 5.682709296695643e+18 | 1261.8149406885866 | 1500.9975986859042 | 3.332884187936323e-13 0.002 | 0.05 | 1.3206873938685026 | 5.682709296695643e+18 | 1261.8149406885866 | 600.9990394743195 | 1.3344859428039857e-13 0.002 | 0.075 | 1.3206537523643644 | 5.682709296695643e+18 | 1261.8149406885866 | 400.999359649544 | 8.903974438857353e-14 0.002 | 0.1 | 1.3206198511578071 | 5.682709296695643e+18 | 1261.8149406885866 | 300.99951973715895 | 6.683531944266162e-14 0.002 | 0.2 | 1.3204897899439867 | 5.682709296695643e+18 | 1261.8149406885866 | 150.99975986857845 | 3.35286820237931e-14 0.002 | 0.5 | 1.3201731005146333 | 5.682709296695643e+18 | 1261.8149406885866 | 60.99990394743123 | 1.3544699572472225e-14 0.002 | 0.75 | 1.31999365892768 | 5.682709296695643e+18 | 1261.8149406885866 | 40.9999359649541 | 9.103814583289815e-15 0.002 | 1 | 1.3198874884067837 | 5.682709296695643e+18 | 1261.8149406885866 | 30.99995197371558 | 6.88337208869862e-15 0.005 | 1e-05 | 1.3982813940010914 | 4.668179943166246e+18 | 1036.5441711993042 | 2999970.9839561055 | 6.661273719190909e-10 0.005 | 2e-05 | 1.3887224316295734 | 4.668179943166246e+18 | 1036.5441711993042 | 1499985.9920229833 | 3.330637969918245e-10 0.005 | 5e-05 | 1.3683630926291095 | 4.668179943166246e+18 | 1036.5441711993042 | 599994.9967733887 | 1.3322565201554253e-10 0.005 | 7.5e-05 | 1.3574022549734979 | 4.668179943166246e+18 | 1036.5441711993042 | 399996.99784555874 | 8.881717535781569e-11 0.005 | 0.0001 | 1.349656526523195 | 4.668179943166246e+18 | 1036.5441711993042 | 299997.9983864001 | 6.661293703000839e-11 0.005 | 0.0002 | 1.3338899892293052 | 4.668179943166246e+18 | 1036.5441711993042 | 149999.49919137143 | 3.330657953690062e-11 0.005 | 0.0005 | 1.3231314606160027 | 4.668179943166246e+18 | 1036.5441711993042 | 60000.39967648298 | 1.332276504150864e-11 0.005 | 0.00075 | 1.3213658709060543 | 4.668179943166246e+18 | 1036.5441711993042 | 40000.59978412077 | 8.88191737582939e-12 0.005 | 0.001 | 1.3207434325090537 | 4.668179943166246e+18 | 1036.5441711993042 | 30000.699838098662 | 6.661493543025068e-12 0.005 | 0.002 | 1.3203109441739103 | 4.668179943166246e+18 | 1036.5441711993042 | 15000.849919034625 | 3.330857793811731e-12 0.005 | 0.005 | 1.3204473592564487 | 4.668179943166246e+18 | 1036.5441711993042 | 6000.939967613371 | 1.332476344287541e-12 0.005 | 0.0075 | 1.3205266147949006 | 4.668179943166246e+18 | 1036.5441711993042 | 4000.9599784081925 | 8.883915777265089e-13 0.005 | 0.01 | 1.3205722511348608 | 4.668179943166246e+18 | 1036.5441711993042 | 3000.96998380598 | 6.663491944460764e-13 0.005 | 0.02 | 1.3206416674457944 | 4.668179943166246e+18 | 1036.5441711993042 | 1500.9849919029339 | 3.3328561952548827e-13 0.005 | 0.05 | 1.3206543062552825 | 4.668179943166246e+18 | 1036.5441711993042 | 600.9939967611581 | 1.334474745731469e-13 0.005 | 0.075 | 1.3206317178276634 | 4.668179943166246e+18 | 1036.5441711993042 | 400.9959978407684 | 8.903899791707212e-14 0.005 | 0.1 | 1.3206034290426698 | 4.668179943166246e+18 | 1036.5441711993042 | 300.9969983805791 | 6.683475958903597e-14 0.005 | 0.2 | 1.3204818942961354 | 4.668179943166246e+18 | 1036.5441711993042 | 150.99849919028716 | 3.352840209697997e-14 0.005 | 0.5 | 1.3201703838941496 | 4.668179943166246e+18 | 1036.5441711993042 | 60.99939967611464 | 1.3544587601746958e-14 0.005 | 0.75 | 1.3199921011576619 | 4.668179943166246e+18 | 1036.5441711993042 | 40.99959978407642 | 9.103739936139648e-15 0.005 | 1 | 1.3198865127506894 | 4.668179943166246e+18 | 1036.5441711993042 | 30.999699838057268 | 6.883316103335983e-15 0.0075 | 1e-05 | 1.4043005938815813 | 3.035757992362053e+18 | 674.0736840620383 | 2999933.4646548154 | 6.661190409606589e-10 0.0075 | 2e-05 | 1.399403263321862 | 3.035757992362053e+18 | 674.0736840620383 | 1499967.2322215145 | 3.3305963147911887e-10 0.0075 | 5e-05 | 1.3868056479529616 | 3.035757992362053e+18 | 674.0736840620383 | 599987.4928811593 | 1.3322398581675705e-10 0.0075 | 7.5e-05 | 1.3783815287747978 | 3.035757992362053e+18 | 674.0736840620383 | 399991.995248878 | 8.881606455821212e-11 0.0075 | 0.0001 | 1.3714017808065482 | 3.035757992362053e+18 | 674.0736840620383 | 299994.2464359732 | 6.661210392965815e-11 0.0075 | 0.0002 | 1.3527834516582864 | 3.035757992362053e+18 | 674.0736840620383 | 149997.62321875116 | 3.3306162987301305e-11 0.0075 | 0.0005 | 1.3320471009308383 | 3.035757992362053e+18 | 674.0736840620383 | 59999.64928699191 | 1.3322598421570555e-11 0.0075 | 0.00075 | 1.3264972569268667 | 3.035757992362053e+18 | 674.0736840620383 | 40000.099524720186 | 8.881806295928426e-12 0.0075 | 0.001 | 1.323914758231772 | 3.035757992362053e+18 | 674.0736840620383 | 30000.32464347602 | 6.661410233083313e-12 0.0075 | 0.002 | 1.3209497744221874 | 3.035757992362053e+18 | 674.0736840620383 | 15000.662321737802 | 3.330816138844073e-12 0.0075 | 0.005 | 1.3203066752602552 | 3.035757992362053e+18 | 674.0736840620383 | 6000.864928690148 | 1.3324596822994801e-12 0.0075 | 0.0075 | 1.3203572891448587 | 3.035757992362053e+18 | 674.0736840620383 | 4000.909952460054 | 8.883804697346184e-13 0.0075 | 0.01 | 1.3204149453606322 | 3.035757992362053e+18 | 674.0736840620383 | 3000.932464345157 | 6.66340863452221e-13 0.0075 | 0.02 | 1.3205391985944273 | 3.035757992362053e+18 | 674.0736840620383 | 1500.9662321724022 | 3.3328145402853387e-13 0.0075 | 0.05 | 1.3206077953595634 | 3.035757992362053e+18 | 674.0736840620383 | 600.9864928689558 | 1.3344580837436743e-13 0.0075 | 0.075 | 1.3206001441008732 | 3.035757992362053e+18 | 674.0736840620383 | 400.99099524596835 | 8.903788711788615e-14 0.0075 | 0.1 | 1.3205796752193628 | 3.035757992362053e+18 | 674.0736840620383 | 300.99324643447414 | 6.68339264896454e-14 0.0075 | 0.2 | 1.320470314359398 | 3.035757992362053e+18 | 674.0736840620383 | 150.99662321723687 | 3.352798554728517e-14 0.0075 | 0.5 | 1.3201663674229878 | 3.035757992362053e+18 | 674.0736840620383 | 60.99864928689451 | 1.3544420981869035e-14 0.0075 | 0.75 | 1.3199897942613736 | 3.035757992362053e+18 | 674.0736840620383 | 40.99909952459634 | 9.103628856221033e-15 0.0075 | 1 | 1.3198850669656146 | 3.035757992362053e+18 | 674.0736840620383 | 30.999324643447252 | 6.883232793397032e-15 0.01 | 1e-05 | 1.4061992061350406 | 2.761018195749752e+18 | 613.0691944660764 | 2999880.93890562 | 6.661073779014304e-10 0.01 | 2e-05 | 1.4035355143415935 | 2.761018195749752e+18 | 613.0691944660764 | 1499940.9693439263 | 3.330537999488406e-10 0.01 | 5e-05 | 1.3956271825129098 | 2.761018195749752e+18 | 613.0691944660764 | 599976.9877288077 | 1.3322165320435347e-10 0.01 | 7.5e-05 | 1.3897600752421233 | 2.761018195749752e+18 | 613.0691944660764 | 399984.9918146977 | 8.881450948343643e-11 0.01 | 0.0001 | 1.3845164268268453 | 2.761018195749752e+18 | 613.0691944660764 | 299988.99385947915 | 6.661093762338569e-11 0.01 | 0.0002 | 1.3683402006556789 | 2.761018195749752e+18 | 613.0691944660764 | 149994.99692907472 | 3.330557983384768e-11 0.01 | 0.0005 | 1.3439837129674193 | 2.761018195749752e+18 | 613.0691944660764 | 59998.59877132549 | 1.3322365160234437e-11 0.01 | 0.00075 | 1.3350937119914896 | 2.761018195749752e+18 | 613.0691944660764 | 39999.399180852306 | 8.881650788350971e-12 0.01 | 0.001 | 1.3301839804868678 | 2.761018195749752e+18 | 613.0691944660764 | 29999.79938562744 | 6.661293602411842e-12 0.01 | 0.002 | 1.3231293127498855 | 2.761018195749752e+18 | 613.0691944660764 | 15000.399692801779 | 3.330757823505732e-12 0.01 | 0.005 | 1.320483075844326 | 2.761018195749752e+18 | 613.0691944660764 | 6000.7598771196235 | 1.3324363561650062e-12 0.01 | 0.0075 | 1.320310037166994 | 2.761018195749752e+18 | 613.0691944660764 | 4000.8399180791366 | 8.883649189781765e-13 0.01 | 0.01 | 1.3203089432056807 | 2.761018195749752e+18 | 613.0691944660764 | 3000.8799385589855 | 6.663292003847822e-13 0.01 | 0.02 | 1.3204273423452135 | 2.761018195749752e+18 | 613.0691944660764 | 1500.9399692794188 | 3.332756224948372e-13 0.01 | 0.05 | 1.3205480387004696 | 2.761018195749752e+18 | 613.0691944660764 | 600.9759877117691 | 1.3344347576089024e-13 0.01 | 0.075 | 1.3205583477737064 | 2.761018195749752e+18 | 613.0691944660764 | 400.9839918078364 | 8.903633204223302e-14 0.01 | 0.1 | 1.3205477786808832 | 2.761018195749752e+18 | 613.0691944660764 | 300.9879938558759 | 6.683276018290572e-14 0.01 | 0.2 | 1.3204544411843475 | 2.761018195749752e+18 | 613.0691944660764 | 150.99399692793608 | 3.352740239391496e-14 0.01 | 0.5 | 1.32016079672088 | 2.761018195749752e+18 | 613.0691944660764 | 60.99759877117414 | 1.3544187720520937e-14 0.01 | 0.75 | 1.3199865871021215 | 2.761018195749752e+18 | 613.0691944660764 | 40.998399180782776 | 9.103473348655639e-15 0.01 | 1 | 1.3198830550759417 | 2.761018195749752e+18 | 613.0691944660764 | 30.998799385587123 | 6.883116162722996e-15 0.02 | 1e-05 | 1.399196782174963 | 6.535273978061908e+18 | 1451.122328535594 | 2999520.807994778 | 6.660274127756111e-10 0.02 | 2e-05 | 1.4033527819989597 | 6.535273978061908e+18 | 1451.122328535594 | 1499760.9038745863 | 3.330138173828404e-10 0.02 | 5e-05 | 1.4040547659101899 | 6.535273978061908e+18 | 1451.122328535594 | 599904.9615295177 | 1.3320566017538786e-10 0.02 | 7.5e-05 | 1.4028711082099645 | 6.535273978061908e+18 | 1451.122328535594 | 399936.9743491899 | 8.880384746427825e-11 0.02 | 0.0001 | 1.4014091897450747 | 6.535273978061908e+18 | 1451.122328535594 | 299952.980761774 | 6.660294110933362e-11 0.02 | 0.0002 | 1.395251376299013 | 6.535273978061908e+18 | 1451.122328535594 | 149976.99038140592 | 3.33015815770845e-11 0.02 | 0.0005 | 1.379670117393505 | 6.535273978061908e+18 | 1451.122328535594 | 59991.39615214057 | 1.3320765857503097e-11 0.02 | 0.00075 | 1.369923208565533 | 6.535273978061908e+18 | 1451.122328535594 | 39994.59743472255 | 8.880584586528639e-12 0.02 | 0.001 | 1.3622921102964674 | 6.535273978061908e+18 | 1451.122328535594 | 29996.198076042776 | 6.6604939510479024e-12 0.02 | 0.002 | 1.343949612469127 | 6.535273978061908e+18 | 1451.122328535594 | 14998.599038009137 | 3.3303579978236935e-12 0.02 | 0.005 | 1.3272089427378695 | 6.535273978061908e+18 | 1451.122328535594 | 6000.039615202972 | 1.3322764258922807e-12 0.02 | 0.0075 | 1.3235324618626603 | 6.535273978061908e+18 | 1451.122328535594 | 4000.359743468308 | 8.8825829879642e-13 0.02 | 0.01 | 1.321992371490607 | 6.535273978061908e+18 | 1451.122328535594 | 3000.5198076010442 | 6.662492352485048e-13 0.02 | 0.02 | 1.3204676186157416 | 6.535273978061908e+18 | 1451.122328535594 | 1500.7599038005399 | 3.3323563992671887e-13 0.02 | 0.05 | 1.320292865976956 | 6.535273978061908e+18 | 1451.122328535594 | 600.9039615201932 | 1.3342748273363752e-13 0.02 | 0.075 | 1.3203431410941742 | 6.535273978061908e+18 | 1451.122328535594 | 400.9359743467918 | 8.902567002406587e-14 0.02 | 0.1 | 1.3203699372606283 | 6.535273978061908e+18 | 1451.122328535594 | 300.95198076009075 | 6.682476366927997e-14 0.02 | 0.2 | 1.3203560056794423 | 6.535273978061908e+18 | 1451.122328535594 | 150.97599038004483 | 3.352340413710238e-14 0.02 | 0.5 | 1.3201242272878075 | 6.535273978061908e+18 | 1451.122328535594 | 60.99039615201791 | 1.3542588417795966e-14 0.02 | 0.75 | 1.3199652973344644 | 6.535273978061908e+18 | 1451.122328535594 | 40.993597434678556 | 9.102407146838977e-15 0.02 | 1 | 1.3198696409051942 | 6.535273978061908e+18 | 1451.122328535594 | 30.99519807600895 | 6.882316511360498e-15 0.05 | 1e-05 | 1.1802626763642516 | 2.451250238981105e+18 | 544.2868908869481 | 2997002.130533861 | 6.654681540338683e-10 0.05 | 2e-05 | 1.2638668839656086 | 2.451250238981105e+18 | 544.2868908869481 | 1498501.5652295854 | 3.327341880309443e-10 0.05 | 5e-05 | 1.3399914336689571 | 2.451250238981105e+18 | 544.2868908869481 | 599401.2260861265 | 1.3309380843787333e-10 0.05 | 7.5e-05 | 1.3607917455561922 | 2.451250238981105e+18 | 544.2868908869481 | 399601.1507228575 | 8.872927963984478e-11 0.05 | 0.0001 | 1.3717399149158809 | 2.451250238981105e+18 | 544.2868908869481 | 299701.1130407243 | 6.654701524071978e-11 0.05 | 0.0002 | 1.3883571964198669 | 2.451250238981105e+18 | 544.2868908869481 | 149851.05651985976 | 3.32736186425508e-11 0.05 | 0.0005 | 1.3960809195491266 | 2.451250238981105e+18 | 544.2868908869481 | 59941.02260788699 | 1.3309580683770636e-11 0.05 | 0.00075 | 1.395825918790349 | 2.451250238981105e+18 | 544.2868908869481 | 39961.015071857575 | 8.873127804033837e-12 0.05 | 0.001 | 1.3944270862080228 | 2.451250238981105e+18 | 544.2868908869481 | 29971.011303891435 | 6.654901364176221e-12 0.05 | 0.002 | 1.3868038656183426 | 2.451250238981105e+18 | 544.2868908869481 | 14986.005651939386 | 3.327561704389167e-12 0.05 | 0.005 | 1.3676102367336624 | 2.451250238981105e+18 | 544.2868908869481 | 5995.002260776849 | 1.3311579085188649e-12 0.05 | 0.0075 | 1.3569104532324432 | 2.451250238981105e+18 | 544.2868908869481 | 3997.0015071830558 | 8.875126205472163e-13 0.05 | 0.01 | 1.349296442539508 | 2.451250238981105e+18 | 544.2868908869481 | 2998.0011303873994 | 6.656899765616674e-13 0.05 | 0.02 | 1.333727803082708 | 2.451250238981105e+18 | 544.2868908869481 | 1499.5005651935342 | 3.3295601058325944e-13 0.05 | 0.05 | 1.3230623528240189 | 2.451250238981105e+18 | 544.2868908869481 | 600.400226077414 | 1.3331563099625887e-13 0.05 | 0.075 | 1.3212917568278946 | 2.451250238981105e+18 | 544.2868908869481 | 400.6001507182712 | 8.895110219914653e-14 0.05 | 0.1 | 1.3206495137172838 | 2.451250238981105e+18 | 544.2868908869481 | 300.70011303870484 | 6.676883780059147e-14 0.05 | 0.2 | 1.3201079927565529 | 2.451250238981105e+18 | 544.2868908869481 | 150.85005651935109 | 3.3495441202757955e-14 0.05 | 0.5 | 1.3199430633401978 | 2.451250238981105e+18 | 544.2868908869481 | 60.94002260774041 | 1.3531403244058195e-14 0.05 | 0.75 | 1.319849098238076 | 2.451250238981105e+18 | 544.2868908869481 | 40.960015071826966 | 9.094950364347147e-15 0.05 | 1 | 1.319793746770628 | 2.451250238981105e+18 | 544.2868908869481 | 30.970011303870212 | 6.876723924491615e-15 0.075 | 1e-05 | 1.0289125857143073 | 7.107621288921808e+19 | 15782.089610553847 | 2993261.3407904925 | 6.646375318531944e-10 0.075 | 2e-05 | 1.0683357422994502 | 7.107621288921808e+19 | 15782.089610553847 | 1496631.1703660372 | 3.3231887694241395e-10 0.075 | 5e-05 | 1.178819448585938 | 7.107621288921808e+19 | 15782.089610553847 | 598653.0681405907 | 1.329276840024353e-10 0.075 | 7.5e-05 | 1.2298615813478657 | 7.107621288921808e+19 | 15782.089610553847 | 399102.37875589466 | 8.861853001549284e-11 0.075 | 0.0001 | 1.2625315302909863 | 7.107621288921808e+19 | 15782.089610553847 | 299327.0340663092 | 6.646395302263502e-11 0.075 | 0.0002 | 1.324126963007039 | 7.107621288921808e+19 | 15782.089610553847 | 149664.01703300458 | 3.323208753358666e-11 0.075 | 0.0005 | 1.3699287056818163 | 7.107621288921808e+19 | 15782.089610553847 | 59866.20681313908 | 1.3292968240183684e-11 0.075 | 0.00075 | 1.3804567145240605 | 7.107621288921808e+19 | 15782.089610553847 | 39911.13787534549 | 8.862052841639543e-12 0.075 | 0.001 | 1.3852806670848943 | 7.107621288921808e+19 | 15782.089610553847 | 29933.603406512757 | 6.6465951423816964e-12 0.075 | 0.002 | 1.3897667849500999 | 7.107621288921808e+19 | 15782.089610553847 | 14967.301703253015 | 3.323408593492564e-12 0.075 | 0.005 | 1.3829509021399886 | 7.107621288921808e+19 | 15782.089610553847 | 5987.5206812987935 | 1.329496664159445e-12 0.075 | 0.0075 | 1.3758254841268973 | 7.107621288921808e+19 | 15782.089610553847 | 3992.0137875321957 | 8.864051243078643e-13 0.075 | 0.01 | 1.3694975238619214 | 7.107621288921808e+19 | 15782.089610553847 | 2994.2603406490966 | 6.648593543821183e-13 0.075 | 0.02 | 1.3518655278680716 | 7.107621288921808e+19 | 15782.089610553847 | 1497.6301703244728 | 3.325406994935049e-13 0.075 | 0.05 | 1.3317271553627905 | 7.107621288921808e+19 | 15782.089610553847 | 599.6520681297815 | 1.3314950656035528e-13 0.075 | 0.075 | 1.3262949589526694 | 7.107621288921808e+19 | 15782.089610553847 | 400.1013787531852 | 8.884035257521132e-14 0.075 | 0.1 | 1.3237555826770966 | 7.107621288921808e+19 | 15782.089610553847 | 300.3260340648877 | 6.668577558263949e-14 0.075 | 0.2 | 1.3207755491672855 | 7.107621288921808e+19 | 15782.089610553847 | 150.66301703244372 | 3.345391009378223e-14 0.075 | 0.5 | 1.3198839335916897 | 7.107621288921808e+19 | 15782.089610553847 | 60.865206812977505 | 1.3514790800467914e-14 0.075 | 0.75 | 1.3197723726186077 | 7.107621288921808e+19 | 15782.089610553847 | 40.91013787531836 | 9.083875401953626e-15 0.075 | 1 | 1.3197346497546034 | 7.107621288921808e+19 | 15782.089610553847 | 30.932603406488738 | 6.8684177026964695e-15 0.1 | 1e-05 | 1.0798981068284321 | 2.7591862966009334e+18 | 612.6624311433145 | 2988038.7413438815 | 6.6347788182239e-10 0.1 | 2e-05 | 1.0344111473051316 | 2.7591862966009334e+18 | 612.6624311433145 | 1494019.8706805534 | 3.3173905193540984e-10 0.1 | 5e-05 | 1.0485665416876455 | 2.7591862966009334e+18 | 612.6624311433145 | 597608.5482320086 | 1.3269575399199786e-10 0.1 | 7.5e-05 | 1.0859812540934977 | 2.7591862966009334e+18 | 612.6624311433145 | 398406.0321522719 | 8.846391000900053e-11 0.1 | 0.0001 | 1.119927128930482 | 2.7591862966009334e+18 | 612.6624311433145 | 298804.774114426 | 6.634798801795094e-11 0.1 | 0.0002 | 1.2092033599794523 | 2.7591862966009334e+18 | 612.6624311433145 | 149402.88705536202 | 3.3174105030866934e-11 0.1 | 0.0005 | 1.3058031834922963 | 2.7591862966009334e+18 | 612.6624311433145 | 59761.75482214557 | 1.3269775239109898e-11 0.1 | 0.00075 | 1.3347573633948888 | 2.7591862966009334e+18 | 612.6624311433145 | 39841.503214742595 | 8.846590840936884e-12 0.1 | 0.001 | 1.3503821090955594 | 2.7591862966009334e+18 | 612.6624311433145 | 29881.377411065852 | 6.634998641855872e-12 0.1 | 0.002 | 1.3741948871332808 | 2.7591862966009334e+18 | 612.6624311433145 | 14941.188705526605 | 3.317610343228995e-12 0.1 | 0.005 | 1.3835828525725806 | 2.7591862966009334e+18 | 612.6624311433145 | 5977.075482208737 | 1.32717736405413e-12 0.1 | 0.0075 | 1.3817104633773385 | 2.7591862966009334e+18 | 612.6624311433145 | 3985.0503214727046 | 8.848589242377757e-13 0.1 | 0.01 | 1.3784840301739398 | 2.7591862966009334e+18 | 612.6624311433145 | 2989.0377411040727 | 6.636997043294618e-13 0.1 | 0.02 | 1.3653609610262425 | 2.7591862966009334e+18 | 612.6624311433145 | 1495.0188705520156 | 3.3196087446718884e-13 0.1 | 0.05 | 1.3428726584962258 | 2.7591862966009334e+18 | 612.6624311433145 | 598.6075482207964 | 1.3291757654982836e-13 0.1 | 0.075 | 1.3344025162997275 | 2.7591862966009334e+18 | 612.6624311433145 | 399.40503214719104 | 8.868573256819246e-14 0.1 | 0.1 | 1.3296947739898703 | 2.7591862966009334e+18 | 612.6624311433145 | 299.8037741103933 | 6.656981057737561e-14 0.1 | 0.2 | 1.3228751118106656 | 2.7591862966009334e+18 | 612.6624311433145 | 150.40188705519594 | 3.3395927591150163e-14 0.1 | 0.5 | 1.320130553011318 | 2.7591862966009334e+18 | 612.6624311433145 | 60.76075482207834 | 1.3491597799415076e-14 0.1 | 0.75 | 1.3198268366427486 | 2.7591862966009334e+18 | 612.6624311433145 | 40.84050321471894 | 9.068413401251739e-15 0.1 | 1 | 1.3197458467857657 | 2.7591862966009334e+18 | 612.6624311433145 | 30.880377411039177 | 6.856821202170055e-15 0.2 | 1e-05 | 1.0348780714321624 | 4.67602569807196e+18 | 1038.286278747682 | 2952676.85155188 | 6.556259649741225e-10 0.2 | 2e-05 | 1.099054270540231 | 4.67602569807196e+18 | 1038.286278747682 | 1476338.9256853214 | 3.2781309348924236e-10 0.2 | 5e-05 | 1.1186681575679016 | 4.67602569807196e+18 | 1038.286278747682 | 590536.1702524232 | 1.3112537061764032e-10 0.2 | 7.5e-05 | 1.1064632913455217 | 4.67602569807196e+18 | 1038.286278747682 | 393691.1135068378 | 8.741698776112146e-11 0.2 | 0.0001 | 1.0918067697066494 | 4.67602569807196e+18 | 1038.286278747682 | 295268.585126299 | 6.556279633114204e-11 0.2 | 0.0002 | 1.0475901909837684 | 4.67602569807196e+18 | 1038.286278747682 | 147634.7925633495 | 3.278150918791789e-11 0.2 | 0.0005 | 1.0237060377513294 | 4.67602569807196e+18 | 1038.286278747682 | 59054.51702521031 | 1.3112736901901359e-11 0.2 | 0.00075 | 1.042861264113044 | 4.67602569807196e+18 | 1038.286278747682 | 39370.01135012432 | 8.741898616132353e-12 0.2 | 0.001 | 1.06754823612224 | 4.67602569807196e+18 | 1038.286278747682 | 29527.75851261451 | 6.556479473255219e-12 0.2 | 0.002 | 1.150251960798737 | 4.67602569807196e+18 | 1038.286278747682 | 14764.379256294498 | 3.2783507589272394e-12 0.2 | 0.005 | 1.261346117268915 | 4.67602569807196e+18 | 1038.286278747682 | 5906.351702518194 | 1.3114735303339383e-12 0.2 | 0.0075 | 1.2986164432111493 | 4.67602569807196e+18 | 1038.286278747682 | 3937.9011350116666 | 8.743897017574979e-13 0.2 | 0.01 | 1.3193653278045643 | 4.67602569807196e+18 | 1038.286278747682 | 2953.6758512586903 | 6.558477874693414e-13 0.2 | 0.02 | 1.3511101649720023 | 4.67602569807196e+18 | 1038.286278747682 | 1477.3379256293329 | 3.280349160371305e-13 0.2 | 0.05 | 1.3613657406180064 | 4.67602569807196e+18 | 1038.286278747682 | 591.5351702517372 | 1.3134719317780812e-13 0.2 | 0.075 | 1.3577244256580758 | 4.67602569807196e+18 | 1038.286278747682 | 394.69011350115574 | 8.763881032017989e-14 0.2 | 0.1 | 1.3532120690730092 | 4.67602569807196e+18 | 1038.286278747682 | 296.2675851258668 | 6.578461889136617e-14 0.2 | 0.2 | 1.3396807959929744 | 4.67602569807196e+18 | 1038.286278747682 | 148.63379256293308 | 3.300333174814553e-14 0.2 | 0.5 | 1.3258785979549306 | 4.67602569807196e+18 | 1038.286278747682 | 60.053517025172994 | 1.3334559462213179e-14 0.2 | 0.75 | 1.322792798552994 | 4.67602569807196e+18 | 1038.286278747682 | 40.36901135011535 | 8.963721176450468e-15 0.2 | 1 | 1.3215357179523244 | 4.67602569807196e+18 | 1038.286278747682 | 30.5267585125865 | 6.778302033569106e-15 0.5 | 1e-05 | 0.2409679449189826 | 986526591194416.2 | 0.21905290718980203 | 2725466.7391109765 | 6.051751853222101e-10 0.5 | 2e-05 | 0.2449083880329806 | 986526591194416.2 | 0.21905290718980203 | 1362733.8714743387 | 3.0258770410947793e-10 0.5 | 5e-05 | 0.26603719418856653 | 986526591194416.2 | 0.21905290718980203 | 545094.1490345697 | 1.2103521496932717e-10 0.5 | 7.5e-05 | 0.30414351709212 | 986526591194416.2 | 0.21905290718980203 | 363396.4327585239 | 8.069021734303215e-11 0.5 | 0.0001 | 0.346645595986304 | 986526591194416.2 | 0.21905290718980203 | 272547.57459431764 | 6.051771852407076e-11 0.5 | 0.0002 | 0.48965105730942426 | 986526591194416.2 | 0.21905290718980203 | 136274.28731608146 | 3.0258970288539514e-11 0.5 | 0.0005 | 0.7000258094500523 | 986526591194416.2 | 0.21905290718980203 | 54510.31493092727 | 1.2103721343176781e-11 0.5 | 0.00075 | 0.7803152404535806 | 986526591194416.2 | 0.21905290718980203 | 36340.543287970344 | 8.069221577138373e-12 0.5 | 0.001 | 0.8291864833507567 | 986526591194416.2 | 0.21905290718980203 | 27255.65746622731 | 6.051971694060423e-12 0.5 | 0.002 | 0.9148557625309567 | 986526591194416.2 | 0.21905290718980203 | 13628.328733301962 | 3.0260968693744864e-12 0.5 | 0.005 | 0.9605054768523519 | 986526591194416.2 | 0.21905290718980203 | 5451.931493366734 | 1.2105719745229523e-12 0.5 | 0.0075 | 0.9664986040888414 | 986526591194416.2 | 0.21905290718980203 | 3634.9543289177077 | 8.071219978850647e-13 0.5 | 0.01 | 0.9715770597539491 | 986526591194416.2 | 0.21905290718980203 | 2726.4657466908498 | 6.053970095656002e-13 0.5 | 0.02 | 1.0046221692212907 | 986526591194416.2 | 0.21905290718980203 | 1363.7328733473096 | 3.028095270856811e-13 0.5 | 0.05 | 1.1047795554227962 | 986526591194416.2 | 0.21905290718980203 | 546.093149339378 | 1.212570375973283e-13 0.5 | 0.075 | 1.1581424516578132 | 986526591194416.2 | 0.21905290718980203 | 364.3954328929862 | 8.091203993320888e-14 0.5 | 0.1 | 1.1943040342674949 | 986526591194416.2 | 0.21905290718980203 | 273.54657466976516 | 6.073954110114358e-14 0.5 | 0.2 | 1.2650403625958069 | 986526591194416.2 | 0.21905290718980203 | 137.27328733490157 | 3.048079285303852e-14 0.5 | 0.5 | 1.3129959855689675 | 986526591194416.2 | 0.21905290718980203 | 55.50931493396519 | 1.2325543904171442e-14 0.5 | 0.75 | 1.3205631535380162 | 986526591194416.2 | 0.21905290718980203 | 37.3395432893108 | 8.29104413775612e-15 0.5 | 1 | 1.3228809465642861 | 986526591194416.2 | 0.21905290718980203 | 28.254657466983357 | 6.273794254548405e-15 0.75 | 1e-05 | 0.22244997019985202 | 3532751603282.348 | 0.0007844284340490999 | 2442727.9074514243 | 5.423945531493999e-10 0.75 | 2e-05 | 0.20530606055190978 | 3532751603282.348 | 0.0007844284340490999 | 1221364.8760184385 | 2.71197481364824e-10 0.75 | 5e-05 | 0.2085533760825017 | 3532751603282.348 | 0.0007844284340490999 | 488546.65174448 | 1.0847914827404992e-10 0.75 | 7.5e-05 | 0.21981047089589112 | 3532751603282.348 | 0.0007844284340490999 | 325698.1161772274 | 7.231950953139941e-11 0.75 | 0.0001 | 0.22929183435193062 | 3532751603282.348 | 0.0007844284340490999 | 244273.84276236538 | 5.423968890968864e-11 0.75 | 0.0002 | 0.25206556710843314 | 3532751603282.348 | 0.0007844284340490999 | 122137.4256038195 | 2.7119956414760505e-11 0.75 | 0.0005 | 0.2911465816898953 | 3532751603282.348 | 0.0007844284340490999 | 48855.57125495263 | 1.0848116017692672e-11 0.75 | 0.00075 | 0.32523987895140793 | 3532751603282.348 | 0.0007844284340490999 | 32570.714320105762 | 7.232151393333944e-12 0.75 | 0.001 | 0.3599082896642641 | 3532751603282.348 | 0.0007844284340490999 | 24428.285796380867 | 5.4241690686531435e-12 0.75 | 0.002 | 0.4758051392903277 | 3532751603282.348 | 0.0007844284340490999 | 12214.64294041651 | 2.7121955660051067e-12 0.75 | 0.005 | 0.6515260879776189 | 3532751603282.348 | 0.0007844284340490999 | 4886.457186301037 | 1.0850114554152939e-12 0.75 | 0.0075 | 0.719735699372577 | 3532751603282.348 | 0.0007844284340490999 | 3257.9714590352833 | 7.234149854785173e-13 0.75 | 0.01 | 0.7617582568840169 | 3532751603282.348 | 0.0007844284340490999 | 2443.7285948394847 | 5.426167503851353e-13 0.75 | 0.02 | 0.8408992984214761 | 3532751603282.348 | 0.0007844284340490999 | 1222.364297842004 | 2.7141939758879106e-13 0.75 | 0.05 | 0.9210906761308859 | 3532751603282.348 | 0.0007844284340490999 | 489.5457192381459 | 1.0870098582097441e-13 0.75 | 0.075 | 0.9624492708631205 | 3532751603282.348 | 0.0007844284340490999 | 326.6971461737777 | 7.254133875229167e-14 0.75 | 0.1 | 0.9968207350967736 | 3532751603282.348 | 0.0007844284340490999 | 245.272859635963 | 5.446151521670006e-14 0.75 | 0.2 | 1.0921914353907196 | 3532751603282.348 | 0.0007844284340490999 | 123.13642982220428 | 2.7341779911750192e-14 0.75 | 0.5 | 1.2113887072395224 | 3532751603282.348 | 0.0007844284340490999 | 49.85457192989512 | 1.1069938727880118e-14 0.75 | 0.75 | 1.2502765485666585 | 3532751603282.348 | 0.0007844284340490999 | 33.56971462008021 | 7.453974020261758e-15 0.75 | 1 | 1.2715113409155707 | 3532751603282.348 | 0.0007844284340490999 | 25.427285965116457 | 5.645991666440077e-15 1.0 | 1e-05 | 0.32017904466341335 | 63914600437.045494 | 1.4191892202985e-05 | 2130369.113210278 | 4.730369680872655e-10 1.0 | 2e-05 | 0.30074023142713785 | 63914600437.045494 | 1.4191892202985e-05 | 1065202.8089999421 | 2.3652253688942574e-10 1.0 | 5e-05 | 0.2698729604691427 | 63914600437.045494 | 1.4191892202985e-05 | 426085.9842798409 | 9.461009404351037e-11 1.0 | 7.5e-05 | 0.25597138718717743 | 63914600437.045494 | 1.4191892202985e-05 | 284058.28740730474 | 6.307361020303598e-11 1.0 | 0.0001 | 0.2465204142749437 | 63914600437.045494 | 1.4191892202985e-05 | 213044.2022636677 | 4.730531572320456e-11 1.0 | 0.0002 | 0.22985421062285608 | 63914600437.045494 | 1.4191892202985e-05 | 106522.77866351408 | 2.365280830385654e-11 1.0 | 0.0005 | 0.23919407237622184 | 63914600437.045494 | 1.4191892202985e-05 | 42609.75407314652 | 9.461266009124563e-12 1.0 | 0.00075 | 0.25539979901779764 | 63914600437.045494 | 1.4191892202985e-05 | 28406.842361027593 | 6.3075860892220155e-12 1.0 | 0.001 | 0.26957943340881735 | 63914600437.045494 | 1.4191892202985e-05 | 21305.38413787569 | 4.730745603670637e-12 1.0 | 0.002 | 0.31063876807449725 | 63914600437.045494 | 1.4191892202985e-05 | 10653.193844260486 | 2.365484218338595e-12 1.0 | 0.005 | 0.39366095954964797 | 63914600437.045494 | 1.4191892202985e-05 | 4261.87796378237 | 9.463270087067532e-13 1.0 | 0.0075 | 0.44699159009612993 | 63914600437.045494 | 1.4191892202985e-05 | 2841.585372310853 | 6.309587013555114e-13 1.0 | 0.01 | 0.49016736300398767 | 63914600437.045494 | 1.4191892202985e-05 | 2131.439052904146 | 4.73274542423884e-13 1.0 | 0.02 | 0.601753244691417 | 63914600437.045494 | 1.4191892202985e-05 | 1066.2195442053355 | 2.367482974564207e-13 1.0 | 0.05 | 0.7405769938099279 | 63914600437.045494 | 1.4191892202985e-05 | 427.0878219429186 | 9.483254669160748e-14 1.0 | 0.075 | 0.7971634633262156 | 63914600437.045494 | 1.4191892202985e-05 | 285.0585485931724 | 6.329571280287381e-14 1.0 | 0.1 | 0.8372760302121763 | 63914600437.045494 | 1.4191892202985e-05 | 214.0439116815894 | 4.7527295805946814e-14 1.0 | 0.2 | 0.9403435142235637 | 63914600437.045494 | 1.4191892202985e-05 | 107.52195601832754 | 2.3874670244856132e-14 1.0 | 0.5 | 1.0888716539372492 | 63914600437.045494 | 1.4191892202985e-05 | 43.60878244993884 | 9.68309487035831e-15 1.0 | 0.75 | 1.1499813978771305 | 63914600437.045494 | 1.4191892202985e-05 | 29.40585497293815 | 6.529411449948819e-15 1.0 | 1 | 1.1884279444711836 | 63914600437.045494 | 1.4191892202985e-05 | 22.304391232070714 | 4.952569739218474e-15 2.0 | 1e-05 | 0.37717527079872715 | 4995157.201350917 | 1.1091477073123895e-09 | 931785.0004932241 | 2.0689783230958805e-10 2.0 | 2e-05 | 0.36554739044403184 | 4995157.201350917 | 1.1091477073123895e-09 | 513815.8738706699 | 1.1409004271782262e-10 2.0 | 5e-05 | 0.3556659370168733 | 4995157.201350917 | 1.1091477073123895e-09 | 219045.86744260634 | 4.86379530967543e-11 2.0 | 7.5e-05 | 0.35350368702102114 | 4995157.201350917 | 1.1091477073123895e-09 | 148197.11674392316 | 3.290637023843316e-11 2.0 | 0.0001 | 0.35281735349795634 | 4995157.201350917 | 1.1091477073123895e-09 | 111978.62460525152 | 2.4864249460521464e-11 2.0 | 0.0002 | 0.35343475255303497 | 4995157.201350917 | 1.1091477073123895e-09 | 56624.48346402158 | 1.257316105985264e-11 2.0 | 0.0005 | 0.355024241545291 | 4995157.201350917 | 1.1091477073123895e-09 | 22805.496043360723 | 5.063837359067396e-12 2.0 | 0.00075 | 0.3537163573870455 | 4995157.201350917 | 1.1091477073123895e-09 | 15227.168148584675 | 3.3811105356795045e-12 2.0 | 0.001 | 0.3515909417690454 | 4995157.201350917 | 1.1091477073123895e-09 | 11429.33503584049 | 2.53782218258902e-12 2.0 | 0.002 | 0.3437425980245205 | 4995157.201350917 | 1.1091477073123895e-09 | 5721.711678609502 | 1.2704752091717846e-12 2.0 | 0.005 | 0.33843782328439914 | 4995157.201350917 | 1.1091477073123895e-09 | 2290.8581499578336 | 5.086726928466753e-13 2.0 | 0.0075 | 0.3434911844399524 | 4995157.201350917 | 1.1091477073123895e-09 | 1527.8054038404043 | 3.392409472980705e-13 2.0 | 0.01 | 0.3514509386543811 | 4995157.201350917 | 1.1091477073123895e-09 | 1146.1915618887642 | 2.545056525279952e-13 2.0 | 0.02 | 0.387954531271237 | 4995157.201350917 | 1.1091477073123895e-09 | 573.6614252411332 | 1.273784245283978e-13 2.0 | 0.05 | 0.4710748508056674 | 4995157.201350917 | 1.1091477073123895e-09 | 230.08032761773444 | 5.1088095446901614e-14 2.0 | 0.075 | 0.5172647100633149 | 4995157.201350917 | 1.1091477073123895e-09 | 153.72255305569766 | 3.413326356131955e-14 2.0 | 0.1 | 0.5526131904880446 | 4995157.201350917 | 1.1091477073123895e-09 | 115.5427903033019 | 2.565565322483241e-14 2.0 | 0.2 | 0.6463104722195474 | 4995157.201350917 | 1.1091477073123895e-09 | 58.272051797845364 | 1.2938994719623534e-14 2.0 | 0.5 | 0.7972439807925047 | 4995157.201350917 | 1.1091477073123895e-09 | 23.908978317115814 | 5.308859644585121e-15 2.0 | 0.75 | 0.8774340306364913 | 4995157.201350917 | 1.1091477073123895e-09 | 16.27267555945542 | 3.613259815672491e-15 2.0 | 1 | 0.9386486865700187 | 4995157.201350917 | 1.1091477073123895e-09 | 12.454515425128461 | 2.7654579571053575e-15 5.0 | 1e-05 | 0.7290963744458004 | 555.6368040064718 | 1.233761546274241e-13 | 554.9108780960489 | 1.232149666954394e-13 5.0 | 2e-05 | 0.7291162264236534 | 555.6368040064718 | 1.233761546274241e-13 | 554.1868499306446 | 1.230542001474976e-13 5.0 | 5e-05 | 0.7291757523892824 | 555.6368040064718 | 1.233761546274241e-13 | 552.0260777777686 | 1.2257441234847923e-13 5.0 | 7.5e-05 | 0.7292253230382069 | 555.6368040064718 | 1.233761546274241e-13 | 550.2382870339221 | 1.2217744305907321e-13 5.0 | 0.0001 | 0.7292748625029511 | 555.6368040064718 | 1.233761546274241e-13 | 548.4620596604188 | 1.2178304135366665e-13 5.0 | 0.0002 | 0.7294727088701934 | 555.6368040064718 | 1.233761546274241e-13 | 541.4705687792582 | 1.2023061852312238e-13 5.0 | 0.0005 | 0.7300632679255016 | 555.6368040064718 | 1.233761546274241e-13 | 521.5279990898892 | 1.1580247851525655e-13 5.0 | 0.00075 | 0.7305520056062914 | 555.6368040064718 | 1.233761546274241e-13 | 505.99987151525244 | 1.1235454156272083e-13 5.0 | 0.001 | 0.7310376814323818 | 555.6368040064718 | 1.233761546274241e-13 | 491.37136147684004 | 1.091063598305997e-13 5.0 | 0.002 | 0.7329502858896298 | 555.6368040064718 | 1.233761546274241e-13 | 440.45228221100956 | 9.780005299187201e-14 5.0 | 0.005 | 0.7384152522426066 | 555.6368040064718 | 1.233761546274241e-13 | 336.07266536594227 | 7.46231222072829e-14 5.0 | 0.0075 | 0.7426859104004214 | 555.6368040064718 | 1.233761546274241e-13 | 280.7085388287963 | 6.232981660332288e-14 5.0 | 0.01 | 0.746729529407029 | 555.6368040064718 | 1.233761546274241e-13 | 241.04573230445396 | 5.3522904398407337e-14 5.0 | 0.02 | 0.761091174552371 | 555.6368040064718 | 1.233761546274241e-13 | 154.1683480263636 | 3.4232249929458633e-14 5.0 | 0.05 | 0.7938460924485605 | 555.6368040064718 | 1.233761546274241e-13 | 74.43525519938765 | 1.6527946833241914e-14 5.0 | 0.075 | 0.8150678569112069 | 555.6368040064718 | 1.233761546274241e-13 | 52.217266828137994 | 1.1594562383118844e-14 5.0 | 0.1 | 0.8333792535955515 | 555.6368040064718 | 1.233761546274241e-13 | 40.32070591377203 | 8.952995214921884e-15 5.0 | 0.2 | 0.891988844835201 | 555.6368040064718 | 1.233761546274241e-13 | 21.382869766048387 | 4.747950869365611e-15 5.0 | 0.5 | 1.0120030229593937 | 555.6368040064718 | 1.233761546274241e-13 | 9.336977975122132 | 2.0732255856797125e-15 5.0 | 0.75 | 1.0827038377848377 | 555.6368040064718 | 1.233761546274241e-13 | 6.585973686347973 | 1.4623799252317878e-15 5.0 | 1 | 1.1385693673981578 | 555.6368040064718 | 1.233761546274241e-13 | 5.20005538287452 | 1.1546442430786553e-15 7.5 | 1e-05 | 1.157972716294107 | 65.4206227436607 | 1.4526296331065658e-14 | 65.40750197167729 | 1.452338293443429e-14 7.5 | 2e-05 | 1.157977684282373 | 65.4206227436607 | 1.4526296331065658e-14 | 65.39438654331285 | 1.4520470724324685e-14 7.5 | 5e-05 | 1.1579925874856083 | 65.4206227436607 | 1.4526296331065658e-14 | 65.35507228730174 | 1.4511741205880777e-14 7.5 | 7.5e-05 | 1.158005005949029 | 65.4206227436607 | 1.4526296331065658e-14 | 65.32234705876327 | 1.450447474543887e-14 7.5 | 0.0001 | 1.1580174236191763 | 65.4206227436607 | 1.4526296331065658e-14 | 65.28965509555222 | 1.449721567138345e-14 7.5 | 0.0002 | 1.158067086367233 | 65.4206227436607 | 1.4526296331065658e-14 | 65.15921888360955 | 1.4468253014234723e-14 7.5 | 0.0005 | 1.1582159984650573 | 65.4206227436607 | 1.4526296331065658e-14 | 64.77106671380922 | 1.4382065919040614e-14 7.5 | 0.00075 | 1.1583400046414274 | 65.4206227436607 | 1.4526296331065658e-14 | 64.45117618545733 | 1.431103595305346e-14 7.5 | 0.001 | 1.158463931527119 | 65.4206227436607 | 1.4526296331065658e-14 | 64.13447892974612 | 1.4240715036028223e-14 7.5 | 0.002 | 1.1589588465897571 | 65.4206227436607 | 1.4526296331065658e-14 | 62.89868505684149 | 1.3966313673750338e-14 7.5 | 0.005 | 1.1604360020625986 | 65.4206227436607 | 1.4526296331065658e-14 | 59.46547650580334 | 1.3203988237409833e-14 7.5 | 0.0075 | 1.1616583144933834 | 65.4206227436607 | 1.4526296331065658e-14 | 56.8825404345771 | 1.2630461217927791e-14 7.5 | 0.01 | 1.162872831495904 | 65.4206227436607 | 1.4526296331065658e-14 | 54.51817052965834 | 1.2105465636493471e-14 7.5 | 0.02 | 1.1676546346371224 | 65.4206227436607 | 1.4526296331065658e-14 | 46.77181528449089 | 1.0385429246471321e-14 7.5 | 0.05 | 1.1813220595206992 | 65.4206227436607 | 1.4526296331065658e-14 | 32.91391372384854 | 7.308356969348515e-15 7.5 | 0.075 | 1.1920256826753837 | 65.4206227436607 | 1.4526296331065658e-14 | 26.484230061148875 | 5.88068040067144e-15 7.5 | 0.1 | 1.2021902002798597 | 65.4206227436607 | 1.4526296331065658e-14 | 22.21088199657382 | 4.931806517965725e-15 7.5 | 0.2 | 1.238514238689644 | 65.4206227436607 | 1.4526296331065658e-14 | 13.695471638717011 | 3.04100558928089e-15 7.5 | 0.5 | 1.3202321940460673 | 65.4206227436607 | 1.4526296331065658e-14 | 6.759169742060504 | 1.5008371749970503e-15 7.5 | 0.75 | 1.3692070811077155 | 65.4206227436607 | 1.4526296331065658e-14 | 4.957375592564131 | 1.1007585049158954e-15 7.5 | 1 | 1.4074503211972995 | 65.4206227436607 | 1.4526296331065658e-14 | 4.014324376831694 | 8.91359070294516e-16 10.0 | 1e-05 | 1.4570968886841438 | 19.175864900543605 | 4.2578973459369795e-15 | 19.17448930484128 | 4.25759190233272e-15 10.0 | 2e-05 | 1.4570993187331154 | 19.175864900543605 | 4.2578973459369795e-15 | 19.173113917340398 | 4.257286504958468e-15 10.0 | 5e-05 | 1.4571066086154154 | 19.175864900543605 | 4.2578973459369795e-15 | 19.168989003573497 | 4.256370590110747e-15 10.0 | 7.5e-05 | 1.4571126832141403 | 19.175864900543605 | 4.2578973459369795e-15 | 19.165553005574086 | 4.255607645292444e-15 10.0 | 0.0001 | 1.4571187575372515 | 19.175864900543605 | 4.2578973459369795e-15 | 19.16211830691505 | 4.254844988985662e-15 10.0 | 0.0002 | 1.4571430520739157 | 19.175864900543605 | 4.2578973459369795e-15 | 19.1483924909525 | 4.251797245602984e-15 10.0 | 0.0005 | 1.4572159092372137 | 19.175864900543605 | 4.2578973459369795e-15 | 19.10733923905199 | 4.242681592503847e-15 10.0 | 0.00075 | 1.4572765932520455 | 19.175864900543605 | 4.2578973459369795e-15 | 19.073269807537834 | 4.235116659043266e-15 10.0 | 0.001 | 1.4573372497498747 | 19.175864900543605 | 4.2578973459369795e-15 | 19.039328340305325 | 4.227580139361048e-15 10.0 | 0.002 | 1.4575796008997441 | 19.175864900543605 | 4.2578973459369795e-15 | 18.904827802739543 | 4.1977150206350496e-15 10.0 | 0.005 | 1.458304024876359 | 19.175864900543605 | 4.2578973459369795e-15 | 18.51309687468338 | 4.110733281477903e-15 10.0 | 0.0075 | 1.4589047145245595 | 19.175864900543605 | 4.2578973459369795e-15 | 18.199513783630792 | 4.041103847913961e-15 10.0 | 0.01 | 1.4595026976623582 | 19.175864900543605 | 4.2578973459369795e-15 | 17.896962960525293 | 3.9739240699277574e-15 10.0 | 0.02 | 1.461867908452495 | 19.175864900543605 | 4.2578973459369795e-15 | 16.78620389115527 | 3.727286011202595e-15 10.0 | 0.05 | 1.4687164362187255 | 19.175864900543605 | 4.2578973459369795e-15 | 14.185809264834253 | 3.1498824137519706e-15 10.0 | 0.075 | 1.4741565641494006 | 19.175864900543605 | 4.2578973459369795e-15 | 12.59424951147128 | 2.7964851571019087e-15 10.0 | 0.1 | 1.479371019681004 | 19.175864900543605 | 4.2578973459369795e-15 | 11.345519722685243 | 2.5192114444927955e-15 10.0 | 0.2 | 1.4982457886579086 | 19.175864900543605 | 4.2578973459369795e-15 | 8.230533331790687 | 1.8275455219797647e-15 10.0 | 0.5 | 1.5410098659533928 | 19.175864900543605 | 4.2578973459369795e-15 | 4.79897439707101 | 1.0655863740429728e-15 10.0 | 0.75 | 1.566357435218295 | 19.175864900543605 | 4.2578973459369795e-15 | 3.7223151344244276 | 8.265199934297368e-16 10.0 | 1 | 1.5858716852811872 | 19.175864900543605 | 4.2578973459369795e-15 | 3.1211613132222813 | 6.930370307017333e-16 20.0 | 1e-05 | 1.7286486756519486 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1428984173511347 | 4.758190324752076e-16 20.0 | 2e-05 | 1.7286487644074315 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1428806003827185 | 4.758150763134946e-16 20.0 | 5e-05 | 1.7286490306649474 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1428271528103546 | 4.758032085684049e-16 20.0 | 7.5e-05 | 1.728649252535983 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1427826169852673 | 4.757933196287184e-16 20.0 | 0.0001 | 1.7286494743977214 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1427380846311563 | 4.757834314597434e-16 20.0 | 0.0002 | 1.7286503617517068 | 2.1429162348750763 | 4.758229887602719e-16 | 2.142559989916348 | 4.757438864891545e-16 20.0 | 0.0005 | 1.7286530229215018 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1420260386871086 | 4.756253254994089e-16 20.0 | 0.00075 | 1.728655239541275 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1415814604086187 | 4.755266092912033e-16 20.0 | 0.001 | 1.7286574552328238 | 2.1429162348750763 | 4.758229887602719e-16 | 2.141137228134385 | 4.754279699113761e-16 20.0 | 0.002 | 1.7286663087282645 | 2.1429162348750763 | 4.758229887602719e-16 | 2.139363751013921 | 4.750341788848192e-16 20.0 | 0.005 | 1.728692780526595 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1340762407576834 | 4.738601157589358e-16 20.0 | 0.0075 | 1.728714739279412 | 2.1429162348750763 | 4.758229887602719e-16 | 2.129707326476959 | 4.72890021913521e-16 20.0 | 0.01 | 1.7287366067447667 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1253719446187818 | 4.719273737616229e-16 20.0 | 0.02 | 1.7288231748695115 | 2.1429162348750763 | 4.758229887602719e-16 | 2.1083581370698803 | 4.681495495861566e-16 20.0 | 0.05 | 1.7290745149199336 | 2.1429162348750763 | 4.758229887602719e-16 | 2.0602694484403035 | 4.574717157180393e-16 20.0 | 0.075 | 1.7292748667352151 | 2.1429162348750763 | 4.758229887602719e-16 | 2.0232719267437416 | 4.49256615629721e-16 20.0 | 0.1 | 1.7294674543189974 | 2.1429162348750763 | 4.758229887602719e-16 | 1.9887693614232316 | 4.4159550714422826e-16 20.0 | 0.2 | 1.7301684663530894 | 2.1429162348750763 | 4.758229887602719e-16 | 1.8712610423214955 | 4.1550341885387875e-16 20.0 | 0.5 | 1.731772516977182 | 2.1429162348750763 | 4.758229887602719e-16 | 1.6422723278704805 | 3.6465771022131232e-16 20.0 | 0.75 | 1.7327288355012884 | 2.1429162348750763 | 4.758229887602719e-16 | 1.526875789694727 | 3.3903453149236083e-16 20.0 | 1 | 1.7334658393718365 | 2.1429162348750763 | 4.758229887602719e-16 | 1.446630124465455 | 3.212164144595808e-16 50.0 | 1e-05 | 1.7391411403978454 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013610614050856 | 2.2234682126700225e-16 50.0 | 2e-05 | 1.7391411404618256 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013610477854813 | 2.223468182428426e-16 50.0 | 5e-05 | 1.7391411406537576 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013610069283048 | 2.2234680917072697e-16 50.0 | 7.5e-05 | 1.739141140813692 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013609728825315 | 2.223468016110467e-16 50.0 | 0.0001 | 1.7391411409736186 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013609388384612 | 2.2234679405174456e-16 50.0 | 0.0002 | 1.7391411416132443 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013608026792107 | 2.2234676381831757e-16 50.0 | 0.0005 | 1.739141143531356 | 1.001361075024962 | 2.223468242912223e-16 | 1.00136039436488 | 2.223466731543233e-16 50.0 | 0.00075 | 1.7391411451289047 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013600542900682 | 2.223465976425461e-16 50.0 | 0.001 | 1.739141146725655 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013597143852402 | 2.2234652216851285e-16 50.0 | 0.002 | 1.7391411531046923 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013583564632107 | 2.223462206492523e-16 50.0 | 0.005 | 1.7391411721656653 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013542989220954 | 2.2234531969413837e-16 50.0 | 0.0075 | 1.7391411879631309 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013509361089248 | 2.2234457299961647e-16 50.0 | 0.01 | 1.7391412036824188 | 1.001361075024962 | 2.223468242912223e-16 | 1.0013475899545612 | 2.223438300040928e-16 50.0 | 0.02 | 1.7391412657893082 | 1.001361075024962 | 2.223468242912223e-16 | 1.001334369473718 | 2.2234089446764705e-16 50.0 | 0.05 | 1.7391414450146108 | 1.001361075024962 | 2.223468242912223e-16 | 1.001296219922522 | 2.2233242356562365e-16 50.0 | 0.075 | 1.7391415867303404 | 1.001361075024962 | 2.223468242912223e-16 | 1.001266056190751 | 2.223257258717195e-16 50.0 | 0.1 | 1.739141722006669 | 1.001361075024962 | 2.223468242912223e-16 | 1.0012372643858902 | 2.223193328067841e-16 50.0 | 0.2 | 1.739142206764682 | 1.001361075024962 | 2.223468242912223e-16 | 1.0011341005544734 | 2.222964258345946e-16 50.0 | 0.5 | 1.7391432733304766 | 1.001361075024962 | 2.223468242912223e-16 | 1.0009071775605212 | 2.222460388080541e-16 50.0 | 0.75 | 1.7391438828572479 | 1.001361075024962 | 2.223468242912223e-16 | 1.0007775303834385 | 2.222172513518391e-16 50.0 | 1 | 1.739144340031251 | 1.001361075024962 | 2.223468242912223e-16 | 1.0006803060255958 | 2.2219566320771284e-16 75.0 | 1e-05 | 1.7391475402804035 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001490198573 | 2.2204463801408666e-16 75.0 | 2e-05 | 1.7391475402804109 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001490183668 | 2.220446380137557e-16 75.0 | 5e-05 | 1.7391475402804306 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001490138963 | 2.2204463801276305e-16 75.0 | 7.5e-05 | 1.739147540280447 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001490101718 | 2.2204463801193603e-16 75.0 | 0.0001 | 1.739147540280464 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001490064465 | 2.2204463801110886e-16 75.0 | 0.0002 | 1.7391475402805299 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001489915489 | 2.220446380078009e-16 75.0 | 0.0005 | 1.7391475402807293 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001489468737 | 2.2204463799788104e-16 75.0 | 0.00075 | 1.7391475402808945 | 1.0000001490213473 | 2.220446380144175e-16 | 1.000000148909665 | 2.2204463798961905e-16 75.0 | 0.001 | 1.7391475402810592 | 1.0000001490213473 | 2.220446380144175e-16 | 1.000000148872475 | 2.220446379813612e-16 75.0 | 0.002 | 1.7391475402817214 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001487238992 | 2.2204463794837075e-16 75.0 | 0.005 | 1.7391475402836962 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001482799472 | 2.220446378497936e-16 75.0 | 0.0075 | 1.7391475402853336 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001479120069 | 2.2204463776809444e-16 75.0 | 0.01 | 1.7391475402869623 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001475458884 | 2.220446376867998e-16 75.0 | 0.02 | 1.7391475402933996 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001460993597 | 2.220446373656059e-16 75.0 | 0.05 | 1.7391475403119736 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001419250921 | 2.220446364387323e-16 75.0 | 0.075 | 1.7391475403266603 | 1.0000001490213473 | 2.220446380144175e-16 | 1.000000138624508 | 2.2204463570585544e-16 75.0 | 0.1 | 1.7391475403406795 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001354739514 | 2.220446350062913e-16 75.0 | 0.2 | 1.739147540390915 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000001241844543 | 2.220446324995194e-16 75.0 | 0.5 | 1.7391475405014307 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000000993475626 | 2.220446269846216e-16 75.0 | 0.75 | 1.7391475405645824 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000000851550528 | 2.2204462383325136e-16 75.0 | 1 | 1.7391475406119468 | 1.0000001490213473 | 2.220446380144175e-16 | 1.0000000745106707 | 2.2204462146972376e-16 100.0 | 1e-05 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.000000000000426 | 2.220446049251259e-16 100.0 | 2e-05 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.000000000000426 | 2.220446049251259e-16 100.0 | 5e-05 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004259 | 2.2204460492512587e-16 100.0 | 7.5e-05 | 1.739147540943495 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004259 | 2.2204460492512587e-16 100.0 | 0.0001 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004259 | 2.2204460492512587e-16 100.0 | 0.0002 | 1.739147540943495 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004259 | 2.2204460492512587e-16 100.0 | 0.0005 | 1.739147540943495 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004257 | 2.220446049251258e-16 100.0 | 0.00075 | 1.7391475409434956 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004259 | 2.2204460492512587e-16 100.0 | 0.001 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004254 | 2.2204460492512577e-16 100.0 | 0.002 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004252 | 2.2204460492512572e-16 100.0 | 0.005 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004237 | 2.220446049251254e-16 100.0 | 0.0075 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004228 | 2.220446049251252e-16 100.0 | 0.01 | 1.739147540943495 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004219 | 2.22044604925125e-16 100.0 | 0.02 | 1.7391475409434956 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004174 | 2.22044604925124e-16 100.0 | 0.05 | 1.739147540943495 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000004057 | 2.220446049251214e-16 100.0 | 0.075 | 1.7391475409434949 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000003963 | 2.220446049251193e-16 100.0 | 0.1 | 1.7391475409434953 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000003872 | 2.220446049251173e-16 100.0 | 0.2 | 1.7391475409434958 | 1.000000000000426 | 2.220446049251259e-16 | 1.000000000000355 | 2.2204460492511014e-16 100.0 | 0.5 | 1.7391475409434958 | 1.000000000000426 | 2.220446049251259e-16 | 1.000000000000284 | 2.2204460492509437e-16 100.0 | 0.75 | 1.7391475409434956 | 1.000000000000426 | 2.220446049251259e-16 | 1.0000000000002436 | 2.220446049250854e-16 100.0 | 1 | 1.7391475409434962 | 1.000000000000426 | 2.220446049251259e-16 | 1.000000000000213 | 2.220446049250786e-16 200.0 | 1e-05 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000007 | 2.2204460492503146e-16 200.0 | 2e-05 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000004 | 2.220446049250314e-16 200.0 | 5e-05 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 7.5e-05 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000007 | 2.2204460492503146e-16 200.0 | 0.0001 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.0002 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.0005 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.00075 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000004 | 2.220446049250314e-16 200.0 | 0.001 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.002 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.005 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000004 | 2.220446049250314e-16 200.0 | 0.0075 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.01 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000004 | 2.220446049250314e-16 200.0 | 0.02 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.05 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.075 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000004 | 2.220446049250314e-16 200.0 | 0.1 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.2 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000004 | 2.220446049250314e-16 200.0 | 0.5 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 0.75 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000002 | 2.2204460492503136e-16 200.0 | 1 | 1.739147540943497 | 1.0000000000000002 | 2.2204460492503136e-16 | 1.0000000000000004 | 2.220446049250314e-16 500.0 | 1e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 2e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 5e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 7.5e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.0001 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.0002 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.0005 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.00075 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.001 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.002 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.005 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.0075 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.01 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.02 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.075 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.1 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.2 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.5 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 0.75 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 500.0 | 1 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 1e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 2e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 5e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 7.5e-05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.0001 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.0002 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.0005 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.00075 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.001 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.002 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.005 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.0075 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.01 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.02 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.05 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.075 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.1 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.2 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.5 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 0.75 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 1000.0 | 1 | 1.739147540943497 | 1.0 | 2.220446049250313e-16 | 1.0 | 2.220446049250313e-16 Results saved in solution.pickle
Out[19]:
======================================================================================================================== Results of radial basis function run: Basis function type : gaussian Shape parameter : 0.75 Regularization parameter : 2e-05 Number of terms in RBF model : 31 RBF Expression: -------------------------- 4.038157591240296 + 237.8499310618983*(-1551.5244569724653*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.5518182560031254)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.2731190901874679)**2)**0.5)**2) -334.71891272365644*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.4300326512391472)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.8610629174234938)**2)**0.5)**2) + 235.58220306279952*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.11035664436222248)**2)**0.5)**2) + 27.886153329442664*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.6822668136056202)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.8414199967160326)**2)**0.5)**2) -1415.9829578480585*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.13520535030959882)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5811256814195772)**2)**0.5)**2) -990.5466782787646*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.908216841775921)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.20132180755786522)**2)**0.5)**2) -343.1090937508468*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.18501228557936236)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.09666279198240106)**2)**0.5)**2) -723.9922784804203*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.22021941796179972)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -1.0)**2)**0.5)**2) -8.525226874282353*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.8266365897830592)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.16264387978612238)**2)**0.5)**2) + 255.07386591870466*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.831288905218487)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.26866011456766936)**2)**0.5)**2) + 528.0399010676083*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.43739443999592303)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.9599385359524552)**2)**0.5)**2) + 432.7400228938864*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.8327482651149929)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.3330572615505772)**2)**0.5)**2) -761.4482838759395*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.17485274298531528)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.3711736045917253)**2)**0.5)**2) + 854.420692794522*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.0009932894634358945)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.2461895993357224)**2)**0.5)**2) + 158.4825871481679*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.8102253955808789)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546)**2)**0.5)**2) -2542.5018400507633*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.6086081776493459)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.6109286075828757)**2)**0.5)**2) -1366.5543839281318*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.1028770983224648)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.38063145554753214)**2)**0.5)**2) + 342.1406582035254*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.0325313822086792)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.9084335925839453)**2)**0.5)**2) + 107.60917336005593*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -1.0)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.046386199083146756)**2)**0.5)**2) + 49.67120678308112*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.9074219145138307)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5830026053427076)**2)**0.5)**2) + 897.301866820119*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.7557460774371298)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.6383110566137417)**2)**0.5)**2) -8.627269377401944*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.5911925356865033)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.005381283583344741)**2)**0.5)**2) -136.91234826515438*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.21031274399063923)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5495609776179707)**2)**0.5)**2) + 735.7407100550239*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.7830308832306688)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.24439265066740792)**2)**0.5)**2) + 541.7309338027707*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.28803140536417743)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.8689715765112118)**2)**0.5)**2) -45.9213997601471*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.993941455353761)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.9026639746834018)**2)**0.5)**2) + 1712.678844540206*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.36343808704287006)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.6057949629205177)**2)**0.5)**2) + 3266.9168170432567*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.3586109792561233)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.3372905989731323)**2)**0.5)**2) -827.2715689735041*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.17758838147122427)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.23089727866378493)**2)**0.5)**2) + 919.032553083629*exp(- (0.75*(((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.04112204875223326)**2 + ((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5088149538788017)**2)**0.5)**2)) -------------------------- Model training errors: ----------------------- Mean Squared Error (MSE) : 0.00044438098274641986 Root Mean Squared Error (RMSE) : 0.02108034588773201 Goodness of fit (R2) : 0.9918215668806268 ========================================================================================================================
Step 5: View model metrics¶
We can view the Root Mean Square Error (RMSE) and $R^2$ values of the RBF fit based on the training data:
In [20]:
print('R2: ', rbf_class.R2, '\nRMSE: ', rbf_class.rmse)
R2: 0.9918215668806268 RMSE: 0.02108034588773201
Step 6 (Optional): Generate Pyomo expression output at any unsampled point¶
Based on the model we trained, we can predict the surrogate value at any previously unsampled point. We do this by calling the function predict_output.
Let us again evaluate the RBF surrogate at the same set of points we considered for the polynomial model:
- $x_{1}=5$, $x_{2}=8$ (true function value: 57.9908)
- $x_{1}=-3$, $x_{2}=10$ (true function value: 4.2461)
- $x_{1}=-2$, $x_{2}=3$ (true function value: 50.8899)
In [21]:
unsampled_points = np.array([ [5, 8], [-3, 10], [-2, 3]])
ys = rbf_class.predict_output(unsampled_points)
print(ys)
[[61.06401281] [12.7142377 ] [50.9750144 ]]
The results from the RBF surrogate are similar to those obtained from the polynomial.
For RBF models, the Pyomo expression is generated by calling generate_expression on the results object, while parity plots may be viewed with the parity_residual_plots method.
Further information about using PySMO's RBF tool and features can be found in the documentation.
Training Kriging models¶
The KrigingModel class trains Kriging from data. For details about Kriging models, users should consult the documentation.
As an example, we will again consider the Brainin function. The same dataset loaded previously will be used.
Step 1: Load the data and import the Kriging tool¶
In [22]:
import numpy as np
brainin_data = np.loadtxt('brainin_30.txt')
from idaes.core.surrogate.pysmo.kriging import KrigingModel
Step 2: Specify the Kriging settings and initialize the KrigingModel class¶
The KrigingModel class takes a number of keyword arguments:
- XY_data : The dataset for Kriging training. training_data is expected to contain xy_data,
with the output values (y) in the last column.
It also takes a number of optional arguments:
- regularization : Boolean variable determining whether regularization is done. Default is True.
- numerical_gradients : Boolean variable which determines whether numerical gradients are used when
solving the max. likelihood optimization problem. Default is True.
- fname : Filename for saving (.pickle extension)
- overwrite : Option determining whether any existing file with the same name supplied in 'fname'
should be overwritten.
For this demonstration, we will train a Kriging model with regularization:
In [23]:
krg_class = KrigingModel(XY_data=brainin_data, overwrite=True)
Warning: solution.pickle already exists; previous file will be overwritten.
Step 3: Extract variable names (optional)¶
Next, we extract Pyomo variable names from the dataset.
In [24]:
vars = krg_class.get_feature_vector()
Step 4: Train the Kriging surrogate¶
Next, we train the RBF surrogate by calling training:
In [25]:
krg_class.training()
Optimizing kriging parameters using L-BFGS-B algorithm... Final results ================ Theta: [5.3960043 0.21786916] Mean: [[399.8740733]] Regularization parameter: 1.000000000001e-06 Results saved in solution.pickle
Out[25]:
======================================================================================================================== Results of Kriging run: Kriging mean : [[399.8740733]] Kriging variance : [[91163.53616269]] Kriging weights : [5.3960043 0.21786916] Regularization parameter : 1.000000000001e-06 Number of terms in Kriging model : 31 Kriging Expression: -------------------- 399.8740733017947 + (-26677.086861269083*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.5518182560031254)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.2731190901874679)**2)) -179468.18427625299*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.4300326512391472)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.8610629174234938)**2)) + 100547.57581360638*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.11035664436222248)**2)) -33262.15678686742*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.6822668136056202)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.8414199967160326)**2)) + 95941.67298658937*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.13520535030959882)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5811256814195772)**2)) + 142684.7207657257*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.908216841775921)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.20132180755786522)**2)) + 13684.313556239009*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.18501228557936236)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.09666279198240106)**2)) + 6889.785981370136*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.22021941796179972)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -1.0)**2)) -121336.95087997802*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.8266365897830592)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.16264387978612238)**2)) -116368.078969853*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.831288905218487)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.26866011456766936)**2)) + 137345.35259931907*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.43739443999592303)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.9599385359524552)**2)) -120650.50247030333*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.8327482651149929)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.3330572615505772)**2)) + 115011.32720454037*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.17485274298531528)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.3711736045917253)**2)) -136216.06318300962*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.0009932894634358945)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.2461895993357224)**2)) + 15840.168997440487*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.8102253955808789)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546)**2)) -16618.30784879136*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.6086081776493459)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.6109286075828757)**2)) -119490.16069413722*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.1028770983224648)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.38063145554753214)**2)) -21048.803165102378*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.0325313822086792)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.9084335925839453)**2)) -20287.55144619965*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -1.0)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.046386199083146756)**2)) -13917.940920275636*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.9074219145138307)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5830026053427076)**2)) + 84236.79319164483*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.7557460774371298)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.6383110566137417)**2)) -5640.99826487992*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.5911925356865033)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.005381283583344741)**2)) -126632.90425059432*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.21031274399063923)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5495609776179707)**2)) + 208601.0077510439*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.7830308832306688)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.24439265066740792)**2)) -31992.625017974526*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.28803140536417743)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.8689715765112118)**2)) + 2004.0939995171502*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.993941455353761)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.9026639746834018)**2)) + 112631.0142006306*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.36343808704287006)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.6057949629205177)**2)) + 11001.494287014008*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.3586109792561233)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.3372905989731323)**2)) -58592.92022290826*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.17758838147122427)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.23089727866378493)**2)) + 101781.91392371058*exp(- (5.396004300925354*((IndexedParam[0] + 4.929217157135412)/14.643030012320084 -0.04112204875223326)**2 + 0.21786916191639816*((IndexedParam[1] -0.22882456869508516)/14.45053220191546 -0.5088149538788017)**2))) -------------------------- Model training errors: ----------------------- Mean Squared Error (MSE) : 0.009252009602248907 Root Mean Squared Error (RMSE) : 0.0961873671655946 Goodness of fit (R2) : 0.999996990151358 ========================================================================================================================
The returned information correspond to the Kriging model parameters.
As we can see, the optimization problem was solved using SciPy's L-BFGS algorithm which makes use of gradients. A different algorithm (Basinhopping) is used when no numerical gradients are computed (when numerical_gradients is set to False). The user should try this.
Step 5: View model metrics¶
We can view the RMSE and $R^2$ values of the Kriging fit based on the training data:
In [26]:
print('R2: ', krg_class.training_R2, '\nRMSE: ', krg_class.training_rmse)
R2: 0.999996990151358 RMSE: 0.0961873671655946
Step 6 (Optional): Generate Pyomo expression output at any unsampled point¶
Again, based on the model we trained, we evaluate the surrogate at a set of off-design points:
- $x_{1}=5$, $x_{2}=8$ (true function value: 57.9908)
- $x_{1}=-3$, $x_{2}=10$ (true function value: 4.2461)
- $x_{1}=-2$, $x_{2}=3$ (true function value: 50.8899)
We do this by calling the function predict_output:
In [27]:
unsampled_points = np.array([ [5, 8], [-3, 10], [-2, 3]])
ys = krg_class.predict_output(unsampled_points)
print(ys)
[[57.41251715] [ 4.38081976] [50.82440456]]
The Kriging model performs very well, predicting all three points fairly accurately.
For Kriging models, the Pyomo expression is generated by calling generate_expression on the results object:
In [28]:
print(krg_class.generate_expression([m.x[1], m.x[2]]))
-26677.086861269083*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.5518182560031254)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.2731190901874679)**2)) - 179468.18427625299*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.4300326512391472)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.8610629174234938)**2)) + 100547.57581360638*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.11035664436222248)**2)) - 33262.15678686742*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.6822668136056202)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.8414199967160326)**2)) + 95941.67298658937*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.13520535030959882)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.5811256814195772)**2)) + 142684.7207657257*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.908216841775921)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.20132180755786522)**2)) + 13684.313556239009*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.18501228557936236)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.09666279198240106)**2)) + 6889.785981370136*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.22021941796179972)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 1.0)**2)) - 121336.95087997802*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.8266365897830592)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.16264387978612238)**2)) - 116368.078969853*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.831288905218487)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.26866011456766936)**2)) + 137345.35259931907*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.43739443999592303)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.9599385359524552)**2)) - 120650.50247030333*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.8327482651149929)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.3330572615505772)**2)) + 115011.32720454037*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.17485274298531528)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.3711736045917253)**2)) - 136216.06318300962*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.0009932894634358945)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.2461895993357224)**2)) + 15840.168997440487*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.8102253955808789)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546)**2)) - 16618.30784879136*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.6086081776493459)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.6109286075828757)**2)) - 119490.16069413722*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.1028770983224648)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.38063145554753214)**2)) - 21048.803165102378*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.0325313822086792)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.9084335925839453)**2)) - 20287.55144619965*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 1.0)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.046386199083146756)**2)) - 13917.940920275636*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.9074219145138307)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.5830026053427076)**2)) + 84236.79319164483*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.7557460774371298)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.6383110566137417)**2)) - 5640.99826487992*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.5911925356865033)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.005381283583344741)**2)) - 126632.90425059432*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.21031274399063923)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.5495609776179707)**2)) + 208601.0077510439*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.7830308832306688)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.24439265066740792)**2)) - 31992.625017974526*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.28803140536417743)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.8689715765112118)**2)) + 2004.0939995171502*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.993941455353761)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.9026639746834018)**2)) + 112631.0142006306*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.36343808704287006)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.6057949629205177)**2)) + 11001.494287014008*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.3586109792561233)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.3372905989731323)**2)) - 58592.92022290826*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.17758838147122427)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.23089727866378493)**2)) + 101781.91392371058*exp(- (5.396004300925354*((x[1] + 4.929217157135412)/14.643030012320084 - 0.04112204875223326)**2 + 0.21786916191639816*((x[2] - 0.22882456869508516)/14.45053220191546 - 0.5088149538788017)**2)) + 399.8740733017947
As we can see, expressing a Kriging model algebraically is pretty complicated!
Parity plots for the Kriging model may be viewed with the parity_residual_plots method.
Further information about using PySMO's Kriging tool and features can be found in the documentation.
Summary¶
PySMO allows IDAES users to sample design spaces and generate different types of surrogate models. Further information about PySMO's capabilities may be found in the documentation.