1D Poisson's Equation

This example is taken from [4]. Consider a simple 1D Poisson’s equation with Dirichlet boundary conditions. The solution is given by

\[u(x)=\sin (2 \pi x)+0.1 \sin (50 \pi x)\]

using ModelingToolkit, IntervalSets, Sophon
using Optimization, OptimizationOptimJL, Zygote
using CairoMakie

@parameters x
@variables u(..)
Dₓ² = Differential(x)^2

f(x) = -4 * π^2 * sin(2 * π * x) - 250 * π^2 * sin(50 * π * x)
eq = Dₓ²(u(x)) ~ f(x)
domain = [x ∈ 0 .. 1]
bcs = [u(0) ~ 0, u(1) ~ 0]

@named poisson = PDESystem(eq, bcs, domain, [x], [u(x)])

\[ \begin{align} \frac{\mathrm{d}}{\mathrm{d}x} \frac{\mathrm{d} u\left( x \right)}{\mathrm{d}x} =& - 39.478 \sin\left( 6.2832 x \right) - 2467.4 \sin\left( 157.08 x \right) \end{align} \]

chain = Siren(1, 16, 32, 16, 1)
pinn = PINN(chain)
sampler = QuasiRandomSampler(200, 1)
strategy = NonAdaptiveTraining(1 , 50)

prob = Sophon.discretize(poisson, pinn, sampler, strategy)
@showprogress res = Optimization.solve(prob, BFGS(); maxiters=2000)

phi = pinn.phi
xs = 0:0.001:1
u_true = @. sin(2 * pi * xs) + 0.1 * sin(50 * pi * xs)
us = phi(xs', res.u)
fig = Figure()
axis = Axis(fig[1, 1])
lines!(xs, u_true; label="Ground Truth")
lines!(xs, vec(us); label="Prediction")
axislegend(axis)
fig

0.0%┣                                         ┫ 1/2.0k [00:20<Inf:Inf, InfGs/it]

3.687249e+06 0.1%┣                             ┫ 2/2.0k [00:20<11:01:30, 20s/it]

3.016728e+06 0.1%┣                             ┫ 3/2.0k [00:20<05:31:49, 10s/it]

2.799223e+06 0.2%┣                              ┫ 4/2.0k [00:20<03:42:11, 7s/it]

2.593727e+06 0.2%┣                              ┫ 5/2.0k [00:20<02:47:45, 5s/it]

2.349256e+06 0.3%┣                              ┫ 6/2.0k [00:20<02:14:32, 4s/it]

1.967595e+06 0.3%┣                              ┫ 7/2.0k [00:20<01:52:32, 3s/it]

1.874130e+06 0.4%┣▏                             ┫ 8/2.0k [00:20<01:36:41, 3s/it]

1.850949e+06 0.4%┣▏                             ┫ 9/2.0k [00:20<01:24:53, 3s/it]

1.734018e+06 0.5%┣▏                            ┫ 10/2.0k [00:21<01:15:57, 2s/it]

1.667616e+06 0.5%┣▏                            ┫ 11/2.0k [00:21<01:08:34, 2s/it]

1.529575e+06 0.6%┣▏                            ┫ 12/2.0k [00:21<01:02:30, 2s/it]

1.452332e+06 0.6%┣▏                               ┫ 13/2.0k [00:21<57:29, 2s/it]

1.396178e+06 0.7%┣▎                               ┫ 14/2.0k [00:21<53:15, 2s/it]

1.318164e+06 0.7%┣▎                               ┫ 15/2.0k [00:21<49:42, 2s/it]

1.246803e+06 0.8%┣▎                               ┫ 16/2.0k [00:21<46:30, 1s/it]

1.182296e+06 0.8%┣▎                               ┫ 17/2.0k [00:21<43:43, 1s/it]

1.110743e+06 0.9%┣▎                               ┫ 18/2.0k [00:21<41:15, 1s/it]

1.051656e+06 0.9%┣▎                               ┫ 19/2.0k [00:21<39:05, 1s/it]

9.753187e+05 1.0%┣▎                               ┫ 20/2.0k [00:21<37:07, 1s/it]

9.544630e+05 1.0%┣▍                               ┫ 21/2.0k [00:21<35:25, 1s/it]

8.294083e+05 1.1%┣▍                               ┫ 22/2.0k [00:22<33:49, 1s/it]

7.532046e+05 1.1%┣▍                               ┫ 23/2.0k [00:22<32:23, 1it/s]

6.503855e+05 1.2%┣▍                               ┫ 24/2.0k [00:22<31:04, 1it/s]

5.712806e+05 1.2%┣▍                               ┫ 25/2.0k [00:22<29:52, 1it/s]

5.229110e+05 1.3%┣▍                               ┫ 26/2.0k [00:22<28:47, 1it/s]

4.810070e+05 1.3%┣▍                               ┫ 27/2.0k [00:22<27:45, 1it/s]

4.231839e+05 1.4%┣▌                               ┫ 28/2.0k [00:22<26:47, 1it/s]

3.748783e+05 1.4%┣▌                               ┫ 29/2.0k [00:22<25:54, 1it/s]

3.325009e+05 1.5%┣▌                               ┫ 30/2.0k [00:22<25:04, 1it/s]

2.932228e+05 1.5%┣▌                               ┫ 31/2.0k [00:22<24:21, 1it/s]

2.296430e+05 1.6%┣▌                               ┫ 32/2.0k [00:22<23:37, 1it/s]

2.077682e+05 1.6%┣▌                               ┫ 33/2.0k [00:22<22:56, 1it/s]

1.803278e+05 1.7%┣▌                               ┫ 34/2.0k [00:22<22:19, 1it/s]

1.591974e+05 1.7%┣▋                               ┫ 35/2.0k [00:23<21:43, 2it/s]

1.418064e+05 1.8%┣▋                               ┫ 36/2.0k [00:23<21:10, 2it/s]

1.232340e+05 1.8%┣▋                               ┫ 37/2.0k [00:23<20:38, 2it/s]

1.188881e+05 1.9%┣▋                               ┫ 38/2.0k [00:23<20:07, 2it/s]

1.077311e+05 1.9%┣▋                               ┫ 39/2.0k [00:23<19:38, 2it/s]

9.826872e+04 2.0%┣▋                               ┫ 40/2.0k [00:23<19:10, 2it/s]

8.826507e+04 2.0%┣▋                               ┫ 41/2.0k [00:23<18:44, 2it/s]

8.187875e+04 2.1%┣▊                               ┫ 42/2.0k [00:23<18:20, 2it/s]

7.322824e+04 2.1%┣▊                               ┫ 43/2.0k [00:23<17:56, 2it/s]

6.631895e+04 2.2%┣▊                               ┫ 44/2.0k [00:23<17:33, 2it/s]

5.893821e+04 2.2%┣▊                               ┫ 45/2.0k [00:23<17:12, 2it/s]

5.671810e+04 2.3%┣▊                               ┫ 46/2.0k [00:23<16:51, 2it/s]

4.758161e+04 2.3%┣▊                               ┫ 47/2.0k [00:23<16:33, 2it/s]

4.065279e+04 2.4%┣▊                               ┫ 48/2.0k [00:23<16:13, 2it/s]

3.830858e+04 2.4%┣▉                               ┫ 49/2.0k [00:23<15:55, 2it/s]

3.204164e+04 2.5%┣▉                               ┫ 50/2.0k [00:24<15:38, 2it/s]

2.704826e+04 2.5%┣▉                               ┫ 51/2.0k [00:24<15:21, 2it/s]

2.329090e+04 2.6%┣▉                               ┫ 52/2.0k [00:24<15:05, 2it/s]

2.060940e+04 2.6%┣▉                               ┫ 53/2.0k [00:24<14:50, 2it/s]

1.863391e+04 2.7%┣▉                               ┫ 54/2.0k [00:24<14:35, 2it/s]

1.445327e+04 2.8%┣█                               ┫ 56/2.0k [00:24<14:06, 2it/s]

1.222592e+04 2.9%┣█                               ┫ 58/2.0k [00:24<13:39, 2it/s]

1.086677e+04 2.9%┣█                               ┫ 59/2.0k [00:24<13:26, 2it/s]

9.947117e+03 3.0%┣█                               ┫ 60/2.0k [00:24<13:14, 2it/s]

9.297021e+03 3.0%┣█                               ┫ 61/2.0k [00:24<13:03, 2it/s]

8.482200e+03 3.1%┣█                               ┫ 62/2.0k [00:24<12:51, 3it/s]

8.136126e+03 3.1%┣█                               ┫ 63/2.0k [00:24<12:40, 3it/s]

7.028370e+03 3.2%┣█                               ┫ 65/2.0k [00:24<12:18, 3it/s]

6.566331e+03 3.3%┣█                               ┫ 66/2.0k [00:24<12:08, 3it/s]

5.969739e+03 3.3%┣█                               ┫ 67/2.0k [00:25<11:59, 3it/s]

5.312190e+03 3.4%┣█                               ┫ 69/2.0k [00:25<11:40, 3it/s]

5.130354e+03 3.5%┣█▏                              ┫ 70/2.0k [00:25<11:31, 3it/s]

4.828151e+03 3.5%┣█▏                              ┫ 71/2.0k [00:25<11:22, 3it/s]

4.451421e+03 3.6%┣█▏                              ┫ 72/2.0k [00:25<11:14, 3it/s]

4.045530e+03 3.7%┣█▏                              ┫ 74/2.0k [00:25<10:58, 3it/s]

3.836022e+03 3.7%┣█▏                              ┫ 75/2.0k [00:25<10:50, 3it/s]

3.399210e+03 3.8%┣█▎                              ┫ 77/2.0k [00:25<10:34, 3it/s]

3.272857e+03 3.9%┣█▎                              ┫ 78/2.0k [00:25<10:27, 3it/s]

3.087247e+03 4.0%┣█▎                              ┫ 80/2.0k [00:25<10:13, 3it/s]

2.927304e+03 4.0%┣█▎                              ┫ 81/2.0k [00:25<10:07, 3it/s]

2.815682e+03 4.1%┣█▎                              ┫ 82/2.0k [00:25<10:00, 3it/s]

2.693424e+03 4.1%┣█▎                              ┫ 83/2.0k [00:25<09:54, 3it/s]

2.569423e+03 4.2%┣█▍                              ┫ 84/2.0k [00:25<09:47, 3it/s]

2.205547e+03 4.3%┣█▍                              ┫ 86/2.0k [00:26<09:35, 3it/s]

1.997215e+03 4.4%┣█▍                              ┫ 88/2.0k [00:26<09:25, 3it/s]

1.731550e+03 4.5%┣█▍                              ┫ 90/2.0k [00:26<09:14, 3it/s]

1.541743e+03 4.6%┣█▌                              ┫ 92/2.0k [00:26<09:03, 4it/s]

1.386163e+03 4.7%┣█▌                              ┫ 94/2.0k [00:26<08:53, 4it/s]

1.301592e+03 4.7%┣█▌                              ┫ 95/2.0k [00:26<08:48, 4it/s]

1.138224e+03 4.8%┣█▌                              ┫ 97/2.0k [00:26<08:39, 4it/s]

1.026856e+03 4.9%┣█▋                              ┫ 98/2.0k [00:26<08:34, 4it/s]

8.500193e+02 5.0%┣█▌                             ┫ 100/2.0k [00:26<08:25, 4it/s]

7.499472e+02 5.1%┣█▋                             ┫ 102/2.0k [00:26<08:16, 4it/s]

6.905464e+02 5.1%┣█▋                             ┫ 103/2.0k [00:26<08:12, 4it/s]

6.086530e+02 5.2%┣█▋                             ┫ 105/2.0k [00:27<08:04, 4it/s]

5.506066e+02 5.3%┣█▋                             ┫ 107/2.0k [00:27<07:56, 4it/s]

4.674990e+02 5.4%┣█▊                             ┫ 109/2.0k [00:27<07:48, 4it/s]

4.396545e+02 5.5%┣█▊                             ┫ 110/2.0k [00:27<07:45, 4it/s]

4.137001e+02 5.5%┣█▊                             ┫ 111/2.0k [00:27<07:41, 4it/s]

3.810266e+02 5.6%┣█▊                             ┫ 113/2.0k [00:27<07:34, 4it/s]

3.502825e+02 5.7%┣█▉                             ┫ 115/2.0k [00:27<07:27, 4it/s]

3.116843e+02 5.8%┣█▉                             ┫ 117/2.0k [00:27<07:20, 4it/s]

2.945479e+02 5.9%┣█▉                             ┫ 118/2.0k [00:27<07:17, 4it/s]

2.607077e+02 6.0%┣█▉                             ┫ 120/2.0k [00:27<07:11, 4it/s]

2.409911e+02 6.1%┣██                             ┫ 122/2.0k [00:27<07:05, 4it/s]

2.254114e+02 6.2%┣██                             ┫ 124/2.0k [00:27<06:58, 4it/s]

2.088430e+02 6.3%┣██                             ┫ 126/2.0k [00:28<06:53, 5it/s]

2.029867e+02 6.3%┣██                             ┫ 127/2.0k [00:28<06:50, 5it/s]

1.918758e+02 6.4%┣██                             ┫ 129/2.0k [00:28<06:44, 5it/s]

1.785412e+02 6.5%┣██                             ┫ 131/2.0k [00:28<06:39, 5it/s]

1.659480e+02 6.6%┣██                             ┫ 133/2.0k [00:28<06:34, 5it/s]

1.526760e+02 6.7%┣██                             ┫ 135/2.0k [00:28<06:29, 5it/s]

1.427205e+02 6.8%┣██▏                            ┫ 137/2.0k [00:28<06:23, 5it/s]

1.319246e+02 6.9%┣██▏                            ┫ 139/2.0k [00:28<06:18, 5it/s]

1.265358e+02 7.0%┣██▏                            ┫ 141/2.0k [00:28<06:14, 5it/s]

1.190565e+02 7.1%┣██▏                            ┫ 143/2.0k [00:28<06:09, 5it/s]

1.082568e+02 7.2%┣██▎                            ┫ 145/2.0k [00:28<06:05, 5it/s]

1.019797e+02 7.3%┣██▎                            ┫ 147/2.0k [00:28<06:00, 5it/s]

9.600458e+01 7.4%┣██▎                            ┫ 149/2.0k [00:28<05:56, 5it/s]

8.992210e+01 7.5%┣██▍                            ┫ 151/2.0k [00:29<05:52, 5it/s]

8.135230e+01 7.6%┣██▍                            ┫ 153/2.0k [00:29<05:47, 5it/s]

7.686698e+01 7.7%┣██▍                            ┫ 154/2.0k [00:29<05:46, 5it/s]

6.904921e+01 7.8%┣██▍                            ┫ 156/2.0k [00:29<05:42, 5it/s]

6.495923e+01 7.9%┣██▌                            ┫ 158/2.0k [00:29<05:38, 5it/s]

6.229605e+01 8.0%┣██▌                            ┫ 160/2.0k [00:29<05:34, 6it/s]

5.974095e+01 8.0%┣██▌                            ┫ 161/2.0k [00:29<05:32, 6it/s]

5.566049e+01 8.1%┣██▌                            ┫ 163/2.0k [00:29<05:29, 6it/s]

5.191041e+01 8.2%┣██▋                            ┫ 165/2.0k [00:29<05:25, 6it/s]

4.858172e+01 8.3%┣██▋                            ┫ 167/2.0k [00:29<05:22, 6it/s]

4.650186e+01 8.4%┣██▋                            ┫ 169/2.0k [00:29<05:18, 6it/s]

4.397840e+01 8.5%┣██▋                            ┫ 171/2.0k [00:29<05:15, 6it/s]

4.252310e+01 8.6%┣██▊                            ┫ 173/2.0k [00:29<05:11, 6it/s]

4.219990e+01 8.7%┣██▊                            ┫ 174/2.0k [00:29<05:10, 6it/s]

4.010396e+01 8.8%┣██▊                            ┫ 176/2.0k [00:29<05:07, 6it/s]

3.886346e+01 8.9%┣██▊                            ┫ 178/2.0k [00:30<05:04, 6it/s]

3.702093e+01 9.0%┣██▉                            ┫ 180/2.0k [00:30<05:01, 6it/s]

3.487924e+01 9.1%┣██▉                            ┫ 182/2.0k [00:30<04:58, 6it/s]

3.348214e+01 9.2%┣██▉                            ┫ 184/2.0k [00:30<04:55, 6it/s]

3.274482e+01 9.2%┣██▉                            ┫ 185/2.0k [00:30<04:54, 6it/s]

3.104629e+01 9.3%┣███                            ┫ 187/2.0k [00:30<04:51, 6it/s]

2.923060e+01 9.4%┣███                            ┫ 189/2.0k [00:30<04:48, 6it/s]

2.735956e+01 9.5%┣███                            ┫ 191/2.0k [00:30<04:45, 6it/s]

2.603638e+01 9.6%┣███                            ┫ 193/2.0k [00:30<04:43, 6it/s]

2.472359e+01 9.7%┣███                            ┫ 195/2.0k [00:30<04:40, 6it/s]

2.328968e+01 9.8%┣███                            ┫ 197/2.0k [00:30<04:38, 6it/s]

2.188993e+01 9.9%┣███                            ┫ 199/2.0k [00:30<04:35, 7it/s]

2.066509e+01 10.0%┣███                           ┫ 201/2.0k [00:30<04:33, 7it/s]

1.922760e+01 10.2%┣███                           ┫ 204/2.0k [00:30<04:29, 7it/s]

1.828871e+01 10.3%┣███                           ┫ 206/2.0k [00:30<04:27, 7it/s]

1.743255e+01 10.4%┣███▏                          ┫ 208/2.0k [00:31<04:24, 7it/s]

1.636769e+01 10.5%┣███▏                          ┫ 210/2.0k [00:31<04:22, 7it/s]

1.556716e+01 10.6%┣███▏                          ┫ 212/2.0k [00:31<04:20, 7it/s]

1.495780e+01 10.7%┣███▏                          ┫ 214/2.0k [00:31<04:18, 7it/s]

1.382506e+01 10.8%┣███▎                          ┫ 216/2.0k [00:31<04:15, 7it/s]

1.352458e+01 10.8%┣███▎                          ┫ 217/2.0k [00:31<04:15, 7it/s]

1.292824e+01 10.9%┣███▎                          ┫ 219/2.0k [00:31<04:12, 7it/s]

1.217463e+01 11.0%┣███▎                          ┫ 221/2.0k [00:31<04:10, 7it/s]

1.131087e+01 11.1%┣███▍                          ┫ 223/2.0k [00:31<04:08, 7it/s]

1.046380e+01 11.2%┣███▍                          ┫ 225/2.0k [00:31<04:06, 7it/s]

9.773248e+00 11.3%┣███▍                          ┫ 227/2.0k [00:31<04:04, 7it/s]

8.958232e+00 11.4%┣███▍                          ┫ 229/2.0k [00:31<04:02, 7it/s]

8.192393e+00 11.6%┣███▌                          ┫ 232/2.0k [00:31<03:59, 7it/s]

7.816585e+00 11.7%┣███▌                          ┫ 234/2.0k [00:31<03:58, 7it/s]

7.352166e+00 11.8%┣███▌                          ┫ 236/2.0k [00:31<03:56, 7it/s]

6.679582e+00 11.9%┣███▋                          ┫ 239/2.0k [00:31<03:53, 8it/s]

6.354195e+00 12.1%┣███▋                          ┫ 242/2.0k [00:32<03:50, 8it/s]

6.145783e+00 12.2%┣███▋                          ┫ 244/2.0k [00:32<03:48, 8it/s]

5.992693e+00 12.3%┣███▊                          ┫ 247/2.0k [00:32<03:46, 8it/s]

5.652736e+00 12.5%┣███▊                          ┫ 250/2.0k [00:32<03:43, 8it/s]

5.365139e+00 12.6%┣███▉                          ┫ 253/2.0k [00:32<03:41, 8it/s]

5.140968e+00 12.7%┣███▉                          ┫ 255/2.0k [00:32<03:39, 8it/s]

4.924108e+00 12.8%┣███▉                          ┫ 257/2.0k [00:32<03:37, 8it/s]

4.709006e+00 13.0%┣████                          ┫ 260/2.0k [00:32<03:35, 8it/s]

4.569934e+00 13.1%┣████                          ┫ 262/2.0k [00:32<03:33, 8it/s]

4.262312e+00 13.2%┣████                          ┫ 265/2.0k [00:32<03:31, 8it/s]

4.074901e+00 13.3%┣████                          ┫ 267/2.0k [00:32<03:30, 8it/s]

3.895033e+00 13.4%┣████                          ┫ 269/2.0k [00:32<03:28, 8it/s]

3.708010e+00 13.6%┣████                          ┫ 272/2.0k [00:32<03:26, 8it/s]

3.517270e+00 13.7%┣████▏                         ┫ 275/2.0k [00:32<03:24, 8it/s]

3.334877e+00 13.9%┣████▏                         ┫ 278/2.0k [00:32<03:22, 9it/s]

3.107648e+00 14.0%┣████▏                         ┫ 281/2.0k [00:33<03:20, 9it/s]

2.935605e+00 14.2%┣████▎                         ┫ 284/2.0k [00:33<03:18, 9it/s]

2.794504e+00 14.3%┣████▎                         ┫ 287/2.0k [00:33<03:16, 9it/s]

2.666167e+00 14.5%┣████▍                         ┫ 290/2.0k [00:33<03:14, 9it/s]

2.484609e+00 14.6%┣████▍                         ┫ 293/2.0k [00:33<03:12, 9it/s]

2.350675e+00 14.8%┣████▍                         ┫ 296/2.0k [00:33<03:10, 9it/s]

2.265988e+00 14.9%┣████▌                         ┫ 299/2.0k [00:33<03:08, 9it/s]

2.190719e+00 15.1%┣████▌                         ┫ 302/2.0k [00:33<03:06, 9it/s]

2.083366e+00 15.2%┣████▋                         ┫ 305/2.0k [00:33<03:04, 9it/s]

2.001454e+00 15.4%┣████▋                         ┫ 308/2.0k [00:33<03:03, 9it/s]

1.897563e+00 15.5%┣████▋                         ┫ 311/2.0k [00:33<03:01, 9it/s]

1.826827e+00 15.7%┣████▊                         ┫ 314/2.0k [00:33<02:59, 9it/s]

1.759040e+00 15.8%┣████▊                         ┫ 316/2.0k [00:33<02:58, 9it/s]

1.676413e+00 15.9%┣████▋                        ┫ 319/2.0k [00:33<02:57, 10it/s]

1.572775e+00 16.1%┣████▊                        ┫ 322/2.0k [00:33<02:55, 10it/s]

1.460022e+00 16.2%┣████▊                        ┫ 325/2.0k [00:33<02:53, 10it/s]

1.366123e+00 16.4%┣████▊                        ┫ 328/2.0k [00:34<02:52, 10it/s]

1.308952e+00 16.5%┣████▉                        ┫ 330/2.0k [00:34<02:51, 10it/s]

1.187019e+00 16.6%┣████▉                        ┫ 333/2.0k [00:34<02:49, 10it/s]

1.128470e+00 16.7%┣████▉                        ┫ 335/2.0k [00:34<02:48, 10it/s]

1.054182e+00 16.9%┣█████                        ┫ 338/2.0k [00:34<02:47, 10it/s]

9.891104e-01 17.0%┣█████                        ┫ 341/2.0k [00:34<02:45, 10it/s]

9.404362e-01 17.2%┣█████                        ┫ 344/2.0k [00:34<02:44, 10it/s]

9.242868e-01 17.3%┣█████                        ┫ 346/2.0k [00:34<02:43, 10it/s]

8.903608e-01 17.4%┣█████                        ┫ 349/2.0k [00:34<02:42, 10it/s]

8.396299e-01 17.6%┣█████▏                       ┫ 353/2.0k [00:34<02:40, 10it/s]

7.943028e-01 17.9%┣█████▏                       ┫ 358/2.0k [00:34<02:37, 10it/s]

7.615860e-01 18.0%┣█████▎                       ┫ 361/2.0k [00:34<02:36, 11it/s]

7.453710e-01 18.1%┣█████▎                       ┫ 363/2.0k [00:34<02:35, 11it/s]

7.067483e-01 18.3%┣█████▎                       ┫ 367/2.0k [00:34<02:34, 11it/s]

6.782401e-01 18.5%┣█████▍                       ┫ 370/2.0k [00:34<02:32, 11it/s]

6.373948e-01 18.7%┣█████▍                       ┫ 375/2.0k [00:35<02:30, 11it/s]

6.086799e-01 18.9%┣█████▌                       ┫ 379/2.0k [00:35<02:28, 11it/s]

5.870873e-01 19.1%┣█████▌                       ┫ 383/2.0k [00:35<02:27, 11it/s]

5.617611e-01 19.3%┣█████▋                       ┫ 387/2.0k [00:35<02:25, 11it/s]

5.304783e-01 19.5%┣█████▊                       ┫ 391/2.0k [00:35<02:24, 11it/s]

4.888558e-01 19.8%┣█████▊                       ┫ 396/2.0k [00:35<02:22, 11it/s]

4.600697e-01 20.0%┣█████▉                       ┫ 401/2.0k [00:35<02:20, 11it/s]

4.316363e-01 20.2%┣█████▉                       ┫ 405/2.0k [00:35<02:18, 12it/s]

4.075348e-01 20.5%┣██████                       ┫ 410/2.0k [00:35<02:16, 12it/s]

3.821375e-01 20.7%┣██████                       ┫ 415/2.0k [00:35<02:14, 12it/s]

3.604446e-01 20.9%┣██████                       ┫ 419/2.0k [00:35<02:13, 12it/s]

3.291269e-01 21.2%┣██████▏                      ┫ 424/2.0k [00:35<02:11, 12it/s]

3.124499e-01 21.4%┣██████▏                      ┫ 429/2.0k [00:35<02:10, 12it/s]

2.973808e-01 21.7%┣██████▎                      ┫ 434/2.0k [00:35<02:08, 12it/s]

2.824897e-01 21.8%┣██████▍                      ┫ 437/2.0k [00:35<02:07, 12it/s]

2.659404e-01 22.1%┣██████▍                      ┫ 442/2.0k [00:35<02:05, 12it/s]

2.500431e-01 22.3%┣██████▌                      ┫ 447/2.0k [00:36<02:04, 13it/s]

2.346858e-01 22.6%┣██████▌                      ┫ 452/2.0k [00:36<02:02, 13it/s]

2.189178e-01 22.9%┣██████▋                      ┫ 458/2.0k [00:36<02:00, 13it/s]

2.031678e-01 23.2%┣██████▊                      ┫ 465/2.0k [00:36<01:58, 13it/s]

1.877473e-01 23.5%┣██████▉                      ┫ 471/2.0k [00:36<01:56, 13it/s]

1.739303e-01 23.8%┣███████                      ┫ 477/2.0k [00:36<01:55, 13it/s]

1.648065e-01 24.1%┣███████                      ┫ 482/2.0k [00:36<01:53, 13it/s]

1.540131e-01 24.4%┣███████                      ┫ 488/2.0k [00:36<01:52, 14it/s]

1.431633e-01 24.7%┣███████▏                     ┫ 495/2.0k [00:36<01:50, 14it/s]

1.311931e-01 25.1%┣███████▎                     ┫ 502/2.0k [00:36<01:48, 14it/s]

1.210777e-01 25.4%┣███████▍                     ┫ 509/2.0k [00:36<01:46, 14it/s]

1.137955e-01 25.7%┣███████▌                     ┫ 514/2.0k [00:36<01:45, 14it/s]

1.053258e-01 26.0%┣███████▌                     ┫ 521/2.0k [00:36<01:43, 14it/s]

9.560892e-02 26.4%┣███████▋                     ┫ 528/2.0k [00:36<01:41, 15it/s]

8.799748e-02 26.7%┣███████▊                     ┫ 535/2.0k [00:36<01:40, 15it/s]

8.215089e-02 27.1%┣███████▉                     ┫ 542/2.0k [00:36<01:38, 15it/s]

8.052627e-02 27.2%┣███████▉                     ┫ 544/2.0k [00:36<01:38, 15it/s]

7.823881e-02 27.3%┣████████                     ┫ 547/2.0k [00:36<01:37, 15it/s]

7.429586e-02 27.7%┣████████                     ┫ 554/2.0k [00:37<01:36, 15it/s]

7.085252e-02 28.0%┣████████▏                    ┫ 560/2.0k [00:37<01:34, 15it/s]

6.686093e-02 28.3%┣████████▏                    ┫ 566/2.0k [00:37<01:33, 15it/s]

6.376571e-02 28.6%┣████████▎                    ┫ 572/2.0k [00:37<01:32, 16it/s]

6.166225e-02 28.8%┣████████▍                    ┫ 576/2.0k [00:37<01:31, 16it/s]

5.750279e-02 29.1%┣████████▌                    ┫ 583/2.0k [00:37<01:30, 16it/s]

5.350779e-02 29.5%┣████████▌                    ┫ 590/2.0k [00:37<01:28, 16it/s]

5.034308e-02 29.8%┣████████▋                    ┫ 597/2.0k [00:37<01:27, 16it/s]

4.775199e-02 30.2%┣████████▊                    ┫ 604/2.0k [00:37<01:26, 16it/s]

4.610843e-02 30.4%┣████████▉                    ┫ 608/2.0k [00:37<01:25, 16it/s]

4.333277e-02 30.7%┣█████████                    ┫ 615/2.0k [00:37<01:24, 17it/s]

4.055534e-02 31.0%┣█████████                    ┫ 621/2.0k [00:37<01:23, 17it/s]

3.786397e-02 31.3%┣█████████                    ┫ 627/2.0k [00:37<01:22, 17it/s]

3.553254e-02 31.6%┣█████████▏                   ┫ 633/2.0k [00:37<01:21, 17it/s]

3.359253e-02 31.9%┣█████████▎                   ┫ 639/2.0k [00:37<01:20, 17it/s]

3.122103e-02 32.3%┣█████████▍                   ┫ 646/2.0k [00:37<01:19, 17it/s]

2.905806e-02 32.6%┣█████████▌                   ┫ 653/2.0k [00:37<01:17, 17it/s]

2.705599e-02 32.9%┣█████████▌                   ┫ 659/2.0k [00:37<01:16, 18it/s]

2.549941e-02 33.2%┣█████████▋                   ┫ 665/2.0k [00:38<01:16, 18it/s]

2.422352e-02 33.5%┣█████████▊                   ┫ 671/2.0k [00:38<01:15, 18it/s]

2.285028e-02 33.9%┣█████████▉                   ┫ 678/2.0k [00:38<01:14, 18it/s]

2.148411e-02 34.2%┣██████████                   ┫ 684/2.0k [00:38<01:13, 18it/s]

2.021005e-02 34.5%┣██████████                   ┫ 690/2.0k [00:38<01:12, 18it/s]

1.917471e-02 34.8%┣██████████                   ┫ 696/2.0k [00:38<01:11, 18it/s]

1.819583e-02 35.1%┣██████████▏                  ┫ 702/2.0k [00:38<01:10, 18it/s]

1.739089e-02 35.4%┣██████████▎                  ┫ 709/2.0k [00:38<01:09, 19it/s]

1.665738e-02 35.8%┣██████████▍                  ┫ 716/2.0k [00:38<01:08, 19it/s]

1.577176e-02 36.1%┣██████████▌                  ┫ 723/2.0k [00:38<01:07, 19it/s]

1.482454e-02 36.5%┣██████████▋                  ┫ 730/2.0k [00:38<01:06, 19it/s]

1.430615e-02 36.7%┣██████████▋                  ┫ 734/2.0k [00:38<01:06, 19it/s]

1.349637e-02 37.0%┣██████████▊                  ┫ 741/2.0k [00:38<01:05, 19it/s]

1.299996e-02 37.3%┣██████████▉                  ┫ 747/2.0k [00:38<01:04, 19it/s]

1.250176e-02 37.6%┣███████████                  ┫ 753/2.0k [00:38<01:04, 20it/s]

1.186662e-02 37.9%┣███████████                  ┫ 759/2.0k [00:38<01:03, 20it/s]

1.143545e-02 38.2%┣███████████                  ┫ 765/2.0k [00:38<01:02, 20it/s]

1.095586e-02 38.6%┣███████████▏                 ┫ 772/2.0k [00:39<01:01, 20it/s]

1.059362e-02 38.9%┣███████████▎                 ┫ 778/2.0k [00:39<01:01, 20it/s]

1.019434e-02 39.2%┣███████████▍                 ┫ 784/2.0k [00:39<01:00, 20it/s]

9.776690e-03 39.5%┣███████████▌                 ┫ 791/2.0k [00:39<00:59, 20it/s]

9.445892e-03 39.8%┣███████████▌                 ┫ 797/2.0k [00:39<00:59, 21it/s]

9.146112e-03 40.2%┣███████████▋                 ┫ 804/2.0k [00:39<00:58, 21it/s]

8.858313e-03 40.5%┣███████████▊                 ┫ 810/2.0k [00:39<00:57, 21it/s]

8.521814e-03 40.8%┣███████████▉                 ┫ 816/2.0k [00:39<00:57, 21it/s]

8.167949e-03 41.1%┣████████████                 ┫ 822/2.0k [00:39<00:56, 21it/s]

7.952144e-03 41.4%┣████████████                 ┫ 828/2.0k [00:39<00:55, 21it/s]

7.675209e-03 41.6%┣████████████                 ┫ 832/2.0k [00:39<00:55, 21it/s]

7.294012e-03 41.9%┣████████████▏                ┫ 839/2.0k [00:39<00:54, 21it/s]

7.024751e-03 42.2%┣████████████▎                ┫ 845/2.0k [00:39<00:54, 22it/s]

6.687432e-03 42.5%┣████████████▍                ┫ 851/2.0k [00:39<00:53, 22it/s]

6.452674e-03 42.8%┣████████████▍                ┫ 857/2.0k [00:39<00:53, 22it/s]

6.354092e-03 43.1%┣████████████▌                ┫ 862/2.0k [00:39<00:52, 22it/s]

6.084379e-03 43.4%┣████████████▋                ┫ 869/2.0k [00:39<00:51, 22it/s]

5.931764e-03 43.8%┣████████████▊                ┫ 877/2.0k [00:39<00:51, 22it/s]

5.713330e-03 44.2%┣████████████▉                ┫ 885/2.0k [00:40<00:50, 22it/s]

5.576138e-03 44.6%┣█████████████                ┫ 893/2.0k [00:40<00:49, 23it/s]

5.560699e-03 44.7%┣█████████████                ┫ 895/2.0k [00:40<00:49, 23it/s]

5.432058e-03 45.1%┣█████████████                ┫ 903/2.0k [00:40<00:48, 23it/s]

5.305968e-03 45.5%┣█████████████▏               ┫ 910/2.0k [00:40<00:48, 23it/s]

5.170576e-03 45.8%┣█████████████▎               ┫ 916/2.0k [00:40<00:47, 23it/s]

4.993048e-03 46.1%┣█████████████▍               ┫ 922/2.0k [00:40<00:47, 23it/s]

4.863774e-03 46.4%┣█████████████▌               ┫ 929/2.0k [00:40<00:46, 23it/s]

4.806003e-03 46.7%┣█████████████▌               ┫ 934/2.0k [00:40<00:46, 23it/s]

4.632655e-03 47.1%┣█████████████▋               ┫ 942/2.0k [00:40<00:45, 24it/s]

4.535106e-03 47.4%┣█████████████▊               ┫ 949/2.0k [00:40<00:44, 24it/s]

4.464539e-03 47.8%┣█████████████▉               ┫ 956/2.0k [00:40<00:44, 24it/s]

4.373546e-03 48.1%┣██████████████               ┫ 963/2.0k [00:40<00:43, 24it/s]

4.308327e-03 48.3%┣██████████████               ┫ 966/2.0k [00:40<00:43, 24it/s]

4.212476e-03 48.7%┣██████████████▏              ┫ 974/2.0k [00:40<00:43, 24it/s]

4.136007e-03 49.0%┣██████████████▏              ┫ 981/2.0k [00:40<00:42, 24it/s]

4.033967e-03 49.4%┣██████████████▎              ┫ 988/2.0k [00:40<00:41, 24it/s]

3.905576e-03 49.7%┣██████████████▍              ┫ 995/2.0k [00:40<00:41, 25it/s]

3.806022e-03 50.1%┣██████████████              ┫ 1.0k/2.0k [00:41<00:40, 25it/s]

3.714740e-03 50.5%┣██████████████▏             ┫ 1.0k/2.0k [00:41<00:40, 25it/s]

3.606534e-03 50.9%┣██████████████▎             ┫ 1.0k/2.0k [00:41<00:39, 25it/s]

3.520669e-03 51.4%┣██████████████▍             ┫ 1.0k/2.0k [00:41<00:39, 25it/s]

3.463226e-03 51.8%┣██████████████▌             ┫ 1.0k/2.0k [00:41<00:38, 25it/s]

3.413995e-03 52.0%┣██████████████▌             ┫ 1.0k/2.0k [00:41<00:38, 25it/s]

3.330557e-03 52.4%┣██████████████▊             ┫ 1.0k/2.0k [00:41<00:37, 26it/s]

3.201924e-03 52.8%┣██████████████▉             ┫ 1.1k/2.0k [00:41<00:37, 26it/s]

3.082455e-03 53.3%┣███████████████             ┫ 1.1k/2.0k [00:41<00:36, 26it/s]

2.986226e-03 53.7%┣███████████████             ┫ 1.1k/2.0k [00:41<00:35, 26it/s]

2.933136e-03 54.0%┣███████████████▏            ┫ 1.1k/2.0k [00:41<00:35, 26it/s]

2.847374e-03 54.5%┣███████████████▎            ┫ 1.1k/2.0k [00:41<00:34, 26it/s]

2.758699e-03 54.9%┣███████████████▍            ┫ 1.1k/2.0k [00:41<00:34, 27it/s]

2.635295e-03 55.4%┣███████████████▌            ┫ 1.1k/2.0k [00:41<00:33, 27it/s]

2.558047e-03 55.7%┣███████████████▋            ┫ 1.1k/2.0k [00:41<00:33, 27it/s]

2.439520e-03 56.1%┣███████████████▊            ┫ 1.1k/2.0k [00:41<00:32, 27it/s]

2.367009e-03 56.6%┣███████████████▉            ┫ 1.1k/2.0k [00:41<00:32, 27it/s]

2.256406e-03 57.1%┣████████████████            ┫ 1.1k/2.0k [00:41<00:31, 28it/s]

2.194602e-03 57.5%┣████████████████            ┫ 1.2k/2.0k [00:42<00:31, 28it/s]

2.112240e-03 58.0%┣████████████████▎           ┫ 1.2k/2.0k [00:42<00:30, 28it/s]

2.014785e-03 58.4%┣████████████████▍           ┫ 1.2k/2.0k [00:42<00:30, 28it/s]

1.939773e-03 58.9%┣████████████████▌           ┫ 1.2k/2.0k [00:42<00:29, 28it/s]

1.895068e-03 59.4%┣████████████████▋           ┫ 1.2k/2.0k [00:42<00:29, 28it/s]

1.878257e-03 59.7%┣████████████████▊           ┫ 1.2k/2.0k [00:42<00:28, 29it/s]

1.835011e-03 60.0%┣████████████████▉           ┫ 1.2k/2.0k [00:42<00:28, 29it/s]

1.791688e-03 60.4%┣█████████████████           ┫ 1.2k/2.0k [00:42<00:28, 29it/s]

1.756118e-03 60.7%┣█████████████████           ┫ 1.2k/2.0k [00:42<00:27, 29it/s]

1.691631e-03 61.1%┣█████████████████           ┫ 1.2k/2.0k [00:42<00:27, 29it/s]

1.661264e-03 61.4%┣█████████████████▏          ┫ 1.2k/2.0k [00:42<00:26, 29it/s]

1.613962e-03 61.8%┣█████████████████▎          ┫ 1.2k/2.0k [00:42<00:26, 29it/s]

1.569850e-03 62.2%┣█████████████████▍          ┫ 1.2k/2.0k [00:42<00:26, 30it/s]

1.542027e-03 62.6%┣█████████████████▌          ┫ 1.3k/2.0k [00:42<00:25, 30it/s]

1.512801e-03 62.9%┣█████████████████▋          ┫ 1.3k/2.0k [00:42<00:25, 30it/s]

1.475900e-03 63.3%┣█████████████████▊          ┫ 1.3k/2.0k [00:42<00:25, 30it/s]

1.452574e-03 63.5%┣█████████████████▉          ┫ 1.3k/2.0k [00:42<00:24, 30it/s]

1.412974e-03 64.0%┣██████████████████          ┫ 1.3k/2.0k [00:42<00:24, 30it/s]

1.387896e-03 64.3%┣██████████████████          ┫ 1.3k/2.0k [00:42<00:24, 30it/s]

1.367458e-03 64.7%┣██████████████████          ┫ 1.3k/2.0k [00:43<00:23, 30it/s]

1.334319e-03 65.0%┣██████████████████▏         ┫ 1.3k/2.0k [00:43<00:23, 31it/s]

1.311897e-03 65.3%┣██████████████████▎         ┫ 1.3k/2.0k [00:43<00:23, 31it/s]

1.284367e-03 65.7%┣██████████████████▍         ┫ 1.3k/2.0k [00:43<00:22, 31it/s]

1.257259e-03 66.2%┣██████████████████▌         ┫ 1.3k/2.0k [00:43<00:22, 31it/s]

1.203717e-03 66.7%┣██████████████████▊         ┫ 1.3k/2.0k [00:43<00:21, 31it/s]

1.159889e-03 67.1%┣██████████████████▉         ┫ 1.3k/2.0k [00:43<00:21, 31it/s]

1.126083e-03 67.4%┣██████████████████▉         ┫ 1.3k/2.0k [00:43<00:21, 31it/s]

1.089076e-03 67.8%┣███████████████████         ┫ 1.4k/2.0k [00:43<00:20, 32it/s]

1.049461e-03 68.3%┣███████████████████▏        ┫ 1.4k/2.0k [00:43<00:20, 32it/s]

9.849604e-04 68.7%┣███████████████████▎        ┫ 1.4k/2.0k [00:43<00:20, 32it/s]

9.417840e-04 69.1%┣███████████████████▍        ┫ 1.4k/2.0k [00:43<00:19, 32it/s]

8.769238e-04 69.5%┣███████████████████▌        ┫ 1.4k/2.0k [00:43<00:19, 32it/s]

8.177395e-04 70.0%┣███████████████████▋        ┫ 1.4k/2.0k [00:43<00:19, 32it/s]

7.719005e-04 70.5%┣███████████████████▊        ┫ 1.4k/2.0k [00:43<00:18, 33it/s]

7.154437e-04 71.0%┣███████████████████▉        ┫ 1.4k/2.0k [00:43<00:18, 33it/s]

6.787786e-04 71.3%┣████████████████████        ┫ 1.4k/2.0k [00:43<00:17, 33it/s]

6.225880e-04 71.8%┣████████████████████        ┫ 1.4k/2.0k [00:43<00:17, 33it/s]

5.769212e-04 72.3%┣████████████████████▎       ┫ 1.4k/2.0k [00:44<00:17, 33it/s]

5.344773e-04 72.8%┣████████████████████▍       ┫ 1.5k/2.0k [00:44<00:16, 33it/s]

5.097298e-04 73.2%┣████████████████████▌       ┫ 1.5k/2.0k [00:44<00:16, 34it/s]

4.930512e-04 73.6%┣████████████████████▋       ┫ 1.5k/2.0k [00:44<00:16, 34it/s]

4.769815e-04 74.1%┣████████████████████▊       ┫ 1.5k/2.0k [00:44<00:15, 34it/s]

4.625357e-04 74.6%┣████████████████████▉       ┫ 1.5k/2.0k [00:44<00:15, 34it/s]

4.398070e-04 75.0%┣█████████████████████       ┫ 1.5k/2.0k [00:44<00:15, 34it/s]

4.155922e-04 75.4%┣█████████████████████       ┫ 1.5k/2.0k [00:44<00:14, 34it/s]

3.859251e-04 75.9%┣█████████████████████▎      ┫ 1.5k/2.0k [00:44<00:14, 35it/s]

3.631654e-04 76.4%┣█████████████████████▍      ┫ 1.5k/2.0k [00:44<00:14, 35it/s]

3.543712e-04 76.9%┣█████████████████████▌      ┫ 1.5k/2.0k [00:44<00:13, 35it/s]

3.438153e-04 77.2%┣█████████████████████▋      ┫ 1.5k/2.0k [00:44<00:13, 35it/s]

3.246423e-04 77.7%┣█████████████████████▊      ┫ 1.6k/2.0k [00:44<00:13, 35it/s]

3.087461e-04 78.2%┣█████████████████████▉      ┫ 1.6k/2.0k [00:44<00:12, 35it/s]

2.925008e-04 78.7%┣██████████████████████      ┫ 1.6k/2.0k [00:44<00:12, 36it/s]

2.834174e-04 79.0%┣██████████████████████      ┫ 1.6k/2.0k [00:44<00:12, 36it/s]

2.746756e-04 79.4%┣██████████████████████▎     ┫ 1.6k/2.0k [00:44<00:12, 36it/s]

2.673456e-04 79.9%┣██████████████████████▍     ┫ 1.6k/2.0k [00:44<00:11, 36it/s]

2.619896e-04 80.4%┣██████████████████████▌     ┫ 1.6k/2.0k [00:44<00:11, 36it/s]

2.566898e-04 80.8%┣██████████████████████▋     ┫ 1.6k/2.0k [00:45<00:11, 36it/s]

2.506119e-04 81.3%┣██████████████████████▊     ┫ 1.6k/2.0k [00:45<00:10, 36it/s]

2.434811e-04 81.8%┣███████████████████████     ┫ 1.6k/2.0k [00:45<00:10, 37it/s]

2.345284e-04 82.3%┣███████████████████████     ┫ 1.6k/2.0k [00:45<00:10, 37it/s]

2.279708e-04 82.7%┣███████████████████████▏    ┫ 1.7k/2.0k [00:45<00:09, 37it/s]

2.252035e-04 83.0%┣███████████████████████▎    ┫ 1.7k/2.0k [00:45<00:09, 37it/s]

2.194848e-04 83.5%┣███████████████████████▍    ┫ 1.7k/2.0k [00:45<00:09, 37it/s]

2.118448e-04 84.0%┣███████████████████████▌    ┫ 1.7k/2.0k [00:45<00:09, 37it/s]

2.030400e-04 84.5%┣███████████████████████▋    ┫ 1.7k/2.0k [00:45<00:08, 38it/s]

2.009093e-04 84.7%┣███████████████████████▊    ┫ 1.7k/2.0k [00:45<00:08, 38it/s]

1.981349e-04 85.2%┣███████████████████████▉    ┫ 1.7k/2.0k [00:45<00:08, 38it/s]

1.947511e-04 85.7%┣████████████████████████    ┫ 1.7k/2.0k [00:45<00:08, 38it/s]

1.910733e-04 86.2%┣████████████████████████▏   ┫ 1.7k/2.0k [00:45<00:07, 38it/s]

1.877767e-04 86.6%┣████████████████████████▎   ┫ 1.7k/2.0k [00:45<00:07, 38it/s]

1.847318e-04 86.9%┣████████████████████████▍   ┫ 1.7k/2.0k [00:45<00:07, 38it/s]

1.796687e-04 87.4%┣████████████████████████▌   ┫ 1.7k/2.0k [00:45<00:07, 39it/s]

1.749959e-04 87.9%┣████████████████████████▋   ┫ 1.8k/2.0k [00:45<00:06, 39it/s]

1.715932e-04 88.4%┣████████████████████████▊   ┫ 1.8k/2.0k [00:45<00:06, 39it/s]

1.667290e-04 88.8%┣████████████████████████▉   ┫ 1.8k/2.0k [00:46<00:06, 39it/s]

1.607875e-04 89.3%┣█████████████████████████   ┫ 1.8k/2.0k [00:46<00:05, 39it/s]

1.573171e-04 89.8%┣█████████████████████████▏  ┫ 1.8k/2.0k [00:46<00:05, 39it/s]

1.522889e-04 90.3%┣█████████████████████████▎  ┫ 1.8k/2.0k [00:46<00:05, 40it/s]

1.495920e-04 90.6%┣█████████████████████████▍  ┫ 1.8k/2.0k [00:46<00:05, 40it/s]

1.471885e-04 91.0%┣█████████████████████████▌  ┫ 1.8k/2.0k [00:46<00:05, 40it/s]

1.451792e-04 91.5%┣█████████████████████████▋  ┫ 1.8k/2.0k [00:46<00:04, 40it/s]

1.434688e-04 92.0%┣█████████████████████████▊  ┫ 1.8k/2.0k [00:46<00:04, 40it/s]

1.407036e-04 92.5%┣██████████████████████████  ┫ 1.9k/2.0k [00:46<00:04, 40it/s]

1.370180e-04 93.0%┣██████████████████████████  ┫ 1.9k/2.0k [00:46<00:03, 40it/s]

1.345419e-04 93.4%┣██████████████████████████▏ ┫ 1.9k/2.0k [00:46<00:03, 41it/s]

1.332631e-04 93.9%┣██████████████████████████▎ ┫ 1.9k/2.0k [00:46<00:03, 41it/s]

1.320952e-04 94.4%┣██████████████████████████▍ ┫ 1.9k/2.0k [00:46<00:03, 41it/s]

1.309526e-04 94.7%┣██████████████████████████▌ ┫ 1.9k/2.0k [00:46<00:03, 41it/s]

1.284015e-04 95.1%┣██████████████████████████▋ ┫ 1.9k/2.0k [00:46<00:02, 41it/s]

1.275730e-04 95.6%┣██████████████████████████▊ ┫ 1.9k/2.0k [00:46<00:02, 41it/s]

1.264767e-04 96.1%┣███████████████████████████ ┫ 1.9k/2.0k [00:46<00:02, 41it/s]

1.257943e-04 96.4%┣███████████████████████████ ┫ 1.9k/2.0k [00:46<00:02, 42it/s]

1.245031e-04 96.9%┣███████████████████████████▏┫ 1.9k/2.0k [00:46<00:02, 42it/s]

1.231696e-04 97.3%┣███████████████████████████▎┫ 1.9k/2.0k [00:47<00:01, 42it/s]

1.228674e-04 97.7%┣███████████████████████████▍┫ 2.0k/2.0k [00:47<00:01, 42it/s]

1.219406e-04 98.1%┣███████████████████████████▌┫ 2.0k/2.0k [00:47<00:01, 42it/s]

1.209161e-04 98.4%┣███████████████████████████▌┫ 2.0k/2.0k [00:47<00:01, 42it/s]

1.205982e-04 98.8%┣███████████████████████████▋┫ 2.0k/2.0k [00:47<00:01, 42it/s]

1.200536e-04 99.2%┣███████████████████████████▊┫ 2.0k/2.0k [00:47<00:00, 42it/s]

1.181568e-04 99.7%┣████████████████████████████┫ 2.0k/2.0k [00:47<00:00, 43it/s]

1.178800e-04 100.0%┣███████████████████████████┫ 2.0k/2.0k [00:47<00:00, 43it/s]

Compute the relative L2 error

using Integrals

u_analytical(x,p) = sin.(2 * pi .* x) + 0.1 * sin.(50 * pi .* x)
error(x,p) = abs2.(vec(phi([x;;],res.u)) .- u_analytical(x,p))

relative_L2_error = solve(IntegralProblem(error,0,1),HCubatureJL(),reltol=1e-3,abstol=1e-3) ./ solve(IntegralProblem((x,p) -> abs2.(u_analytical(x,p)),0, 1),HCubatureJL(),reltol=1e-3,abstol=1e-3)
1-element Vector{Float64}:
 2.724438778516777e-9