Public
Authored by Frederik Schnack

CCG residual version

Initial mesh:

image

Output:

level =  0
nⱼ =  36 which took 0.052092790603637695 s
error =  3.589524310678518 which took 0.020485877990722656 s
tolerance =  0.08460000000000001
mesh refined which took 0.015832185745239258 s
---------- level 0 took 0.20894956588745117 s ------------

level =  1
nⱼ =  48 which took 0.0019307136535644531 s
mⱼ =  5 which took 0.014891862869262695 s
δⱼ =  0.013126074526166941
error =  2.772628774795587 which took 0.009831666946411133 s
mesh refined which took 0.0042133331298828125 s
---------- level 1 took 0.032209157943725586 s ------------

level =  2
nⱼ =  69 which took 0.0019652843475341797 s
mⱼ =  4 which took 0.010799169540405273 s
δⱼ =  0.027156920856268
error =  2.1776651748689506 which took 0.013390779495239258 s
mesh refined which took 0.006398916244506836 s
---------- level 2 took 0.03415656089782715 s ------------

level =  3
nⱼ =  115 which took 0.002709627151489258 s
mⱼ =  3 which took 0.017399311065673828 s
δⱼ =  0.03799331167547584
error =  1.6620959229209324 which took 0.029681921005249023 s
mesh refined which took 0.015787839889526367 s
---------- level 3 took 0.06606364250183105 s ------------

level =  4
nⱼ =  134 which took 0.0035123825073242188 s
mⱼ =  1 which took 0.024500370025634766 s
δⱼ =  0.02662708612021174
error =  1.5104245691979308 which took 0.04893946647644043 s
mesh refined which took 0.05667877197265625 s
---------- level 4 took 0.136016845703125 s ------------

level =  5
nⱼ =  170 which took 0.003367185592651367 s
mⱼ =  1 which took 0.023004770278930664 s
δⱼ =  0.038702546241846415
error =  1.329638816176021 which took 0.044661521911621094 s
mesh refined which took 0.02322697639465332 s
---------- level 5 took 0.09552335739135742 s ------------

level =  6
nⱼ =  218 which took 0.003875255584716797 s
mⱼ =  1 which took 0.039270639419555664 s
δⱼ =  0.03967120202416681
error =  1.1619372485306598 which took 0.05426478385925293 s
mesh refined which took 0.03534579277038574 s
---------- level 6 took 0.13682270050048828 s ------------

level =  7
nⱼ =  248 which took 0.005029916763305664 s
mⱼ =  1 which took 0.03787827491760254 s
δⱼ =  0.02895274864323608
error =  1.0871979522547435 which took 0.07219934463500977 s
mesh refined which took 0.03248715400695801 s
---------- level 7 took 0.15065431594848633 s ------------

level =  8
nⱼ =  294 which took 0.003549814224243164 s
mⱼ =  1 which took 0.036223649978637695 s
δⱼ =  0.0377883777211573
error =  1.0060648280920081 which took 0.07451844215393066 s
mesh refined which took 0.03809976577758789 s
---------- level 8 took 0.15430068969726562 s ------------

level =  9
nⱼ =  453 which took 0.007724761962890625 s
mⱼ =  1 which took 0.05041074752807617 s
δⱼ =  0.03667712716833925
error =  0.8034173790199557 which took 0.1120150089263916 s
mesh refined which took 0.04465770721435547 s
---------- level 9 took 0.21564579010009766 s ------------

level =  10
nⱼ =  542 which took 0.004900932312011719 s
mⱼ =  1 which took 0.06983494758605957 s
δⱼ =  0.026825649506906833
error =  0.7148055643585794 which took 0.12760639190673828 s
mesh refined which took 0.052507877349853516 s
---------- level 10 took 0.25748300552368164 s ------------

level =  11
nⱼ =  623 which took 0.00558018684387207 s
mⱼ =  1 which took 0.08188509941101074 s
δⱼ =  0.02443915323142596
error =  0.6695617347996469 which took 0.13857626914978027 s
mesh refined which took 0.06211352348327637 s
---------- level 11 took 0.28870654106140137 s ------------

level =  12
nⱼ =  819 which took 0.004020214080810547 s
mⱼ =  1 which took 0.1167607307434082 s
δⱼ =  0.023389341397738683
error =  0.5850432062725879 which took 0.17141461372375488 s
mesh refined which took 0.07935547828674316 s
---------- level 12 took 0.372891902923584 s ------------

level =  13
nⱼ =  961 which took 0.007401227951049805 s
mⱼ =  1 which took 0.1323988437652588 s
δⱼ =  0.022023946819350118
error =  0.5418286837849083 which took 0.19699382781982422 s
mesh refined which took 0.09448409080505371 s
---------- level 13 took 0.4319474697113037 s ------------

level =  14
nⱼ =  1260 which took 0.009451150894165039 s
mⱼ =  1 which took 0.18529176712036133 s
δⱼ =  0.019844210303101704
error =  0.47602764650560203 which took 0.24820971488952637 s
mesh refined which took 0.12563371658325195 s
---------- level 14 took 0.5691120624542236 s ------------

level =  15
nⱼ =  1707 which took 0.00891733169555664 s
mⱼ =  1 which took 0.2194375991821289 s
δⱼ =  0.017131453718458457
error =  0.4026382684575586 which took 0.3180503845214844 s
mesh refined which took 0.15901923179626465 s
---------- level 15 took 0.7074482440948486 s ------------

level =  16
nⱼ =  2098 which took 0.011401653289794922 s
mⱼ =  1 which took 0.23757624626159668 s
δⱼ =  0.015224487967946813
error =  0.3583536053260954 which took 0.3885939121246338 s
mesh refined which took 0.19638395309448242 s
---------- level 16 took 0.8347175121307373 s ------------

level =  17
nⱼ =  2283 which took 0.008605480194091797 s
mⱼ =  1 which took 0.2676551342010498 s
δⱼ =  0.012606383025640792
error =  0.3434232544791432 which took 0.4224700927734375 s
mesh refined which took 0.2239847183227539 s
---------- level 17 took 0.92364501953125 s ------------

level =  18
nⱼ =  2880 which took 0.011539936065673828 s
mⱼ =  1 which took 0.33236265182495117 s
δⱼ =  0.012566140558205202
error =  0.30749909012906385 which took 0.5286545753479004 s
mesh refined which took 0.2676844596862793 s
---------- level 18 took 1.1408066749572754 s ------------

level =  19
nⱼ =  3378 which took 0.012983560562133789 s
mⱼ =  1 which took 0.40049099922180176 s
δⱼ =  0.010748811889933655
error =  0.2847999212607402 which took 0.608201265335083 s
mesh refined which took 0.33404040336608887 s
---------- level 19 took 1.3562936782836914 s ------------

level =  20
nⱼ =  4357 which took 0.021976232528686523 s
mⱼ =  1 which took 0.5444016456604004 s
δⱼ =  0.011210744882890606
error =  0.25291155739219223 which took 0.7919211387634277 s
mesh refined which took 0.41776442527770996 s
---------- level 20 took 1.77695631980896 s ------------

level =  21
nⱼ =  5359 which took 0.01965928077697754 s
mⱼ =  1 which took 0.6854193210601807 s
δⱼ =  0.009339699023588263
error =  0.2286139427842272 which took 0.9583978652954102 s
mesh refined which took 0.5291705131530762 s
---------- level 21 took 2.194326162338257 s ------------

level =  22
nⱼ =  7296 which took 0.03328537940979004 s
mⱼ =  1 which took 0.9944047927856445 s
δⱼ =  0.008770277524762008
error =  0.19232029286044142 which took 1.3120157718658447 s
mesh refined which took 0.6808440685272217 s
---------- level 22 took 3.022439479827881 s ------------

level =  23
nⱼ =  8084 which took 0.03216195106506348 s
mⱼ =  1 which took 1.1401050090789795 s
δⱼ =  0.007120794745375671
error =  0.18094299266157793 which took 1.4295341968536377 s
mesh refined which took 0.769000768661499 s
---------- level 23 took 3.372718572616577 s ------------

level =  24
nⱼ =  8832 which took 0.03711867332458496 s
mⱼ =  1 which took 1.2743003368377686 s
δⱼ =  0.006466421239269172
error =  0.17341261391749752 which took 1.6140389442443848 s
mesh refined which took 0.8399348258972168 s
---------- level 24 took 3.766042947769165 s ------------

level =  25
nⱼ =  11493 which took 0.03813481330871582 s
mⱼ =  1 which took 1.745192050933838 s
δⱼ =  0.006614765739078151
error =  0.15353383851730865 which took 1.9197642803192139 s
mesh refined which took 1.052826166152954 s
---------- level 25 took 4.757462501525879 s ------------

level =  26
nⱼ =  13866 which took 0.04078507423400879 s
mⱼ =  1 which took 2.1436712741851807 s
δⱼ =  0.005650080151082483
error =  0.14019517289959332 which took 2.3767683506011963 s
mesh refined which took 1.2984693050384521 s
---------- level 26 took 5.861567497253418 s ------------

level =  27
nⱼ =  16844 which took 0.06576681137084961 s
mⱼ =  1 which took 2.8098318576812744 s
δⱼ =  0.005241454367270634
error =  0.1274075120266847 which took 3.178680896759033 s
mesh refined which took 1.7507026195526123 s
---------- level 27 took 7.805772304534912 s ------------

level =  28
nⱼ =  20090 which took 0.07349634170532227 s
mⱼ =  1 which took 3.734689712524414 s
δⱼ =  0.0047383528810985625

error =  0.11734961173473375 which took 3.534965991973877 s
mesh refined which took 1.911214828491211 s
---------- level 28 took 9.2551748752594 s ------------

level =  29
nⱼ =  26556 which took 0.08678126335144043 s
mⱼ =  1 which took 5.31338906288147 s
δⱼ =  0.004493402384120442
error =  0.10125171495706485 which took 4.598304748535156 s
mesh refined which took 2.576777696609497 s
---------- level 29 took 12.575908899307251 s ------------

level =  30
nⱼ =  31738 which took 0.10872888565063477 s
mⱼ =  1 which took 7.259253740310669 s
δⱼ =  0.0039092161992360175
error =  0.09169465407698195 which took 5.518657922744751 s
mesh refined which took 3.020479440689087 s
---------- level 30 took 15.908429384231567 s ------------

level =  31
nⱼ =  33780 which took 0.08083510398864746 s
mⱼ =  1 which took 8.01367473602295 s
δⱼ =  0.0033552513132664315
error =  0.08878583037880998 which took 5.970690727233887 s
mesh refined which took 3.3554673194885254 s
---------- level 31 took 17.42255210876465 s ------------

level =  32
nⱼ =  39209 which took 0.09499597549438477 s
mⱼ =  1 which took 9.896912097930908 s
δⱼ =  0.003241478575012659
error =  0.08299916685634527 which took 6.982182264328003 s
--------------------------------------------------------
tolerance reached, which took 32 levels and 

image

image

image

image

Final mesh: image

snippet2_v2.ipynb 8.57 KB
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment