11:17:38 [INFO] Required matrix memory: 0.69 MB
11:17:38 [INFO] mesh is not defined, creating evenly distributed mesh of length 51
11:17:38 [INFO] system discretized
11:17:38 [INFO] Calculating jacobian
11:17:38 [INFO] Jacobian calculation time:
11:17:38 [INFO] Elapsed 28 ms
11:17:38 [INFO] (" ms ", 28.0)
11:17:38 [INFO] Bandwidth calculated:
11:17:38 [INFO] (kl, ku) = Some((8, 3))
11:17:38 [INFO] Jacobian lambdified
11:17:38 [INFO] Residuals vector lambdified
11:17:38 [INFO] 
 
 ____________END OF GENERATE BVP ________________________________________________________________
11:17:38 [INFO] 
 
 solving system of equations with Newton-Raphson method! 
 

11:17:38 [INFO] (6, 50)
11:17:38 [INFO] norm of the initial residual = 1.053874698251816
11:17:38 [INFO] 
 number of iterations with old jac 0 is higher then threshold 3
11:17:38 [INFO] 
 
 JACOBIAN (RE)CALCULATED! 
 

11:17:38 [INFO] jacobian recalculation time: 
11:17:38 [INFO] Elapsed 1 ms
11:17:38 [INFO] 
 
 residual norm = 1.053874698251816 
11:17:38 [INFO] 
 
 fboundary  = 1
11:17:38 [INFO] 
 
 Damping coefficient = 1
11:17:38 [INFO] 
 
 residual norm = 0.12362378412905398 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 27.638180592981612
11:17:38 [INFO] Elapsed 15 ms
11:17:38 [INFO] damping coefficient decreased to 0.7071067811865475
11:17:38 [INFO] 
 
 Damping coefficient = 0.7071067811865475
11:17:38 [INFO] 
 
 residual norm = 0.3139806382940293 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 14.4860302601236
11:17:38 [INFO] Elapsed 22 ms
11:17:38 [INFO] damping coefficient decreased to 0.24999999999999997
11:17:38 [INFO] 
 
 damped_step number 2 
11:17:38 [INFO] 
 
 Damping coefficient = 0.24999999999999997
11:17:38 [INFO] 
 
 residual norm = 0.7904267435012543 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 7.467004373543234
11:17:38 [INFO] Elapsed 31 ms
11:17:38 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:38 [INFO] 
 
  step norm =  7.467004373543234, weight norm = 7.467004373543234, convergence condition = 0.000003227550885450139
11:17:38 [INFO] 
 
 end of iteration 1 with jac age 1 
 

11:17:38 [INFO] 
 
 residual norm = 0.7904267435012543 
11:17:38 [INFO] 
 
 fboundary  = 1
11:17:38 [INFO] 
 
 Damping coefficient = 1
11:17:38 [INFO] 
 
 residual norm = 0.1212836614551693 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 26.88023256193245
11:17:38 [INFO] Elapsed 14 ms
11:17:38 [INFO] damping coefficient decreased to 0.7071067811865475
11:17:38 [INFO] 
 
 Damping coefficient = 0.7071067811865475
11:17:38 [INFO] 
 
 residual norm = 0.24179089551110527 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 16.683515918205064
11:17:38 [INFO] Elapsed 22 ms
11:17:38 [INFO] damping coefficient decreased to 0.24999999999999997
11:17:38 [INFO] 
 
 damped_step number 2 
11:17:38 [INFO] 
 
 Damping coefficient = 0.24999999999999997
11:17:38 [INFO] 
 
 residual norm = 0.5931011198463517 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 8.120439831696817
11:17:38 [INFO] Elapsed 30 ms
11:17:38 [INFO] damping coefficient decreased to 0.04419417382415921
11:17:38 [INFO] 
 
 damped_step number 3 
11:17:38 [INFO] 
 
 Damping coefficient = 0.04419417382415921
11:17:38 [INFO] 
 
 residual norm = 0.755510135460278 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 7.419349564188868
11:17:38 [INFO] Elapsed 38 ms
11:17:38 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:38 [INFO] 
 
  step norm =  7.419349564188868, weight norm = 7.419349564188868, convergence condition = 0.000003716399837335049
11:17:38 [INFO] 
 
 end of iteration 2 with jac age 2 
 

11:17:38 [INFO] 
 
 residual norm = 0.755510135460278 
11:17:38 [INFO] 
 
 fboundary  = 1
11:17:38 [INFO] 
 
 Damping coefficient = 1
11:17:38 [INFO] 
 
 residual norm = 0.11942924648441665 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 26.345768922734443
11:17:38 [INFO] Elapsed 14 ms
11:17:38 [INFO] damping coefficient decreased to 0.7071067811865475
11:17:38 [INFO] 
 
 Damping coefficient = 0.7071067811865475
11:17:38 [INFO] 
 
 residual norm = 0.23227485932326217 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 16.76225421696539
11:17:38 [INFO] Elapsed 23 ms
11:17:38 [INFO] damping coefficient decreased to 0.24999999999999997
11:17:38 [INFO] 
 
 damped_step number 2 
11:17:38 [INFO] 
 
 Damping coefficient = 0.24999999999999997
11:17:38 [INFO] 
 
 residual norm = 0.5669910698055807 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 8.330239429041328
11:17:38 [INFO] Elapsed 32 ms
11:17:38 [INFO] damping coefficient decreased to 0.04419417382415921
11:17:38 [INFO] 
 
 damped_step number 3 
11:17:38 [INFO] 
 
 Damping coefficient = 0.04419417382415921
11:17:38 [INFO] 
 
 residual norm = 0.7221438712650066 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 7.419565836572473
11:17:38 [INFO] Elapsed 40 ms
11:17:38 [INFO] damping coefficient decreased to 0.0039062499999999987
11:17:38 [INFO] 
 
 damped_step number 4 
11:17:38 [INFO] 
 
 Damping coefficient = 0.0039062499999999987
11:17:38 [INFO] 
 
 residual norm = 0.7525606001666887 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 7.417072493594616
11:17:38 [INFO] Elapsed 48 ms
11:17:38 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:38 [INFO] 
 
  step norm =  7.417072493594616, weight norm = 7.417072493594616, convergence condition = 0.000003759440863932329
11:17:38 [INFO] 
 
 end of iteration 3 with jac age 3 
 

11:17:38 [INFO] 
 
 residual norm = 0.7525606001666887 
11:17:38 [INFO] 
 
 fboundary  = 1
11:17:38 [INFO] 
 
 Damping coefficient = 1
11:17:38 [INFO] 
 
 residual norm = 0.11923155133890083 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 26.289250944885012
11:17:38 [INFO] Elapsed 16 ms
11:17:38 [INFO] damping coefficient decreased to 0.7071067811865475
11:17:38 [INFO] 
 
 Damping coefficient = 0.7071067811865475
11:17:38 [INFO] 
 
 residual norm = 0.23146690231791797 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 16.762658283727536
11:17:38 [INFO] Elapsed 26 ms
11:17:38 [INFO] damping coefficient decreased to 0.24999999999999997
11:17:38 [INFO] 
 
 damped_step number 2 
11:17:38 [INFO] 
 
 Damping coefficient = 0.24999999999999997
11:17:38 [INFO] 
 
 residual norm = 0.5647859912434534 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 8.348376837354072
11:17:38 [INFO] Elapsed 33 ms
11:17:38 [INFO] damping coefficient decreased to 0.04419417382415921
11:17:38 [INFO] 
 
 damped_step number 3 
11:17:38 [INFO] 
 
 Damping coefficient = 0.04419417382415921
11:17:38 [INFO] 
 
 residual norm = 0.7193253869220041 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 7.42129197328372
11:17:38 [INFO] Elapsed 42 ms
11:17:38 [INFO] damping coefficient decreased to 0.0039062499999999987
11:17:38 [INFO] 
 
 damped_step number 4 
11:17:38 [INFO] 
 
 Damping coefficient = 0.0039062499999999987
11:17:38 [INFO] 
 
 residual norm = 0.7496226413683598 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 7.4151471742269734
11:17:38 [INFO] Elapsed 50 ms
11:17:38 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:38 [INFO] 
 
  step norm =  7.4151471742269734, weight norm = 7.4151471742269734, convergence condition = 0.0000038024746886135056
11:17:38 [INFO] 
 
 end of iteration 4 with jac age 4 
 

11:17:38 [INFO] 
 number of iterations with old jac 4 is higher then threshold 3
11:17:38 [INFO] 
 
 JACOBIAN (RE)CALCULATED! 
 

11:17:38 [INFO] jacobian recalculation time: 
11:17:38 [INFO] Elapsed 2 ms
11:17:38 [INFO] 
 
 residual norm = 0.7496226413683598 
11:17:38 [INFO] 
 
 fboundary  = 1
11:17:38 [INFO] 
 
 Damping coefficient = 1
11:17:38 [INFO] 
 
 residual norm = 0.2339960015397812 
11:17:38 [INFO] 
 
 L2 norm of undamped step = 58.42984685467341
11:17:38 [INFO] Elapsed 15 ms
11:17:38 [INFO] damping coefficient decreased to 0.7071067811865475
11:17:38 [INFO] 
 
 Damping coefficient = 0.7071067811865475
11:17:38 [INFO] 
 
 residual norm = 0.24936008304583648 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 24.659988095107472
11:17:39 [INFO] Elapsed 22 ms
11:17:39 [INFO] damping coefficient decreased to 0.24999999999999997
11:17:39 [INFO] 
 
 damped_step number 2 
11:17:39 [INFO] 
 
 Damping coefficient = 0.24999999999999997
11:17:39 [INFO] 
 
 residual norm = 0.562279068916019 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 10.686118224250317
11:17:39 [INFO] Elapsed 31 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:39 [INFO] 
 
  step norm =  10.686118224250317, weight norm = 10.686118224250317, convergence condition = 0.000009839069448880256
11:17:39 [INFO] 
 
 end of iteration 5 with jac age 1 
 

11:17:39 [INFO] 
 
 residual norm = 0.562279068916019 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 0.10503769818422413 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 22.486632412870694
11:17:39 [INFO] Elapsed 14 ms
11:17:39 [INFO] damping coefficient decreased to 0.7071067811865475
11:17:39 [INFO] 
 
 Damping coefficient = 0.7071067811865475
11:17:39 [INFO] 
 
 residual norm = 0.17788768500148705 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 11.10002961106509
11:17:39 [INFO] Elapsed 22 ms
11:17:39 [INFO] damping coefficient decreased to 0.24999999999999997
11:17:39 [INFO] 
 
 damped_step number 2 
11:17:39 [INFO] 
 
 Damping coefficient = 0.24999999999999997
11:17:39 [INFO] 
 
 residual norm = 0.42233105274026833 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 4.883357341449434
11:17:39 [INFO] Elapsed 31 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:39 [INFO] 
 
  step norm =  4.883357341449434, weight norm = 4.883357341449434, convergence condition = 0.00001362435311162393
11:17:39 [INFO] 
 
 end of iteration 6 with jac age 2 
 

11:17:39 [INFO] 
 
 residual norm = 0.42233105274026833 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 0.04416120250147185 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 9.558076870864493
11:17:39 [INFO] Elapsed 84 ms
11:17:39 [INFO] damping coefficient decreased to 0.7071067811865475
11:17:39 [INFO] 
 
 Damping coefficient = 0.7071067811865475
11:17:39 [INFO] 
 
 residual norm = 0.1265073792211707 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 6.8791244185728315
11:17:39 [INFO] Elapsed 102 ms
11:17:39 [INFO] damping coefficient decreased to 0.24999999999999997
11:17:39 [INFO] 
 
 damped_step number 2 
11:17:39 [INFO] 
 
 Damping coefficient = 0.24999999999999997
11:17:39 [INFO] 
 
 residual norm = 0.31683576418894366 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 4.846368550722914
11:17:39 [INFO] Elapsed 110 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:39 [INFO] 
 
  step norm =  4.846368550722914, weight norm = 4.846368550722914, convergence condition = 0.00001522885085156274
11:17:39 [INFO] 
 
 end of iteration 7 with jac age 3 
 

11:17:39 [INFO] 
 
 residual norm = 0.31683576418894366 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 0.030233010680436857 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 1.2785666444290296
11:17:39 [INFO] Elapsed 17 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:39 [INFO] 
 
  step norm =  1.2785666444290296, weight norm = 1.2785666444290296, convergence condition = 0.000021795138806164317
11:17:39 [INFO] 
 
 end of iteration 8 with jac age 4 
 

11:17:39 [INFO] 
 number of iterations with old jac 4 is higher then threshold 3
11:17:39 [INFO] 
 
 JACOBIAN (RE)CALCULATED! 
 

11:17:39 [INFO] jacobian recalculation time: 
11:17:39 [INFO] Elapsed 3 ms
11:17:39 [INFO] 
 
 residual norm = 0.030233010680436857 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 0.0008956682166862503 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 0.025559064652228864
11:17:39 [INFO] Elapsed 15 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:39 [INFO] 
 
  step norm =  0.025559064652228864, weight norm = 0.025559064652228864, convergence condition = 0.00002170006315321213
11:17:39 [INFO] 
 
 end of iteration 9 with jac age 1 
 

11:17:39 [INFO] 
 
 residual norm = 0.0008956682166862503 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 6.621107936382199e-5 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 0.0018053357074719468
11:17:39 [INFO] Elapsed 13 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:39 [INFO] 
 
  step norm =  0.0018053357074719468, weight norm = 0.0018053357074719468, convergence condition = 0.0000216830427703902
11:17:39 [INFO] 
 
 end of iteration 10 with jac age 2 
 

11:17:39 [INFO] 
 
 residual norm = 6.621107936382199e-5 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 4.726669240769536e-6 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 0.00013264594346661966
11:17:39 [INFO] Elapsed 14 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has not converged yet)
11:17:39 [INFO] 
 
  step norm =  0.00013264594346661966, weight norm = 0.00013264594346661966, convergence condition = 0.000021683768913490485
11:17:39 [INFO] 
 
 end of iteration 11 with jac age 3 
 

11:17:39 [INFO] 
 
 residual norm = 4.726669240769536e-6 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 3.4364979564713294e-7 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 0.000009552112397943027
11:17:39 [INFO] Elapsed 13 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has converged)
11:17:39 [INFO] 
 
 step norm =  0.000009552112397943027, weight norm = 0.000009552112397943027, convergence condition = 0.000021683700975247254
11:17:39 [INFO] 
 
 Solution has converged, breaking the loop!
11:17:39 [INFO] residual norm of the solution = 0.00000034364979564713294
11:17:39 [INFO] 
 
 solution found for the current grid 300
11:17:39 [INFO] solving with new grid!
11:17:39 [INFO] ================GRID REFINEMENT===================
11:17:39 [INFO] Constructing full solution
11:17:39 [INFO] old mesh len 51
11:17:39 [INFO] shape of solution 6, 51
11:17:39 [INFO] y= 0.3333333333333333, 2.3470503235374265
11:17:39 [INFO] delta 0.2013716990204093 (2.3470503235374265, 0.3333333333333333) for component j: 0
11:17:39 [INFO] gamma 0.6692802691756459 for component j: 0
11:17:39 [INFO] 
 
 for row 0 found intervals to be refined: [] of length 0 
 

11:17:39 [INFO] y= 0.4695455924857138, 1e-10
11:17:39 [INFO] delta 0.046954559238571385 (0.4695455924857138, 0.0000000001) for component j: 1
11:17:39 [INFO] gamma 0.1048081614517704 for component j: 1
11:17:39 [INFO] tau 0.010944161036031519 (0.4586014314496823, 0.4695455924857138)> delta 0.046954559238571385 for i 1
11:17:39 [INFO] eta 0.1373826022154967  (-0.6845906540170726, -0.547208051801576)> gamma 0.1048081614517704 for i 1
11:17:39 [INFO]  conditions vaiolation at index 1, => N = 1
11:17:39 [INFO] tau 0.013691813080341453 (0.44490961836934084, 0.4586014314496823)> delta 0.046954559238571385 for i 2
11:17:39 [INFO] eta 0.12352098013892765  (-0.8081116341560003, -0.6845906540170726)> gamma 0.1048081614517704 for i 2
11:17:39 [INFO]  conditions vaiolation at index 2, => N = 1
11:17:39 [INFO] 
 
 for row 1 found intervals to be refined: [1, 0] of length 2 
 

11:17:39 [INFO] y= 0.999, 0.014004576928736778
11:17:39 [INFO] delta 0.09849954230712632 (0.999, 0.014004576928736778) for component j: 2
11:17:39 [INFO] gamma 0.3278465372031076 for component j: 2
11:17:39 [INFO] 
 
 for row 2 found intervals to be refined: [1, 0] of length 2 
 

11:17:39 [INFO] y= -0.0009452756755774057, 1e-10
11:17:39 [INFO] delta 0.00009452757755774058 (0.0000000001, -0.0009452756755774057) for component j: 3
11:17:39 [INFO] gamma 0.00021112753812835012 for component j: 3
11:17:39 [INFO] tau 0.0000222801383290385 (-0.0009229955372483672, -0.0009452756755774057)> delta 0.00009452757755774058 for i 1
11:17:39 [INFO] eta 0.00027445795453644975  (0.0013884648709883748, 0.001114006916451925)> gamma 0.00021112753812835012 for i 1
11:17:39 [INFO]  conditions vaiolation at index 1, => N = 1
11:17:39 [INFO] tau 0.000027769297419767496 (-0.0008952262398285997, -0.0009229955372483672)> delta 0.00009452757755774058 for i 2
11:17:39 [INFO] eta 0.0002466093649638136  (0.0016350742359521884, 0.0013884648709883748)> gamma 0.00021112753812835012 for i 2
11:17:39 [INFO]  conditions vaiolation at index 2, => N = 1
11:17:39 [INFO] 
 
 for row 3 found intervals to be refined: [1, 0] of length 2 
 

11:17:39 [INFO] y= 0.001, 0.9871891104669605
11:17:39 [INFO] delta 0.09861891104669605 (0.9871891104669605, 0.001) for component j: 4
11:17:39 [INFO] gamma 0.32824361358807086 for component j: 4
11:17:39 [INFO] 
 
 for row 4 found intervals to be refined: [1, 0] of length 2 
 

11:17:39 [INFO] y= 0.0009464207608832945, 1e-10
11:17:39 [INFO] delta 0.00009464206608832945 (0.0009464207608832945, 0.0000000001) for component j: 5
11:17:39 [INFO] gamma 0.00021138324867667365 for component j: 5
11:17:39 [INFO] tau 0.000022307123264442004 (0.0009241136376188525, 0.0009464207608832945)> delta 0.00009464206608832945 for i 1
11:17:39 [INFO] eta 0.0002747903687472177  (-0.0013901465319693179, -0.0011153561632221002)> gamma 0.00021138324867667365 for i 1
11:17:39 [INFO]  conditions vaiolation at index 1, => N = 1
11:17:39 [INFO] tau 0.000027802930639386357 (0.0008963107069794662, 0.0009241136376188525)> delta 0.00009464206608832945 for i 2
11:17:39 [INFO] eta 0.0002469080499028567  (-0.0016370545818721746, -0.0013901465319693179)> gamma 0.00021138324867667365 for i 2
11:17:39 [INFO]  conditions vaiolation at index 2, => N = 1
11:17:39 [INFO] 
 
 for row 5 found intervals to be refined: [1, 0] of length 2 
 

11:17:39 [INFO] total new points to add: 2
11:17:39 [INFO] 
 
 points added: 0.01 at index 0 
11:17:39 [INFO] 
 
 column added: 
  ┌                        ┐
  │     0.4004112751170067 │
  │     0.4640735119676981 │
  │     0.9661423297610856 │
  │ -0.0009341356064128865 │
  │    0.03389747326817707 │
  │  0.0009352671992510736 │
  └                        ┘

 at index 1 
11:17:39 [INFO] 
 
 points added: 0.03 at index 1 
11:17:39 [INFO] 
 
 column added: 
  ┌                        ┐
  │     0.5330037071077776 │
  │    0.45175552490951154 │
  │     0.9012014442459806 │
  │ -0.0009091108885384835 │
  │    0.09891702684865711 │
  │  0.0009102121722991594 │
  └                        ┘

 at index 3 
11:17:39 [INFO] created new grid of length 53
11:17:39 [INFO] 
 
 new_initial_guess of shape(6, 53)
11:17:39 [INFO] Extracting unknown variables from full solution
11:17:39 [INFO] Sorted biased BC tuples: [(0, 0, 0.3333333333333333), (2, 0, 0.999), (4, 0, 0.001), (313, 1, 1e-10), (315, 1, 1e-10), (317, 1, 1e-10)] and the length of full_vector 318
11:17:39 [INFO] found BC position 0 with value 0.3333333333333333
11:17:39 [INFO] found BC position 2 with value 0.999
11:17:39 [INFO] found BC position 4 with value 0.001
11:17:39 [INFO] found BC position 313 with value 0.0000000001
11:17:39 [INFO] found BC position 315 with value 0.0000000001
11:17:39 [INFO] found BC position 317 with value 0.0000000001
11:17:39 [INFO] ================GRID REFINEMENT ENDED===================
11:17:39 [INFO] 
 
 grid refinement counter = 1 
 

11:17:39 [INFO] 
 
 number of marked intervals is 2, new grid is needed 
 

11:17:39 [INFO] new guess of shape 6 52
11:17:39 [INFO] new mesh length 53
11:17:39 [INFO] Required matrix memory: 0.74 MB
11:17:39 [INFO] mesh with n_steps = 53 is defined directly
11:17:39 [INFO] system discretized
11:17:39 [INFO] Bandwidth provided: Some((8, 3))
11:17:39 [INFO] Calculating jacobian
11:17:39 [INFO] Jacobian calculation time:
11:17:39 [INFO] Elapsed 17 ms
11:17:39 [INFO] (" ms ", 17.0)
11:17:39 [INFO] Jacobian lambdified
11:17:39 [INFO] Residuals vector lambdified
11:17:39 [INFO] 
 
 ____________END OF GENERATE BVP ________________________________________________________________
11:17:39 [INFO] 
 
 solving system of equations with Newton-Raphson method! 
 

11:17:39 [INFO] (6, 52)
11:17:39 [INFO] norm of the initial residual = 0.0017950504778655163
11:17:39 [INFO] 
 number of iterations with old jac 0 is higher then threshold 3
11:17:39 [INFO] 
 
 JACOBIAN (RE)CALCULATED! 
 

11:17:39 [INFO] jacobian recalculation time: 
11:17:39 [INFO] Elapsed 2 ms
11:17:39 [INFO] 
 
 residual norm = 0.0017950504778655163 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 4.217537533699929e-8 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 0.0000008674406450817753
11:17:39 [INFO] Elapsed 16 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has converged)
11:17:39 [INFO] 
 
 step norm =  0.0000008674406450817753, weight norm = 0.0000008674406450817753, convergence condition = 0.0000229000711937569
11:17:39 [INFO] 
 
 Solution has converged, breaking the loop!
11:17:39 [INFO] residual norm of the solution = 0.00000004217537533699929
11:17:39 [INFO] 
 
 solution found for the current grid 312
11:17:39 [INFO] solving with new grid!
11:17:39 [INFO] ================GRID REFINEMENT===================
11:17:39 [INFO] Constructing full solution
11:17:39 [INFO] old mesh len 53
11:17:39 [INFO] shape of solution 6, 53
11:17:39 [INFO] y= 0.3333333333333333, 2.3470443670530887
11:17:39 [INFO] delta 0.20137110337197553 (2.3470443670530887, 0.3333333333333333) for component j: 0
11:17:39 [INFO] gamma 0.6713693541801391 for component j: 0
11:17:39 [INFO] 
 
 for row 0 found intervals to be refined: [] of length 0 
 

11:17:39 [INFO] y= 0.47100892140252737, 1e-10
11:17:39 [INFO] delta 0.04710089213025274 (0.47100892140252737, 0.0000000001) for component j: 1
11:17:39 [INFO] gamma 0.10480102626775943 for component j: 1
11:17:39 [INFO] 
 
 for row 1 found intervals to be refined: [] of length 0 
 

11:17:39 [INFO] y= 0.999, 0.014017591025621688
11:17:39 [INFO] delta 0.09849824089743783 (0.999, 0.014017591025621688) for component j: 2
11:17:39 [INFO] gamma 0.3288636663213531 for component j: 2
11:17:39 [INFO] 
 
 for row 2 found intervals to be refined: [] of length 0 
 

11:17:39 [INFO] y= -0.0009482037743574734, 1e-10
11:17:39 [INFO] delta 0.00009482038743574735 (0.0000000001, -0.0009482037743574734) for component j: 3
11:17:39 [INFO] gamma 0.00021111405044902206 for component j: 3
11:17:39 [INFO] 
 
 for row 3 found intervals to be refined: [] of length 0 
 

11:17:39 [INFO] y= 0.001, 0.9871760806084011
11:17:39 [INFO] delta 0.09861760806084012 (0.9871760806084011, 0.001) for component j: 4
11:17:39 [INFO] gamma 0.3292619746187732 for component j: 4
11:17:39 [INFO] 
 
 for row 4 found intervals to be refined: [] of length 0 
 

11:17:39 [INFO] y= 0.0009493524060777725, 1e-10
11:17:39 [INFO] delta 0.00009493523060777726 (0.0009493524060777725, 0.0000000001) for component j: 5
11:17:39 [INFO] gamma 0.00021136974466152414 for component j: 5
11:17:39 [INFO] 
 
 for row 5 found intervals to be refined: [] of length 0 
 

11:17:39 [INFO] total new points to add: 0
11:17:39 [INFO] created new grid of length 53
11:17:39 [INFO] 
 
 new_initial_guess of shape(6, 53)
11:17:39 [INFO] Extracting unknown variables from full solution
11:17:39 [INFO] Sorted biased BC tuples: [(0, 0, 0.3333333333333333), (2, 0, 0.999), (4, 0, 0.001), (313, 1, 1e-10), (315, 1, 1e-10), (317, 1, 1e-10)] and the length of full_vector 318
11:17:39 [INFO] found BC position 0 with value 0.3333333333333333
11:17:39 [INFO] found BC position 2 with value 0.999
11:17:39 [INFO] found BC position 4 with value 0.001
11:17:39 [INFO] found BC position 313 with value 0.0000000001
11:17:39 [INFO] found BC position 315 with value 0.0000000001
11:17:39 [INFO] found BC position 317 with value 0.0000000001
11:17:39 [INFO] ================GRID REFINEMENT ENDED===================
11:17:39 [INFO] 
 
 grid refinement counter = 2 
 

11:17:39 [INFO] 
 
 number of marked intervals is 0, no new grid is needed 
 

11:17:39 [INFO] no new grid needed - returning to main loop
11:17:39 [INFO] 
 
 end of iteration 1 with jac age 1 
 

11:17:39 [INFO] 
 
 JACOBIAN (RE)CALCULATED! 
 

11:17:39 [INFO] jacobian recalculation time: 
11:17:39 [INFO] Elapsed 1 ms
11:17:39 [INFO] 
 
 residual norm = 4.217537533699929e-8 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 1.9895358264676457e-15 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 0.000000000000043379290257735305
11:17:39 [INFO] Elapsed 16 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has converged)
11:17:39 [INFO] 
 
 step norm =  0.000000000000043379290257735305, weight norm = 0.000000000000043379290257735305, convergence condition = 0.00002290007055957458
11:17:39 [INFO] 
 
 Solution has converged, breaking the loop!
11:17:39 [INFO] residual norm of the solution = 0.0000000000000019895358264676457
11:17:39 [INFO] 
 
 solution found for the current grid 312
11:17:39 [INFO] returning the result
11:17:39 [INFO] 
 
 end of iteration 12 with jac age 1 
 

11:17:39 [INFO] 
 
 JACOBIAN (RE)CALCULATED! 
 

11:17:39 [INFO] jacobian recalculation time: 
11:17:39 [INFO] Elapsed 1 ms
11:17:39 [INFO] 
 
 residual norm = 4.217537533699929e-8 
11:17:39 [INFO] 
 
 fboundary  = 1
11:17:39 [INFO] 
 
 Damping coefficient = 1
11:17:39 [INFO] 
 
 residual norm = 1.9895358264676457e-15 
11:17:39 [INFO] 
 
 L2 norm of undamped step = 0.000000000000043379290257735305
11:17:39 [INFO] Elapsed 15 ms
11:17:39 [INFO] 
 
  Damping coefficient found (solution has converged)
11:17:39 [INFO] 
 
 step norm =  0.000000000000043379290257735305, weight norm = 0.000000000000043379290257735305, convergence condition = 0.00002290007055957458
11:17:39 [INFO] 
 
 Solution has converged, breaking the loop!
11:17:39 [INFO] residual norm of the solution = 0.0000000000000019895358264676457
11:17:39 [INFO] 
 
 solution found for the current grid 312
11:17:39 [INFO] returning the result
11:17:39 [INFO] Elapsed 641 ms
11:17:39 [INFO] Constructing full solution
11:17:39 [INFO] matrix of results has shape (53, 6)
11:17:39 [INFO] length of x mesh : 52
11:17:39 [INFO] number of Ys: 6
11:17:39 [INFO] Matrix memory usage: 0.74 MB
11:17:39 [INFO] 
 
 CALC STATISTICS 
 
 ╭────────────────────────────────────┬───────╮
│ number of grid refinements         │ 2     │
├────────────────────────────────────┼───────┤
│ number of jacobians recalculations │ 6     │
├────────────────────────────────────┼───────┤
│ number of grid points              │ 53    │
├────────────────────────────────────┼───────┤
│ jacobian memory, MB                │ 0     │
├────────────────────────────────────┼───────┤
│ number of jacobian elements        │ 97344 │
├────────────────────────────────────┼───────┤
│ length of y vector                 │ 312   │
├────────────────────────────────────┼───────┤
│ number of iterations               │ 15    │
├────────────────────────────────────┼───────┤
│ time elapsed, s                    │ 0     │
├────────────────────────────────────┼───────┤
│ number of solving linear systems   │ 49    │
╰────────────────────────────────────┴───────╯
11:17:39 [INFO] 
 
 NODES ADDED PER REFINEMENT 
 

11:17:39 [INFO] 
 ╭───┬───╮
│ 0 │ 2 │
├───┼───┤
│ 1 │ 0 │
╰───┴───╯ 

11:17:39 [INFO] Elapsed 1 s
11:17:39 [INFO] Elapsed 14 ms
11:17:39 [INFO] Elapsed 34 ms
11:17:39 [INFO] Elapsed 290 ms
11:17:39 [INFO] Elapsed 599 ms
11:17:39 [INFO] Elapsed 49 ms
11:17:39 [INFO] 
 
 TIMER DATA 
 
 ╭───────────────────────────────┬─────────────╮
│ Jacobian (%,  ms )            │ 1.154, 14   │
├───────────────────────────────┼─────────────┤
│ time elapsed,  s              │ 1           │
├───────────────────────────────┼─────────────┤
│ Linear System (%,  ms )       │ 23.867, 290 │
├───────────────────────────────┼─────────────┤
│ Symbolic Operations (%,  ms ) │ 49.331, 599 │
├───────────────────────────────┼─────────────┤
│ Function (%,  ms )            │ 2.872, 34   │
├───────────────────────────────┼─────────────┤
│ other %                       │ 18.691      │
├───────────────────────────────┼─────────────┤
│ Grid Refinement (%,  ms )     │ 4.084, 49   │
╰───────────────────────────────┴─────────────╯
11:17:39 [INFO]  
 
 Program ended
