In [1]:
from pycalphad.model import TestModel
import pycalphad
import pycalphad.variables as v
from pycalphad.core.halton import halton
import numpy as np
import itertools
rng = np.random.RandomState(1769)
In [2]:
def sample_quasi(comp_count, pdof=10):
"""
Sample 'pdof * (sum(comp_count) - len(comp_count))' points in
composition space for the sublattice configuration specified
by 'comp_count'. Points are sampled quasi-randomly from a Halton sequence.
A Halton sequence is like a uniform random distribution, but the
result will always be the same for a given 'comp_count' and 'pdof'.
Note: For systems with only one component, only one point will be
returned, regardless of 'pdof'. This is because the degrees of freedom
are zero for that case.
Parameters
----------
comp_count : list
Number of components in each sublattice.
pdof : int
Number of points to sample per degree of freedom.
Returns
-------
ndarray of generated points satisfying the mass balance.
Examples
--------
>>> comps = [8,1] # 8 components in sublattice 1; only 1 in sublattice 2
>>> pts = point_sample(comps, pdof=20) # 7 d.o.f, returns a 140x7 ndarray
"""
# Generate Halton sequence with appropriate dimensions and size
pts = halton(sum(comp_count),
pdof * (sum(comp_count) - len(comp_count)), scramble=True)
# Convert low-discrepancy sequence to normalized exponential
# This will be uniformly distributed over the simplices
pts = -np.log(pts)
cur_idx = 0
for ctx in comp_count:
end_idx = cur_idx + ctx
pts[:, cur_idx:end_idx] /= pts[:, cur_idx:end_idx].sum(axis=1)[:, None]
cur_idx = end_idx
if len(pts) == 0:
pts = np.atleast_2d([1] * len(comp_count))
return pts
def sample_pseudo(comp_count, pdof=10):
"""
Sample 'pdof * (sum(comp_count) - len(comp_count))' points in
composition space for the sublattice configuration specified
by 'comp_count'. Points are sampled pseudo-randomly from a symmetric
Dirichlet distribution.
Note: For systems with only one component, only one point will be
returned, regardless of 'pdof'. This is because the degrees of freedom
are zero for that case.
Parameters
----------
comp_count : list
Number of components in each sublattice.
pdof : int
Number of points to sample per degree of freedom.
Returns
-------
ndarray of generated points satisfying the mass balance.
Examples
--------
>>> comps = [8,1] # 8 components in sublattice 1; only 1 in sublattice 2
>>> pts = point_sample(comps, pdof=20) # 7 d.o.f, returns a 140x7 ndarray
"""
# Generate Halton sequence with appropriate dimensions and size
pts = np.empty((pdof * (sum(comp_count) - len(comp_count)), sum(comp_count), ))
cur_idx = 0
for ctx in comp_count:
end_idx = cur_idx + ctx
pts[:, cur_idx:end_idx] = rng.dirichlet(np.ones(ctx, dtype=np.int),
size=(pdof * (sum(comp_count) - len(comp_count))))
cur_idx = end_idx
if len(pts) == 0:
pts = np.atleast_2d([1] * len(comp_count))
return pts
def sample_baduniform(comp_count, pdof=10):
"""
Sample 'pdof * (sum(comp_count) - len(comp_count))' points in
composition space for the sublattice configuration specified
by 'comp_count'. Points are sampled pseudo-randomly from a uniform
distribution and then normalized to unity.
Note: For systems with only one component, only one point will be
returned, regardless of 'pdof'. This is because the degrees of freedom
are zero for that case.
Parameters
----------
comp_count : list
Number of components in each sublattice.
pdof : int
Number of points to sample per degree of freedom.
Returns
-------
ndarray of generated points satisfying the mass balance.
Examples
--------
>>> comps = [8,1] # 8 components in sublattice 1; only 1 in sublattice 2
>>> pts = point_sample(comps, pdof=20) # 7 d.o.f, returns a 140x7 ndarray
"""
# Generate Halton sequence with appropriate dimensions and size
pts = np.empty((pdof * (sum(comp_count) - len(comp_count)), sum(comp_count), ))
cur_idx = 0
for ctx in comp_count:
end_idx = cur_idx + ctx
pts[:, cur_idx:end_idx] = rng.uniform(size=(pdof * (sum(comp_count) - len(comp_count)), ctx))
pts[:, cur_idx:end_idx] /= pts[:, cur_idx:end_idx].sum(axis=1)[:, None]
cur_idx = end_idx
if len(pts) == 0:
pts = np.atleast_2d([1] * len(comp_count))
return pts
def sample_fixed_grid(comp_count, pdof=10):
"""
Sample _roughly_ 'pdof * (sum(comp_count) - len(comp_count))' points in
composition space for the sublattice configuration specified
by 'comp_count'. Points are generated on a uniform grid by rejection sampling.
Note: For systems with only one component, only one point will be
returned, regardless of 'pdof'. This is because the degrees of freedom
are zero for that case.
Parameters
----------
comp_count : list
Number of components in each sublattice.
pdof : int
Number of points to sample per degree of freedom.
Returns
-------
ndarray of generated points satisfying the mass balance.
"""
print(result_simplex.shape)
for ctx in comp_count:
n = ctx # number of components
points_per_edge = 10
x = np.linspace(1e-12, 1-1e-12, points_per_edge)
result = np.array(np.meshgrid(*itertools.repeat(x, n-1)))
result_simplex = result[:, result.sum(axis=0) <= 1]
# Add dependent component
temp = np.empty((result_simplex.shape[0]+1, result_simplex.shape[1]))
temp[:-1, :] = result_simplex
temp[-1, :] = 1 - result_simplex.sum(axis=0)
result_simplex = temp
def test_opt(num_components, num_samples, sampler, solutions):
result_array = np.empty((num_samples))
comps = ['X{}'.format(str(i+1).zfill(len(str(num_components)))) for i in range(num_components)]
print('{} COMPONENTS'.format(num_components))
for idx in range(num_samples):
print('TRIAL {}/{}'.format(idx+1, num_samples))
solution = solutions[idx]
tm = TestModel(None, comps, 'TEST', solution=solution, kmax=2)
dbf = pycalphad.Database()
dbf.add_phase('TEST', {}, [1])
dbf.add_phase_constituents('TEST', [comps])
conds = {v.P:101325, v.T:300}
conds.update({v.X(i):sol for i, sol in list(zip(comps, solution))[:-1]})
grid_opts = {'sampler': sampler, 'pdens': 100*num_components}
eq = pycalphad.equilibrium(dbf, comps, 'TEST', conds, model=tm, verbose=False, grid_opts=grid_opts)
if eq.GM < 1e-2:
result_array[idx] = eq.attrs['solve_iterations']
else:
print(eq.GM)
print(eq.X.values)
print(solution)
result_array[idx] = np.nan
return result_array
In [3]:
%%time
comp_iters_pseudo = []
comp_iters_quasi = []
comp_iters_normunity = []
max_comps = 10
num_trials = 20
for num_comps in range(2, max_comps+1):
solutions = [rng.dirichlet(np.ones(num_comps, dtype=np.int)) for i in range(num_trials)]
from pycalphad.core.lower_convex_hull import lower_convex_hull, _initialize_array
print('PSEUDO')
results_pseudo = test_opt(num_comps, num_trials, sample_pseudo, solutions)
print('QUASI')
results_quasi = test_opt(num_comps, num_trials, sample_quasi, solutions)
print('NORMUNITY')
results_normunity = test_opt(num_comps, num_trials, sample_baduniform, solutions)
comp_iters_pseudo.append(results_pseudo)
comp_iters_quasi.append(results_quasi)
comp_iters_normunity.append(results_normunity)
comp_iters_pseudo = np.array(comp_iters_pseudo, dtype=np.float)
comp_iters_quasi = np.array(comp_iters_quasi, dtype=np.float)
comp_iters_normunity = np.array(comp_iters_normunity, dtype=np.float)
PSEUDO
2 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
2 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
2 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
PSEUDO
3 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
3 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
3 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
PSEUDO
4 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
4 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
4 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
PSEUDO
5 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
5 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
5 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1)>
array([[[[[[ 7238.68588873]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.6361
* X_X2 (X_X2) float64 0.0401
* X_X3 (X_X3) float64 0.02001
* X_X4 (X_X4) float64 0.2639
[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 6.08307823e-01 1.34597549e-01 7.48176512e-06 2.46601915e-01
1.04852314e-02]
[ 5.07932112e-01 3.97586318e-06 9.81424797e-02 2.84523475e-01
1.09397957e-01]
[ 3.44334551e-01 5.76482064e-02 3.19692179e-05 5.52583078e-01
4.54021957e-02]
[ 4.53327847e-01 2.71505425e-04 2.61954220e-02 3.65522236e-01
1.54682989e-01]]]]]]]]
[ 0.63614401 0.04009614 0.02001309 0.26387137 0.0398754 ]
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
PSEUDO
6 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>
array([[[[[[[ 14325.51248918]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.005255
* X_X2 (X_X2) float64 0.01646
* X_X3 (X_X3) float64 0.2695
* X_X4 (X_X4) float64 0.4737
* X_X5 (X_X5) float64 0.0004482
[[[[[[[[[ 2.28994549e-02 2.63686386e-02 1.41746895e-01
5.49870809e-01 9.80551357e-07 2.59113222e-01]
[ 2.52882907e-03 1.33567198e-01 4.90304457e-01
3.24080139e-01 2.07881999e-03 4.74405570e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12]
[ 8.37532226e-03 4.41192924e-02 4.00108276e-01
4.77598536e-01 1.94258016e-02 5.03727722e-02]
[ 2.11603941e-03 2.23942894e-02 3.30205564e-01
2.35238421e-01 3.94165622e-04 4.09651520e-01]]]]]]]]]
[ 5.25474112e-03 1.64646957e-02 2.69451682e-01 4.73661002e-01
4.48199795e-04 2.34719679e-01]
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>
array([[[[[[[ 12555.60915059]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.336
* X_X2 (X_X2) float64 0.1077
* X_X3 (X_X3) float64 0.00534
* X_X4 (X_X4) float64 0.09086
* X_X5 (X_X5) float64 0.4575
[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.66080538e-01 4.09157822e-01 2.62032915e-05
4.10971526e-02 2.79313777e-01 4.32450698e-03]
[ 4.82810025e-01 3.89000668e-02 4.55702050e-02
2.24143598e-02 4.09892565e-01 4.12778463e-04]
[ 5.14111255e-01 5.31863628e-02 3.53283697e-04
2.13229225e-01 2.19082119e-01 3.77545277e-05]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12]
[ 4.36960260e-01 1.69089519e-02 1.78538181e-07
1.55286086e-01 3.76765658e-01 1.40788660e-02]]]]]]]]]
[ 0.3360238 0.10770806 0.0053402 0.09086331 0.45751585 0.00254878]
TRIAL 20/20
QUASI
6 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>
array([[[[[[[ 23315.20975776]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.005255
* X_X2 (X_X2) float64 0.01646
* X_X3 (X_X3) float64 0.2695
* X_X4 (X_X4) float64 0.4737
* X_X5 (X_X5) float64 0.0004482
[[[[[[[[[ 4.04723993e-02 1.17201777e-03 1.02783905e-01
6.81601795e-01 1.56040473e-06 1.73968322e-01]
[ 1.22040661e-02 4.83666910e-02 2.82503668e-01
2.93110740e-01 1.21864049e-04 3.63692972e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12]
[ 1.49666711e-02 7.24832569e-03 2.85433043e-01
2.78381683e-01 6.88639087e-03 4.07083886e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00]]]]]]]]]
[ 5.25474112e-03 1.64646957e-02 2.69451682e-01 4.73661002e-01
4.48199795e-04 2.34719679e-01]
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>
array([[[[[[[ 9395.60280102]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.336
* X_X2 (X_X2) float64 0.1077
* X_X3 (X_X3) float64 0.00534
* X_X4 (X_X4) float64 0.09086
* X_X5 (X_X5) float64 0.4575
[[[[[[[[[ 6.77486496e-01 3.45970907e-02 1.69615539e-02
2.36129985e-01 3.48244002e-02 4.74495647e-07]
[ 5.48283355e-02 4.19898050e-01 5.63290984e-03
1.91002461e-01 3.28637368e-01 8.74960982e-07]
[ 5.36472193e-01 2.81966942e-01 1.73428686e-02
7.60516103e-03 1.56611263e-01 1.57277379e-06]
[ 2.46457641e-01 1.48773099e-02 2.12470196e-02
1.29319472e-01 5.88097310e-01 1.24760401e-06]
[ 2.64258268e-01 3.45666704e-03 3.82536937e-03
4.30175287e-03 7.24156040e-01 1.90222818e-06]
[ 6.01352316e-01 1.31506394e-01 5.58706897e-10
1.35043948e-01 1.23548989e-01 8.54835207e-03]]]]]]]]]
[ 0.3360238 0.10770806 0.0053402 0.09086331 0.45751585 0.00254878]
TRIAL 20/20
NORMUNITY
6 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>
array([[[[[[[ 13818.1806527]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.005255
* X_X2 (X_X2) float64 0.01646
* X_X3 (X_X3) float64 0.2695
* X_X4 (X_X4) float64 0.4737
* X_X5 (X_X5) float64 0.0004482
[[[[[[[[[ 4.61554052e-02 8.27797832e-02 2.27583028e-01
3.70677232e-01 6.75748450e-04 2.72128803e-01]
[ 1.52135241e-03 3.08921218e-01 2.90501438e-01
2.26290076e-01 2.94584845e-07 1.72765621e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12]
[ nan nan nan
nan nan nan]
[ nan nan nan
nan nan nan]]]]]]]]]
[ 5.25474112e-03 1.64646957e-02 2.69451682e-01 4.73661002e-01
4.48199795e-04 2.34719679e-01]
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>
array([[[[[[[ 12488.49977685]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.5454
* X_X2 (X_X2) float64 0.005528
* X_X3 (X_X3) float64 0.135
* X_X4 (X_X4) float64 0.08867
* X_X5 (X_X5) float64 0.0082
[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 3.72980817e-01 4.65234656e-02 1.31283865e-01
9.55051871e-02 2.48863598e-03 3.51218029e-01]
[ 4.56802832e-01 1.49080445e-03 3.94107230e-01
4.37233356e-02 3.29358153e-06 1.03872504e-01]
[ 3.58038337e-01 1.11561921e-05 2.66068523e-02
2.91603860e-01 2.30479354e-02 3.00691859e-01]
[ 3.24653212e-01 2.41257088e-05 2.07002562e-01
9.34844005e-02 4.24244391e-02 3.32411260e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00]]]]]]]]]
[ 0.54544083 0.00552811 0.13496374 0.0886669 0.00820048 0.21719994]
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>
array([[[[[[[ 21452.94040026]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.336
* X_X2 (X_X2) float64 0.1077
* X_X3 (X_X3) float64 0.00534
* X_X4 (X_X4) float64 0.09086
* X_X5 (X_X5) float64 0.4575
[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.25256905e-01 4.23827402e-01 2.56755900e-02
5.16846078e-02 2.73535480e-01 2.00150336e-05]
[ 5.27543048e-01 2.21567822e-02 8.42526024e-02
5.54984906e-02 3.10532561e-01 1.65154555e-05]
[ 1.08764846e-01 2.69558227e-01 2.12049496e-03
2.93867309e-01 3.25688992e-01 1.30132014e-07]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12]
[ 2.57116944e-01 1.70649815e-01 2.50484573e-08
1.67493091e-01 3.85881083e-01 1.88590413e-02]]]]]]]]]
[ 0.3360238 0.10770806 0.0053402 0.09086331 0.45751585 0.00254878]
TRIAL 20/20
PSEUDO
7 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>
array([[[[[[[[ 12143.2016423]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.2864
* X_X2 (X_X2) float64 0.01361
* X_X3 (X_X3) float64 0.3705
* X_X4 (X_X4) float64 0.1148
* X_X5 (X_X5) float64 0.07925
* X_X6 (X_X6) float64 0.1344
[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 4.06795681e-01 8.15709375e-02 3.53371408e-01
1.07040370e-02 2.33356673e-03 1.42483583e-01
2.74078647e-03]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 3.20983587e-01 2.11604841e-02 1.56648979e-01
7.38654361e-02 2.98633773e-01 1.28670610e-01
3.71309076e-05]
[ 3.41347449e-01 4.21348088e-02 2.56336192e-01
5.33628045e-02 2.00197490e-01 1.02216075e-01
4.40518132e-03]
[ 8.19245733e-02 6.94143385e-03 4.06373279e-01
4.20412062e-02 4.10129593e-02 4.21576857e-01
1.29691339e-04]
[ 3.26122702e-01 1.27546656e-02 3.10233766e-01
2.25030638e-01 1.77543621e-02 1.05679822e-01
2.42404372e-03]]]]]]]]]]
[ 0.2863812 0.0136059 0.37047733 0.114804 0.07925465 0.13435309
0.00112385]
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
7 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>
array([[[[[[[[ 20768.12436207]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.2956
* X_X2 (X_X2) float64 0.000142
* X_X3 (X_X3) float64 0.06756
* X_X4 (X_X4) float64 0.03174
* X_X5 (X_X5) float64 0.1409
* X_X6 (X_X6) float64 0.2146
[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.82319327e-01 1.65483781e-03 1.06801155e-01
1.03120283e-02 1.06629194e-01 2.97884072e-01
2.94399385e-01]
[ 4.01622708e-01 4.33715605e-06 2.96283160e-01
1.31752327e-02 7.12994466e-02 2.01019883e-01
1.65952320e-02]
[ 1.41232782e-01 5.85016252e-05 4.75439332e-02
7.55396312e-02 2.08764150e-01 2.62837737e-01
2.64023266e-01]
[ 1.30011937e-01 2.14423163e-05 4.96361670e-02
3.78632117e-02 4.34443001e-01 8.79283090e-02
2.60095933e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00]]]]]]]]]]
[ 2.95591933e-01 1.41958677e-04 6.75577666e-02 3.17416602e-02
1.40863083e-01 2.14621331e-01 2.49482267e-01]
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
7 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>
array([[[[[[[[ 30129.37048082]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.2956
* X_X2 (X_X2) float64 0.000142
* X_X3 (X_X3) float64 0.06756
* X_X4 (X_X4) float64 0.03174
* X_X5 (X_X5) float64 0.1409
* X_X6 (X_X6) float64 0.2146
[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 2.53562387e-01 1.34072367e-03 4.46139593e-04
3.50878202e-02 3.31969453e-01 2.32949790e-01
1.44643687e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 2.37993183e-01 8.00158467e-05 7.41135674e-02
7.85937575e-02 1.33041239e-01 2.50885301e-01
2.25292937e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00]]]]]]]]]]
[ 2.95591933e-01 1.41958677e-04 6.75577666e-02 3.17416602e-02
1.40863083e-01 2.14621331e-01 2.49482267e-01]
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>
array([[[[[[[[ 15543.65030879]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.0001097
* X_X2 (X_X2) float64 0.2581
* X_X3 (X_X3) float64 0.2186
* X_X4 (X_X4) float64 0.04642
* X_X5 (X_X5) float64 0.2661
* X_X6 (X_X6) float64 0.1149
[[[[[[[[[[ 3.11675739e-03 3.05939450e-01 1.60983103e-01
7.21524643e-02 2.77349720e-01 1.90515173e-02
1.61406988e-01]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 7.78576033e-07 2.58993420e-01 3.08435205e-01
1.40567349e-01 1.35289094e-01 5.48925683e-02
1.01821586e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12]
[ 6.15393783e-07 2.57676965e-01 1.97439094e-01
5.88892870e-02 2.12131437e-01 2.12841783e-01
6.10208189e-02]
[ 9.59530652e-06 2.02220894e-01 1.20179140e-01
2.85381378e-02 2.53040728e-01 1.78094773e-01
2.17916732e-01]]]]]]]]]]
[ 1.09720560e-04 2.58123592e-01 2.18574574e-01 4.64156446e-02
2.66086284e-01 1.14907390e-01 9.57827936e-02]
TRIAL 20/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>
array([[[[[[[[ 9276.12132551]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.5347
* X_X2 (X_X2) float64 0.2542
* X_X3 (X_X3) float64 0.08581
* X_X4 (X_X4) float64 0.05228
* X_X5 (X_X5) float64 0.01141
* X_X6 (X_X6) float64 0.03852
[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 3.68784363e-01 2.99669911e-01 2.09561728e-01
6.16606033e-02 5.35216108e-07 3.26033780e-02
2.77194820e-02]
[ 4.01529562e-01 3.44737716e-01 3.45576770e-02
1.18296457e-01 3.42961326e-05 9.38142043e-02
7.03008768e-03]
[ 4.64047990e-01 2.68104176e-01 2.48725422e-02
1.01905726e-01 1.67509930e-02 2.40809625e-02
1.00237611e-01]
[ 4.17028951e-01 2.89627588e-01 1.46138096e-01
1.06610025e-04 8.62343173e-02 4.03585411e-02
2.05058962e-02]
[ 3.15252340e-01 3.25542784e-01 1.60642141e-01
3.87306275e-02 5.72606478e-08 9.67421065e-02
6.30899434e-02]]]]]]]]]]
[ 0.53468409 0.25423942 0.08580726 0.05227689 0.01140856 0.03851943
0.02306436]
PSEUDO
8 COMPONENTS
TRIAL 1/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>
array([[[[[[[[[ 24925.26079685]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.132
* X_X2 (X_X2) float64 0.3108
* X_X3 (X_X3) float64 0.001656
* X_X4 (X_X4) float64 0.007117
* X_X5 (X_X5) float64 0.05303
* X_X6 (X_X6) float64 0.1749
* X_X7 (X_X7) float64 0.01697
[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 2.45350891e-01 2.53024277e-03 2.19646487e-02
6.64647951e-06 1.84315590e-01 2.60247963e-01
5.57802357e-04 2.85026216e-01]
[ 2.08126209e-01 1.61910448e-01 1.31610565e-07
1.52374550e-02 2.59663396e-02 1.70396324e-01
2.38899459e-02 3.94473147e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12]
[ 4.56930588e-02 2.99722910e-01 5.12090565e-04
1.41072319e-05 2.07422880e-01 1.57307182e-01
9.25090928e-02 1.96818679e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00]]]]]]]]]]]
[ 0.13198901 0.31075186 0.00165636 0.00711656 0.05302571 0.17494918
0.01697014 0.30354119]
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>
array([[[[[[[[[ 29020.59743032]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.01202
* X_X2 (X_X2) float64 0.002629
* X_X3 (X_X3) float64 0.06178
* X_X4 (X_X4) float64 0.1926
* X_X5 (X_X5) float64 0.001547
* X_X6 (X_X6) float64 0.2561
* X_X7 (X_X7) float64 0.2868
[[[[[[[[[[[ 3.06804350e-02 9.66445553e-06 2.54806365e-01
2.71455875e-02 4.90583095e-03 1.84318101e-02
2.22672374e-01 4.41347933e-01]
[ 2.63840442e-02 8.12518019e-03 1.12363698e-01
2.21027081e-01 1.56394501e-06 1.77915799e-01
2.34147510e-01 2.20035124e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 4.66985013e-03 1.90365211e-05 1.93950589e-02
6.56178124e-01 5.26194990e-03 1.46536685e-01
4.30737026e-02 1.24865593e-01]
[ 8.44995295e-03 6.80495372e-05 1.03343849e-02
2.44006807e-01 1.95509978e-02 2.26380340e-01
3.28080674e-01 1.63128793e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00]]]]]]]]]]]
[ 0.01202487 0.00262912 0.06178035 0.19263393 0.00154735 0.25611985
0.28678588 0.18647865]
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
8 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>
array([[[[[[[[[ 26711.88091663]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.01202
* X_X2 (X_X2) float64 0.002629
* X_X3 (X_X3) float64 0.06178
* X_X4 (X_X4) float64 0.1926
* X_X5 (X_X5) float64 0.001547
* X_X6 (X_X6) float64 0.2561
* X_X7 (X_X7) float64 0.2868
[[[[[[[[[[[ 2.68793576e-02 2.61567032e-04 9.13722020e-02
1.01270066e-01 1.87975902e-03 2.42202336e-01
1.08268326e-01 4.27866386e-01]
[ 1.30322758e-03 2.04288838e-02 4.76188312e-02
3.78046249e-02 1.40714227e-07 1.86153776e-01
6.09296668e-01 9.73938478e-02]
[ 1.17118219e-02 3.86411293e-04 4.18077159e-01
4.95621115e-02 1.32493090e-02 1.76350675e-01
2.82255743e-01 4.84067699e-02]
[ 2.41304820e-02 6.55383823e-04 6.14534678e-02
4.07421770e-01 8.11797282e-07 1.04007343e-01
1.04706854e-01 2.97623887e-01]
[ 1.10913058e-05 6.86044693e-03 2.31760827e-01
3.73398133e-02 2.20013913e-02 2.25373639e-01
3.11472468e-01 1.65180323e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00]]]]]]]]]]]
[ 0.01202487 0.00262912 0.06178035 0.19263393 0.00154735 0.25611985
0.28678588 0.18647865]
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
8 COMPONENTS
TRIAL 1/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>
array([[[[[[[[[ 36582.43384545]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.132
* X_X2 (X_X2) float64 0.3108
* X_X3 (X_X3) float64 0.001656
* X_X4 (X_X4) float64 0.007117
* X_X5 (X_X5) float64 0.05303
* X_X6 (X_X6) float64 0.1749
* X_X7 (X_X7) float64 0.01697
[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 3.10938021e-01 1.53445960e-01 2.44556899e-02
1.01461956e-02 1.80289338e-03 2.62122842e-01
3.28575053e-02 2.04230892e-01]
[ 2.18353347e-01 2.13403942e-01 4.95389178e-04
7.20272596e-02 1.92592930e-02 2.54290171e-01
6.83001625e-02 1.53870435e-01]
[ 1.14242449e-01 1.99545072e-01 3.14120203e-04
5.41035957e-03 2.26866191e-01 2.47855708e-01
2.27734152e-02 1.82992684e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12]
[ 2.21314783e-01 1.93722590e-01 1.27473530e-03
3.61336483e-03 1.35766167e-02 2.42844156e-01
1.86729123e-01 1.36924631e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00]]]]]]]]]]]
[ 0.13198901 0.31075186 0.00165636 0.00711656 0.05302571 0.17494918
0.01697014 0.30354119]
TRIAL 2/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>
array([[[[[[[[[ 11926.98627586]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.00885
* X_X2 (X_X2) float64 0.04943
* X_X3 (X_X3) float64 0.1276
* X_X4 (X_X4) float64 0.01819
* X_X5 (X_X5) float64 0.22
* X_X6 (X_X6) float64 0.2321
* X_X7 (X_X7) float64 0.01136
[[[[[[[[[[[ 1.16752330e-02 1.04250710e-02 1.20763895e-01
9.25891505e-02 1.44203206e-01 3.34224302e-01
4.55614035e-05 2.86073581e-01]
[ 7.11727074e-03 1.07663566e-01 3.76451511e-01
3.95537060e-03 6.86824701e-02 1.95410549e-01
1.69328577e-07 2.40719093e-01]
[ 2.07324560e-02 1.04925040e-02 1.80631466e-01
1.52932000e-03 8.64077363e-02 3.87146177e-01
5.24730335e-02 2.60587307e-01]
[ 1.11282438e-01 1.33940705e-02 1.84401208e-01
1.18986031e-02 2.23817579e-01 2.73018199e-01
1.30848454e-05 1.82174818e-01]
[ 7.38646546e-02 1.24286333e-01 1.09249905e-01
1.45298742e-07 2.30374058e-01 2.40698386e-01
2.80464999e-04 2.21246052e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12]
[ 1.37591238e-06 5.44590777e-02 8.71072338e-02
2.55573236e-02 3.63978451e-01 9.67690412e-02
1.06368041e-02 3.61490693e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00]]]]]]]]]]]
[ 0.00885012 0.04943079 0.12760323 0.01819249 0.22004347 0.23211777
0.01136232 0.3323998 ]
TRIAL 3/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>
array([[[[[[[[[ 24547.35806647]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.04442
* X_X2 (X_X2) float64 0.3364
* X_X3 (X_X3) float64 0.004386
* X_X4 (X_X4) float64 0.03402
* X_X5 (X_X5) float64 0.491
* X_X6 (X_X6) float64 0.0195
* X_X7 (X_X7) float64 0.008923
[[[[[[[[[[[ 2.02814262e-01 2.67469671e-01 7.72620840e-03
2.85150903e-02 2.10278425e-01 4.87925142e-03
7.52373872e-02 2.03079704e-01]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 1.76930236e-02 2.95818610e-01 2.28377433e-01
1.77384750e-05 2.94607456e-01 4.51802180e-02
6.04574040e-03 1.12259781e-01]
[ 1.03108584e-01 3.50225910e-01 2.84071408e-03
1.21792349e-01 3.16964852e-01 4.45465079e-02
1.06145239e-05 6.05104687e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 7.58496988e-02 1.66031071e-01 1.76688893e-03
2.49963341e-04 4.67322178e-01 6.49287863e-02
6.45353627e-02 1.59316051e-01]
[ 9.81651754e-02 2.40820809e-01 3.90078223e-02
3.19462279e-02 3.24115032e-01 8.38429704e-05
1.73096624e-01 9.27644660e-02]
[ 1.02327975e-01 1.46922057e-01 1.57648926e-03
7.53443205e-04 3.16702839e-01 1.08870542e-05
9.33689101e-03 4.22369419e-01]]]]]]]]]]]
[ 0.04442347 0.33642228 0.00438618 0.03401742 0.49100443 0.01949908
0.00892316 0.06132397]
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>
array([[[[[[[[[ 44246.02970379]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.01202
* X_X2 (X_X2) float64 0.002629
* X_X3 (X_X3) float64 0.06178
* X_X4 (X_X4) float64 0.1926
* X_X5 (X_X5) float64 0.001547
* X_X6 (X_X6) float64 0.2561
* X_X7 (X_X7) float64 0.2868
[[[[[[[[[[[ 5.56878341e-02 3.54251350e-03 1.46870307e-02
2.94141587e-01 4.03335091e-04 9.17200204e-02
2.84088988e-01 2.55728691e-01]
[ 6.89594046e-03 3.17844302e-02 3.41487518e-01
2.61129143e-01 2.82067120e-04 1.68831556e-01
4.73453276e-02 1.42244018e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12]
[ 1.06276494e-04 4.45265368e-03 2.45366581e-01
2.30643063e-01 3.28399307e-02 1.66489557e-01
1.00770205e-01 2.19331733e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00]]]]]]]]]]]
[ 0.01202487 0.00262912 0.06178035 0.19263393 0.00154735 0.25611985
0.28678588 0.18647865]
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
PSEUDO
9 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 22351.08264421]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.00283
* X_X2 (X_X2) float64 0.158
* X_X3 (X_X3) float64 0.006395
* X_X4 (X_X4) float64 0.009869
* X_X5 (X_X5) float64 0.2794
* X_X6 (X_X6) float64 0.2289
* X_X7 (X_X7) float64 0.06501
* X_X8 (X_X8) float64 0.1835
[[[[[[[[[[[[ 1.88635980e-05 4.14667245e-01 1.01248930e-03
1.15237908e-02 1.65610570e-01 2.28220611e-01
1.15068280e-01 3.20651772e-02 3.18129734e-02]
[ 7.47166899e-05 2.02856552e-01 1.89631513e-02
2.88676023e-03 2.96893520e-01 1.58372861e-01
1.66552429e-02 1.23812493e-01 1.79484703e-01]
[ 7.48148142e-03 7.06868878e-02 2.54706622e-03
9.20269885e-02 1.65583645e-01 4.49637066e-01
1.47156951e-02 1.56861469e-01 4.04597013e-02]
[ 7.85949870e-03 8.28920592e-03 6.85691087e-06
4.38913540e-02 2.68776444e-01 2.74402939e-01
9.43394301e-02 1.10110475e-01 1.92323797e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.05863462e-04 2.48613957e-02 1.57773917e-02
2.66265783e-03 3.23541866e-01 2.17030867e-01
1.61565756e-01 1.52883514e-01 1.01470688e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12]
[ 1.27827640e-02 2.63674075e-01 9.71426803e-08
9.84036312e-03 2.00336310e-02 4.57602593e-01
8.81016931e-02 1.07389786e-01 4.05749970e-02]]]]]]]]]]]]
[ 0.00282968 0.15800187 0.00639466 0.00986923 0.27939811 0.22891086
0.06500936 0.18352608 0.06606016]
TRIAL 15/20
TRIAL 16/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 27853.15998024]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.2817
* X_X2 (X_X2) float64 0.0007679
* X_X3 (X_X3) float64 0.3302
* X_X4 (X_X4) float64 0.2655
* X_X5 (X_X5) float64 0.01434
* X_X6 (X_X6) float64 0.009814
* X_X7 (X_X7) float64 0.008701
* X_X8 (X_X8) float64 0.02087
[[[[[[[[[[[[ 4.53072223e-01 1.00000000e-12 3.10862158e-01
3.90503127e-02 2.16292438e-02 2.73578118e-02
2.85883642e-02 1.97978570e-02 9.96420292e-02]
[ 7.61935055e-01 2.34086371e-03 1.19631493e-01
4.97718074e-02 1.32086261e-02 1.10246063e-05
8.66218111e-03 1.52039705e-02 2.92349777e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 4.63165792e-02 1.41015257e-04 3.26997956e-01
4.61739831e-01 1.02728291e-01 1.90861723e-02
1.73606934e-02 3.14107671e-06 2.56263207e-02]
[ 2.49153940e-01 3.16829138e-02 3.66476705e-01
1.89254098e-01 3.44386537e-02 4.00141249e-04
4.45125538e-02 2.84115954e-02 5.56693994e-02]
[ 2.34915758e-02 6.82461048e-08 5.24923059e-01
2.02424382e-01 7.41610300e-03 1.02710726e-02
2.19432630e-01 5.58269889e-03 6.45841078e-03]
[ 2.64126330e-01 1.95286539e-06 5.56801498e-01
5.82959408e-02 6.76872063e-03 2.47907623e-02
2.02202073e-03 6.10893330e-02 2.61034417e-02]
[ 1.99322244e-01 3.28286823e-06 4.46886686e-01
8.44958981e-02 3.33925162e-02 8.00680450e-03
1.01966388e-02 2.11664072e-02 1.96529522e-01]]]]]]]]]]]]
[ 0.28171799 0.00076793 0.33020705 0.26546957 0.01433725 0.00981422
0.00870094 0.02086778 0.06811727]
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
9 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 27130.52965317]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.2817
* X_X2 (X_X2) float64 0.0007679
* X_X3 (X_X3) float64 0.3302
* X_X4 (X_X4) float64 0.2655
* X_X5 (X_X5) float64 0.01434
* X_X6 (X_X6) float64 0.009814
* X_X7 (X_X7) float64 0.008701
* X_X8 (X_X8) float64 0.02087
[[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 5.44747482e-01 5.11107711e-02 7.75277005e-02
1.59300939e-01 2.54185534e-02 1.49869532e-03
5.35822023e-03 5.19889092e-02 8.30487293e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.33487631e-01 1.98865335e-06 4.44565558e-02
3.46186651e-01 6.23758072e-02 6.73565508e-03
1.62227144e-03 1.01811750e-02 2.94952264e-01]
[ 1.00010211e-01 1.27030860e-03 2.47213449e-01
4.76712724e-01 3.99958761e-02 6.40957662e-05
5.92502536e-02 9.50530523e-03 6.59777766e-02]
[ 1.06305319e-01 4.78529776e-03 8.71312188e-02
4.65190409e-01 1.44192214e-02 1.20718239e-04
1.13341776e-02 3.14251950e-02 2.79288443e-01]
[ 3.05843814e-01 1.32689315e-05 4.00680731e-01
1.03029987e-01 7.90990919e-03 7.74854915e-03
8.61569512e-03 1.00345461e-01 6.58125860e-02]
[ 5.80904825e-01 9.43147913e-07 1.09710141e-01
2.20955384e-01 3.39797375e-03 2.36307198e-02
3.45458874e-03 3.49580559e-02 2.29873677e-02]]]]]]]]]]]]
[ 0.28171799 0.00076793 0.33020705 0.26546957 0.01433725 0.00981422
0.00870094 0.02086778 0.06811727]
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
9 COMPONENTS
TRIAL 1/20
TRIAL 2/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 18845.35210356]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.02335
* X_X2 (X_X2) float64 0.06201
* X_X3 (X_X3) float64 0.0399
* X_X4 (X_X4) float64 0.03426
* X_X5 (X_X5) float64 0.3044
* X_X6 (X_X6) float64 0.0172
* X_X7 (X_X7) float64 0.38
* X_X8 (X_X8) float64 0.1372
[[[[[[[[[[[[ 1.38470880e-01 3.75910811e-02 8.29710176e-02
7.85944496e-04 3.05692290e-01 1.05376609e-01
2.83576918e-01 4.54897921e-02 4.54674600e-05]
[ 6.38642361e-03 7.49551392e-02 3.04119255e-01
3.69701112e-07 2.05764674e-01 1.56615833e-02
1.99463695e-01 1.82134667e-01 1.15141929e-02]
[ 5.75022704e-02 4.62050137e-03 9.79584711e-02
5.08338298e-02 2.46979785e-01 2.49614713e-02
2.54960957e-01 2.62182413e-01 3.01692083e-07]
[ 2.31628775e-03 3.64822978e-02 4.94493563e-03
2.13758097e-01 3.37032976e-01 4.51031308e-02
3.13586045e-01 3.99255322e-02 6.85069769e-03]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.45838153e-03 1.07952168e-02 4.25684339e-02
6.33233694e-02 2.10801668e-01 2.46123724e-01
2.79600566e-01 1.44328527e-01 1.13394214e-07]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12]
[ 5.51169099e-09 2.24522430e-01 2.38209963e-03
3.09627411e-02 3.36091529e-01 1.24883447e-02
2.53679182e-01 1.38679511e-01 1.19415693e-03]
[ 7.32878038e-06 8.58043598e-02 3.22480630e-02
6.46635289e-03 3.12293445e-01 1.29966531e-02
3.05240926e-01 5.56875890e-02 1.89255283e-01]]]]]]]]]]]]
[ 0.02334898 0.06201314 0.03989883 0.03426188 0.3043589 0.01719865
0.37995435 0.13717176 0.00179351]
TRIAL 3/20
TRIAL 4/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 17229.16218313]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.005729
* X_X2 (X_X2) float64 0.2758
* X_X3 (X_X3) float64 0.04265
* X_X4 (X_X4) float64 0.01342
* X_X5 (X_X5) float64 0.04879
* X_X6 (X_X6) float64 0.06658
* X_X7 (X_X7) float64 0.3424
* X_X8 (X_X8) float64 0.02674
[[[[[[[[[[[[ 8.69841068e-02 1.71778179e-01 1.33808650e-06
4.20061211e-03 1.75310765e-01 5.58322147e-02
2.87204061e-01 2.20819570e-03 2.16480527e-01]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 4.41070407e-03 2.59354668e-01 2.45190260e-08
6.93565693e-05 3.53542177e-02 1.80890893e-01
2.24746879e-01 1.21827671e-01 1.73345586e-01]
[ 1.89479075e-02 2.05710293e-01 2.99295827e-05
6.36906645e-02 1.78370853e-01 1.96414509e-02
1.96446024e-01 6.09502329e-02 2.56212643e-01]
[ 1.65397598e-05 2.46260170e-01 6.15121840e-02
4.40088784e-03 1.50386842e-01 1.39016888e-01
2.38657416e-01 1.59323769e-02 1.43816696e-01]
[ 9.55388533e-08 1.41062047e-01 9.87435994e-02
1.04178342e-02 2.51380054e-02 5.16631863e-02
3.85285871e-01 8.26732641e-03 2.79422035e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12]
[ 4.80530670e-03 2.01155838e-01 2.15616077e-04
1.38203502e-01 3.89232302e-02 1.57996886e-01
2.22134688e-01 6.47996457e-02 1.71765287e-01]
[ 8.43277284e-06 9.12344747e-02 3.77618278e-02
9.25555555e-03 1.53410813e-01 1.18325447e-01
2.52802517e-01 1.82089349e-02 3.18991997e-01]]]]]]]]]]]]
[ 0.00572916 0.27580921 0.04265009 0.01342425 0.04878618 0.06658208
0.34240152 0.02674375 0.17787376]
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 21237.55386294]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.2367
* X_X2 (X_X2) float64 0.02297
* X_X3 (X_X3) float64 0.006571
* X_X4 (X_X4) float64 0.3339
* X_X5 (X_X5) float64 0.08391
* X_X6 (X_X6) float64 0.00423
* X_X7 (X_X7) float64 0.05765
* X_X8 (X_X8) float64 0.1325
[[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.61280345e-01 2.50884548e-03 2.50971286e-05
1.22001364e-01 1.61454770e-02 5.21866265e-03
1.36395347e-01 2.41875877e-01 2.14548984e-01]
[ 1.99043674e-01 8.46513687e-02 1.78488420e-02
2.10263033e-01 7.67601052e-02 1.02006476e-06
3.45780183e-02 1.75158191e-01 2.01695748e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.10435544e-01 5.41510293e-02 3.42579749e-03
3.41523389e-03 2.97605171e-01 3.46558805e-05
8.09757987e-03 1.66758243e-01 3.56076746e-01]
[ 2.00706994e-01 1.16943625e-02 1.79307717e-07
1.95186127e-01 1.92135919e-01 2.08540779e-02
8.84845103e-02 1.52427984e-01 1.38509845e-01]
[ 2.58389865e-01 5.83934714e-02 7.50950020e-08
1.98174023e-01 1.17143369e-01 8.16756596e-03
1.34095877e-01 1.19114427e-01 1.06521327e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12]
[ 3.18023778e-01 2.18592598e-02 3.02135104e-02
1.97055195e-01 1.46568197e-01 9.60574129e-07
5.23903405e-02 1.07938613e-01 1.25950146e-01]]]]]]]]]]]]
[ 0.23665472 0.0229676 0.00657122 0.33394364 0.08390696 0.00422985
0.05765071 0.13248851 0.12158678]
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 22714.1544021]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.04242
* X_X2 (X_X2) float64 0.007688
* X_X3 (X_X3) float64 0.002395
* X_X4 (X_X4) float64 0.03897
* X_X5 (X_X5) float64 0.05962
* X_X6 (X_X6) float64 0.2164
* X_X7 (X_X7) float64 0.3567
* X_X8 (X_X8) float64 0.08278
[[[[[[[[[[[[ 1.55860083e-01 2.24433685e-02 1.93698328e-04
5.17020701e-03 5.56081772e-02 2.14183690e-01
1.65223056e-01 8.52643689e-02 2.96053351e-01]
[ 3.78208281e-02 1.61055479e-02 2.97884298e-04
8.11677031e-03 1.27880673e-01 2.98259518e-01
2.02556959e-01 1.37018104e-01 1.71943716e-01]
[ 2.37544097e-02 2.32386898e-06 4.02441712e-02
7.89375554e-02 3.36242735e-02 3.55310292e-01
3.76335808e-01 5.51274875e-02 3.66636784e-02]
[ 5.84899745e-03 7.85509743e-04 4.43242957e-03
2.17887089e-01 6.78627470e-02 2.52146182e-01
2.59422149e-01 5.64763332e-03 1.85967263e-01]
[ 1.10867089e-02 7.09962104e-07 4.77729651e-03
9.38480165e-02 1.05196278e-01 3.71429744e-01
1.35005463e-01 2.01465267e-01 7.71905164e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12]
[ 5.16645352e-02 4.67254542e-02 6.44319867e-03
9.19397437e-02 1.26993232e-05 1.16450644e-01
2.39934723e-01 2.15500589e-01 2.31328413e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00]]]]]]]]]]]]
[ 0.04242133 0.00768829 0.00239485 0.03896605 0.05961641 0.2164378
0.35667422 0.08278424 0.19301679]
TRIAL 14/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 30379.06864119]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.00283
* X_X2 (X_X2) float64 0.158
* X_X3 (X_X3) float64 0.006395
* X_X4 (X_X4) float64 0.009869
* X_X5 (X_X5) float64 0.2794
* X_X6 (X_X6) float64 0.2289
* X_X7 (X_X7) float64 0.06501
* X_X8 (X_X8) float64 0.1835
[[[[[[[[[[[[ 3.61163653e-02 1.45635676e-01 8.76563458e-04
1.19819003e-03 2.76769221e-01 9.04968428e-02
1.05361720e-01 3.39366111e-01 4.17931040e-03]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.20017521e-07 2.36206741e-01 5.40237248e-02
9.61529998e-03 1.28720063e-01 1.73024791e-01
9.41448647e-02 2.50784268e-01 5.34800277e-02]
[ 4.18321526e-03 1.78722848e-01 3.95401471e-06
7.10895770e-02 1.00289014e-01 2.56274323e-01
7.15545066e-02 2.17225868e-01 1.00656695e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 9.86002829e-07 1.65909542e-01 9.66264355e-04
7.92454894e-04 1.34336649e-01 1.68737104e-01
2.09153510e-01 1.78762304e-01 1.41341186e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12]
[ 2.69187535e-04 2.21965048e-01 1.22727952e-02
1.19947527e-02 1.84964679e-01 1.83535729e-01
5.10774946e-02 1.64619637e-01 1.69300678e-01]]]]]]]]]]]]
[ 0.00282968 0.15800187 0.00639466 0.00986923 0.27939811 0.22891086
0.06500936 0.18352608 0.06606016]
TRIAL 15/20
TRIAL 16/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>
array([[[[[[[[[[ 48780.99231019]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X1 (X_X1) float64 0.2817
* X_X2 (X_X2) float64 0.0007679
* X_X3 (X_X3) float64 0.3302
* X_X4 (X_X4) float64 0.2655
* X_X5 (X_X5) float64 0.01434
* X_X6 (X_X6) float64 0.009814
* X_X7 (X_X7) float64 0.008701
* X_X8 (X_X8) float64 0.02087
[[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 2.42921399e-01 6.83067510e-02 2.41270201e-01
1.50118177e-01 1.15174153e-02 3.33431800e-02
1.30536663e-02 1.77151406e-02 2.21754069e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12]
[ 3.51726388e-01 1.61249531e-05 6.10922172e-02
8.60933635e-02 9.07572905e-02 3.31300194e-03
4.17503580e-02 1.40716148e-02 3.51179641e-01]
[ 1.28061724e-01 3.44103099e-04 1.43617125e-01
2.67312792e-01 4.79024662e-02 5.88413942e-02
9.94614714e-04 1.20403769e-01 2.32522011e-01]
[ 2.33027583e-01 5.22866946e-05 1.19297223e-01
2.06664912e-01 9.61852251e-03 9.91220718e-02
7.12543029e-02 5.29279094e-02 2.08035189e-01]
[ 2.34747950e-01 6.81532920e-04 2.52842441e-01
1.13472497e-01 8.64418639e-03 4.62828716e-02
1.19118499e-01 1.57041647e-01 6.71683747e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00]]]]]]]]]]]]
[ 0.28171799 0.00076793 0.33020705 0.26546957 0.01433725 0.00981422
0.00870094 0.02086778 0.06811727]
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
PSEUDO
10 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
QUASI
10 COMPONENTS
TRIAL 1/20
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
TRIAL 8/20
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
NORMUNITY
10 COMPONENTS
TRIAL 1/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>
array([[[[[[[[[[[ 32992.27139106]]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X01 (X_X01) float64 0.06489
* X_X02 (X_X02) float64 0.02209
* X_X03 (X_X03) float64 0.08005
* X_X04 (X_X04) float64 0.02539
* X_X05 (X_X05) float64 0.002807
* X_X06 (X_X06) float64 0.2206
* X_X07 (X_X07) float64 0.3263
* X_X08 (X_X08) float64 0.1224
* X_X09 (X_X09) float64 0.005118
[[[[[[[[[[[[[ 9.28319273e-02 1.12761786e-01 2.15646733e-01
1.32299314e-02 7.36246031e-03 7.59389544e-02
1.91644900e-01 2.33460356e-01 1.86335353e-05
5.71043187e-02]
[ 8.25166714e-02 2.76327247e-03 1.58166329e-01
2.21489913e-02 1.51720532e-04 2.36911395e-01
2.09283602e-01 2.25419923e-01 2.54921061e-02
3.71459893e-02]
[ 1.94918877e-01 2.22781758e-03 2.38002395e-01
4.95664475e-03 7.31740448e-04 8.57676219e-03
1.84525940e-01 1.41064936e-01 2.58138390e-03
2.22413503e-01]
[ 6.59858807e-02 1.40062627e-01 1.24390384e-01
9.23242504e-02 3.21058098e-05 8.41272960e-02
1.84474408e-01 1.58008662e-01 9.28513967e-04
1.49665873e-01]
[ 2.48505716e-02 7.28459165e-02 2.08625197e-01
8.84015698e-03 3.67285303e-02 8.49977447e-02
2.54332356e-01 1.73203619e-01 8.09777923e-06
1.35567810e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 8.94638922e-02 2.56407669e-02 4.34686382e-02
7.63292427e-02 9.67864954e-04 1.03662969e-01
2.43247000e-01 2.86652030e-01 1.34947787e-02
1.17072817e-01]
[ 2.14050967e-01 9.25548789e-04 1.52001262e-02
6.15170509e-02 2.24791634e-05 2.21769992e-01
2.09324243e-01 6.76362252e-02 1.12105510e-02
1.98342816e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00]]]]]]]]]]]]]
[ 0.06488928 0.02209467 0.08005051 0.02538987 0.00280749 0.22061855
0.32629228 0.12243823 0.00511819 0.13030094]
TRIAL 2/20
TRIAL 3/20
TRIAL 4/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>
array([[[[[[[[[[[ 35906.97331219]]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X01 (X_X01) float64 0.2537
* X_X02 (X_X02) float64 0.2332
* X_X03 (X_X03) float64 0.1363
* X_X04 (X_X04) float64 0.04557
* X_X05 (X_X05) float64 0.01003
* X_X06 (X_X06) float64 0.1446
* X_X07 (X_X07) float64 0.009234
* X_X08 (X_X08) float64 0.1441
* X_X09 (X_X09) float64 0.01678
[[[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.66598660e-01 1.48727186e-01 1.98077174e-01
1.94907352e-01 2.08502521e-02 2.43630074e-01
3.22203084e-07 2.34995689e-02 9.68890280e-04
2.74052044e-03]
[ 2.47661566e-01 1.00401392e-01 6.97649547e-02
7.15243160e-02 2.28845287e-06 2.48195992e-01
3.80687841e-02 2.11423076e-01 1.15366080e-02
1.42102276e-03]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.90025154e-01 7.34768819e-02 2.01056539e-01
9.85015056e-03 4.59820764e-02 1.79292945e-01
1.17500302e-06 2.35286548e-01 1.32448158e-02
5.17837139e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e+00 1.00000000e-12
1.00000000e-12]
[ 2.05680919e-01 2.05086444e-01 1.44232066e-02
1.06670381e-01 2.17649849e-02 2.19287502e-01
1.34422111e-03 1.47097996e-01 7.86266993e-02
1.76457982e-05]
[ 1.20587086e-01 7.79626369e-02 2.56469127e-01
6.42765056e-02 5.67070499e-02 2.31781352e-01
1.75491315e-03 1.38989660e-01 5.14469198e-02
2.47493634e-05]]]]]]]]]]]]]
[ 0.253678 0.23320191 0.13628781 0.04556506 0.0100317 0.14460173
0.00923435 0.14410382 0.01677506 0.00652055]
TRIAL 5/20
TRIAL 6/20
TRIAL 7/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>
array([[[[[[[[[[[ 44334.2276809]]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X01 (X_X01) float64 0.5539
* X_X02 (X_X02) float64 0.001489
* X_X03 (X_X03) float64 0.006258
* X_X04 (X_X04) float64 0.01409
* X_X05 (X_X05) float64 0.02235
* X_X06 (X_X06) float64 0.1907
* X_X07 (X_X07) float64 0.1355
* X_X08 (X_X08) float64 0.03399
* X_X09 (X_X09) float64 0.004718
[[[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.29381233e-01 6.17042257e-04 1.44450715e-02
1.65719829e-02 9.31488358e-03 3.03914733e-01
2.15848090e-01 1.78035377e-01 4.04707504e-02
9.14008363e-02]
[ 2.83992648e-01 5.50126245e-04 1.78171230e-02
2.29944091e-02 1.35887742e-01 1.87106553e-01
6.23213480e-02 2.15302755e-01 3.18445918e-02
4.21827038e-02]
[ 1.26216164e-01 2.04161309e-02 2.50789842e-03
9.59769815e-02 8.78495894e-02 1.95847912e-01
1.08664660e-01 1.64018417e-01 1.35283133e-06
1.98500894e-01]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 8.19885197e-02 1.99736376e-05 2.28176529e-03
9.07380841e-02 2.12067943e-01 1.55349900e-01
2.31126501e-01 2.65467904e-03 4.05685712e-02
1.83204063e-01]
[ 1.04825050e-01 9.41733983e-03 6.49080302e-02
4.84212721e-02 6.99447259e-03 2.20980034e-01
1.74831367e-01 1.67767689e-01 1.17702312e-06
2.01853568e-01]
[ 1.30633062e-01 7.99857971e-03 8.70160487e-03
1.44581187e-01 1.28871506e-01 1.17551674e-01
1.40014939e-01 1.58358571e-01 1.15513617e-07
1.63288763e-01]]]]]]]]]]]]]
[ 0.55389931 0.00148919 0.00625847 0.01409256 0.02235222 0.19074092
0.13547748 0.03398572 0.00471823 0.03698591]
TRIAL 8/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>
array([[[[[[[[[[[ 25161.85918291]]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X01 (X_X01) float64 0.09572
* X_X02 (X_X02) float64 0.02795
* X_X03 (X_X03) float64 0.1509
* X_X04 (X_X04) float64 0.1646
* X_X05 (X_X05) float64 0.03749
* X_X06 (X_X06) float64 0.3075
* X_X07 (X_X07) float64 0.006549
* X_X08 (X_X08) float64 0.04803
* X_X09 (X_X09) float64 0.1564
[[[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.74857248e-01 2.96173558e-02 1.74856487e-01
1.77947578e-01 4.61011568e-02 1.96937302e-01
4.81331395e-09 1.42543234e-01 4.82827745e-02
8.85685981e-03]
[ 1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.11826974e-01 2.28554216e-02 1.67839497e-01
1.84592309e-01 1.29616678e-01 1.38526357e-01
9.52754526e-03 5.08795955e-02 1.84335343e-01
2.80282346e-07]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.80131744e-02 6.27719523e-02 2.72634297e-01
1.31059257e-01 7.07206622e-02 2.81895541e-01
1.96623695e-02 8.80558396e-03 1.34411590e-01
2.55731459e-05]
[ 1.06760933e-01 4.33499537e-03 1.68161042e-01
2.08222986e-01 5.79010217e-03 1.92143339e-01
1.79480472e-02 2.03455406e-01 9.27956882e-02
3.87461275e-04]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12]
[ 1.89372742e-01 2.07566903e-02 1.92335927e-01
2.14661048e-01 2.96712689e-02 1.80872445e-01
2.60699766e-06 9.22777721e-03 1.36420161e-01
2.66793332e-02]]]]]]]]]]]]]
[ 0.09571528 0.0279492 0.15093519 0.16460131 0.03749186 0.30750794
0.0065485 0.04802951 0.15641141 0.00480979]
TRIAL 9/20
TRIAL 10/20
TRIAL 11/20
TRIAL 12/20
TRIAL 13/20
TRIAL 14/20
TRIAL 15/20
TRIAL 16/20
<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>
array([[[[[[[[[[[ 15143.63158122]]]]]]]]]]])
Coordinates:
* P (P) float64 1.013e+05
* T (T) float64 300.0
* X_X01 (X_X01) float64 0.1467
* X_X02 (X_X02) float64 0.08028
* X_X03 (X_X03) float64 0.05512
* X_X04 (X_X04) float64 0.06812
* X_X05 (X_X05) float64 0.004017
* X_X06 (X_X06) float64 0.07684
* X_X07 (X_X07) float64 0.1361
* X_X08 (X_X08) float64 0.01094
* X_X09 (X_X09) float64 0.3856
[[[[[[[[[[[[[ 1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.61789072e-01 1.45742581e-02 4.92647815e-03
2.08922862e-01 5.79231218e-05 9.30901223e-02
1.01954225e-01 1.03521159e-02 2.37300591e-01
1.67032353e-01]
[ 1.82555591e-01 1.73741629e-01 1.39450608e-01
5.72234304e-02 8.75693696e-06 2.10634692e-02
1.24848504e-01 1.81913744e-03 2.77958720e-01
2.13301534e-02]
[ 2.08639885e-01 1.33930760e-01 7.03464244e-03
1.46408215e-01 1.22733032e-04 9.48127269e-02
1.30640584e-01 1.36160838e-02 2.02274441e-01
6.25199282e-02]
[ 9.19589496e-02 6.65839887e-02 1.51409610e-02
2.40449827e-02 5.20503261e-08 1.33658365e-01
2.71054421e-01 3.58197781e-02 2.10785919e-01
1.50952583e-01]
[ 2.02930234e-01 8.94445084e-03 5.50923689e-02
1.23589063e-01 2.72600754e-05 2.04430013e-01
1.45437736e-01 4.47213236e-02 2.04248070e-01
1.05794809e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e+00 1.00000000e-12 1.00000000e-12
1.00000000e-12]
[ 1.85130462e-01 1.02335599e-01 3.25421169e-02
1.68911543e-02 4.35370755e-04 8.15001958e-02
2.03382202e-01 6.64794763e-02 2.67806634e-01
4.34967891e-02]
[ 1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e-12
1.00000000e-12 1.00000000e-12 1.00000000e+00
1.00000000e-12]
[ 2.03913693e-01 7.87279928e-02 1.82021907e-02
5.39568482e-02 2.70794604e-02 2.07898770e-01
2.12264418e-01 2.81805390e-04 1.91656541e-01
6.01828052e-03]]]]]]]]]]]]]
[ 0.14668039 0.08028181 0.0551166 0.068121 0.00401697 0.07684048
0.13608608 0.01094241 0.38563009 0.03628417]
TRIAL 17/20
TRIAL 18/20
TRIAL 19/20
TRIAL 20/20
CPU times: user 1h 11min 7s, sys: 1h 6min 39s, total: 2h 17min 47s
Wall time: 1h 13min 4s
In [4]:
%matplotlib inline
import matplotlib
import matplotlib.pyplot as plt
matplotlib.style.use('bmh')
ax = plt.figure()
plt.errorbar(np.arange(2, max_comps+1), np.nanmean(comp_iters_pseudo, axis=1),
yerr=np.nanstd(comp_iters_pseudo, axis=1), label='Pseudo')
plt.errorbar(np.arange(2, max_comps+1), np.nanmean(comp_iters_quasi, axis=1),
yerr=np.nanstd(comp_iters_quasi, axis=1), label='Quasi')
plt.errorbar(np.arange(2, max_comps+1), np.nanmean(comp_iters_normunity, axis=1),
yerr=np.nanstd(comp_iters_quasi, axis=1), label='NormUnity')
plt.ylim((0,None))
plt.xticks(list(range(1, max_comps+2)))
plt.xlabel('Number of Components')
plt.ylabel('Number of Iterations')
plt.legend(loc='best')
plt.show()
In [5]:
ind = np.arange(2, num_comps+1) # the x locations for the groups
width = 0.5 # the width of the bars
fig, ax = plt.subplots()
rects1 = ax.bar(2*ind, 100*np.invert(np.isnan(comp_iters_quasi)).sum(1)/comp_iters_quasi.shape[-1],
width, color='r', label='Quasi')
rects2 = ax.bar(2*ind + width, 100*np.invert(np.isnan(comp_iters_pseudo)).sum(1)/comp_iters_pseudo.shape[-1],
width, color='b', label='Pseudo')
rects3 = ax.bar(2*ind + 2*width, 100*np.invert(np.isnan(comp_iters_normunity)).sum(1)/comp_iters_normunity.shape[-1],
width, color='g', label='NormUnity')
# add some text for labels, title and axes ticks
ax.set_xticks(2*ind + 1.5*width)
ax.set_xticklabels(tuple(ind))
ax.set_xlabel('Number of Components')
ax.set_ylabel('% Global Minima Found')
plt.legend(loc=3)
plt.show()
In [ ]:
Content source: richardotis/pycalphad-sandbox
Similar notebooks: