Dimension adaptive sampling tutorial applied to the fusion example
Here, briefly describe the concept behind dimension-adaptive sparse grids, starting from a standard Stochastic Collocation (SC) campaign. Following this, a dimension adaptive EasyVVUQ script using the fusion example is presented. We will assume you are familiar with the basics of EasyVVUQ.
The explanation below is a copy of that in dimension_adaptive_tutorial.ipynb and describes a simpler 2D case. The fusion example has 10 varying quantities.
Standard SC
In a standard EasyVVUQ Campaign, a Stochastic Collocation sampler object might be created via:
.. code:: python
sampler = uq.sampling.SCSampler(vary=vary, polynomial_order=2)
Here the specified polynomial_order, and the number of inputs in vary, determine the number of samples, which increases exponentially fast with an increasing amount of inputs. This is the so-called curse of dimensionality.
Basically, by setting polynomial_order=2 we create a sampling plan through a single tensor product of one-dimensional quadrature nodes with order 3 for every input. It is this tensor product construction that leads to the exponential rise in cost. So if we have 2 inputs x1 and x2, and our one-dimensional quadrature rule of order 2 produces 5 points, we obtain a total of 25 points in the (x1, x2) domain. Likewise, if vary contains 3 inputs, we would need to evaluate the
computational model 125 times, and 10 inputs would require 5**10 = 9765625 model evaluations. For this reason, a standard SC campaign is rarely used beyond 6 or 7 inputs.
Sparse SC
Sparse grids on the other hand, do not create a single tensor product, but build the sampling plan from the ground up by using a linear combination of tensor products involving 1D quadrature rules of different orders.
For two inputs, we might for instance consider using 1D quadrature rules of order [0, 0], [0, 1] and [1, 0], where:
* [0, 0]: a single point in the 2D domain (x1, x2)
* [0, 1]: a line of 3 points with constant x1
* [1, 0]: a line of 3 points with constant x2
In the case of sparse grids it is common to select a nested quadrature rule. This means that the quadrature rule of order p contains all points of the same rule of order p-1. When taking the linear combinations, a nested rule ensures that many points will conincide, which yields efficient sampling plans, especially in higher dimensions. If our nested 1D rule of order 1 and 2 generates the points [0.5] and [0, 0.5, 1] we obtain a sampling plan consisting of
* [0, 0]: [0.5, 0.5]
* [0, 1]: [0.5, 0.0], [0.5, 0.5], [0.5, 1.0]
* [1, 0]: [0.0, 0.5], [0.5, 0.5], [1.0, 0.5],
which gives a total of 5 unique points, compared to a corresponding standard SC campaign with [1, 1], which would generate 9 unique points ([0, 0.5, 1] x [0, 0.5, 1.0]). Note that sparse grids do not circumvent the curse of dimensionality, although they can postpone its effect to higher dimensions.
Dimension-adaptive SC
What we described above is an isotropic sparse grid, since the multi indices [0, 0], [1, 0], [0,1] result in a sampling plan where both inputs end up with the same number of samples. However, in practice model parameters are rarely equally important. The idea behind dimension-adaptive sampling is to build the sampling plan in an iterative fashion, find out which (combination of) parameters are important as we go, and then place more samples along those directions. This results in a
anisotropic sampling plan, where the important inputs get relatively high number of samples. To find out which directions are important we need an appropriate error measure, and we need to split the quadrature order multi indices in an accepted and an admissible set. The accepted set is initialized to [0, 0] in 2D, i.e. we start with just a single code evaluation. Without going into detail, we can think of the admissible set as the candidate refinement directions, from which we must add
a single entry to the accepted set at every iteration.
In our 2D example, at the 1st iteration the candidate set consists of [1, 0] and [0, 1]. That is, we can either refine only x1 or only x2. We must select the multi index which generates the highest error when added to the accepted set. There are a variety of error measures, the two main ones in EasyVVUQ are:
the hierarchical surplus error, and
a variance-based error.
Roughly speaking, the surplus is an interpolation based error, which measures the difference between the code output and the corresponding SC polynomial surrogate, when evaluated at new sample locations. The variance-based error selects the direction in which the variance in the output changes the most. For more information we refer to the references below.
Assume that [1, 0] generated the highest error, and so it is added to the accepted set, now consisting of [0, 0] and [1, 0]. This means that x1 has more points than x2. Also, adding a multi index to the accepted set means that the admissible set changes. In this case, since [1, 0] has been accepted, [2, 0] has become admissible. Note that the new entry [2, 0] also requires new evaluations of the code, and so a new ensemble must be submitted. Again, if we use a
nested rule, the grid of [2, 0] will have a partial overlap with the accepted points, so we only have to evaluate the code at the new points, not all points of [2, 0].
Thus, the admissible set now consists of [0, 1] and [2, 0]. Hence, we now have to option of refining x1 again (to second order), or refining x2 to first order. Assume the latter happens. As both x1 and x2 have been refined to 1st order, [1, 1] has become admissible. If accepted, this multi index results in a simultaneous refinement of both x1 and x2. Note that [1, 1] represents a tensor product, and that therefore it is not the same as [1, 0] and
[0, 1] taken together. We added this example to show that the algoritmn is not limited to one-at-a-time refinement.
To conclude, every time a multi index is accepted, new indices become admissible, and the cycle repeats.
References
Our description of the method here was rather limited, so for more information and applications of this (and similar) methods, see the following references:
Gerstner and M. Griebel. “Dimension–adaptive tensor–product quadrature.” Computing 71.1 (2003): 65-87.
Edeling , H. Arabnejad , R. Sinclair, D. Suleimenova, K. Gopalakrishnan, B. Bosak, D. Groen, I. Mahmood, D. Crommelin, and Peter V Coveney, “The Impact of Uncertainty on Predictions of the CovidSim Epidemiological Code”, Nature Computational Science, 1 (2), 2021.
Loukrezis, U. Römer, and H. De Gersem. “Assessing the performance of Leja and Clenshaw-Curtis collocation for computational electromagnetics with random input data”. International Journal for Uncertainty Quantification , 9(1), 2019.
J.D. Jakeman, M.S. Eldred, G. Geraci, and A. Gorodetsky. “Adaptive multi-index collocation for uncertainty quantification and sensitivity analysis”. Numerical Methods in Engineering , 121(6):1314-1343, 2020.
[1]:
# import packages that we will use
%matplotlib inline
import os
import easyvvuq as uq
import chaospy as cp
import pickle
import time
import numpy as np
import pandas as pd
import matplotlib
if not os.getenv("DISPLAY"): matplotlib.use('Agg')
import matplotlib.pylab as plt
from IPython.display import display
%matplotlib inline
#%matplotlib notebook
[2]:
# we need fipy -- install if not already available
try:
import fipy
except ModuleNotFoundError:
! pip install future
! pip install fipy
import fipy
[3]:
# routine to write out (if needed) the fusion .template file
def write_template(params):
str = ""
first = True
for k in params.keys():
if first:
str += '{"%s": "$%s"' % (k,k) ; first = False
else:
str += ', "%s": "$%s"' % (k,k)
str += '}'
print(str, file=open('fusion.template','w'))
[4]:
# define parameters of the fusion model
def define_params():
return {
"Qe_tot": {"type": "float", "min": 1.0e6, "max": 50.0e6, "default": 2e6},
"H0": {"type": "float", "min": 0.00, "max": 1.0, "default": 0},
"Hw": {"type": "float", "min": 0.01, "max": 100.0, "default": 0.1},
"Te_bc": {"type": "float", "min": 10.0, "max": 1000.0, "default": 100},
"chi": {"type": "float", "min": 0.01, "max": 100.0, "default": 1},
"a0": {"type": "float", "min": 0.2, "max": 10.0, "default": 1},
"R0": {"type": "float", "min": 0.5, "max": 20.0, "default": 3},
"E0": {"type": "float", "min": 1.0, "max": 10.0, "default": 1.5},
"b_pos": {"type": "float", "min": 0.95, "max": 0.99, "default": 0.98},
"b_height": {"type": "float", "min": 3e19, "max": 10e19, "default": 6e19},
"b_sol": {"type": "float", "min": 2e18, "max": 3e19, "default": 2e19},
"b_width": {"type": "float", "min": 0.005, "max": 0.025, "default": 0.01},
"b_slope": {"type": "float", "min": 0.0, "max": 0.05, "default": 0.01},
"nr": {"type": "integer", "min": 10, "max": 1000, "default": 100},
"dt": {"type": "float", "min": 1e-3, "max": 1e3, "default": 100},
"out_file": {"type": "string", "default": "output.csv"}
}
[5]:
# define varying quantities
def define_vary():
vary_all = {
"Qe_tot": cp.Uniform(1.8e6, 2.2e6),
"H0": cp.Uniform(0.0, 0.2),
"Hw": cp.Uniform(0.1, 0.5),
"chi": cp.Uniform(0.8, 1.2),
"Te_bc": cp.Uniform(80.0, 120.0),
"a0": cp.Uniform(0.9, 1.1),
"R0": cp.Uniform(2.7, 3.3),
"E0": cp.Uniform(1.4, 1.6),
"b_pos": cp.Uniform(0.95, 0.99),
"b_height": cp.Uniform(5e19, 7e19),
"b_sol": cp.Uniform(1e19, 3e19),
"b_width": cp.Uniform(0.015, 0.025),
"b_slope": cp.Uniform(0.005, 0.020)
}
vary_2 = {
"Qe_tot": cp.Uniform(1.8e6, 2.2e6),
"Te_bc": cp.Uniform(80.0, 120.0)
}
vary_5 = {
"Qe_tot": cp.Uniform(1.8e6, 2.2e6),
"H0": cp.Uniform(0.0, 0.2),
"Hw": cp.Uniform(0.1, 0.5),
"chi": cp.Uniform(0.8, 1.2),
"Te_bc": cp.Uniform(80.0, 120.0)
}
vary_10 = {
"Qe_tot": cp.Uniform(1.8e6, 2.2e6),
"H0": cp.Uniform(0.0, 0.2),
"Hw": cp.Uniform(0.1, 0.5),
"chi": cp.Uniform(0.8, 1.2),
"Te_bc": cp.Uniform(80.0, 120.0),
"b_pos": cp.Uniform(0.95, 0.99),
"b_height": cp.Uniform(5e19, 7e19),
"b_sol": cp.Uniform(1e19, 3e19),
"b_width": cp.Uniform(0.015, 0.025),
"b_slope": cp.Uniform(0.005, 0.020)
}
return vary_10
[6]:
# define a model to run the fusion code directly from python, expecting a dictionary and returning a dictionary
def run_fusion_model(input):
import json
import fusion
qois = ["te", "ne", "rho", "rho_norm"]
del input['out_file']
return {q: v for q,v in zip(qois, [t.tolist() for t in fusion.solve_Te(**input, plots=False, output=False)])}
[7]:
# routines for plotting the results
def plot_Te(results, title=None):
# plot the calculated Te: mean, with std deviation, 1, 10, 90 and 99%
plt.figure()
rho = results.describe('rho', 'mean')
plt.plot(rho, results.describe('te', 'mean'), 'b-', label='Mean')
plt.plot(rho, results.describe('te', 'mean')-results.describe('te', 'std'), 'b--', label='+1 std deviation')
plt.plot(rho, results.describe('te', 'mean')+results.describe('te', 'std'), 'b--')
plt.fill_between(rho, results.describe('te', 'mean')-results.describe('te', 'std'), results.describe('te', 'mean')+results.describe('te', 'std'), color='b', alpha=0.2)
try:
plt.plot(rho, results.describe('te', '10%'), 'b:', label='10 and 90 percentiles')
plt.plot(rho, results.describe('te', '90%'), 'b:')
plt.fill_between(rho, results.describe('te', '10%'), results.describe('te', '90%'), color='b', alpha=0.1)
plt.fill_between(rho, results.describe('te', '1%'), results.describe('te', '99%'), color='b', alpha=0.05)
except:
print('Problem with some of the percentiles')
plt.legend(loc=0)
plt.xlabel('rho [$m$]')
plt.ylabel('Te [$eV$]')
if not title is None: plt.title(title)
plt.savefig('Te.png')
plt.savefig('Te.pdf')
def plot_ne(results, title=None):
# plot the calculated ne: mean, with std deviation, 1, 10, 90 and 99%
plt.figure()
rho = results.describe('rho', 'mean')
plt.plot(rho, results.describe('ne', 'mean'), 'b-', label='Mean')
plt.plot(rho, results.describe('ne', 'mean')-results.describe('ne', 'std'), 'b--', label='+1 std deviation')
plt.plot(rho, results.describe('ne', 'mean')+results.describe('ne', 'std'), 'b--')
plt.fill_between(rho, results.describe('ne', 'mean')-results.describe('ne', 'std'), results.describe('ne', 'mean')+results.describe('ne', 'std'), color='b', alpha=0.2)
try:
plt.plot(rho, results.describe('ne', '10%'), 'b:', label='10 and 90 percentiles')
plt.plot(rho, results.describe('ne', '90%'), 'b:')
plt.fill_between(rho, results.describe('ne', '10%'), results.describe('ne', '90%'), color='b', alpha=0.1)
plt.fill_between(rho, results.describe('ne', '1%'), results.describe('ne', '99%'), color='b', alpha=0.05)
except:
print('Problem with some of the percentiles')
plt.legend(loc=0)
plt.xlabel('rho [$m$]')
plt.ylabel('ne [$m^{-3}$]')
if not title is None: plt.title(title)
plt.savefig('ne.png')
plt.savefig('ne.pdf')
def plot_sobols_first(results, title=None, field='te'):
# plot the first Sobol results
plt.figure()
rho = results.describe('rho', 'mean')
for k in results.sobols_first()[field].keys(): plt.plot(rho, results.sobols_first()[field][k], label=k)
plt.legend(loc=0)
plt.xlabel('rho [$m$]')
plt.ylabel('sobols_first')
if not title is None: plt.title(field + ': ' + title)
plt.savefig('sobols_first_%s.png' % (field))
plt.savefig('sobols_first_%s.pdf' % (field))
def plot_sobols_second(results, title=None, field='te'):
# plot the second Sobol results
plt.figure()
rho = results.describe('rho', 'mean')
for k1 in results.sobols_second()[field].keys():
for k2 in results.sobols_second()[field][k1].keys():
plt.plot(rho, results.sobols_second()[field][k1][k2], label=k1+'/'+k2)
plt.legend(loc=0, ncol=2)
plt.xlabel('rho [$m$]')
plt.ylabel('sobols_second')
if not title is None: plt.title(field + ': ' + title)
plt.savefig('sobols_second_%s.png' % (field))
plt.savefig('sobols_second_%s.pdf' % (field))
def plot_sobols_total(results, title=None, field='te'):
# plot the total Sobol results
plt.figure()
rho = results.describe('rho', 'mean')
for k in results.sobols_total()[field].keys(): plt.plot(rho, results.sobols_total()[field][k], label=k)
plt.legend(loc=0)
plt.xlabel('rho [$m$]')
plt.ylabel('sobols_total')
if not title is None: plt.title(field + ': ' + title)
plt.savefig('sobols_total_%s.png' % (field))
plt.savefig('sobols_total_%s.pdf' % (field))
def plot_distribution(results, results_df, title=None):
te_dist = results.raw_data['output_distributions']['te']
rho_norm = results.describe('rho_norm', 'mean')
for i in [np.maximum(0, int(i-1))
for i in np.linspace(0,1,5) * rho_norm.shape]:
plt.figure()
pdf_raw_samples = cp.GaussianKDE(results_df.te[i])
pdf_kde_samples = cp.GaussianKDE(te_dist.samples[i])
plt.hist(results_df.te[i], density=True, bins=50, label='histogram of raw samples', alpha=0.25)
if hasattr(te_dist, 'samples'):
plt.hist(te_dist.samples[i], density=True, bins=50, label='histogram of kde samples', alpha=0.25)
plt.plot(np.linspace(pdf_raw_samples.lower, pdf_raw_samples.upper), pdf_raw_samples.pdf(np.linspace(pdf_raw_samples.lower, pdf_raw_samples.upper)), label='PDF (raw samples)')
plt.plot(np.linspace(pdf_kde_samples.lower, pdf_kde_samples.upper), pdf_kde_samples.pdf(np.linspace(pdf_kde_samples.lower, pdf_kde_samples.upper)), label='PDF (kde samples)')
plt.legend(loc=0)
plt.xlabel('Te [$eV$]')
if title is None:
plt.title('Distributions for rho_norm = %0.4f' % (rho_norm[i]))
else:
plt.title('%s\nDistributions for rho_norm = %0.4f' % (title, rho_norm[i]))
plt.savefig('distribution_function_rho_norm=%0.4f.png' % (rho_norm[i]))
plt.savefig('distribution_function_rho_norm=%0.4f.pdf' % (rho_norm[i]))
[8]:
# prepare the test data
if __name__ == '__main__':
test_campaign = uq.Campaign(name='fusion_pce.')
# Add the app (automatically set as current app)
test_campaign.add_app(name="fusion", params=define_params(),
actions=uq.actions.Actions(uq.actions.ExecutePython(run_fusion_model)))
# Associate a sampler with the campaign
test_campaign.set_sampler(uq.sampling.quasirandom.LHCSampler(vary=define_vary(), count=100))
# Perform the actions
test_campaign.execute(nsamples=1000).collate(progress_bar=True)
# Collate the results
test_df = test_campaign.get_collation_result()
100%|███████████████████████████████████████| 1000/1000 [00:19<00:00, 51.40it/s]
[9]:
# calculate the SC surrogates
def test_surrogate():
test_points = np.array(test_df[test_campaign.get_active_sampler().vary.get_keys()])
test_results = test_df['te'].values
test_predictions = np.array([analysis.surrogate('te', tp) for tp in np.array(test_points)])
frms = np.sqrt(((test_predictions - test_results)**2).mean(axis=0)) / test_results.mean(axis=0)
return frms.mean(), frms
[10]:
# routine to run a SC campaign
#def run_sc_case(sc_order=2, local=True, dask=True, batch_size=os.cpu_count(), use_files=True):
if __name__ == '__main__':
sc_order=2; local=True; dask=False; batch_size=os.cpu_count(); use_files=True
if dask:
if local:
print('Running locally')
import multiprocessing.popen_spawn_posix
from dask.distributed import Client, LocalCluster
cluster = LocalCluster(threads_per_worker=1)
client = Client(cluster) # processes=True, threads_per_worker=1)
else:
print('Running using SLURM')
from dask.distributed import Client
from dask_jobqueue import SLURMCluster
cluster = SLURMCluster(
job_extra=['--qos=p.tok.openmp.2h', '--mail-type=end', '--mail-user=dpc@rzg.mpg.de', '-t 2:00:00'],
queue='p.tok.openmp',
cores=8,
memory='8 GB',
processes=8)
cluster.scale(32)
print(cluster)
print(cluster.job_script())
client = Client(cluster)
print(client)
else:
import concurrent.futures
# client = concurrent.futures.ProcessPoolExecutor(max_workers=batch_size)
client = concurrent.futures.ThreadPoolExecutor(max_workers=batch_size)
# client = None
times = np.zeros(7)
time_start = time.time()
time_start_whole = time_start
# Set up a fresh campaign called "fusion_sc."
my_campaign = uq.Campaign(name='fusion_sc_adaptive.')
# Define parameter space
params = define_params()
# Create an encoder and decoder for sc test app
if use_files:
encoder = uq.encoders.GenericEncoder(template_fname='fusion.template',
delimiter='$',
target_filename='fusion_in.json')
decoder = uq.decoders.SimpleCSV(target_filename="output.csv",
output_columns=["te", "ne", "rho", "rho_norm"])
execute = uq.actions.ExecuteLocal('python3 %s/fusion_model.py fusion_in.json' % (os.getcwd()))
actions = uq.actions.Actions(uq.actions.CreateRunDirectory('/tmp'),
uq.actions.Encode(encoder), execute, uq.actions.Decode(decoder))
else:
actions = uq.actions.Actions(uq.actions.ExecutePython(run_fusion_model))
# Add the app (automatically set as current app)
my_campaign.add_app(name="fusion", params=params, actions=actions)
time_end = time.time()
times[1] = time_end-time_start
print('Time for phase 1 = %.3f' % (times[1]))
time_start = time.time()
# Associate a sampler with the campaign
sampler = uq.sampling.SCSampler(vary=define_vary(), polynomial_order=1,
quadrature_rule="C",
sparse=True, growth=True,
midpoint_level1=True,
dimension_adaptive=True)
my_campaign.set_sampler(sampler)
my_campaign.draw_samples()
print('Number of samples = %s' % my_campaign.get_active_sampler().count)
time_end = time.time()
times[2] = time_end-time_start
print('Time for phase 2 = %.3f' % (times[2]))
time_start = time.time()
# Perform the actions
my_campaign.execute(pool=client).collate(progress_bar=True)
if dask:
client.close()
client.shutdown()
time_end = time.time()
times[3] = time_end-time_start
print('Time for phase 3 = %.3f' % (times[3]))
time_start = time.time()
# Collate the results
results_df = my_campaign.get_collation_result()
time_end = time.time()
times[4] = time_end-time_start
print('Time for phase 4 = %.3f' % (times[4]))
time_start = time.time()
# Post-processing analysis
qoi_cols = ["te", "ne", "rho", "rho_norm"]
results = my_campaign.analyse(qoi_cols=qoi_cols)
time_end = time.time()
times[5] = time_end-time_start
print('Time for phase 5 = %.3f' % (times[5]))
time_start = time.time()
# Save the results
pickle.dump(results, open('fusion_results.pickle','bw'))
time_end = time.time()
times[6] = time_end-time_start
print('Time for phase 6 = %.3f' % (times[6]))
times[0] = time_end - time_start_whole
#return results_df, results, times, sc_order, my_campaign.get_active_sampler().count
Time for phase 1 = 0.007
Number of samples = 1
Time for phase 2 = 0.012
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 1/1 [00:03<00:00, 3.24s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Time for phase 3 = 3.248
Time for phase 4 = 0.007
Time for phase 5 = 0.110
Time for phase 6 = 0.001
[11]:
analysis = uq.analysis.SCAnalysis(sampler=sampler, qoi_cols=qoi_cols)
my_campaign.apply_analysis(analysis)
[12]:
S = []
[13]:
frms_mean, frms = test_surrogate()
S.append([frms_mean, frms])
[14]:
analysis.l_norm
[14]:
array([[1, 1, 1, 1, 1, 1, 1, 1, 1, 1]])
A standard SC (or PCE) campaign would be over at this point. Except we have thus far only sampled a single point in the stochastic domain. To show this, we define the following function to plot 2D slices of the accepted points in the 20 dimensional input space. The analysis.l_norm array contains the accepted multi indices.
[15]:
def plot_te(analysis):
te_mean, te_std = analysis.get_moments('te')[0], np.sqrt(analysis.get_moments('te')[1])
rho_mean, rho_std = analysis.get_moments('rho')[0], np.sqrt(analysis.get_moments('rho')[1])
plt.figure()
rho = rho_mean
plt.plot(rho, te_mean, 'b-', label='Mean')
plt.plot(rho, te_mean-te_std, 'b--', label='+1 std deviation')
plt.plot(rho, te_mean+te_std, 'b--')
plt.fill_between(rho, te_mean-te_std, te_mean+te_std, color='b', alpha=0.2)
plt.legend(loc=0)
plt.xlabel('rho [$m$]')
plt.ylabel('Te [$eV$]')
[16]:
plot_te(analysis)
[17]:
def plot_grid_2D():
labels = list(my_campaign.get_active_sampler().vary.get_keys())
values = list(my_campaign.get_active_sampler().vary.get_values())
L = (len(labels)+1)//2
C = int(np.ceil(np.sqrt((10+1)//2)))
R = int(np.ceil(L / C))
fig = plt.figure(figsize=[12,12/C*R])
ax=[]
ic=0
for i in range(L-1):
xd = values[ic*2].upper[0] - values[ic*2].lower[0]
yd = values[ic*2+1].upper[0] - values[ic*2+1].lower[0]
ax.append(fig.add_subplot(R, C, ic+1,
xlim=[values[ic*2].lower[0] - xd/10, values[ic*2].upper[0] + xd/10],
ylim=[values[ic*2+1].lower[0] - yd/10, values[ic*2+1].upper[0] + yd/10],
xlabel=labels[ic*2], ylabel=labels[ic*2+1])
)
ic += 1
xd = values[ic*2].upper[0] - values[ic*2].lower[0]
yd = values[ic*2+1].upper[0] - values[ic*2+1].lower[0]
ax.append(fig.add_subplot(R, C, ic+1,
xlim=[values[-2].lower[0] - xd/10, values[-2].upper[0] + xd/10],
ylim=[values[-1].lower[0] - yd/10, values[-1].upper[0] + yd/10],
xlabel=labels[-2], ylabel=labels[-1])
)
accepted_grid = sampler.generate_grid(analysis.l_norm)
ic=0
for i in range(L-1):
ax[i].plot(accepted_grid[:,ic*2], accepted_grid[:,ic*2+1], 'o', alpha=0.25)
ic += 1
ax[-1].plot(accepted_grid[:,-2], accepted_grid[:,-1], 'o', alpha=0.25)
plt.tight_layout()
[18]:
plot_grid_2D()
To refine the sampling plan, we need to:
Compute the candidate directions of the admissible set. This is done in the
look_aheadsubroutine.Run the ensemble of the new points. This is done exactly the same as before.
Accept the direction with the highest error. This is done in the
adapt_dimensionsubroutine.
[19]:
def refine_sampling_plan(number_of_refinements):
"""
Refine the sampling plan.
Parameters
----------
number_of_refinements (int)
The number of refinement iterations that must be performed.
Returns
-------
None. The new accepted indices are stored in analysis.l_norm and the admissible indices
in sampler.admissible_idx.
"""
for i in range(number_of_refinements):
# compute the admissible indices
sampler.look_ahead(analysis.l_norm)
# run the ensemble
my_campaign.execute().collate(progress_bar=True)
# accept one of the multi indices of the new admissible set
data_frame = my_campaign.get_collation_result()
analysis.adapt_dimension('te', data_frame)
Note that the subroutine above uses the surplus error by default. To select the variance-based error use analysis.adapt_dimension('f', data_frame, method='var') instead.
[20]:
# refine the sampling plan once and then do the analysis to see the results.
refine_sampling_plan(1)
my_campaign.apply_analysis(analysis)
5%|██████████████▏ | 1/20 [00:04<01:24, 4.46s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:08<00:00, 2.39it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
[21]:
frms_mean, frms = test_surrogate()
S.append([frms_mean, frms])
[22]:
analysis.l_norm
[22]:
array([[1, 1, 1, 1, 1, 1, 1, 1, 1, 1],
[1, 1, 2, 1, 1, 1, 1, 1, 1, 1]])
[23]:
plot_te(analysis)
[24]:
# plot the 2D slices again. Note that the most important input (Hw) got refined.
plot_grid_2D()
[25]:
# repeat
refine_sampling_plan(1)
my_campaign.apply_analysis(analysis)
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 2/2 [00:02<00:00, 1.38s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
[26]:
frms_mean, frms = test_surrogate()
S.append([frms_mean, frms])
[27]:
plot_te(analysis)
[28]:
# Now chi got refined.
plot_grid_2D()
[29]:
# again
refine_sampling_plan(1)
my_campaign.apply_analysis(analysis)
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 6/6 [00:02<00:00, 2.00it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
[30]:
frms_mean, frms = test_surrogate()
S.append([frms_mean, frms])
[31]:
plot_te(analysis)
[32]:
# now H0 got refined to first order
plot_grid_2D()
[33]:
# perform 100 refinements
for i in range(100):
print('%s/%s' % (i,100))
refine_sampling_plan(1)
my_campaign.apply_analysis(analysis)
frms_mean, frms = test_surrogate()
S.append([frms_mean, frms])
0/100
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 10/10 [00:03<00:00, 2.55it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
1/100
7%|████████████████████▏ | 1/14 [01:10<15:11, 70.09s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 14/14 [02:23<00:00, 10.24s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
2/100
0it [00:00, ?it/s]
3/100
6%|███████████████▋ | 1/18 [01:16<21:40, 76.48s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.141e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 18/18 [02:30<00:00, 8.38s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.859e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
4/100
0it [00:00, ?it/s]
5/100
0it [00:00, ?it/s]
6/100
0it [00:00, ?it/s]
7/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [01:14<00:00, 9.30s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
8/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [01:14<00:00, 9.35s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
9/100
0it [00:00, ?it/s]
10/100
5%|████████████▊ | 1/22 [01:16<26:38, 76.12s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 22/22 [02:31<00:00, 6.87s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
11/100
6%|█████████████████▋ | 1/16 [01:15<18:54, 75.63s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [02:30<00:00, 9.43s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
12/100
0it [00:00, ?it/s]
13/100
4%|██████████▉ | 1/26 [01:17<32:11, 77.24s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
50%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ | 13/26 [02:38<02:19, 10.69s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 26/26 [03:56<00:00, 9.11s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
14/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [01:15<00:00, 9.45s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
15/100
5%|██▏ | 1/20 [01:16<24:08, 76.24s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|███████████████████████████████████████████| 20/20 [02:30<00:00, 7.52s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
16/100
0it [00:00, ?it/s]
17/100
6%|█████████████████▋ | 1/16 [01:15<18:56, 75.77s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [02:28<00:00, 9.29s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
18/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 4/4 [01:13<00:00, 18.35s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
19/100
0it [00:00, ?it/s]
20/100
4%|███████████▊ | 1/24 [01:16<29:26, 76.80s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 24/24 [02:31<00:00, 6.32s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
21/100
0it [00:00, ?it/s]
22/100
3%|█████████▍ | 1/30 [01:15<36:27, 75.42s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
67%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ | 20/30 [01:27<00:28, 2.85s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 30/30 [01:30<00:00, 3.02s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
23/100
0it [00:00, ?it/s]
24/100
3%|████████▎ | 1/34 [00:07<04:19, 7.87s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
41%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ | 14/34 [00:14<00:18, 1.08it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
68%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▊ | 23/34 [00:14<00:04, 2.53it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 34/34 [00:19<00:00, 1.70it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
25/100
24%|███████████████████████████████████████████████████████████████████ | 9/38 [00:05<00:13, 2.20it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
58%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▎ | 22/38 [00:11<00:06, 2.60it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
71%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▎ | 27/38 [00:17<00:06, 1.62it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 38/38 [00:19<00:00, 1.93it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
26/100
0it [00:00, ?it/s]
27/100
0it [00:00, ?it/s]
28/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 4/4 [00:03<00:00, 1.09it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
29/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 4/4 [00:03<00:00, 1.14it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
30/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [00:05<00:00, 1.40it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
31/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [00:03<00:00, 2.02it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
32/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [00:04<00:00, 1.79it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
33/100
0it [00:00, ?it/s]
34/100
6%|█████████████████▋ | 1/16 [00:05<01:18, 5.21s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:08<00:00, 1.97it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
35/100
6%|█████████████████▋ | 1/16 [00:05<01:18, 5.25s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:09<00:00, 1.70it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
36/100
0it [00:00, ?it/s]
37/100
0it [00:00, ?it/s]
38/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [00:04<00:00, 1.99it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
39/100
0it [00:00, ?it/s]
40/100
0it [00:00, ?it/s]
41/100
6%|█████████████████▋ | 1/16 [00:04<01:10, 4.68s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:07<00:00, 2.16it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
42/100
0it [00:00, ?it/s]
43/100
0it [00:00, ?it/s]
44/100
0it [00:00, ?it/s]
45/100
6%|█████████████████▋ | 1/16 [00:05<01:17, 5.18s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:07<00:00, 2.02it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
46/100
0it [00:00, ?it/s]
47/100
4%|██████████ | 1/28 [00:05<02:15, 5.03s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
46%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▉ | 13/28 [00:10<00:10, 1.49it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 28/28 [00:12<00:00, 2.21it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
48/100
6%|█████████████████▋ | 1/16 [00:05<01:17, 5.17s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:07<00:00, 2.02it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
49/100
0it [00:00, ?it/s]
50/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [00:05<00:00, 1.51it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
51/100
3%|███████▊ | 1/36 [00:05<02:56, 5.03s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
36%|█████████████████████████████████████████████████████████████████████████████████████████████████████▊ | 13/36 [00:09<00:14, 1.54it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 36/36 [00:15<00:00, 2.39it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
52/100
4%|███████████▊ | 1/24 [00:05<02:14, 5.83s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 24/24 [00:11<00:00, 2.18it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
53/100
0it [00:00, ?it/s]
54/100
0it [00:00, ?it/s]
55/100
3%|████████▊ | 1/32 [00:10<05:29, 10.64s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
66%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ | 21/32 [00:15<00:05, 1.94it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 32/32 [00:19<00:00, 1.66it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
56/100
0it [00:00, ?it/s]
57/100
6%|█████████████████▋ | 1/16 [00:05<01:20, 5.37s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:08<00:00, 1.94it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
58/100
0it [00:00, ?it/s]
59/100
3%|███████▊ | 1/36 [00:05<03:08, 5.38s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
64%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▏ | 23/36 [00:10<00:04, 3.07it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 36/36 [00:15<00:00, 2.35it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
60/100
3%|████████▊ | 1/32 [00:04<02:33, 4.95s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
41%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████▌ | 13/32 [00:10<00:13, 1.44it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 32/32 [00:14<00:00, 2.22it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
61/100
0it [00:00, ?it/s]
62/100
2%|███████ | 1/40 [00:06<04:00, 6.18s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
50%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ | 20/40 [00:11<00:08, 2.44it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
88%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▊ | 35/40 [00:16<00:01, 3.01it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 40/40 [00:19<00:00, 2.07it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
63/100
22%|███████████████████████████████████████████████████████████████▋ | 9/40 [00:05<00:13, 2.26it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
45%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▉ | 18/40 [00:11<00:11, 1.95it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
75%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▌ | 30/40 [00:17<00:04, 2.25it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 40/40 [00:19<00:00, 2.04it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
64/100
0it [00:00, ?it/s]
65/100
0it [00:00, ?it/s]
66/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [00:03<00:00, 2.03it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
67/100
4%|███████████▊ | 1/24 [00:05<02:06, 5.50s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 24/24 [00:10<00:00, 2.20it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
68/100
0it [00:00, ?it/s]
69/100
4%|██████████ | 1/28 [00:06<02:42, 6.01s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
46%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▉ | 13/28 [00:11<00:11, 1.32it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 28/28 [00:14<00:00, 1.94it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
70/100
0it [00:00, ?it/s]
71/100
0it [00:00, ?it/s]
72/100
0it [00:00, ?it/s]
73/100
6%|█████████████████▋ | 1/16 [00:04<01:13, 4.91s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:07<00:00, 2.12it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
74/100
0it [00:00, ?it/s]
75/100
0it [00:00, ?it/s]
76/100
2%|███████ | 1/40 [00:04<03:04, 4.74s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
52%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ | 21/40 [00:10<00:06, 2.74it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
62%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▎ | 25/40 [00:15<00:08, 1.70it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 40/40 [00:17<00:00, 2.24it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
77/100
0it [00:00, ?it/s]
78/100
0it [00:00, ?it/s]
79/100
2%|██████▍ | 1/44 [00:04<03:27, 4.82s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
52%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▍ | 23/44 [00:09<00:06, 3.25it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
64%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▍ | 28/44 [00:14<00:07, 2.15it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 44/44 [00:18<00:00, 2.44it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
80/100
6%|█████████████████▋ | 1/16 [00:05<01:19, 5.33s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:08<00:00, 1.99it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
81/100
0it [00:00, ?it/s]
82/100
6%|█████████████████▋ | 1/16 [00:05<01:21, 5.46s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
62%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▎ | 10/16 [00:05<00:02, 2.78it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:08<00:00, 1.88it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
83/100
21%|██████████████████████████████████████████████████████████▊ | 10/48 [00:05<00:15, 2.49it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
33%|██████████████████████████████████████████████████████████████████████████████████████████████ | 16/48 [00:10<00:19, 1.68it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
60%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▍ | 29/48 [00:15<00:09, 1.95it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 48/48 [00:20<00:00, 2.38it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
84/100
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 4/4 [00:03<00:00, 1.12it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
85/100
0it [00:00, ?it/s]
86/100
6%|█████████████████▋ | 1/16 [00:06<01:40, 6.71s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:09<00:00, 1.67it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
87/100
0it [00:00, ?it/s]
88/100
6%|█████████████████▋ | 1/16 [00:08<02:02, 8.16s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:10<00:00, 1.47it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
89/100
2%|█████ | 1/56 [00:04<04:01, 4.39s/it]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
23%|█████████████████████████████████████████████████████████████████▍ | 13/56 [00:09<00:26, 1.64it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
45%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▉ | 25/56 [00:13<00:16, 1.85it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
66%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▎ | 37/56 [00:18<00:08, 2.15it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 56/56 [00:22<00:00, 2.51it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
90/100
6%|█████████████████▋ | 1/16 [00:05<01:16, 5.12s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:07<00:00, 2.03it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
91/100
0it [00:00, ?it/s]
92/100
6%|█████████████████▋ | 1/16 [00:05<01:23, 5.59s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:08<00:00, 1.95it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
93/100
3%|████████▊ | 1/32 [00:04<02:34, 4.99s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
41%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████▌ | 13/32 [00:10<00:12, 1.47it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 32/32 [00:13<00:00, 2.36it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
94/100
0it [00:00, ?it/s]
95/100
2%|██████▍ | 1/44 [00:04<03:26, 4.80s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
30%|███████████████████████████████████████████████████████████████████████████████████▎ | 13/44 [00:09<00:19, 1.60it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
57%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▏ | 25/44 [00:14<00:09, 1.94it/s]
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 44/44 [00:18<00:00, 2.42it/s]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
96/100
6%|█████████████████▋ | 1/16 [01:20<20:00, 80.01s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [02:35<00:00, 9.73s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
97/100
0it [00:00, ?it/s]
98/100
2%|█████▍ | 1/52 [01:16<1:05:11, 76.69s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
25%|██████████████████████████████████████████████████████████████████████▌ | 13/52 [02:33<06:40, 10.26s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
48%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▌ | 25/52 [03:49<03:36, 8.01s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
71%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████▋ | 37/52 [05:05<01:48, 7.24s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 1.800e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.200e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 52/52 [06:20<00:00, 7.31s/it]
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
Volume = 88.8264396098042 m^3
Heating power = 2.000e+06 W
99/100
0it [00:00, ?it/s]
[34]:
plot_te(analysis)
[35]:
# Plot the slices again.
plot_grid_2D()
Post processing
There are a number of post-processing step we can take. Below we show the ‘adaptation table’ which displays which multi indices were refined at every iteration.
[36]:
analysis.adaptation_table()
plt.savefig('Adaptation_table_DASC.png')
plt.savefig('Adaptation_table_DASC.pdf')
<Figure size 640x480 with 0 Axes>
[37]:
analysis.adaptation_table?
Signature: analysis.adaptation_table(**kwargs)
Docstring:
Plots a color-coded table of the quadrature-order refinement.
Shows in what order the parameters were refined, and unlike
adaptation_histogram, this also shows higher-order refinements.
Parameters
----------
**kwargs: can contain kwarg 'order' to specify the order in which
the variables on the x axis are plotted (e.g. in order of decreasing
1st order Sobol index).
Returns
-------
None.
File: /Volumes/UserData/dpc/GIT/EasyVVUQ/env_3.12/lib/python3.12/site-packages/easyvvuq/analysis/sc_analysis.py
Type: method
We can also make a histogram which visualises the adaptation. This displays only a first-order information, i.e. only the maximum quadrature order per input. It therefore does not display that certain inputs were refined simultaneously:
[38]:
analysis.adaptation_histogram()
To get a list of the error magnitudes associated to the multi indices that were selected use:
[39]:
plt.semilogy(analysis.get_adaptation_errors())
[39]:
[<matplotlib.lines.Line2D at 0x17a867d40>]
This shows a nice decrease of the error.
To compute the mean and variance of the code output we use:
[40]:
analysis.get_uncertainty_amplification('te')
-----------------
Mean CV input = 21.3848 %
Mean CV output = 18.5262 %
Uncertainty amplification factor = 0.1853/0.2138 = 0.8663
-----------------
[40]:
0.8663242668726164
[41]:
plot_te(analysis)
We can see that the estimates are fair, although not yet fully converged in this case. Note however that these results are computed only with the accepted set of multi indices. At the end, we can merge the accepted and admissible set (thereby using all samples), and recompute the results:
[42]:
analysis.merge_accepted_and_admissible()
df = my_campaign.get_collation_result()
results = analysis.analyse(df)
[43]:
frms_mean, frms = test_surrogate()
S.append([frms_mean, frms])
[44]:
plot_te(analysis)
[45]:
# Plot the slices again.
plot_grid_2D()
plt.savefig('Grid_DASC.png')
plt.savefig('Grid_DASC.pdf')
This improved our estimates. Note however, that if we would refine again from this point, the new admissble set will be very large, since we added all previous admissible indices to the accepted set. This opens up a wide range of possible new candidate directions, making the corresponding ensemble very large.
Thus if we are still not happy about the result, we first have to undo the merging via analysis.undo_merge(), before refining again.
[46]:
# This will undo the merge, and reproduce the old results
#analysis.undo_merge()
#df = campaign.get_collation_result()
#results = analysis.analyse(df)
#print('Mean = %.4e' % results.describe('f', 'mean'))
#print('Standard deviation = %.4e' % results.describe('f', 'std'))
[47]:
analysis.plot_stat_convergence()
We can also display the Sobol sensitivity indices via:
[48]:
sobols = []
# retrieve the Sobol indices from the results object
params = list(sampler.vary.get_keys())
for param in params:
sobols.append(results._get_sobols_first('te', param))
sobols = np.array(sobols).mean(axis=1)
# make a bar chart
fig = plt.figure()
ax = fig.add_subplot(111, title='First-order Sobol indices (mean over rho)')
ax.bar(range(len(sobols)), height=np.array(sobols).flatten())
ax.set_xticks(range(len(sobols)))
ax.set_xticklabels(params)
plt.xticks(rotation=90)
plt.tight_layout()
plt.savefig('Average_Sobol_DASC.png')
plt.savefig('Average_Sobol_DASC.pdf')
[49]:
dict(zip(params, sobols))
[49]:
{'Qe_tot': 0.0826519003065851,
'H0': 0.019762879974652064,
'Hw': 0.21007714324357513,
'chi': 0.34090854895953654,
'Te_bc': 0.025553700331106294,
'b_pos': 0.06080314397550473,
'b_height': 0.20208219034487726,
'b_sol': 0.030533812380133626,
'b_width': 0.000842415509072117,
'b_slope': 0.0017000632610276914}
[50]:
te_mean, te_std = analysis.get_moments('te')[0], np.sqrt(analysis.get_moments('te')[1])
rho_mean, rho_std = analysis.get_moments('rho')[0], np.sqrt(analysis.get_moments('rho')[1])
We can also plot the Sobol first order indices as a function of rho
[51]:
# plot the Sobol first order indices as a function of rho
sobols = []
# retrieve the Sobol indices from the results object
params = list(sampler.vary.get_keys())
for param in params:
sobols.append(results._get_sobols_first('te', param))
fig = plt.figure()
plt.plot(rho_mean, np.array(sobols).T, label=params)
plt.legend(loc=0)
plt.xlabel('rho [m]')
plt.ylabel('First order Sobol indices');
[52]:
# plot the performance of the SC surrogates
if __name__ == '__main__':
for i in range(len(S)):
plt.semilogy(rho_mean, S[i][1])
plt.xlabel('rho [m]') ; plt.ylabel('fractional RMS for predicted Te') #; plt.legend(loc=0, ncol=5)
plt.title('Performance of the SC surrogate')
plt.show()
# plt.savefig('SC_surrogate.png')
# plt.savefig('SC_surrogate.pdf')
[ ]:
[53]:
# plot the convergence of the surrogate based on 1000 random points
if __name__ == '__main__':
plt.figure()
plt.semilogy([s[0] for s in S], 'o-')
plt.xlabel('element number')
plt.ylabel('fractional RMS error for the SC surrogate')
plt.legend(loc=0)
plt.savefig('Convergence_DASC_surrogate.png')
plt.savefig('Convergence_DASC_surrogate.pdf')
/var/folders/6f/rn14629n60j16dc99dtk7bs4000ctx/T/ipykernel_27273/1120173477.py:7: UserWarning: No artists with labels found to put in legend. Note that artists whose label start with an underscore are ignored when legend() is called with no argument.
plt.legend(loc=0)
[ ]: