From def88b4aa52dd56ccc1b9109a8030d849c3d4d8e Mon Sep 17 00:00:00 2001 From: Laura Prieto Saavedra <40216050+lpsaavedra@users.noreply.github.com> Date: Mon, 29 Jan 2024 10:31:52 -0500 Subject: [PATCH] Fix read checkpoint for the matrix-free application (#1007) Description of the problem The restart option was not working in parallel for the matrix-free application. Description of the solution Use locally relevant vectors for the present and previous solutions regardless of the type of vector. In addition, use a locally relevant vector for mesh refinement (in the init_temporary_vector() call). How Has This Been Tested? All existing tests pass. A new test tgv_restart_bdf1 is added to the MF application in serial and in parallel to verify restart feature. --- .../lethe-fluid-matrix-free/CMakeLists.txt | 9 + .../tgv_restart_bdf1.mpirun=1.output | 461 ++++++++++++++++++ .../tgv_restart_bdf1.mpirun=2.output | 461 ++++++++++++++++++ .../tgv_restart_bdf1.prm | 204 ++++++++ .../tgv_restart_bdf1/L2Error_FD.checkpoint | 1 + .../tgv_restart_bdf1/enstrophy.checkpoint | 1 + .../kinetic_energy.checkpoint | 1 + .../tgv_restart_bdf1/restart.pvdhandler | 2 + .../restart.simulationcontrol | 8 + .../tgv_restart_bdf1/restart.triangulation | Bin 0 -> 12544 bytes .../restart.triangulation.info | 2 + .../restart.triangulation_fixed.data | Bin 0 -> 397320 bytes source/solvers/navier_stokes_base.cc | 17 +- 13 files changed, 1157 insertions(+), 10 deletions(-) create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=1.output create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=2.output create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.prm create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/L2Error_FD.checkpoint create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/enstrophy.checkpoint create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/kinetic_energy.checkpoint create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.pvdhandler create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.simulationcontrol create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation.info create mode 100644 applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation_fixed.data diff --git a/applications_tests/lethe-fluid-matrix-free/CMakeLists.txt b/applications_tests/lethe-fluid-matrix-free/CMakeLists.txt index 719955bf51..77b57199b6 100644 --- a/applications_tests/lethe-fluid-matrix-free/CMakeLists.txt +++ b/applications_tests/lethe-fluid-matrix-free/CMakeLists.txt @@ -4,6 +4,15 @@ string(TOLOWER ${CMAKE_BUILD_TYPE} _build_type) file(COPY cylinder_structured.msh DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/cylinder_kelly_estimator.${_build_type}") +file(COPY tgv_restart_bdf1/restart.pvdhandler DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") +file(COPY tgv_restart_bdf1/restart.simulationcontrol DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") +file(COPY tgv_restart_bdf1/restart.triangulation DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") +file(COPY tgv_restart_bdf1/restart.triangulation_fixed.data DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") +file(COPY tgv_restart_bdf1/restart.triangulation.info DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") +file(COPY tgv_restart_bdf1/enstrophy.checkpoint DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") +file(COPY tgv_restart_bdf1/kinetic_energy.checkpoint DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") +file(COPY tgv_restart_bdf1/L2Error_FD.checkpoint DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/tgv_restart_bdf1.${_build_type}") + deal_ii_pickup_tests() if(CMAKE_BUILD_TYPE STREQUAL "Debug") diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=1.output b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=1.output new file mode 100644 index 0000000000..5806bdb792 --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=1.output @@ -0,0 +1,461 @@ +Running on 1 MPI rank(s)... + Number of active cells: 1 + Number of degrees of freedom: 12 + Volume of triangulation: 39.4784 +************************ +---> Simulation Restart +************************ + Number of active cells: 1024 + Number of degrees of freedom: 3267 + Volume of triangulation: 39.4784 + +******************************************************************************** +Transient iteration: 51 Time: 0.51 Time step: 0.01 CFL: 0.0166875 +******************************************************************************** +Enstrophy : 0.0660159 +Kinetic energy : 0.0329552 +L2 error velocity : 0.0110063 + +******************************************************************************** +Transient iteration: 52 Time: 0.52 Time step: 0.01 CFL: 0.0163593 +******************************************************************************** +Enstrophy : 0.0634444 +Kinetic energy : 0.0316715 +L2 error velocity : 0.0110007 + +******************************************************************************** +Transient iteration: 53 Time: 0.53 Time step: 0.01 CFL: 0.0160375 +******************************************************************************** +Enstrophy : 0.0609731 +Kinetic energy : 0.0304378 +L2 error velocity : 0.0109909 + +******************************************************************************** +Transient iteration: 54 Time: 0.54 Time step: 0.01 CFL: 0.0157221 +******************************************************************************** +Enstrophy : 0.058598 +Kinetic energy : 0.0292522 +L2 error velocity : 0.0109773 + +******************************************************************************** +Transient iteration: 55 Time: 0.55 Time step: 0.01 CFL: 0.0154128 +******************************************************************************** +Enstrophy : 0.0563155 +Kinetic energy : 0.0281128 +L2 error velocity : 0.0109599 + +******************************************************************************** +Transient iteration: 56 Time: 0.56 Time step: 0.01 CFL: 0.0151096 +******************************************************************************** +Enstrophy : 0.0541219 +Kinetic energy : 0.0270177 +L2 error velocity : 0.0109389 + +******************************************************************************** +Transient iteration: 57 Time: 0.57 Time step: 0.01 CFL: 0.0148124 +******************************************************************************** +Enstrophy : 0.0520137 +Kinetic energy : 0.0259653 +L2 error velocity : 0.0109145 + +******************************************************************************** +Transient iteration: 58 Time: 0.58 Time step: 0.01 CFL: 0.0145211 +******************************************************************************** +Enstrophy : 0.0499876 +Kinetic energy : 0.0249539 +L2 error velocity : 0.0108868 + +******************************************************************************** +Transient iteration: 59 Time: 0.59 Time step: 0.01 CFL: 0.0142355 +******************************************************************************** +Enstrophy : 0.0480405 +Kinetic energy : 0.0239819 +L2 error velocity : 0.0108559 + +******************************************************************************** +Transient iteration: 60 Time: 0.6 Time step: 0.01 CFL: 0.0139555 +******************************************************************************** +Enstrophy : 0.0461692 +Kinetic energy : 0.0230477 +L2 error velocity : 0.0108221 + +******************************************************************************* +Transient iteration: 61 Time: 0.61 Time step: 0.01 CFL: 0.013681 +******************************************************************************* +Enstrophy : 0.0443708 +Kinetic energy : 0.0221499 +L2 error velocity : 0.0107853 + +******************************************************************************** +Transient iteration: 62 Time: 0.62 Time step: 0.01 CFL: 0.0134119 +******************************************************************************** +Enstrophy : 0.0426424 +Kinetic energy : 0.0212872 +L2 error velocity : 0.0107458 + +******************************************************************************* +Transient iteration: 63 Time: 0.63 Time step: 0.01 CFL: 0.013148 +******************************************************************************* +Enstrophy : 0.0409814 +Kinetic energy : 0.020458 +L2 error velocity : 0.0107036 + +******************************************************************************** +Transient iteration: 64 Time: 0.64 Time step: 0.01 CFL: 0.0128894 +******************************************************************************** +Enstrophy : 0.0393851 +Kinetic energy : 0.0196611 +L2 error velocity : 0.0106589 + +******************************************************************************** +Transient iteration: 65 Time: 0.65 Time step: 0.01 CFL: 0.0126359 +******************************************************************************** +Enstrophy : 0.0378509 +Kinetic energy : 0.0188952 +L2 error velocity : 0.0106118 + +******************************************************************************** +Transient iteration: 66 Time: 0.66 Time step: 0.01 CFL: 0.0123874 +******************************************************************************** +Enstrophy : 0.0363765 +Kinetic energy : 0.0181592 +L2 error velocity : 0.0105624 + +******************************************************************************** +Transient iteration: 67 Time: 0.67 Time step: 0.01 CFL: 0.0121437 +******************************************************************************** +Enstrophy : 0.0349596 +Kinetic energy : 0.0174519 +L2 error velocity : 0.0105109 + +******************************************************************************** +Transient iteration: 68 Time: 0.68 Time step: 0.01 CFL: 0.0119048 +******************************************************************************** +Enstrophy : 0.0335978 +Kinetic energy : 0.0167721 +L2 error velocity : 0.0104572 + +******************************************************************************** +Transient iteration: 69 Time: 0.69 Time step: 0.01 CFL: 0.0116707 +******************************************************************************** +Enstrophy : 0.0322891 +Kinetic energy : 0.0161188 +L2 error velocity : 0.0104016 + +******************************************************************************** +Transient iteration: 70 Time: 0.7 Time step: 0.01 CFL: 0.0114411 +******************************************************************************** +Enstrophy : 0.0310314 +Kinetic energy : 0.0154909 +L2 error velocity : 0.0103441 + +******************************************************************************** +Transient iteration: 71 Time: 0.71 Time step: 0.01 CFL: 0.0112161 +******************************************************************************** +Enstrophy : 0.0298226 +Kinetic energy : 0.0148875 +L2 error velocity : 0.0102849 + +******************************************************************************** +Transient iteration: 72 Time: 0.72 Time step: 0.01 CFL: 0.0109955 +******************************************************************************** +Enstrophy : 0.0286609 +Kinetic energy : 0.0143076 +L2 error velocity : 0.0102239 + +******************************************************************************** +Transient iteration: 73 Time: 0.73 Time step: 0.01 CFL: 0.0107792 +******************************************************************************** +Enstrophy : 0.0275445 +Kinetic energy : 0.0137503 +L2 error velocity : 0.0101614 + +******************************************************************************** +Transient iteration: 74 Time: 0.74 Time step: 0.01 CFL: 0.0105672 +******************************************************************************** +Enstrophy : 0.0264716 +Kinetic energy : 0.0132147 +L2 error velocity : 0.0100973 + +******************************************************************************** +Transient iteration: 75 Time: 0.75 Time step: 0.01 CFL: 0.0103593 +******************************************************************************** +Enstrophy : 0.0254405 +Kinetic energy : 0.0126999 +L2 error velocity : 0.0100318 + +******************************************************************************** +Transient iteration: 76 Time: 0.76 Time step: 0.01 CFL: 0.0101555 +******************************************************************************** +Enstrophy : 0.0244495 +Kinetic energy : 0.0122052 +L2 error velocity : 0.00996494 + +********************************************************************************* +Transient iteration: 77 Time: 0.77 Time step: 0.01 CFL: 0.00995578 +********************************************************************************* +Enstrophy : 0.0234971 +Kinetic energy : 0.0117298 +L2 error velocity : 0.00989683 + +********************************************************************************* +Transient iteration: 78 Time: 0.78 Time step: 0.01 CFL: 0.00975995 +********************************************************************************* +Enstrophy : 0.0225819 +Kinetic energy : 0.0112729 +L2 error velocity : 0.00982753 + +********************************************************************************* +Transient iteration: 79 Time: 0.79 Time step: 0.01 CFL: 0.00956798 +********************************************************************************* +Enstrophy : 0.0217022 +Kinetic energy : 0.0108338 +L2 error velocity : 0.0097571 + +********************************************************************************* +Transient iteration: 80 Time: 0.8 Time step: 0.01 CFL: 0.00937978 +********************************************************************************* +Enstrophy : 0.0208569 +Kinetic energy : 0.0104118 +L2 error velocity : 0.00968563 + +********************************************************************************* +Transient iteration: 81 Time: 0.81 Time step: 0.01 CFL: 0.00919528 +********************************************************************************* +Enstrophy : 0.0200445 +Kinetic energy : 0.0100062 +L2 error velocity : 0.00961319 + +********************************************************************************* +Transient iteration: 82 Time: 0.82 Time step: 0.01 CFL: 0.00901441 +********************************************************************************* +Enstrophy : 0.0192637 +Kinetic energy : 0.00961645 +L2 error velocity : 0.00953983 + +******************************************************************************** +Transient iteration: 83 Time: 0.83 Time step: 0.01 CFL: 0.0088371 +******************************************************************************** +Enstrophy : 0.0185133 +Kinetic energy : 0.00924186 +L2 error velocity : 0.00946562 + +********************************************************************************* +Transient iteration: 84 Time: 0.84 Time step: 0.01 CFL: 0.00866328 +********************************************************************************* +Enstrophy : 0.0177922 +Kinetic energy : 0.00888187 +L2 error velocity : 0.00939062 + +********************************************************************************* +Transient iteration: 85 Time: 0.85 Time step: 0.01 CFL: 0.00849287 +********************************************************************************* +Enstrophy : 0.0170991 +Kinetic energy : 0.0085359 +L2 error velocity : 0.0093149 + +********************************************************************************* +Transient iteration: 86 Time: 0.86 Time step: 0.01 CFL: 0.00832582 +********************************************************************************* +Enstrophy : 0.0164331 +Kinetic energy : 0.00820341 +L2 error velocity : 0.00923852 + +********************************************************************************* +Transient iteration: 87 Time: 0.87 Time step: 0.01 CFL: 0.00816205 +********************************************************************************* +Enstrophy : 0.015793 +Kinetic energy : 0.00788386 +L2 error velocity : 0.00916152 + +********************************************************************************* +Transient iteration: 88 Time: 0.88 Time step: 0.01 CFL: 0.00800151 +********************************************************************************* +Enstrophy : 0.0151778 +Kinetic energy : 0.00757677 +L2 error velocity : 0.00908396 + +********************************************************************************* +Transient iteration: 89 Time: 0.89 Time step: 0.01 CFL: 0.00784412 +********************************************************************************* +Enstrophy : 0.0145866 +Kinetic energy : 0.00728163 +L2 error velocity : 0.00900589 + +********************************************************************************* +Transient iteration: 90 Time: 0.9 Time step: 0.01 CFL: 0.00768983 +********************************************************************************* +Enstrophy : 0.0140184 +Kinetic energy : 0.006998 +L2 error velocity : 0.00892737 + +********************************************************************************* +Transient iteration: 91 Time: 0.91 Time step: 0.01 CFL: 0.00753857 +********************************************************************************* +Enstrophy : 0.0134723 +Kinetic energy : 0.00672541 +L2 error velocity : 0.00884844 + +********************************************************************************* +Transient iteration: 92 Time: 0.92 Time step: 0.01 CFL: 0.00739029 +********************************************************************************* +Enstrophy : 0.0129476 +Kinetic energy : 0.00646344 +L2 error velocity : 0.00876915 + +********************************************************************************* +Transient iteration: 93 Time: 0.93 Time step: 0.01 CFL: 0.00724493 +********************************************************************************* +Enstrophy : 0.0124432 +Kinetic energy : 0.00621167 +L2 error velocity : 0.00868954 + +********************************************************************************* +Transient iteration: 94 Time: 0.94 Time step: 0.01 CFL: 0.00710242 +********************************************************************************* +Enstrophy : 0.0119585 +Kinetic energy : 0.00596971 +L2 error velocity : 0.00860966 + +********************************************************************************* +Transient iteration: 95 Time: 0.95 Time step: 0.01 CFL: 0.00696272 +********************************************************************************* +Enstrophy : 0.0114927 +Kinetic energy : 0.00573717 +L2 error velocity : 0.00852955 + +********************************************************************************* +Transient iteration: 96 Time: 0.96 Time step: 0.01 CFL: 0.00682576 +********************************************************************************* +Enstrophy : 0.011045 +Kinetic energy : 0.0055137 +L2 error velocity : 0.00844925 + +******************************************************************************** +Transient iteration: 97 Time: 0.97 Time step: 0.01 CFL: 0.0066915 +******************************************************************************** +Enstrophy : 0.0106148 +Kinetic energy : 0.00529892 +L2 error velocity : 0.00836879 + +********************************************************************************* +Transient iteration: 98 Time: 0.98 Time step: 0.01 CFL: 0.00655988 +********************************************************************************* +Enstrophy : 0.0102013 +Kinetic energy : 0.00509252 +L2 error velocity : 0.00828822 + +********************************************************************************* +Transient iteration: 99 Time: 0.99 Time step: 0.01 CFL: 0.00643085 +********************************************************************************* +Enstrophy : 0.00980397 +Kinetic energy : 0.00489415 +L2 error velocity : 0.00820757 + +********************************************************************************* +Transient iteration: 100 Time: 1 Time step: 0.01 CFL: 0.00630436 +********************************************************************************* +Enstrophy : 0.00942208 +Kinetic energy : 0.00470351 +L2 error velocity : 0.00812688 + time error_velocity +0.0100 9.702604e-04 +0.0200 1.373517e-03 +0.0300 1.840106e-03 +0.0400 2.317500e-03 +0.0500 2.788763e-03 +0.0600 3.247543e-03 +0.0700 3.691204e-03 +0.0800 4.118644e-03 +0.0900 4.529473e-03 +0.1000 4.923664e-03 +0.1100 5.301385e-03 +0.1200 5.662912e-03 +0.1300 6.008589e-03 +0.1400 6.338792e-03 +0.1500 6.653919e-03 +0.1600 6.954380e-03 +0.1700 7.240586e-03 +0.1800 7.512949e-03 +0.1900 7.771878e-03 +0.2000 8.017774e-03 +0.2100 8.251036e-03 +0.2200 8.472051e-03 +0.2300 8.681201e-03 +0.2400 8.878859e-03 +0.2500 9.065390e-03 +0.2600 9.241148e-03 +0.2700 9.406480e-03 +0.2800 9.561726e-03 +0.2900 9.707214e-03 +0.3000 9.843266e-03 +0.3100 9.970194e-03 +0.3200 1.008830e-02 +0.3300 1.019789e-02 +0.3400 1.029924e-02 +0.3500 1.039264e-02 +0.3600 1.047835e-02 +0.3700 1.055665e-02 +0.3800 1.062779e-02 +0.3900 1.069203e-02 +0.4000 1.074960e-02 +0.4100 1.080075e-02 +0.4200 1.084570e-02 +0.4300 1.088467e-02 +0.4400 1.091790e-02 +0.4500 1.094558e-02 +0.4600 1.096792e-02 +0.4700 1.098512e-02 +0.4800 1.099738e-02 +0.4900 1.100489e-02 +0.5000 1.100781e-02 +0.5100 1.100635e-02 +0.5200 1.100065e-02 +0.5300 1.099090e-02 +0.5400 1.097726e-02 +0.5500 1.095987e-02 +0.5600 1.093890e-02 +0.5700 1.091449e-02 +0.5800 1.088679e-02 +0.5900 1.085593e-02 +0.6000 1.082205e-02 +0.6100 1.078529e-02 +0.6200 1.074576e-02 +0.6300 1.070359e-02 +0.6400 1.065891e-02 +0.6500 1.061182e-02 +0.6600 1.056244e-02 +0.6700 1.051088e-02 +0.6800 1.045724e-02 +0.6900 1.040163e-02 +0.7000 1.034414e-02 +0.7100 1.028487e-02 +0.7200 1.022391e-02 +0.7300 1.016135e-02 +0.7400 1.009728e-02 +0.7500 1.003178e-02 +0.7600 9.964941e-03 +0.7700 9.896829e-03 +0.7800 9.827526e-03 +0.7900 9.757104e-03 +0.8000 9.685634e-03 +0.8100 9.613187e-03 +0.8200 9.539827e-03 +0.8300 9.465619e-03 +0.8400 9.390625e-03 +0.8500 9.314905e-03 +0.8600 9.238516e-03 +0.8700 9.161516e-03 +0.8800 9.083957e-03 +0.8900 9.005891e-03 +0.9000 8.927369e-03 +0.9100 8.848439e-03 +0.9200 8.769148e-03 +0.9300 8.689541e-03 +0.9400 8.609660e-03 +0.9500 8.529549e-03 +0.9600 8.449246e-03 +0.9700 8.368792e-03 +0.9800 8.288222e-03 +0.9900 8.207573e-03 +1.0000 8.126880e-03 \ No newline at end of file diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=2.output b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=2.output new file mode 100644 index 0000000000..8256a90da3 --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.mpirun=2.output @@ -0,0 +1,461 @@ +Running on 2 MPI rank(s)... + Number of active cells: 1 + Number of degrees of freedom: 12 + Volume of triangulation: 39.4784 +************************ +---> Simulation Restart +************************ + Number of active cells: 1024 + Number of degrees of freedom: 3267 + Volume of triangulation: 39.4784 + +******************************************************************************** +Transient iteration: 51 Time: 0.51 Time step: 0.01 CFL: 0.0166875 +******************************************************************************** +Enstrophy : 0.0660159 +Kinetic energy : 0.0329552 +L2 error velocity : 0.0110063 + +******************************************************************************** +Transient iteration: 52 Time: 0.52 Time step: 0.01 CFL: 0.0163593 +******************************************************************************** +Enstrophy : 0.0634444 +Kinetic energy : 0.0316715 +L2 error velocity : 0.0110007 + +******************************************************************************** +Transient iteration: 53 Time: 0.53 Time step: 0.01 CFL: 0.0160375 +******************************************************************************** +Enstrophy : 0.0609731 +Kinetic energy : 0.0304378 +L2 error velocity : 0.0109909 + +******************************************************************************** +Transient iteration: 54 Time: 0.54 Time step: 0.01 CFL: 0.0157221 +******************************************************************************** +Enstrophy : 0.058598 +Kinetic energy : 0.0292522 +L2 error velocity : 0.0109773 + +******************************************************************************** +Transient iteration: 55 Time: 0.55 Time step: 0.01 CFL: 0.0154128 +******************************************************************************** +Enstrophy : 0.0563155 +Kinetic energy : 0.0281128 +L2 error velocity : 0.0109599 + +******************************************************************************** +Transient iteration: 56 Time: 0.56 Time step: 0.01 CFL: 0.0151096 +******************************************************************************** +Enstrophy : 0.0541219 +Kinetic energy : 0.0270177 +L2 error velocity : 0.0109389 + +******************************************************************************** +Transient iteration: 57 Time: 0.57 Time step: 0.01 CFL: 0.0148124 +******************************************************************************** +Enstrophy : 0.0520137 +Kinetic energy : 0.0259653 +L2 error velocity : 0.0109145 + +******************************************************************************** +Transient iteration: 58 Time: 0.58 Time step: 0.01 CFL: 0.0145211 +******************************************************************************** +Enstrophy : 0.0499876 +Kinetic energy : 0.0249539 +L2 error velocity : 0.0108868 + +******************************************************************************** +Transient iteration: 59 Time: 0.59 Time step: 0.01 CFL: 0.0142355 +******************************************************************************** +Enstrophy : 0.0480405 +Kinetic energy : 0.0239819 +L2 error velocity : 0.0108559 + +******************************************************************************** +Transient iteration: 60 Time: 0.6 Time step: 0.01 CFL: 0.0139555 +******************************************************************************** +Enstrophy : 0.0461692 +Kinetic energy : 0.0230477 +L2 error velocity : 0.0108221 + +******************************************************************************* +Transient iteration: 61 Time: 0.61 Time step: 0.01 CFL: 0.013681 +******************************************************************************* +Enstrophy : 0.0443708 +Kinetic energy : 0.0221499 +L2 error velocity : 0.0107853 + +******************************************************************************** +Transient iteration: 62 Time: 0.62 Time step: 0.01 CFL: 0.0134119 +******************************************************************************** +Enstrophy : 0.0426424 +Kinetic energy : 0.0212872 +L2 error velocity : 0.0107458 + +******************************************************************************* +Transient iteration: 63 Time: 0.63 Time step: 0.01 CFL: 0.013148 +******************************************************************************* +Enstrophy : 0.0409814 +Kinetic energy : 0.020458 +L2 error velocity : 0.0107036 + +******************************************************************************** +Transient iteration: 64 Time: 0.64 Time step: 0.01 CFL: 0.0128894 +******************************************************************************** +Enstrophy : 0.0393851 +Kinetic energy : 0.0196611 +L2 error velocity : 0.0106589 + +******************************************************************************** +Transient iteration: 65 Time: 0.65 Time step: 0.01 CFL: 0.0126359 +******************************************************************************** +Enstrophy : 0.0378509 +Kinetic energy : 0.0188952 +L2 error velocity : 0.0106118 + +******************************************************************************** +Transient iteration: 66 Time: 0.66 Time step: 0.01 CFL: 0.0123874 +******************************************************************************** +Enstrophy : 0.0363765 +Kinetic energy : 0.0181592 +L2 error velocity : 0.0105624 + +******************************************************************************** +Transient iteration: 67 Time: 0.67 Time step: 0.01 CFL: 0.0121437 +******************************************************************************** +Enstrophy : 0.0349596 +Kinetic energy : 0.0174519 +L2 error velocity : 0.0105109 + +******************************************************************************** +Transient iteration: 68 Time: 0.68 Time step: 0.01 CFL: 0.0119048 +******************************************************************************** +Enstrophy : 0.0335978 +Kinetic energy : 0.0167721 +L2 error velocity : 0.0104572 + +******************************************************************************** +Transient iteration: 69 Time: 0.69 Time step: 0.01 CFL: 0.0116707 +******************************************************************************** +Enstrophy : 0.0322891 +Kinetic energy : 0.0161188 +L2 error velocity : 0.0104016 + +******************************************************************************** +Transient iteration: 70 Time: 0.7 Time step: 0.01 CFL: 0.0114411 +******************************************************************************** +Enstrophy : 0.0310314 +Kinetic energy : 0.0154909 +L2 error velocity : 0.0103441 + +******************************************************************************** +Transient iteration: 71 Time: 0.71 Time step: 0.01 CFL: 0.0112161 +******************************************************************************** +Enstrophy : 0.0298226 +Kinetic energy : 0.0148875 +L2 error velocity : 0.0102849 + +******************************************************************************** +Transient iteration: 72 Time: 0.72 Time step: 0.01 CFL: 0.0109955 +******************************************************************************** +Enstrophy : 0.0286609 +Kinetic energy : 0.0143076 +L2 error velocity : 0.0102239 + +******************************************************************************** +Transient iteration: 73 Time: 0.73 Time step: 0.01 CFL: 0.0107792 +******************************************************************************** +Enstrophy : 0.0275445 +Kinetic energy : 0.0137503 +L2 error velocity : 0.0101614 + +******************************************************************************** +Transient iteration: 74 Time: 0.74 Time step: 0.01 CFL: 0.0105672 +******************************************************************************** +Enstrophy : 0.0264716 +Kinetic energy : 0.0132147 +L2 error velocity : 0.0100973 + +******************************************************************************** +Transient iteration: 75 Time: 0.75 Time step: 0.01 CFL: 0.0103593 +******************************************************************************** +Enstrophy : 0.0254405 +Kinetic energy : 0.0126999 +L2 error velocity : 0.0100318 + +******************************************************************************** +Transient iteration: 76 Time: 0.76 Time step: 0.01 CFL: 0.0101555 +******************************************************************************** +Enstrophy : 0.0244495 +Kinetic energy : 0.0122052 +L2 error velocity : 0.00996494 + +********************************************************************************* +Transient iteration: 77 Time: 0.77 Time step: 0.01 CFL: 0.00995578 +********************************************************************************* +Enstrophy : 0.0234971 +Kinetic energy : 0.0117298 +L2 error velocity : 0.00989683 + +********************************************************************************* +Transient iteration: 78 Time: 0.78 Time step: 0.01 CFL: 0.00975995 +********************************************************************************* +Enstrophy : 0.0225819 +Kinetic energy : 0.0112729 +L2 error velocity : 0.00982753 + +********************************************************************************* +Transient iteration: 79 Time: 0.79 Time step: 0.01 CFL: 0.00956798 +********************************************************************************* +Enstrophy : 0.0217022 +Kinetic energy : 0.0108338 +L2 error velocity : 0.0097571 + +********************************************************************************* +Transient iteration: 80 Time: 0.8 Time step: 0.01 CFL: 0.00937978 +********************************************************************************* +Enstrophy : 0.0208569 +Kinetic energy : 0.0104118 +L2 error velocity : 0.00968563 + +********************************************************************************* +Transient iteration: 81 Time: 0.81 Time step: 0.01 CFL: 0.00919528 +********************************************************************************* +Enstrophy : 0.0200445 +Kinetic energy : 0.0100062 +L2 error velocity : 0.00961319 + +********************************************************************************* +Transient iteration: 82 Time: 0.82 Time step: 0.01 CFL: 0.00901441 +********************************************************************************* +Enstrophy : 0.0192637 +Kinetic energy : 0.00961645 +L2 error velocity : 0.00953983 + +******************************************************************************** +Transient iteration: 83 Time: 0.83 Time step: 0.01 CFL: 0.0088371 +******************************************************************************** +Enstrophy : 0.0185133 +Kinetic energy : 0.00924186 +L2 error velocity : 0.00946562 + +********************************************************************************* +Transient iteration: 84 Time: 0.84 Time step: 0.01 CFL: 0.00866328 +********************************************************************************* +Enstrophy : 0.0177922 +Kinetic energy : 0.00888187 +L2 error velocity : 0.00939062 + +********************************************************************************* +Transient iteration: 85 Time: 0.85 Time step: 0.01 CFL: 0.00849287 +********************************************************************************* +Enstrophy : 0.0170991 +Kinetic energy : 0.0085359 +L2 error velocity : 0.0093149 + +********************************************************************************* +Transient iteration: 86 Time: 0.86 Time step: 0.01 CFL: 0.00832582 +********************************************************************************* +Enstrophy : 0.0164331 +Kinetic energy : 0.00820341 +L2 error velocity : 0.00923852 + +********************************************************************************* +Transient iteration: 87 Time: 0.87 Time step: 0.01 CFL: 0.00816205 +********************************************************************************* +Enstrophy : 0.015793 +Kinetic energy : 0.00788386 +L2 error velocity : 0.00916152 + +********************************************************************************* +Transient iteration: 88 Time: 0.88 Time step: 0.01 CFL: 0.00800151 +********************************************************************************* +Enstrophy : 0.0151778 +Kinetic energy : 0.00757677 +L2 error velocity : 0.00908396 + +********************************************************************************* +Transient iteration: 89 Time: 0.89 Time step: 0.01 CFL: 0.00784412 +********************************************************************************* +Enstrophy : 0.0145866 +Kinetic energy : 0.00728163 +L2 error velocity : 0.00900589 + +********************************************************************************* +Transient iteration: 90 Time: 0.9 Time step: 0.01 CFL: 0.00768983 +********************************************************************************* +Enstrophy : 0.0140184 +Kinetic energy : 0.006998 +L2 error velocity : 0.00892737 + +********************************************************************************* +Transient iteration: 91 Time: 0.91 Time step: 0.01 CFL: 0.00753857 +********************************************************************************* +Enstrophy : 0.0134723 +Kinetic energy : 0.00672541 +L2 error velocity : 0.00884844 + +********************************************************************************* +Transient iteration: 92 Time: 0.92 Time step: 0.01 CFL: 0.00739029 +********************************************************************************* +Enstrophy : 0.0129476 +Kinetic energy : 0.00646344 +L2 error velocity : 0.00876915 + +********************************************************************************* +Transient iteration: 93 Time: 0.93 Time step: 0.01 CFL: 0.00724493 +********************************************************************************* +Enstrophy : 0.0124432 +Kinetic energy : 0.00621167 +L2 error velocity : 0.00868954 + +********************************************************************************* +Transient iteration: 94 Time: 0.94 Time step: 0.01 CFL: 0.00710242 +********************************************************************************* +Enstrophy : 0.0119585 +Kinetic energy : 0.00596971 +L2 error velocity : 0.00860966 + +********************************************************************************* +Transient iteration: 95 Time: 0.95 Time step: 0.01 CFL: 0.00696272 +********************************************************************************* +Enstrophy : 0.0114927 +Kinetic energy : 0.00573717 +L2 error velocity : 0.00852955 + +********************************************************************************* +Transient iteration: 96 Time: 0.96 Time step: 0.01 CFL: 0.00682576 +********************************************************************************* +Enstrophy : 0.011045 +Kinetic energy : 0.0055137 +L2 error velocity : 0.00844925 + +******************************************************************************** +Transient iteration: 97 Time: 0.97 Time step: 0.01 CFL: 0.0066915 +******************************************************************************** +Enstrophy : 0.0106148 +Kinetic energy : 0.00529892 +L2 error velocity : 0.00836879 + +********************************************************************************* +Transient iteration: 98 Time: 0.98 Time step: 0.01 CFL: 0.00655988 +********************************************************************************* +Enstrophy : 0.0102013 +Kinetic energy : 0.00509252 +L2 error velocity : 0.00828822 + +********************************************************************************* +Transient iteration: 99 Time: 0.99 Time step: 0.01 CFL: 0.00643085 +********************************************************************************* +Enstrophy : 0.00980397 +Kinetic energy : 0.00489415 +L2 error velocity : 0.00820757 + +********************************************************************************* +Transient iteration: 100 Time: 1 Time step: 0.01 CFL: 0.00630436 +********************************************************************************* +Enstrophy : 0.00942208 +Kinetic energy : 0.00470351 +L2 error velocity : 0.00812688 + time error_velocity +0.0100 9.702604e-04 +0.0200 1.373517e-03 +0.0300 1.840106e-03 +0.0400 2.317500e-03 +0.0500 2.788763e-03 +0.0600 3.247543e-03 +0.0700 3.691204e-03 +0.0800 4.118644e-03 +0.0900 4.529473e-03 +0.1000 4.923664e-03 +0.1100 5.301385e-03 +0.1200 5.662912e-03 +0.1300 6.008589e-03 +0.1400 6.338792e-03 +0.1500 6.653919e-03 +0.1600 6.954380e-03 +0.1700 7.240586e-03 +0.1800 7.512949e-03 +0.1900 7.771878e-03 +0.2000 8.017774e-03 +0.2100 8.251036e-03 +0.2200 8.472051e-03 +0.2300 8.681201e-03 +0.2400 8.878859e-03 +0.2500 9.065390e-03 +0.2600 9.241148e-03 +0.2700 9.406480e-03 +0.2800 9.561726e-03 +0.2900 9.707214e-03 +0.3000 9.843266e-03 +0.3100 9.970194e-03 +0.3200 1.008830e-02 +0.3300 1.019789e-02 +0.3400 1.029924e-02 +0.3500 1.039264e-02 +0.3600 1.047835e-02 +0.3700 1.055665e-02 +0.3800 1.062779e-02 +0.3900 1.069203e-02 +0.4000 1.074960e-02 +0.4100 1.080075e-02 +0.4200 1.084570e-02 +0.4300 1.088467e-02 +0.4400 1.091790e-02 +0.4500 1.094558e-02 +0.4600 1.096792e-02 +0.4700 1.098512e-02 +0.4800 1.099738e-02 +0.4900 1.100489e-02 +0.5000 1.100781e-02 +0.5100 1.100635e-02 +0.5200 1.100065e-02 +0.5300 1.099090e-02 +0.5400 1.097726e-02 +0.5500 1.095987e-02 +0.5600 1.093890e-02 +0.5700 1.091449e-02 +0.5800 1.088679e-02 +0.5900 1.085593e-02 +0.6000 1.082205e-02 +0.6100 1.078529e-02 +0.6200 1.074576e-02 +0.6300 1.070359e-02 +0.6400 1.065891e-02 +0.6500 1.061182e-02 +0.6600 1.056244e-02 +0.6700 1.051088e-02 +0.6800 1.045724e-02 +0.6900 1.040163e-02 +0.7000 1.034414e-02 +0.7100 1.028487e-02 +0.7200 1.022391e-02 +0.7300 1.016135e-02 +0.7400 1.009728e-02 +0.7500 1.003178e-02 +0.7600 9.964941e-03 +0.7700 9.896829e-03 +0.7800 9.827526e-03 +0.7900 9.757104e-03 +0.8000 9.685634e-03 +0.8100 9.613187e-03 +0.8200 9.539827e-03 +0.8300 9.465619e-03 +0.8400 9.390625e-03 +0.8500 9.314905e-03 +0.8600 9.238516e-03 +0.8700 9.161516e-03 +0.8800 9.083957e-03 +0.8900 9.005891e-03 +0.9000 8.927369e-03 +0.9100 8.848439e-03 +0.9200 8.769148e-03 +0.9300 8.689541e-03 +0.9400 8.609660e-03 +0.9500 8.529549e-03 +0.9600 8.449246e-03 +0.9700 8.368792e-03 +0.9800 8.288222e-03 +0.9900 8.207573e-03 +1.0000 8.126880e-03 \ No newline at end of file diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.prm b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.prm new file mode 100644 index 0000000000..ca4c5d36a7 --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1.prm @@ -0,0 +1,204 @@ +# Listing of Parameters +#---------------------- + +set dimension = 2 + +#--------------------------------------------------- +# Simulation Control +#--------------------------------------------------- + +subsection simulation control + set method = bdf1 + set time step = 0.01 # Time step + set number mesh adapt = 0 # If steady, nb mesh adaptation + set time end = 1.00 # End time of simulation + set output frequency = 0 # Frequency of simulation output + set output path = ../ +end + +#--------------------------------------------------- +# FEM +#--------------------------------------------------- + +subsection FEM + set velocity order = 1 + set pressure order = 1 +end + +#--------------------------------------------------- +# Timer +#--------------------------------------------------- + +subsection timer + set type = none # +end + +#--------------------------------------------------- +# Initial condition +#--------------------------------------------------- + +subsection initial conditions + set type = L2projection + subsection uvwp + set Function expression = cos(x)*sin(y); -sin(x)*cos(y); -1./4*(cos(2*x)+cos(2*y)); + end +end + +#--------------------------------------------------- +# Physical Properties +#--------------------------------------------------- + +subsection physical properties + set number of fluids = 1 + subsection fluid 0 + set kinematic viscosity = 1.000 + end +end + +#--------------------------------------------------- +# Analytical Solution +#--------------------------------------------------- + +subsection analytical solution + set verbosity = verbose + set enable = true + subsection uvwp + set Function constants = viscosity=1 + set Function expression = exp(-2*viscosity*t)* cos(x) * sin(y) ; -sin(x) * cos(y) * exp(-2*viscosity*t); 0 + end +end + +#-------------------------------------------------- +# Restart +#-------------------------------------------------- + +subsection restart + # Enable checkpointing + set checkpoint = false + # Prefix for the filename of checkpoints + set filename = restart + # Frequency for checkpointing + set frequency = 1 + # Enable restart from checkpoint + set restart = true +end + +#--------------------------------------------------- +# Post-Processing +#--------------------------------------------------- + +subsection post-processing + set verbosity = verbose + set calculate enstrophy = true + set calculate kinetic energy = true +end + +#--------------------------------------------------- +# Mesh +#--------------------------------------------------- + +subsection mesh + set type = dealii + set grid type = hyper_cube + set grid arguments = 0 : 6.28318530718 : true + set initial refinement = 5 # initial mesh refinement +end + +#--------------------------------------------------- +# Mesh Adaptation Control +#--------------------------------------------------- + +subsection mesh adaptation + set type = none +end + +#--------------------------------------------------- +# Boundary Conditions +#--------------------------------------------------- + +subsection boundary conditions + set number = 2 + subsection bc 0 + set type = periodic + set id = 0 + set periodic_id = 1 + set periodic_direction = 0 + end + subsection bc 1 + set type = periodic + set id = 2 + set periodic_id = 3 + set periodic_direction = 1 + end +end + +#--------------------------------------------------- +# Non-Linear Solver Control +#--------------------------------------------------- + +subsection non-linear solver + subsection fluid dynamics + set verbosity = quiet + set tolerance = 1e-6 + set max iterations = 5 + end +end + +#--------------------------------------------------- +# Linear Solver Control +#--------------------------------------------------- + +subsection linear solver + subsection fluid dynamics + set verbosity = quiet + set method = gmres + set max iters = 5000 + set relative residual = 1e-4 + set minimum residual = 1e-9 + set preconditioner = ilu + set ilu preconditioner fill = 0 + set ilu preconditioner absolute tolerance = 1e-10 + set ilu preconditioner relative tolerance = 1.00 + end +end + +#--------------------------------------------------- +# Linear Solver Control +#--------------------------------------------------- + +subsection linear solver + subsection fluid dynamics + set method = gmres + set max iters = 5000 + set relative residual = 1e-4 + set minimum residual = 1e-7 + set preconditioner = gcmg + set verbosity = quiet + + #MG parameters + set mg verbosity = quiet + set mg min level = -1 + set mg level min cells = 16 + + #smoother + set mg smoother iterations = 10 + set mg smoother eig estimation = true + + # Eigenvalue estimation parameters + set eig estimation degree = 3 + set eig estimation smoothing range = 5 + set eig estimation cg n iterations = 20 + set eig estimation verbosity = quiet + + #coarse-grid solver + set mg coarse grid max iterations = 2000 + set mg coarse grid tolerance = 1e-7 + set mg coarse grid reduce = 1e-4 + set mg coarse grid max krylov vectors = 30 + set mg coarse grid preconditioner = ilu + + set ilu preconditioner fill = 1 + set ilu preconditioner absolute tolerance = 1e-10 + set ilu preconditioner relative tolerance = 1.00 + end +end \ No newline at end of file diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/L2Error_FD.checkpoint b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/L2Error_FD.checkpoint new file mode 100644 index 0000000000..7baf4765e2 --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/L2Error_FD.checkpoint @@ -0,0 +1 @@ +0 0 0 0 2 0 4 time 14 error_velocity 0 0 2 0 0 0 14 error_velocity 0 0 0 0 50 0 0 0 100 9.70260438756851947e-04 100 1.37351724439845491e-03 100 1.84010607302818440e-03 100 2.31749994121434625e-03 100 2.78876322761135520e-03 100 3.24754270524175775e-03 100 3.69120352527292075e-03 100 4.11864350666947716e-03 100 4.52947302130056548e-03 100 4.92366418756171736e-03 100 5.30138467716634483e-03 100 5.66291238323648246e-03 100 6.00858869287675037e-03 100 6.33879153008128950e-03 100 6.65391912906595270e-03 100 6.95437992264262127e-03 100 7.24058606285793438e-03 100 7.51294917727212980e-03 100 7.77187754425011975e-03 100 8.01777419334964020e-03 100 8.25103562307126157e-03 100 8.47205093909883522e-03 100 8.68120128407522638e-03 100 8.87885947266102427e-03 100 9.06538977306088449e-03 100 9.24114779421839234e-03 100 9.40648044994206629e-03 100 9.56172597941684979e-03 100 9.70721400924101641e-03 100 9.84326564608971506e-03 100 9.97019359194816109e-03 100 1.00883022758919815e-02 100 1.01978879978666098e-02 100 1.02992390810181089e-02 100 1.03926360299361770e-02 100 1.04783516927750304e-02 100 1.05566514256825311e-02 100 1.06277932583167135e-02 100 1.06920280595025063e-02 100 1.07495997022842796e-02 100 1.08007452277957966e-02 100 1.08456950075020674e-02 100 1.08846729034549842e-02 100 1.09178964262987305e-02 100 1.09455768908119956e-02 100 1.09679195688370847e-02 100 1.09851238394703850e-02 100 1.09973833364404670e-02 100 1.10048860926087438e-02 100 1.10078146815611386e-02 14 error_velocity 1 c 4 0 0 6 4 time 50 0 100 1.00000000000000002e-02 100 2.00000000000000004e-02 100 2.99999999999999989e-02 100 4.00000000000000008e-02 100 5.00000000000000028e-02 100 6.00000000000000047e-02 100 7.00000000000000067e-02 100 8.00000000000000017e-02 100 8.99999999999999967e-02 100 9.99999999999999917e-02 100 1.09999999999999987e-01 100 1.19999999999999982e-01 100 1.29999999999999977e-01 100 1.39999999999999986e-01 100 1.49999999999999994e-01 100 1.60000000000000003e-01 100 1.70000000000000012e-01 100 1.80000000000000021e-01 100 1.90000000000000030e-01 100 2.00000000000000039e-01 100 2.10000000000000048e-01 100 2.20000000000000057e-01 100 2.30000000000000066e-01 100 2.40000000000000074e-01 100 2.50000000000000056e-01 100 2.60000000000000064e-01 100 2.70000000000000073e-01 100 2.80000000000000082e-01 100 2.90000000000000091e-01 100 3.00000000000000100e-01 100 3.10000000000000109e-01 100 3.20000000000000118e-01 100 3.30000000000000127e-01 100 3.40000000000000135e-01 100 3.50000000000000144e-01 100 3.60000000000000153e-01 100 3.70000000000000162e-01 100 3.80000000000000171e-01 100 3.90000000000000180e-01 100 4.00000000000000189e-01 100 4.10000000000000198e-01 100 4.20000000000000207e-01 100 4.30000000000000215e-01 100 4.40000000000000224e-01 100 4.50000000000000233e-01 100 4.60000000000000242e-01 100 4.70000000000000251e-01 100 4.80000000000000260e-01 100 4.90000000000000269e-01 100 5.00000000000000000e-01 4 time 1 c 4 0 0 6 0 0 0 0 0 0 0 0 0 0 0 diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/enstrophy.checkpoint b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/enstrophy.checkpoint new file mode 100644 index 0000000000..c23d996adb --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/enstrophy.checkpoint @@ -0,0 +1 @@ +0 0 0 0 2 0 4 time 9 enstrophy 0 0 2 0 0 0 9 enstrophy 0 0 0 0 51 0 0 0 100 5.00800056059643395e-01 100 4.81292453184982871e-01 100 4.62544734156999260e-01 100 4.44527299704690182e-01 100 4.27211702324144738e-01 100 4.10570602690655428e-01 100 3.94577726493181846e-01 100 3.79207822943402106e-01 100 3.64436624901343820e-01 100 3.50240810554676485e-01 100 3.36597966591067610e-01 100 3.23486552805436200e-01 100 3.10885868086191564e-01 100 2.98776017726692045e-01 100 2.87137882010295609e-01 100 2.75953086019365912e-01 100 2.65203970620546070e-01 100 2.54873564580438028e-01 100 2.44945557767675820e-01 100 2.35404275399022533e-01 100 2.26234653288830939e-01 100 2.17422214062759800e-01 100 2.08953044298186313e-01 100 2.00813772555184233e-01 100 1.92991548263393459e-01 100 1.85474021431407649e-01 100 1.78249323146645405e-01 100 1.71306046834909798e-01 100 1.64633230250013057e-01 100 1.58220338165055374e-01 100 1.52057245737990759e-01 100 1.46134222525220264e-01 100 1.40441917117969300e-01 100 1.34971342377177811e-01 100 1.29713861243581335e-01 100 1.24661173100587941e-01 100 1.19805300668398645e-01 100 1.15138577408682025e-01 100 1.10653635419910273e-01 100 1.06343393804241890e-01 100 1.02201047487581465e-01 100 9.82200564751545685e-02 100 9.43941355256398690e-02 100 9.07172442275475249e-02 100 8.71835774621714382e-02 100 8.37875562380576699e-02 100 8.05238188825218026e-02 100 7.73872125762883750e-02 100 7.43727852179072585e-02 100 7.14757776050800875e-02 100 6.86916159205671339e-02 9 enstrophy 1 c 12 0 0 14 4 time 51 0 100 0.00000000000000000e+00 100 1.00000000000000002e-02 100 2.00000000000000004e-02 100 2.99999999999999989e-02 100 4.00000000000000008e-02 100 5.00000000000000028e-02 100 6.00000000000000047e-02 100 7.00000000000000067e-02 100 8.00000000000000017e-02 100 8.99999999999999967e-02 100 9.99999999999999917e-02 100 1.09999999999999987e-01 100 1.19999999999999982e-01 100 1.29999999999999977e-01 100 1.39999999999999986e-01 100 1.49999999999999994e-01 100 1.60000000000000003e-01 100 1.70000000000000012e-01 100 1.80000000000000021e-01 100 1.90000000000000030e-01 100 2.00000000000000039e-01 100 2.10000000000000048e-01 100 2.20000000000000057e-01 100 2.30000000000000066e-01 100 2.40000000000000074e-01 100 2.50000000000000056e-01 100 2.60000000000000064e-01 100 2.70000000000000073e-01 100 2.80000000000000082e-01 100 2.90000000000000091e-01 100 3.00000000000000100e-01 100 3.10000000000000109e-01 100 3.20000000000000118e-01 100 3.30000000000000127e-01 100 3.40000000000000135e-01 100 3.50000000000000144e-01 100 3.60000000000000153e-01 100 3.70000000000000162e-01 100 3.80000000000000171e-01 100 3.90000000000000180e-01 100 4.00000000000000189e-01 100 4.10000000000000198e-01 100 4.20000000000000207e-01 100 4.30000000000000215e-01 100 4.40000000000000224e-01 100 4.50000000000000233e-01 100 4.60000000000000242e-01 100 4.70000000000000251e-01 100 4.80000000000000260e-01 100 4.90000000000000269e-01 100 5.00000000000000000e-01 4 time 1 c 12 0 0 14 0 0 0 0 0 0 0 0 0 0 0 diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/kinetic_energy.checkpoint b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/kinetic_energy.checkpoint new file mode 100644 index 0000000000..646b1f5078 --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/kinetic_energy.checkpoint @@ -0,0 +1 @@ +0 0 0 0 2 0 4 time 14 kinetic-energy 0 0 2 0 0 0 14 kinetic-energy 0 0 0 0 51 0 0 0 100 2.49999992602278648e-01 100 2.40261774191911220e-01 100 2.30902890185604948e-01 100 2.21908565118993956e-01 100 2.13264597968605696e-01 100 2.04957340910372487e-01 100 1.96973677776323297e-01 100 1.89301003342272928e-01 100 1.81927203422487666e-01 100 1.74840635739944933e-01 100 1.68030111541929122e-01 100 1.61484877931928866e-01 100 1.55194600889924122e-01 100 1.49149348954229294e-01 100 1.43339577539106894e-01 100 1.37756113863395757e-01 100 1.32390142466314364e-01 100 1.27233191287579483e-01 100 1.22277118289848358e-01 100 1.17514098602333875e-01 100 1.12936612165310174e-01 100 1.08537431855967204e-01 100 1.04309612076881122e-01 100 1.00246477789047023e-01 100 9.63416139721800652e-02 100 9.25888554956175447e-02 100 8.89822773838300912e-02 100 8.55161854611709293e-02 100 8.21851073610761068e-02 100 7.89837838855261404e-02 100 7.59071607011192057e-02 100 7.29503803586381838e-02 100 7.01087746235089238e-02 100 6.73778571050381719e-02 100 6.47533161727949780e-02 100 6.22310081489358957e-02 100 5.98069507657364163e-02 100 5.74773168779894733e-02 100 5.52384284203424891e-02 100 5.30867506000306175e-02 100 5.10188863158367312e-02 100 4.90315707944623119e-02 100 4.71216664358406687e-02 100 4.52861578592531691e-02 100 4.35221471424239587e-02 100 4.18268492460759922e-02 100 4.01975876167232557e-02 100 3.86317899607558635e-02 100 3.71269841831440434e-02 100 3.56807944843488614e-02 100 3.42909376092767615e-02 14 kinetic-energy 1 c 12 0 0 14 4 time 51 0 100 0.00000000000000000e+00 100 1.00000000000000002e-02 100 2.00000000000000004e-02 100 2.99999999999999989e-02 100 4.00000000000000008e-02 100 5.00000000000000028e-02 100 6.00000000000000047e-02 100 7.00000000000000067e-02 100 8.00000000000000017e-02 100 8.99999999999999967e-02 100 9.99999999999999917e-02 100 1.09999999999999987e-01 100 1.19999999999999982e-01 100 1.29999999999999977e-01 100 1.39999999999999986e-01 100 1.49999999999999994e-01 100 1.60000000000000003e-01 100 1.70000000000000012e-01 100 1.80000000000000021e-01 100 1.90000000000000030e-01 100 2.00000000000000039e-01 100 2.10000000000000048e-01 100 2.20000000000000057e-01 100 2.30000000000000066e-01 100 2.40000000000000074e-01 100 2.50000000000000056e-01 100 2.60000000000000064e-01 100 2.70000000000000073e-01 100 2.80000000000000082e-01 100 2.90000000000000091e-01 100 3.00000000000000100e-01 100 3.10000000000000109e-01 100 3.20000000000000118e-01 100 3.30000000000000127e-01 100 3.40000000000000135e-01 100 3.50000000000000144e-01 100 3.60000000000000153e-01 100 3.70000000000000162e-01 100 3.80000000000000171e-01 100 3.90000000000000180e-01 100 4.00000000000000189e-01 100 4.10000000000000198e-01 100 4.20000000000000207e-01 100 4.30000000000000215e-01 100 4.40000000000000224e-01 100 4.50000000000000233e-01 100 4.60000000000000242e-01 100 4.70000000000000251e-01 100 4.80000000000000260e-01 100 4.90000000000000269e-01 100 5.00000000000000000e-01 4 time 1 c 12 0 0 14 0 0 0 0 0 0 0 0 0 0 0 diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.pvdhandler b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.pvdhandler new file mode 100644 index 0000000000..61fef595f1 --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.pvdhandler @@ -0,0 +1,2 @@ +0 +Time File diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.simulationcontrol b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.simulationcontrol new file mode 100644 index 0000000000..d00487ca31 --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.simulationcontrol @@ -0,0 +1,8 @@ +Simulation control +dt_0 0.01 +dt_1 0.01 +dt_2 0.01 +dt_3 0.01 +CFL 0.0166875 +Time 0.5 +Iter 50 diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation new file mode 100644 index 0000000000000000000000000000000000000000..c6fdeaf63b3be00ed2d39e3746d802307e0e8d82 GIT binary patch literal 12544 zcma*rVM~yA6o&D~YT-g6BqBuQAtWS3L_|cwZ?pG(3jPwm=QWrgpYwnG`tj-4e&)=Z;qlYYeD?1PJ#$_rkBht;JvjecM$4I4_gO z3%Ads2j{;GkL!`C2j^w-xOcz|y#qLJhCJ>akm()3d6_(}#|-u0yczPi9+`S@UM7zh zZePb9aQl+k+cu4><1z^0+<7vVz5 zKs`9`2J*NbnR;+uCXd^L8QKHRn<0o`93%=lgImF-vIUCyl;R!u1BUGoR`Vti^b;uF*t9AJiZ*6tHF7hJg&zK z_29f2^0*$EdT?GQkN3sC0onu3`v%D4_8`+9a9$>l_r<;e>cM&60C`-GOg%U+lgI7B z4DA8u&5+0KL8d+6yi6X~V}^Qg-VAwMk4!x{FO$dDi`~!7;Jh2iWF80SW%9TlGt`6gX2|1uWa`0rnLKU}ZlFEj zyc@{l_8`+9a9$>l>v041;Jh2i<9cN3!FicHZVzT?4>)gzJZ=v%?E&Xy^0*!|)PwV8 z$m4os>cM%LJl+@k2Iw8YdEWqe+&duCJAm^tdAu+74Nwox`v%D4dSvRsd6_)!9WX=h z0M45sk9!AXdIxY`CXeeeLp?ZehCHrErXHM^$>V*oZ-DlI^S%M{xIM_U2b`D5<9)Gj zfO>G=H$WcOBU2B~%j9u;FhhI5c{Aj3dyr`lI4_gO^_Za^oHs)r*CSI8&dcO+?|>WV z9l&`vkjK3PGQ9&hFO$dhxPf|b-VNk&Ju>y+yi6YV4w#{L0O!q+$GrnGy#qKelgIU# zp&pz!Lmt;7QxDF|&KA$Mu+@9-KEr9@isN56;Wv@xpB%zYl=(--gHUM*jWayi6W1+|EZ2 z&VL^s*CSI8&dcQS;bJp?7@Rjl9v_X&$H94-Jg&zK_29f2^0*$EdT?GQj~8wiV-Gn0 zV|d&iWZDDH%jEIG?Q-lKP`4Y$Aj~3Adf$f%$LD=nLMt?4b+44ZXl2Ak*NpgW%BszV)LI2&YK~R zPel>oG$;IB$kLu1BUGoR`Vt_TUEE1J1jFJZ=v%?E&Xy^0*#1P!G<#fjq89 zrXHM^$>a85hW3E-X2|3AAk!XjUM7$0F+)8#Z-zXsN2VT}m&xOW+x2({aQZJuGxQGNyczPicR;3h0Ow`$xE?dqgY#y{<9cN3 z!FicHUbx+lJ>dMq@VGt5vW9F&(VYPkHh17Wa`0rnLKU}W@ryMZ-zW>4>Iim z=VkJ^9y8Q~^Jd87dSvRsd6_)!9dHA^132#n^0;?Ergs46W%9TlH&74GyMa8eN2VT} zm&xPa0WcM%LJZ=wepgrKc8_47KAk!Xj zUM7$0aRc?>yc@{ldSvRsd6_(J4`yf&IB$kLZVxi;0q14%xE?dqgY#y{<9cN3!FicH z-WU6R==T9|-ZwxVpO4JHuI0cc{Akk<;YwO z&dcO+J!Yr}=gpAE^~ltN^D=q7FZK=49&p|_KpwXTnf8G5GI_i&_6<-E&ie+)<9cN3 z!FicHZVzT?4>)gzJZ=v%?E&Xy^0*!|)PwV8$m4os>cM%LJicD+er^Wm-9R4Sj?CTQ zyi6X~;|A)%c{h;9^~ltN^D=pSzu5c_gY#y{k28-9R4KBU2B~%j9u;FhhI5c{Aj3dyr`lI4_gO z^_Za^oHs)r*CSI8&dcQSzSuWF?*PvG2FT;y0h!(boR`VteX(zVdT`!1Kpxj4QxDF| zl+k+cu4><1z^0+<7vVz5Ks`9`2J*NbnR;+u pCXd^L8QKHRn<04)y>5 literal 0 HcmV?d00001 diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation.info b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation.info new file mode 100644 index 0000000000..acc59ac17e --- /dev/null +++ b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation.info @@ -0,0 +1,2 @@ +version nproc n_attached_fixed_size_objs n_attached_variable_size_objs n_coarse_cells +5 1 1 0 1 diff --git a/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation_fixed.data b/applications_tests/lethe-fluid-matrix-free/tgv_restart_bdf1/restart.triangulation_fixed.data new file mode 100644 index 0000000000000000000000000000000000000000..290ccf3d383ccd63fbb95a8c696888a694e2bae4 GIT binary patch literal 397320 zcmcG1d05R|_x_PN846M65+PH{kUG1Rp^RmY5{g1(%8-(I44E>Nu?(eoQk=6Jh>8Z9 z6bVI{HIYicwa?S}{`Iau_I3R}&mZq|J-FCI{;a*uT5I2TE0IXFzXkjU`>SNl zm1D|pj6}1nJWk*EZM5V-`N>Jw782VLL+3BFNECP9s*oGv7eV{Ge!L!??LzM8_sfTW zCbNfA*I)48KYpFPj{W=d7dpj>Cqhue0ydEiC_kAEuhlWV)zMj}0% zpHUls8x0Ijm{nT4jOPCP5PrWgQLOrFm(R?`2)aDD;kkZ=1#wV1a}oZTbc#8WHuqnm zxW$`lH9h$E{Z{qZ*X~*bIo!9+%EYS?)Zuy6a1|G4BI=*};7?_tk*IdqQcPw#}IY|7}oKnPv1ZJnJEq*U$Jn+VuVjh!k7Uj&yP8Cuox08S zd$k9}`oRryZ{M+S;%~c2!{2*N0PPs{ z-F6$nu(lufZtu;qaPmhzFsJ>EBs$l&i_NdKdt_(8%WL0NA3ngDkjK`G)X)CP6!#gM zS{evvvp&H5xJ`K?*%_&@K1es5N=;QZ?}IZCIzBJ2&F-nsfiod>n(C9cz!^EG2I;V~ zF@ET6wb(3)G`=pMCOQ^Q)7!OQ-4)J&m)EDg-*195A(d~}(i}LWO%o109}H(xb0+QV z@Eb|wtSf_Hs^})EYCZ!jhhV(Mus*m#(R)L20x4^s16!nD1#s^#k(@vV%=amjnq=)IE%<9}WBA=ufx57?d&QWHYf4(WDj zPE#tqH|72-+j*CW(DAXH({=pNk@KT+#cdnLcKh5uozBW}(CX6W9Kp`!he4pm{j;IO zH6={#)&4ZP!hf9Fv)4WZ%Q+{VR*Rg@<%wVH{66*B*>t)jH|EmWMW+dtb9OxaYqX&)a&t#bFwqQu+BJBWtZWdCw>`oqPtab7FA9QUVHq`4T52<+VD$t z<-9<0!6r)0@>mvK{wF9oxwStLIzEQ=N~0Ffp19?UA157c{aY`ae*Up+XL7)Gf??g$ zCVR7gQ6RbVbo0^aC0R5ubEvU~@=bza-RAtNT@yAGh+Q1dv{>wwO`mGKe$)BlHG*ON zIQMb5uR{>A3TtrvJ2so%{Cwh@cGnxSGvMVlmUBEZhj&jtkVx5$zjvN$N_V+2>{xEj<^2)2+yvA})w==Q%ninWBn6}PP z&n=by%uF9-7~>~9172QZSkG12ta0c?9MQgVac=KnnRILKxiO~Yx4FL?-uS${#<0Fz zzvarXh4G}wy~so_J(HT$>h${6>ozy9F|4mIyR&OYaXhKc&onX_mqq<=UHn!)^_J`m zc<*7_6Rejux_7_{FE_eFvGpd0L-IMZ`<$%dv>5)!_x(!M%`*6sUNyBzKVq2u$` zIS=2)b%*YYVYU!|t-4iOa&Mz3`=EAsE(EI_>gZ;z*Sp`y`#dk(Lfje2~_`JM6G-+8IxTA^Rkvon?yYj{HU%NQ7yLxu) z`Uln}S)_W-weF(}18LuhI`6dLPUhvc-%%$8xTDGVHa?bNaA)ryWSqtBYL`2alNHBh zlk;urr>WZq(PwwQJke>slX-cKVJ&U_V6xS)Oft`LLu6rg96eFEpkNx@2SUf^(VnxzT+<4?F-^*PpR&STDT8*d5vN1@x*D9?Kt3QwL!Y{O+2k%d_CrS^L@Zu zhhRA;>SNW6nTjby>rDNN7K$k}&{n&3n>V+J(D8Ze94zNtKl9@BKDQJ?&WvlX@0LQF zyceA*ZhedJ);U2_k)Sh@pHE&f<>IS|f(m)BU%DZBYuea@)o#Pa%{Ri`|Y zDSdHH&9V6mczKOseK?5@9A5N<__TObSW@_e_Pd;vrEL*GgpSY4YYgi}cY;h~=RYO? zBDcL;J^v{+R`ZIGXhjfSUSnAAntUfO=EYMo$hDiH{fnox&N%$n08xbOJHUDpUv=!) z|IW`5YwYT^A~W|T9h+f4rvLuC1nWuMzhNFxAL>tnpF6c5x8)^uZn<4!(Yrt*bbPEQ z5fMK0Yh7rr_(%Wv4dD}W=+xLF{dPaSL$IENb);Q@iFyF}n9%)>X2DB(&FbZ@;xT~) z>q+c4Ni>f;l_y?t!!Ps7%^Z5lab{R<>pKMNN%XU;zi{_p013ZccSu>BLx&zXz3|8G zyRtLDu+F`3NPAjYrugV8>B%9Ad33t3M$!6tp#;M^pzrjxe|Gzlue}3r#rfq@mnpet zUoHwGLdVCjeqj1A`HfGO_j;Vz%dqsW8Z|27)VeL?$Ae$o1MatQZ#!yBKM*H}+NWAuqd z{UV}B&H`t@rQt8=(8_Jj;o=Z(USmCpuJiP-ukRgAGzK;hJLQ*D%V+(W$f#i18SwHN z!+M8lq002CAkxaEuD{K_@28QubfV-zK&R&W0LwYM)=wMRWRfV}QL=t@m_t65 zwD=ge@p3r9a*m6B=ZC*Ld6669P7BUd=26v$L;rS!GY~pHmUHx91;0+LP86rNT{pk@ zO+MB2w@UuCA)H`2CubPZC;ncfaCz{XjgI;Bre^8=>u?5G&iPB)8%jfx#LqiKpW3>p zfDT_G3i~}NoM1U;Nz=c;1BNHa>6d;lgG=(MdH%|$t(wmO!&+qQ>uCP(g=BSzXFlDuebn+)CXR>eOu!h=nm#5umC7o!(TjyX{$9uJn zSO~T0w(OskHc*>3k~HrZTG6sx#LH`}CqX2--yNW5_sBS7`yc4pt?%Kwzaci72pyl7 z*H}-YvLRttF!byi6q0W$LC@~^nTyd{@aN~{HP(}u_Hm`}2k6;N7;2x_3VL=IPF~z_ zGUmTC;N>-jwO%W$$Br;h9y&K|PfwU9N2>gIJM7whB6NISUSnAMhH7_If_ZXs+@d+7 zV4gf9c73#mx@+289cdgZI|JT21jAaC@KLg*C{Z%GOU%+F=<$_~e*JRL zi&!FbeBL?-!&;PZDL-jNvgD8@{d5z0e6vUNtM2zHmhjd&IMzRURX0461pO|dr=Z7I zv$)%jiH-lA0WYtyoMUlnn)kQ4fuz^vof9o!$18pYabD{kPlS%o%WEv>cxhDV4TT*q zr+yx5%wWgs{n1_FWl8abm)BU%8D04Ka39$5Quy;or3!Yu_MKNh8(14J`wn<{jbSZH zXw);YeI{A=Q*Egk?B!Wr5Z!BQ{D24@pO@Dd)*`zq&u2ZJOVZsElFMN)FK^@NuyflV z5MExRSohyK#k%CVWb3eQx)!jPmz<%K6dUOth(A@muJqCd<_AK@=iQTdd+CG=%oE6y-H#{l2+I{0kG|G} z%@y|Ov`!jb@RE$P@Cf?cGk}^eS2VRMC= zr4#k~i*v}~5sDAZ4g}Dwp)C_@n&$_+dlDGd@vBztnm_jiIa!jgGtE4T*31|`zYKcU zLdWOblfbavV;2+T8vcSP^m}@3M|c$7Y_)av1L$4z?nz)+J9P|O>ZbgXjPdl*P3jd* z^~)5MrOmx-US50d*?9_PIz)VV!}RGNGsUTWZ420JCt-5>tQ5aoa_-CM#wK@P`Yt&9 zKt9ZHgpSY4>+nZ!*-VEtc+A>lcOgqWVO`!vHrq+koO?1uDUaM&8Kg0V75oS8%^xVvq=4RRA=YJ+%WV4;%Z3{J&J@Uxd6DGkv0e)1n?yaY7^9+ZV z*BI6sUM}0>mgbNdmUreatPG-#b%{BBpjRPud|qB-ShpYQl04^e4!L)JLsrS2dsI5( z^nwWJRq*l}!`imH!ORw;a!HGAHm=jV1yi>lRl2Ipy$aqv2`uLvelze}FOMu@n(%Dt z&CL=T*AQAi2+lz0_`G`(Sk6gTnB4Z|i!5^AXw0W_MJb){F+BePoB{8i1eSAd=I@?n zG$fk@6vizd7a*nKS;xBcZ$1OwJqZl!OCMia64-6A486LkAMCdDNs=6Fg5UZ=$LHOX zz_30aGz9Pc3Yf$CVtX}KR@rD1cr5YzlFywVYfv`$#}XC?6&OD-~IA* zbFYGz*I3S3)oJ+KMQ|r~mCj4Bg*$np&de+Gp=T#_d|qB-IVXD7xiLH8PQF}f-QOGT zxDQC~Md^jlOlW zqyhF#51L=;s0;lfq2u%N8pAqS*K_n%*f&iwQ+awB_Dx?$4>vmi{UTmoV^}}CdsuNL z?3*e+s>uHU`=%}vTeYxi?icaaIqpwdIl_FIq}jh}Q!^=1EQ;?N&*sd3_qVpcTA4?B znT*QKRq~?WPJ~3OzJgwNMB)%k{B)=UVHrC(- zo$r09k7DzEz+2~FSnu9jeK{R!(?zitsR5J%)KQFRwAI?L4k5E`#50 zqII^zcSkfB#h+TNx4U_s%v}i*l2nSuaovps}h=r8{AMz+2~FIme{Lyl@Egi#9J$PK$zm(W&Y+iwv6Q z2fTF-hPCdp`H!67E!40XCMnzJ2GTXX9qTtk?^@{iymbzS^~Y#K0|$5uHDcVL;4pX# zHBh;8%6I5p^VT^S)>HppAKM+?Lf!qS+x05Eg<9vbU|dvl@0yp_Sk5s?_5D)}{Z!GU zx(*-KDTpIt?sXpxGaRAg^YR+YIp1&B1SLX0RW!)y+d=ZI(6 z-4BI+sz_(V1*Mq^;>eu`+^06raCmu*VO_W=b>Yh5EF#)+`QeBjZN#fRcBm=C{6OgV zyu8M+e%#Sh)LuUuv~tvsx~jDiH(HoPT!r}oFRwAI%Z8OaTnldnMRRgA6xG{^pD7Go zp42=);N6qJdJ;bCW^}fI_h_xHFZ#5D_h>_UHJr)UiX%eD=iQUQdJy5?D{7M|R(D_3$3;W5EQ)Xn2n{yj?^4*3NOV84mBB1cvq6 zzYDw6EtQh#XD=A_1+9b@>V{s5ogNaQfR62pyl7*H}+tL~%TmJZGS}P?&$LHlWhIQGc zz{)P5UzDCC(k}u1q8s6nPv_4E-Xb{!<28o$obVN*Eudf2eoMkM6VNXj_WIzYWA~}- zJK)`uz;e#fp6!;Hfqqe9-&Fr8pkMT$?&+=>O%ft>eBM0?Ea%)au5dO7{i3Q!?bHm= zFM9MUF?h*5DdF9dz|JPpTqT%*e$n~IfoZKkzvz0XPS~Jesq74R_arc^MJ1o~_wCM= z`1CSZ69}492b1HxE8u@a==i*Q5*XGZ&!{nPcjihaj9g-s1)5XuJ#vpJDN@3_CxK#Z zyl!BeTdt(ELG5rjXigO>I;H5X`R_a67Piu-~__UYb@vdobcndJ7@~Jg?Ikm1~diB zJu5S}RV2v1172QZSc`s6ch)`fMlz}6_Hl}!135Y(>61!&0ueesFRwAIMWgiI2W!2N z^sp@ZnhiRTi(6{$pU^Im@bVhPIx8?+=kaSv%zlZb8|Xl8?Q3k@G&NCn2E26+){_X^ zr1HxX^o!gp8}su)zv!mp*=1I69|#?vx6Z+O5_T4fn>rjwBkrqz?JEHNqT{8NUzfvu zz+2~FJ&8Ul*UL>ozo@tk{nn~i7#(lxY&pF7KH#l$Fs#)!rZ}wvZOxz-vo&)-Tl4Ox z8M7ZlE)qIEZ=Hi--R*C$^D97G^TgSHAwxl1)A@K}j}?%MckYRPCXEAaO{d)! zRZ~D)Q^);~`HtpX#LH{6CvjIC@d5O`MElpJIgV*3-aGd3s-JMqLdWOjHP(}uc9eYW zpqNfX3Wq-{i`$9+^&J1KE8NMvyvBMG@28d@`vLl1qCQ_HDNJoAcDWq8-LUyi=H+!T zg7x2q8QsQ$_N>VM`iI|^?ZnbK0g>I{J`g%SFRwAIkCz9YY69(9(X<}b2}|3F%~~0+ zn+W#-FRwAI2km{A5eV9|qVHYy2RpPAJ1(-C(qh_w_W^I6gXJ8yr-Z}@x+EZtb)y0pJk8i>&6G5YV|KB0{-IVHI{R(Wj~pl<_B+{@A)p#uxTq+T72LRI|HHP^YR+YIjhe9Bx@Ss&u`r; z*Jf2)@rXiKgGx99US4B4M`=l#Z4Kzbi`M)OjdX4+p66@rq1Su{yu8M+7L6L_W4ijY zL}a|A<$%$x#c`hxmA2`XM1+pd%WDj4QNoIBpMMn+(Sof~Fa6fyK~HwANOwpgyu3!S zPT07+$2o`}xM@${@vX&Hug66h@fV`!2zRxeYgS8;v!q#&kOArk#>O@z6BRYE|o zf-31o$@D7R;=ZNG^(ru|Z|t%(aw*G{8qcbc=@*r>7@N*^+%c?;x3^gb`bD(P@|8@# z$RzaVX{KK!bbJhJ=S^C1)jnC$mc19t^o!oQFSTbo?iki)>?mz?YLuDzwB$C0{TT%x-?d%Ut}$wb6Kum#La6& z6Jg=kPzrh#^sE1KnO?;J-Q-hDuR`$nh}H*&HKvJRy7b;o(5s*?Z}*hxRXlQ6mFZRR z^ZI{lA`E#Abt9aC+RI{@Ud72H9BVE?nMW$z`*4D2tdv}KAoEv9HF+DqaYeSe!&(6Z!s}Iw&6FNSYa~@vw zeNrEqE4}B`nY}-w2C*X)*}F6BYzB9~36T zdUka8%y^leU8;GvHB8Trm)BU%d7?j?>DkfQUCfP{o?WeJ%}KeQ9WSpjte+VfGXDc= z=R1EM^FN5*?HA1a4}^}-%WDkln<|r;{{fv|b$$}_KX}?OPKEg&@bVhNdc?!s%>RJe zOdc*~{s+^>t}K=NA7I)~*xie1A~^N=#Plj?u3wW(ucGsd$t#&&h2ZfK?I-N+#WWFm zOdY-s^eSlh{RuL?in!%&=}fOeB;u`e{nqc%P|5FdO|DmgVU6V+ z4C@oQU5r4#h>joHW){;gax+~%fawm)989Gj~auei6O?e%4;5Uo^+f`90Gw5;{JHHO6ZU>l~*GOuvYp@_QrG zFRD>hxX<*9czKOsJ^t$aAkZ(OmGi#J^oy!Lm3Nfu7hyREyL+*mbK9wb>Df`+i;rb` zc7Bp$MoiC6==fO9!R}t{Y~B_?}TD>!v|AEl)d3lXt-Eyuk^FN@wR?L3S z{11Lle-+C74|sWvVZBlHBJ)3><8FPC`5!#}`s2IY{{ZbtyvsE#EXt8S-g|dE^UJi{ z9Cn!bpkY0U?InwzgI^}S6Y%X4^UFLprZR!~WeOc1>q)fgG;P)Js9dSZis8%;jm{d= zwX!q#pkY0UmAx+x0KZInZub#&=9k$u;C53w_@iMxiD2=Eai-_;q#j;n)0kgohpj~p z%m)qYNxb>6gZX9B#Un#@GQUida~3V+ewjBAto14mKUwxOQ@Zfb9p=?b+v;o#U>>~~ z)^pbRM}U7aZFz40VdkGaCH&fc=ASHdd<^T0x?S@M{j#Jp9+k{y{>dW~Lq7Y0M=yqT z=SQc5z(1K@U#t9w`6p+r)>y*)lQFE_<_D@7T4YOYZ{KBJy;L*4l^OHs#jt+$)`0mZ zQ_tAjxy(PgrMpKrxqmV*ud%y2D7QcJ%cNuUZw+RCnJuSE1~I=(q2u%N8oR6W>yDL! zUnbQhjWWN?y$`x6EeDS}US4B8i8Yr7Grvq~_S?jW`DL07Rc6OuOl}9MAld4X#?h?GFCQyu8M+)*N8Q z{F7;B=A|O$pM0bD=9O~)WGv^Dj=pU8(mYX`zuYE21a?)TY;S#y%61`G&UrRcd|wQE zMRxTM?EacW(9U_}@X`lvMCkZf&QYJZFmPOBqO{>&-o=@)t8#zG;;$nsEC`l!F7}wY zx81b}I#U@GH^a$9S54$SK%HJ=3Fm#|;&bfTL zi*gj~747qyuCfOHu4_Adh4#QPvfUO8Ymr@J`^&>ij6^y|ey!?NCD}9P$1dY*3kim` z;~f1wc2=`2OD(pUCyFnhrG@OwgpQA4EfV=2GTHgoNOacKJ7QLqq~3n>&X2XrD2DZn zt05QISq*fw`d!eFD4uF{W-&W63~Q06(F4O1#YUn#O;(cKRg&Vb-m01#Tq%b2F-!M7 z?5te3gnr%~n(vw^$VI(@H85bXRaC3YXb5AI~4ggi5n*(?L8qyW+WL{olImfR1oU_B>F3xK8JwqMt<}Q)eqN=rf zWZwZVuQ9Bxt7psVr-q4k2Ni#2ir4v%_h)sK(D8YBjbVLa{UKTX)FpR8r%=Bvv8Rf2 z1*@ZYd5vM6>pxpoKlOGu8Q$3RpdjX# zx#Li3f4N^KZ=Hi-U4F3@^G~K%{JdW<|KwX0zgsZ>WTE5p);So~4s)h3|71GP(W)2o zPyXIxSff4Ero43yhIMXyMjH4h(}wzSmdro-V@t;|a{pvrUSmCp33YRrUnbqva;!S@ z%N##wh8y$C6goaHud$xQka@n$FOwQopE%6?G7s!ZJjeVpd3lZXB;4)&nO`P#nYMK` z^UG|xepXkxUnVcFF|1Eanauo?>BkPS(ab;jV1Z{g^G_B!J}<8^tW9h~SAu^s-D#Pm z$NZDuI_V!`{>i+&#;~?-t;772Y0JkR&zXO+IN2>!?w`zC=U_SKKd#d>@~$LFnc zFs$`rlx4lW3{5M=#>iaptL!O%S&xsm&cU!w>HSpJ+uO9`d%p6yJn>H*i}kF>$Lb

*Z_JJbb`V*nO?=en`z-puR_GTCxK`p3~ESXdKFZ4Pn}G!Vusy`#IT;+y`Je8QIW<#G1D&!Y|8Lv`b9#=$FRm`mKfHjOqM?b{UZAH!1z3-U!+pq zdLq*=;@y+Lunz5*%=C+B(cT#zOus0o?DToLei1jX5lw{XLFbAuEJ_K(B(<$FG&? zRd~gkzm)4$@bVhN`hLe@4WM5{m(-TX^otsIsJAu){UV{`V_0Ln#;|spt6l~AMYMk9 z8JT|3l-;VCOuvYi*BI6ns@Y7xh^G0@HD~%oC%Pnk(cv$x zE-*d2i)(AoF+Dq><6}7oo84n)vuxLSre{Zczq_W+^z3R|t=q!%?0EMiu$RJ4RNudu`5*KhNjg z&HNAOUgv>>nE!!Z_|K2b|A2Q-0>irZt)PZT6FNRGud$r-V7_=a=-JUW|CIJHJ-fmaVM$ESj+fV1&UtpJ zf$7=NyFW|@Gd;U~O%e0udUm|L#<1Q}w1N2_Q2R>X_ssu5v!{6{=6@h`d|qB-Sm&OV zF#iMk`HoT=^FKHj`B~wVpCG5-T<+_CHx^FR2qKIi}Z4>0W~?C!-h5hlzy zCI-C<`hMw{a;8_Y$w#{_)2k3XKBE1E-MyG5LQG%fY|yKqC3(6Bm|jI_tHrmNUWJIa z&OtO0tov9E0=)`q^~EQP=~e7&IiAS%Dln|EoP%KiLg!I*$(t7s3@f_=@Z1){&c`ou2;d!YYgkwT79~Mei0RQ zDe{f^-C#6m;l@~|UnF#V3~P+n7}n7jw|Rqp5fvS-c{2A)gOTmep)O3nh?mzG*5b&( zOwcc)A|;nG?lC<(+WVkW5!17i#!s?mdUis`$8rvK z_hM&r|HyTwXGiyJT7HP>*=e86_{H?>cf? zYpl+}u-;nj$NUed>1d}k=6?{oNUf3i9|#>E!y2n|Fs!3r?PLB2wBf2%IP*WK%g=kl z{115R91QE?;7`o|fR=l#?7{pGd|v5KkozC-@*2xI!{4h%f}R}}JynX+Z>lzGUDC;! z>DdV#pO@EI&RKo9=Qq%^qar&i{Y@3sM!p?mUobs8US4B4XJ$+3anQ4)qCYW{pJZ1X zEx+QtTdrrv%WDkl-=0Q`!2f`XPSu{6-M`$(&pE9R^FI(eJ}<8^tXIcWGXDek^AGJj zwdj*k@3`m~=6}G;YYb~+OG^p(A5hV>@jDvMd@?$tpX?y_Kfrnt*z6wbNq9e2W`3FU z#i(z#%rA4*=JYz|mnn37tS5oZ?y;T(aT&$@GU?QDb1Rr%W{Z_mbD3Wz@16wKld#Tm zVt$!)PS;l}nO~;mQU6YIzf25kY-WjJt@e{J|73dUr{zfIpS=0r+}X@OS?Kr}*4WGv z!`k_?7xPc1z8NMrnSb&DrJd%?Kbdz=0>e7s(gEh5Onvg^th5YwWJ}8|cscGHK7hw=|ev=F&6EG?`x}FR!tl zgzZ1V{4!~p$G|}5m-+NxmPqcG$;)dDYt5D$nSU}ZrjM+dfAUeAWiOb2ve5B)d5vLh zdU!qaPo`B9ymFa;^06sH=Q00eUS4BZCoLVz{FCX7AU8MWpM3V?!BKMmWGv@kvwJM( zBz4<+-fYWD>hY&cH6p~Ho|(2U!v^$=gpQBp9Bg)v<(zg_X^Eo?UeX@Vbi+UQ44~!H zckC?%{UY8y2`uOAtZ5w6A9Rb<6r_n};E$GZ%KBl=?z^&?K87{6CxK!8{a~z2Ys&Vh z#N)nKMMC3JiYYiv&f!@B+9`7*7ktjQy81cl~`cPjVq$~32V_arc^6<3VU zx^gpzZY8yU)}G1}2fsgZlW9(|wF+KdV>u__YtpO~zg(KC=G#sLo~55_P50-6rl8RA zd3lZHoTvrUPG%_O(ffN(51$8~r9Y-POf~~eL0(>CIp>a^g@&?69(@@da^EMwj}Gu| zbhk|pm7M`EuQ99_n{PTbL@|%*Du28bU79JLy?9*?(}5H^J}<8^tc~B~$n+nV_G~ga z=aVIV**(vh=|J-G8pHZ@s|AS#-}9(<7w;N((`>Ovui@vI4kY^ycDa2ZnanRU=JeqI^ULI|b1bJY{f<89 zR%mU+s;I*ZPxUm>&X}AD8{w$0Ske=f5bP<($_? zqE^CvK=Rx*zQTP#x~FZ4WcPv4@jqx9FzkrP@cx+@Sv5w1898qlc738-GhP#s@5tbj z>9t1NRSphkyspuTVtxo@;pT&s0Cpdcdy|*Ug!_Q3jk)!Z-3Nzb zD=q+bM6|S1iSMBrqplhkJ~8YjO?fy1cuhoSC&pgbP;2yR?*3rL>x=IinI8ff*@})> z8ImN`Sh6*Y<(!np&u`&AAZG?Sr@(zcvTgNx$?pRUYs|A9!+LJ%-OKR%fW%#DYO%w^ zmspmjQ}+8n==d1cm}fhNwS8?^H~4)(eq2fG2)_@AXL6k%`+b06jd`|XSpOS)@ICxK zAUhr0BjEP|QT=_oTK@Zhm)D>2%9tkw2~OFrR9~7YJ<%Pd6flOWemj&h0zkBTXVbGERaglebA6VNXrS@8o{zlaQew0{!o z7xD6XMQJ+oq#zwnPd^m^9u+Cy(l4`~#Cp%@0njfZ1`|88ei6B7({;PNU&PC64C@@5 zx2`ZhAUT82FNOI5d9CH1$>s+_$LHlWhPC3fH%Tx*AeIBB=fnJfykDei$>s;VyvDE= z?+KUz^8@01T+1Ei2gG;s_TlpR0hV)2Qhoox?_MIBRM+9dIt6h=%)RdHcQ2N6PGp#} zei7-t=fOGX7mwVmW7ewjS#j z5f#(?mC!FD)33~Uw-fZPv795GVaI;=5|Pe`3raH;#F0A>xU=8ASk4(arWNZKk-QC$ zUP8Z!Ow_$K&%iiV_WJcPJLm&i1ZP^8zbK@;^j4# zbBZtDDTaA6IT9N52t3XEskJbIk4ygn2S~bL@OyHc#f| zHI{R~ygwijP0!2V{Q>Ef+OHjZ ze;{;xUS4BZd!4v43f>=(mVpM(;r#(Q=XzFr0Q9bTd5vM+n61j*9}x5S&TruT0a=(o zyWaBkec2iC);XJ>4q!7&($_3^$1Cv5*+pGPvzevV%EEInPbQ)+JJ~#$R1F{Ugw2zM zj?Y`?Y;;Z70RA}Se6j0f8BRrnh5QVbY^-LWRu=XrdL5i zG^G8QUWL%{d3jCmtz{lXWUlq7dy(K-q*FI-3)4hs&{`f3dKKhPvXeUKRgl?QU$dEB z1uw5_pKWFyMP!s>k8o}9EUK=&p3gK9mTBKK1HB5e<9qN-(5oOvcg(#n*Q?;=HHP)L z_b*IAzli*Lx;X~)i-?OWdBXIIgpSY4YYgkWTLvpYzlek`chv{|B4Roud>7L%;^j4l zb)Q`>8$rK_Z2mJO1N4iClV60tT)&97&cSky;u)u}uy0Dne9yImeN*Br@vmb0rb5T( zt#h!P<9h8-8`w7`TSn`Dfqhdl`*2ML+c)K{bFiHA-Zo(>?3xi+iZpR$>fn(b$1iI?-e>eZ=Hi-efNy&5_q3X-Uj&(hWE*& zx25A-_CA@n&cU$OeA(j|yiX>_v~J9V_sL}2l-=p4WB#j6d3lZHoR0agTY#P&an|nL z3-s*B!&Gf&re`N~d|qB-ImgFds~+_1$eWOjM?lYx{B$)u#PsZVd5z_qwwH&hfu0@7 z+dj_`^z6vWD*IM)Jv&}rV_4@W&(H+_1LCo#%^UDPApMjEr!fBmq2u%N8pC>1%dwNd z|A3s_9s2?N56Bp;pxey<~XLEc<BC85{NMG;*-p*U#)d8rG9|KehZAysIIi zK3^s&Ol>E2xg5Kly{o}`5}kfcEr<7L-6MCzBlInRW0!ne3H}GMB$k=H)fklc*j3Cgkd*qsc&1mu%WJGBvA?J~9`q_mqu15LULsjP z<@;l%U&PC64C|ivOPV%&C6e-PSEhr05jlS+hRF4cu$*%(`^j8*cTYt3e3xk0v=u8Y zK5&PRbBK$^7UihtMPDlH-~MC`UPi(VJ$*55m2mME+2Du|1PPK zHvLY7v%9fP?dKRcGq&T!%WEv>{QUds73kTKR`3h-^+9@YH-BpD6v_1LczKQGoZF+!UW1+;ng9Ov&o*#3$DLl0Dc7^(ef%%U-U)3D6elGSTmLg3RY0#o==i+6#;_J$h)(~y;g_U(!ImG{SCdGg)0e#y zK(B(A*C^J9LmX{2e@P;=X3TSIN+KO%o5W31K~n|(;;nPAp2V@!^E-i_9ceT*y##u8 zWT>8n7t^y7IzDfmgY_i7b#WOFdUj;Z%2T~T&yHy4d#+`AcD!{C)|1#F9`qIT?8x*V z@06n9@8@UMDaiHgcyUdtJAu9_8FOrrG3c9;F>f*sF?~~^FK;M+u>FxFeeN%F4UDGkSz9}!Sv7UsVLz{2l zmq}iBuyO;xOmejChn>tXQ|S1-yvBMG23cp^z%P>+kjb~eFO!H~JPXo-KR++8v7SVm zCjHCcmq{`%loo(rCRzK@+u>w$&yJVZ7}h=xN0Y%nnM8JQ+7JH8WaAXe0_LAAbbMZ3 zV^}}Gy`ls7ClmWEwN|5Gjxg%GUoi7e=H)ep^|xVt?}C3a32ykUPhgHPt;(RK+&`JO z&cSj{MM~Bp(7Pr>S8lltde@|*&-!gl?^@{iymbzibHdi&UkiHI#CreB^o;#!wZ}o1h0hK2{g;@*2xIek+%q_y)UOy_^yZEnu!3Y}(ObE$ny+9iNxiSkAfd zuZ`}|qCi3yB(F7txpJyO=euRFSnI~V?(WE*{148efn9R!tpa0W<$W6LRBXBOl1lfOk*gh4-j&;HOO1r?=6m z1%Ks>=g#{x&z05At&71Ahm12Vss%qB=${6@WPUh8$LHOXh#wW+75tRRxIqtZTY$gv z;KGlcx`5}(58Z_~zz>IXO}+IE{BX#%o&~1N4~KV8V&JegS+zXvJck4BaSi+Y1U8X;%yeq{b=yn7ND*7H(HD)^(3QC*9fKN=aadO>fw zKN>Hu^%j=C15azRYuKt03E*)ZGcY)jc`{0lnj8hcO!9P8C+3$)%x_NW$^0^fj?c^M z{?eYaz|)$9s=4k^x{@WGWue;64?G#C{5dcO{4z<%qzVi0%Oq`T*R*GTnY_HV|M12I zJgrGitAG(3z~lOBj>#V8$=GQ4;57JUk{dr9m|rGozp3xr%4WYzUS4BZn`@6e3I53> zrmyiC@J}Y!a>n0d{>ehe=jAnqb?AGWVDL{S&c^eYe=<2S_GKFLPv+${hIP~4xDVi; zOt$BKV*bfQbYbNzxqmY6o&=V2bja}O;Fn3dCnV*7UnZ$|Q}UVlWeOdicTWP#IkuPD zd4pdjnY;Ke^UEYVUo70t{4#m>B(R*L?7C|?_+^r*9sV)DOft301ZBBjChwjEhIQOp z%ak?HtJu4J{?Xsit9adE;)(#+Z4o*?@16vP^-aqw`y-%NQT{#lb_eKHbp3hD@D}X0 z@a{=qShv=An_CXO3fs0C&ssvSqM@6i>&i3tWOWfQud$r7B!2YZt~M}-ov};158TZi z9Hy!T!JRC0d|qB-IVa8gOqpVxlz9Hl3@U@3-T5!uf7-#F%*$&m=Uh(F%F%?LUF3LW zu^!yb?>Z-bni%6JI|E)`V^|+6bMAKt`b85vUz^_>`bBHQ1Fx9EzNyghd3lXt{Vl@z z@-XNZZSpNP)q{Rf--ve#zhK{#m)989ntE4v?0|mJq3fAq6X+K?wz!>BKJ}LD40!7t zh5dV8f@eGF=@x#=+#*rhYQp4q%(FeYe>A@i6FhU?>Aohz>NqUpSR8# zbgA+jc(#)Ttw$AvH6}`J)>g*F?o}3kT-?fX2GvTdsJ}bs~ zfM+}DrnY>R&HW_l2(|JM=GmSzU8hqN%#n*+!6f!b8)_`G!vhIPr<;ye4GCY{-H*`Kjcqpnn{8leoeDQ}&F zVIADCYxhm4NvEIumD?U_)X?zmpMEyiro6ls4ISD3+VR#>Q9}RzMrv;)C!mt>_MLQ^TzZRrA413H<+UhJ$>Vth{P~SjgLY}Zm5lI_7JjN-M#n9g zsHWPGC|*5g{6+SiwA)}BJRH6cZeAZ=P^4M_@%wuvs9r0Pj2W}iMRS8IUDY^c{DYVz zarvZ%J?uL%`Jz0c6uuAEv*YD8hIQHPkSfa`nc_hI4=M0j^K5ZfzoJd-H{2J=MHTkjjhz8+ zorC3^1p~~UhQR#biYW2&9OxSvc`vzZ3;iOYl?*y z{*(^=BHlU&%Q+cSiw$;0j3?b>m}vFXse4)?gEqXBQBc>zo-g#S zuiMvd`wPA6#*iJaw41$qd3lZHoVBB#)Gr3_>gNX+Phfr3I{`TdSwB_i_`JNva!#Al z$@W3uUG1#qyc@dNxmN(%#FQ^;K(HM6Q&3So884!+M~y zM)7#?O|I>5fcYqst$S83U~dJ5j?c?$4C_wQex>MwZ}Q^w-ONXsJiNBip1l?1Huqf9n!F?%BSPv+f|z2 zAFy_nHM~>ab~jjC3A7S;_arc^b+j8o7xs=O?GNfy+rT^Ji8o&#e$lLzz{_i_Ct=@9 zbNwec1GAcbn+@R%()v|a7Qq{Fq2u%N8tX~Ss%V@ z{XlbycTWPvx_8#&B^iGt7Z!I&JO!E}lgYx56G3x|{Z8iPHI{QOUpo@(x!aetd{{I; z8g%W(-0(SK#WV$xyBa%7-uS${#&S+SyN#PlmHkLml2-Xc(6#&d@cze+pee}9Yb@vZ zy>8SSf7_23=6rtq&pnUa-?1>tNj*w-1{f9?uQ9AeqW9eoe5sRsIhQ)4A7~CL#q3%9 zts;sD9iNxi7}g@YH=ZWNbrSbfEjkx82fMc1;`p#6itzFp#k%0z0O!m)$%6rIcb|ji zU}!)gWC2^38wkrp0XbBpdYb@%$=Qn$?@bVh#NjyH@W%DlZBH9|Gy^`sCl|3)i zWcpr0$LHlW){`jzHK<@Wco7|0kixu)$kFe|Ix~GQUS4B83FXOQ!;Ha;s7=Bh=1oM% z#)7MIPaUE*_sLHF9pt0JjY^KUfXItRmAG)>bk)b@|$vVHo!K+q%~z2cDL9ndHjIzDfm zgJCVQ57Z?#e@(&Eww zO{U>5bbMZ3V_1tEr*vNnfBs%SD`Z;!&!1bpW*YvyyhgD${;o8xs9K_V%7^(N5XIk} zKg#_NFinJv@!RHuUIi6JxXx|iq#*uad@7J>A7GjY7I`1UpjSape35=-dKDvn*0*7L z6@tf4Sds1XuRqJEL1!Xr2X!Q&%Z9~jn{CW5c_)~%pdK?|;yYcaiwgO;mf zm|g`puMtgzJFVP{L9c@P83)PqDy$7$U&!?;czKOsU4MUb1Lzmg{Bx;;n0`^mI)yf7 zpkE|(d<<)h*BI7y$#zwsUqqv?xXJX3+Bz#_G5sQ5USn9>>vnIk#Ve88Oq(dvFVay} zv5@N*VL4~`d-X`rv!kM?N^$y4)kdvLIyo~vJ1plEf6n;@dUiCm$IB9?XLn)TfmKY; zPU!ep&RKo9=Qq%^qar&i{Y@3sM!p?mUobs8>}+nOKivd+cGTP|;jSs7s>VPFsw08 zMhxq*yB;$C1KKa?%yZ^{5c?soiTNK09UsFQ^JK)ZPWyD7`5(}U)3-)3{{tVbTWQSy z0K*#dWW=!EX3%;Z_#e;(jqd8q|G?F-hoRj6fS1=;&iS$T+&a*+qj#rof5r6d>b`n^ zVtRH$$LHlWmUDC#pVO$$E(-h)=*SGGQs#fqvio8c=6@h`d|qB-STD{Tt_%JLw8}3po%tUG zWvQ4m{{voLV_5GqeZ%|@XnkE!5%WK2+tjU%RkQyAru~H7y_hC~!jL~UpjSbyx-FFH zRrFc*b{EsD5IjDj{e<1Um?nZ$Il~3?D(G;x7@1y0^u2;$rdJ{2t#c4fgvNg!A)r@5 zduoiA=~dW&G!2vMRbW_SIS0e~cF_}M&@ZB&je5xRi+*)f^=A4-LdVCj#&QmZb$;>S z=b&Fid(Rps(=X~9WM;thi+JlC4C~zn`$9m!h#Ef=%k+zuCOthT*DvDcHKK`d`}Nqf zpjSb=NFT`bD(-FE>BRIZgpdEfwLUPcF-?Sh7qiBJUIktETPoA5SY3H3nCVrBczKOz zA}qE~KL&ag^vs9db8pGP_XiN&|7ttrJX2|r5eoqQ`%k+zcj*nrD z@fyQ=rp}R0pkGAS?kbe&7gZT(MKk>(US4BZJ3kWN1N|cU{P<^?e$kWQT1B~j5tehX zyBEtjcdBYSgPt9IuJhp{)3clQX2u1kXD4)gEazZ%FLpK=nH>#4&yG%P93|7UE1vP( zmFd~>);U@cjcItRnLPu@-|@IRn!=X-Bq{s&XO z^*+M<4}^}7VU5)}7}jsE*DeMB1KRV~D4G93%i!<|=6}Fj=U`ahI9z)i{10ee$Vi$0 zfpMl%p4|U{m)BU%DUs%W2R%D#(6^sV&u(az|30Q?Cv<#XUSm0@_TX}N(6gf_i@M76 z>|EC>idUj*N)7Q)O?09*NVeRnIK@0p3sJL6!YvzAo z{Us!m`5y=!pO@Dd*0ZW+nSuWSJ!!d5=6|3%eO);7Kj7sxhIK^Xw0!VCpwdnKWc~-+ zukHOQ_dmdT5&_S<)`4Fp6?wYps$Y{Dg=8D|WPX`gPvT(PD(06-Yg*NKGr!EQGyE!; zU#8IUv7Us5?%I3cmkEFVkiw7J4~*Q_*SunWnOIMv!{_(RFOzQW(8Z4VWqNKs@rn6m zVm*nqi*lG>CKZV~*-X6?YgBt;iYN2S#Cj4wn&HeZlXfp}$Yg$*hVPV@$o(=gtTE4a z4C|I@D}2E}nJQGN{@o7Vz0i+&#_sAKOOCw-zf5YiLO+@LW$qmC z=Nj|N6goaHud%y2c*`&5mr2|27^lhnG820=4rYFtyu8MG5+k1-?*M+8H0oVwAoI&K z=yI@iQS%!RUS4BZU#!X72L8#^^!2gD%s=_f;IaA4KUwJbyu8M+{@|XI4F1WqE`D+j z^H1(?(!KJ2PsIcMeVZ;l?0`Lz3u5k_rqc#+Ebu0?Q{YZkmkv}A!JVx0bM(=C*!dMY zJ}OjqYsUl$X~S){zz;*`6sCZ4JrqW|l1_^)x`fZ^~QeU_FV>;_ctSFO!}){%#rb z%lueyrIGn%3LT%f&cS*T!Q^~f@XMse&SsyPUuJmQA+MNUCU2dC^(1u5GNyrFCcUoq z#Gd(Oil(1ck^5!x);So~$Hza~1^&r2^HSfD%s+X=!^?}9f3nc=dFvbu>)Cx&?ty^|i{F8a>91QEpx7*u*e==SA@P!`pPu}R`)lKf7%*$)6Co%IjX$5|n z^zxhvb>^3ud$^rD^UD-EJ}hJ$Zbr%dlX>eLEaxP&)$%ce zei5nN+ASi)pU?|lg*NaWP3ZW%bq zj?Y`?U|6>>P~He{pr~!kQQzTVxzgC4W?Jkm6mOk_VIB9Vn!SOd!CvY~KThRI^|$ni zk-vrF9Wq2u%N8p}EVmS#P$h8?eX=>@gX zu;Zn%?o_E3Xh`w$8p}B=mUf+@2s>Wx-OFdF!;V+ZBCB+rW(_G`UjIMZ&NQyZulxT8 z8A=ElGiRPfhIGy@b14xrl?<6F$`DD0Oqn7>Br=my5=lB+ArUH#(m;c+j7eqYzt(j- zzeoS|U_ZFNPwvT61C=bQrM7V+ZmheKa@WU1H08M-6CNYvMq)(#8bAh(Fu2;Ey~{sVnh| zW@Os7Mt+gp^(n2XJxfaK;XChH0l$dH=eS7xBG>e}y^&uedQO7UIwEu7HQ*QV$p%Ie zzo>!Ht-Tt4k+81`CPIg`ckcnOf?xajM&eZr&zZ6tc@?tPCs-ep))W(AWzIP-;8pOK z`}<40iUV_fN|9Hg5cM^|MEI{?&tbr;;F$xyOS}rd43j7guR_$RL`OI?x~*h(cWkj@a%X?y>k-Ju0@XFP~_RkU7zYX)ZRVyY_pE=P(7zw<2et2XU7MxeJAnk&YWvfADG7xqMCj0mXiz-n|qP!PjkQBJe8s?`y*(UPY5GE?<#X zA$xs-{Y1TcDJDY6YhQEVRq%rzuOwcDRbG=D$g5C@u5$<`f?mPuDBxA_+G|oJUPZ#3 z&G$9D3QB9LCs11dJzBIJ_(eSB-6x4()TVgR0pu6SU7ymL>N%9wb%wOBHoH+E-{v(& z;ui&m?DYxO2ct~}zUMPosWcyn+abKDB-rF=2_mkzWPiam0n$o&-Qpc?W zU>|vz&mxImRF=Qp1^1IheNAcIZFJFptzjQ|Q1Ls7U)1?j$zjcYGSzdacQ4g*dKi3+ z1fCuDdURFd+0A(sYk)jEx$9FshkEx?&*o&Z?nU6)@o?{%63@=!>DyJvvlCtCP(3GS zOjrQ$?D)(M!zG^GlD})mYj}2)*3>$O(z<9;Knut};Iqo#N%;r1)4xYy{(;=}DXpn> z4yAR&7n@c?{sE7h5F+ItRIO+|0P_z-*Ey8drvf}9A^(6+Z$Ho-^AD1@8otuxABg&z z>N&^r?(ZB9dpDt=~Di|V5_~;asNQ> z`l7z3w6;sHKTU6amimobs+51Q$t-dq?jMNyn$o(eh0>CVP@tvKr;CPialD`(Q^`%)&n=xd<^-? zeEO^Awo-m_@YOz={A5vIQ}1e(e~X@wm&r%Bx+mpjzFE?zJ?3S~U0>AK)VrEjdA)_a zOupNxRS4!~4*ubzm<`WB)YsHZ!enjJ=a84l`%IZ5@g)pLfn`}5KT_(gw|L%uozzv#q@y`N@7KAPP1sh&ga z-BUeh(P`Ur-@PC4jPH#;&H{eX%{jeiB||=%=s5|h=e%w(b4owp7cD=1BP^k1IIDV= zayao{J{qMpbxwlPdTR|M6UgM^+WYKJR=r49FY2Cl0kgT}u1{%Aos*!n-e?fn6f(KE zRguAWi^vT17w5``dXUW}dQO7UI>&mkC1i5(h{qm391mrxwlb0#F&Y@-!hx<=Xg}h9@|M2kRn3wr(jq`QP%apsm=sJg*N%&pb(+To2dF1w2 zQeI|CjfgLpmnpi=p=J`@f4MNo%jCny)>)5vnN36UD>QkTqU#(=>x`ivBOyPT_p9c# z6Z4Z5_ggq)ezM&4Mb|l$);rU)Lm@wzms*xf`N{fk&JD->WYKjFrFB^5dT+>2=56cy zWMY2ugYz|hY4VdreND|IZfK2uG%hz?9dM@(zMI+7>WBB?J6Z1fqQ0hP5+$X#yPkq~ zale*(@!iaPFTekWd6}ZVre+eaCv0kW1K!1zS4ZQ!nT>04eyQf2Eb41YYtOQ>+(n<0 z)xjTw@V&q$AMiVi?*qB(i~5?<+D%`pj^l|Gb)s30l%G7@baW}^CyV-;(z@c~Lq`tp z0iywD@x8!4yliBqc^`(e+Ac-K=p|7x`OZv-oKOVsPqu|Mp$ zD-@#Z98&AG71@yC&X+daGVi=aochjmk4nsPm%F~`I)~CaxcQ<0$Z+TOCl=M${vD?- zy>hT$U&wM7UFT3*?>=Va4H@pdOT&c~;ZgDGiDf@`6#dI`7xgvOb9NaIuUlaWe|y_S z_$-)JWN|J&1G(#q`kLxFRd+|**!_u9&uCtB44wsB9pvWJ2cChbuc@B%)NA3z-#6pc z&HG1efoH)sTv@5sJOfc*Q(7y6Vk2u8{fl)vLo9(!Ho$VGikYQ52E&9MNeZ)W)4(^nPB#OE*D7`m(m_740r z>$i45ua`kmsP^b;7PY!h3+5sufZ2USRfqPiEo!K=%5I)SW+eW2K=j4|DtuT)(8Ya8qJB{u*4{-&tihAa=CX*c6lo9|5D zh`tVIV=U3vyC!&U!k)8wW!PJxVHd@QRhS#|Y-x)9eneR^z6#B~m(wfUhPE-<6u+Zs071v3lw>evHAoK?tOpVGR^ zxXqvRPe9Jo_KI=v{ma(xnphvddnv7(Jaw^z83q4%V(J8#S+Irs-{|A4g3?<1Pr3C< z%T$$HkC+v>&e2-nz8lv$l-37yzCVE(1wYZq*AHeE>};a%+wy<23Q=EAudpe?nZ&Q5 zv$pBLGg^9fu(Sq!_+-Qm_%`GB+j(AvZ#8CEVZR8!(d4c#>g&DE+ih?rp_gG-KP)gs zwXxAMX$_kHGWF{~m>20<&ku%gHD;pShT%7wsIN25+PLCOqLrt4uLkgp^hPc*$7l0( z+rkX^HsgLHwKu`H8r#$S+E>jtny9ZSt!=-bKVSs=2l^#n_Q5^^OKIp9iTej~*BA9Q zrS+4(d{E|qIDRLq;}zIPVD153({TSl)Yp{OPIafefPDl0B7So*>?5#^A!&A+{R67! zShhIjQa?9cv~Esh(r}_vIDXH{gGMD<{J~0(*LIua{>3fYSPG(Qb3NU#a?R zR?-CSTe6ucEO&iTUsF9NtJb-e^ME=2Bx=(G*xP3nj~{);J$+GMQ#~i}u+Mnd%ja{W z>pX(JefEBhPZ!OezNoJ$ty8ttkvXZ+eB;DJBY<1M*rVD{kyjyieNkUiT3gJ1(BIwe z7C&S8cMfnX*rOiXwji%U)Yp{OgIa8U44eu+^K;q?;8w6_G4BtNr2pC70)g?$8OGhJ1#**_3n=TKVTeHPTEHJmTH`(krfIOD<^H=57} z&$h^2Uv!;AX+8KxRWO`s;cJuLyoNI_tYwYIKW0HM5?$v|T92D(q8tF{i}(ny3^?P$ z-tC%Pt~uKx>T9RZjU*<*x81WJy2Cqro6#WYOe6CiglBzui*836!`WXJm$xJm&-}_= zU)0xi=GT{)2<1BC+quCz+qKRz=}hCtpOrpifmcz{Wc_kD`^#>4@9K+Zenoxl=S5M=8G`Z`G`kK;ubo-oPa1V|9`>s0<_tIE<2kqK;A5GNPl-5CytnnTi zFC81U4(_F~QuEK*n)_&?>l~`*gnWIPRt?UZj%fO5ESx)K`MrjJ#`C6f*B4#qP(7#J z$~ro5&Xk8ty09P4ow7+Gc^mM&spvX~>N!*T*S`$sO!>ZUt>?hGQ)XgU)J${URCJv~ zX}!y;Ix^Jwf3D8efvLtizZkI!S!#0E7hUI2T069!w!;DLCx;!q`4gCG%vyKrK4hth zu5&1@$IehL@P+%y`!?_vz*J+e&0gqgSZboarg~1u@OSW4L+&TQU2(Q?N5%rY zBQAG+QD0L%XOsT1^L2q|w>7!(R=6w9LUt@2CfyMy-lde*RL>djsg+a%cy``jCTGH3 zai%w`_j=78aZz7WT5mpe{{`e8@P((Y8bS^NtI^lx5#}GrU0>AKl-8G9tS(#t`3J8f z1BO8k0xK;EQi z9^SL0v~HjAsuM7l`Olgikh#qKSIi%VtYu2;OO3OJd>)>vvK=;3x@TE4Gj})Mv!t{x z9r!sM7|VRX;ZDe0W($h%tkJNRMSV@pB=jQ`y@4&u=kEUI0jyb8e_A_hWX#H4U)0ys zOk&ua_osj@%R6)%eHmD@>_4T!Xk^Ta`kI_c)QIX2Y*`+$h1mmZmN^_-T1Uf}74G?@<{+?{TRNm-{(;=}MSV?a{iMmTXOMfqtG*A!90V2-{oEGw4}^V9 zX#G9e&>eCQ_<$w~%t2r?FSPHY$v>cajzO@~AeUirs-A5+OZV=Nv5W}3cTe@4S4Ayy zAoqZOk7$ZH2yDe|x2u?cAa{ML=QO<^Kc|0XoGPolqjc~7-tpSE@!mbvbL^gN9tyb! zynCKE<{+?JTMpgA`~#}z#AduU?-v%YN%LLCUmA%y2yEl$=B+XRfYMr_XZqntXt{EQ+8eVgnC-Pm z%P_M-?)sG0ijDC{?4tiuE*f|Ya}e0enDH5ye?V!i@ST6AKRL`lgTv;8mh!>+bMz*Lb zjTHXfh}IYNHKnyek)T%|R;g4x%h<9Q)=_bxPp;m787b^*QtOQ~8!z}&sZ6-=W?IR| zOqRFJp`TY4Fbm-y(RB_rlQ^=4&4k(!TpR!)xF#kaA`l9O`Y9_Jp(8)KD zd%#~#Os)kv2rR3ce;Vc=h^}*}Z?|1L+RlaC1HR^Q9n3*sl{k!4bLvc zK_(YZ7+RhW8C}ft`kw&I=90U<=sJhe`fS@v?IDwkw>mWqGrHJ;n?)Tkn@b_M&Y`q^ z__WO>$mHTDrdR1eMi(>hZl%yN)BopZ7veCO6-G z0CO|ht1c^EV_v4*^+ne?RL?n@Q_H~)c-JpRdt+`U3k|m_#=K0?bq>{YCf5D<7;-YX z+IbS@X0qB%ZcWkTWs0stD6JLkK5Es>u2e3&_URjBow74;%9>)vsoeEN*Ey8digsn) z4=SsaelsmF>y&++ay%O|PDR%_wALdEYo4i6uFs9ctW#Ffxp)I+oZ`Aj)YnwciS;jj z=sq!+-`MJOZxNjFN=sPt$8G;DCU<>NUsFA2x=rD(pqsRY*KUga%J_*EQPjlDqFwME9#_s3={P=skK?PNo~HAD_1Lq z#4dqzdFHeH=Eq-)k)DC*If;+kqHsTf&-95kznGP->YQoY1@{$37il5Kj$a$rCk(iD zY<0=e~|7UkEMSc?S0tQuiO&(NLAnvzvc71YA4z{B|8H4bM*W zoCKxyf+@#cO(7qx_m}L&bD1>_yfG75 z%cAEbD6Lbij+o7Yd^9Js6(@nY%<5a&dH)}4S+MIztmyAXw7#gX&GOQ5Plww#|7(8eNQx@0N#-)#+nHB< z!zm8(GJVb@R?cb^$i5}i(?`avsIL!b*-4p<`}Ou$dcq#h27_LQaBpYm)}hFjxB6rLfvB%3trs@m_Y-mtxW}RwlOYFzWvr+dq{%-JJtsl+oQ3_HV(tO| zYEo@Dnmw1ZP{;9lx=7 z-wl{Y!9SwDrh1Okr~~6n;GH~oMgO6YO~icD#u#8Gk=*q~eNFY8wA_k$)8U=G$?~BE zWD~JT$|*fDlStIpRL`+Yy#B=j-pO{;mxn_(5o>(s@@`Egk*KdJt!o{3+0q=&nqw4(^u>s~qqQ0iIc8OVQ z%;3D~zR}K?rowsCkqvIv8T)TuB)Sf{ldF_6+lTH6C|+Y8r;7Ogt$AcX6LNcef|z9&o!4m%Bp_ z0vr3KMik~Bh^})sI&YRT+k13wy>G~scvbGW`x4(>CvPd{9`OD}8-pMRfo=EoO3>sV zh^})et+!P8U?vwom0WQXGP;=e&bfOqn@jHcqU#(=>+3P;|J{bQ>HVIWZjjN%z8|~v z7PGlT*Ey8dE{Bg|CKvyZVm}Zvy4dP0=bM^rE>T}!3S5eN{9IvnuIDV{idfTy{kz~^ zf7`6Ukl=an?}O6kJdUZ1oxa?B*+tmXm%F~GuWLS3;~qa(EdQ&r-B=OZbithh-0N>{ zSm>CY2Y>%^`;cqpm9fPxr*mvzPhZs6#@$kCQHq&J) z{QL1W_v&=1icPJ2P)qOMp1!EBDXqWT2}rckV1l{>pKw*28J=i9KKEb41Y>*c!Em!5ot`pe(ww%$+~Yn58ve$~JI zWYKjF)pMTDPn^*g=0$Z@9M}oDnXF!&_A@arQ||hr>l~`*+!$01b254DgUJgaHW2tg>p7Y-TPr; zFl!se7cYU?wcPbZ*Ey8dKFd}t83A{p9Ot_kb%(o9=g8P2jlaNM zs4+{{Y-pj;ODxo|C%s)g9k6Q2(W)dU2bISoaQnJ`DZ0 zha>81O6%03?hYgJQ{bKRzUGX$^4RU$+Fz~>`v-E@7xguzwfXRrdv6A$f<_K&R`K%K zgs3}VhhYCe)Yp{OMdv&8w}kv;#mkq9&bP{AJ#Q|Gi~6^JAbL)Mnn~Q6@YcT$+((P4 zJW$^f?xW@Qs(mo4Q#6yizUVm#Y9{fexv_yM+($bzBy{jbxR2(g?dE!DOf(ZcCqd05 zJe|DjUWfZ={giXx|AhN!<>4K4YA%bGW)-66Bq*)b2ZFK=0xKb?+Lk3VfR)g~ZDC2h z`nQ?f^+nG~P+E_krN4PQuoCi3y*G3QR>I?##qpKnZZpwy5|q{s_szQo04t#=VR+0M zU?tqx`hIfaf!oqE5cM@RlenQ8%i6#l@rTL*4Zho}S(EIPy*ls=6@DQrT-aJ4$nZ;*VIg6+mU5n#&AdcLfxO2Yr`FJOUo^O_x?QtQD0M9 zZ+TmzvnlY4_HA4Ib20FXcK_I-HKIz*DDP=k*ScFDF$1UGUuI8IGtqMrRL_~|KCVm) zc-JNq@4Oif{GyXHOQ!S+RZGu6^qd5xwW3{+>s_l?O6}hz5k)5$BMY7{U9?Go~b4Faa%lBzQ6wWW8Y-+`#FDXkTXN4Ff) z(qcoumgo-z4&;FkE;sd)?l8IQi~5?<8n(0I_@HO8mQP0R{|g*QlO4t@hSrH=qQ0iJ z{#NzEFD3Su%ZB7&;6RREuySm*F~G-%e?-?g)J$T*`K0rsfM4WXt<$s-z%TM}-@m~M z-Uo8m7hUI2Gl^+B-IAUHzo=+ot)&*gFB+d28a4yo2cqj7Y9>(?dGWCp@Qc#Mwhv*z zFIv)9718K^dh`Ahtp^;M-Dv=@ zHGf8WxH$q_({Vve)GX*lqU#(=>m6U)m(>Ec=JXJsTu)$YM)xmBTl23MiTawFNp!EZ za%nK|y%gUTPA$wTi(O@u>{klUS?>CxzNTgpqjC;RN&&u?!vCm2&!n>0hv)V7>cTr& z)YsHZBBi5^&sE@iDQ@=lwTmf>_4<^QIN;wqS=85*)~CWB4bcJitit!~qlJgx#&)~= z-`hs;K9IY^e9|G?KQD0M9-~Lf|k_WJ7747tE z6?b_XYu_hy_}|h0dmo6dbEuwETPONmZOA{^|7)XiFz~LMHfgQBaA_QqyT0f;hw3@& z?5^!d2Hy3K-t|(d0q?q9qj?vbhsQC|bq>{YqS6|ry94ifLj0sop}@QD8f|YLQuV)Q zAiB<>v{r09yktqkirBqxk|#Y1guVQRx>us^#xc3;i>`Ahtrd6H2EM)XCAQAGpJt1H zN3g{@ea*^s;+g0=ht#_7mVWh?eTj8u*^ZUKDBor{#(TizcxhcE>T9a!Y%3~#ek?G8 zD-@SJmcPo2z1v|}AAAOK*BA9Q)pNpbPJTZ4H~jr(UD}s)$%_qp5fuFyo`I;Zsh(4_ z+q$oDz=Kzqx$n3x=rD(pvHI`)Kzz`1ju^(GLnzu3ZzjFlk;q6ZJKz^__v^zy4A{{aMHN)eKe|{@6eM z_O<{04k#wVQ>*Q%SIpzoyPsQ0yo%M$Plh9}f?^`*+}-{jcood2$7qRHv1khqM_z^O z^$8|I!Ox1#t19Bu1J8_=cojN79z^3gB8rLdpjORsz^h=Vj$D^`6}1WrqmfrZF%hhO z8rE2SC0@N{#9@h7VfMMtdE`|n2qr@B-3||cSHWt&@sfBIU0%#?uHjWsTHkW>-k4LI ztj-)4Eb)utRNL<28Fxx+tyPe!u;DoKv;QOUiw1w%uo3x1a@VJ{ez)q>$%{u))Pct?JmXGjeZ~7`Z{Qa(i~8Fneo;c*9p=a{qO|^;6)^C?@Kp8nH3uYq(Y-Sx zx8NCfO6&LjPosff#1@+PNc^G}!y2yE@QZ|fO)wGG)qbZ~3(w&7hk?hDS26W^fEsxf zvezeAAC%S<6CwFvjxO*jSXg%^@hZ9+-an4K3SnOpOa#xoh^xS>V7~*K6d|u-mErXN z<5h_In$p@Z)72jMMa-xCio`E!w*P1f@{8oIPiam0n$p_Ub7&Uui`WW&MdBBA)tV;p zi$r}*X?>wycTeCKu_1--62ItJ)}`JWei7Ak@>=hi7M_)^F5Rm=7kPH)eq}eqduLS7 z@z@gG40v{|XR{KCXXkl0QsUXkU7zYXX@k%Aj||IDADwhb;@M407*GT6ol(z5C#Plv z@a)(i-#+okvkRQJSK`@GJtrq5|NBbcO!boTJrd7ua$q+PZ?)sG0`*z+5Tji0Z_V2R{*~+YutN#qlKcKXJd%ANg$Uk7i?*+zS{z0{g&2ljR zfYRD_%(wk-oU_#x<5nYEnGM}!tj7EUO6%}dn`|Kefc3gjDCHj%^a{VD$v+VFHPv&T z_BZJZJUgaia(+GX>=yL>B=PLzt}p6qs^|Qv-Qz6q?AUAPDv4+JtXIP^$g>mmHPv&v zp8H`8JUjN@p-kf0C4JfjX$7IuH=3xgDXonL-@F0&2kgUli&2<=uw~DPXPAE=cYRS` zQ(Dj4dfyoG513`mn@*U2u*0Y34$MCg^);n+tJ>aAA^(7NogCqc`3HL4r-f_s4=DB% z_3ovZ2t6y!4g#-&9rVqWconCMzev0a+3OSRC+gixF%gvKoaX?q0{DB$60hP@dX~hi zPzbMc{{KvbyT;DXfLFoJbPJMr6+{18Pu1`$D6OfULuox`T$_8qFJhOYHx(eis7aP> zAo7dku1{%A^&Cp;{fD}C1bz|oSHG3`MQ6ubbVPoU=sJhe`bpnQmw{izwl~{Z7x_hf zEACd)@QZ|fO)wFB@<%=eUInW<6eaO0Jl1LbMP7yc_5c5@4@zr_iLl~|$#vjWus?=D z60c%eeyYT)P$&d_O)wF*<=7bjuY$b^t2GFD6<@ZQ57qE0M14(ZozUFB82Cjjz0sdY z|BzJvEYs%M@)`gMv9s|FKjrlm<2>C@L7S$Y${320bQ(AATc=rPMMa_pc&RL^O`oaO<~j_vcxlz4W@4YoAU@a!n9sdWyewZr(I zn18@(GaD_;KQKsYzaH}skiqQ0hj&cnuAR|3zD1y>KTMV_6(i-hgSvy;2NsIRG> zQ-5jKF2J*6RX)8PkZ0F3ZQ2{;*@^m^>N$UU+&TTP^Q{5$~C9AV*7!_lFa=pk4=P*B+(mIIiu7LbxR#8>88}pNU zMEFd={A9W7Q(DKbd*G=KOi{mRvuHf#Cx<3#y*m#1$&}W!YJC|C`N^zXsQE<9PY!V& zIs)^PDXsgioYJc?k29pRCrry=f|1NI@d6_Kdz4JrN%luOP ze6@LyQ77taY9>)VahfgUWwN@KKawyn)Anog!J52GQD0M9=j?O$gZyNcU*q+D%ujAn zKj%K?C(B)5)Yp{O6BZp94Ef3IP-V+-%un9Bc3=tSCyV-;()x9Jk6g%4X1QtWUt@l9 z)SZ2cH2KL?&*^k}#kA8Faq8k`U9zvhS(W^N;%3)Vmouv86!vaD`-gc1+cv@XmI}@n zUCy`}Y#-yn%jRU-E4yi4VN!tRL|*rrXnyjDqiikchQzCIIB|k&Ue!y(|MfgIRhLoHChbkinbNb zX&edv88PqnmF^w;OLJ37YsJgOCw4W=Q!2*UZa8rIt8&KHH~WSKO<|PQ7meFy;j^;o z#aDl}j8omuOwyx?i{J-Qdjc@u)64H^Dnu?)sv>rh3k+ zlz?vo;a%+hzWJ3c@NRB3xxld@ypu(JP4ygyjsr5a;az;JEYAB3yqimwt!nZ4-#b~< z*Ob=B=NA>>`sx1ksVxtFPFA%Z_^=q)QF7N8^);pSh>(`j`sur2%aJQjq^OoXyxcPo z)={FqrnEk=zt9!ePvr$mc)De(%5})rC%BHn-+|~lhnh)js2`CId6~>HFUcD7G7DCC ze8;>@x$BFrbEuhwk;5m<%VZr+g!*A#=I0Jmw_{$W=sJg*Nn9!VhIyH6_J3WRFfa4` zw$tr2d6}Z?97^jE{%&0%KbdVzaOsKp$$K94>4^Esa@QAK=TKUAQa-s1`N^zCdG1@x zPhPzGzgn|lZ7RCXp|oy$NM8l{$!xjDy62dm-2UCOKbrhxQD0LtiG#iVVqPZmyYwUx z^D@`^8ZXDZOu6fe`kI0mb<>FuPLqXEejkD`N_=9 zZb=a4Cojl*r;GW?qQ0iI-q9rj^OKp~u=5)*KiS6jd%h+=S#+I4^_l~`*tPh-N@EqnuE_?ku--P>UMql24O@eum=sJh$IbU41 zc{G4I(ax&xY%SbJTVFQ!`#q~e(mPpnokMB8Zo08F+gotLMB6PlU6t+G`83Y>yiML$%UkkFFlf_(az^l-8!fuF`C8;PGu2at>vx>eraJ31@t` zE)w-M)pN`b9cnoN&US^ht!z1WQZW1d$o0luIO8REeNkUiJ!ed!-LA24wrf*P;hM9+ zkQ(*Z->3te@e=hl)pM4vdbQ@qz*JULTyp(5Fr>7)oi6Cw__Xv4M14(ZJ;iUfbRKV- z=ly%b4rHqI0)l7axjeb+i~5?$6n1-s>&Hb9th^rnC<9StOmu z8(izigm9N^m3L>e6g-!QzXK{eg4(;Mmk8ymC`8Xm5KM$snVz44SHWIBSY3*|ioZweP1Eoy zD6OeIOG@h%O)g~vzX)5^uLY#qPaP4Q#Je|abN%cnFzf;bVptVdwb4rG4d)t z#?P9Byb9Ut6RZzPYxv>+f0+pDy=ow@f~~QBFYzj-C)H?;yb6V|umAr{1lw@~kypXS zs16QCUPWB~_i-9tg{ZG7tt&s@>jeBF<~!)>FXR_h+qt|B@{8oIPiam0n$kMz`a0wn zF`v)pz97G-c!bY0AE8uE+0E?&3N@QbLPL+#yDJ*R{B;-kQ` zV-rfp=^@Xq<=^}x$g`8XKGk!my?g4}6pp)%JUgZ{!`TFRcCj~{RwK_&^qd6Mb8dw+ z(gmI!8|BmOBl7IVCH1ba;n`7IQ!@!l>%&ctVg3QTeb?+0<{tzF9rnij1G(!{T2nI# zO6#*54`BWQ3!0ebiTMZnsy(=l`3IusBq*&d2M@}I`~x;9bX_3kADHK)yJ+$cM14*5 zoL=ACAkU7aygtwfd3IB0J03)yo!s?BeNFY8zdr{e&yIQQtNs*ucD~t{mB_Ob^)=OV z4jdYXJUceFY*`rc>}C(OIjZ5=iTaw-`tzvin18_7#tRQH|KP9lts=}nkh{L9uPLoB zMZCxS1Lmr~qXFh046Qdk8uJfCeNAayc6}A*AF%LuV?#0jz`)b0fhPZeVn0#uUW$n@ zdD56t;8n0dH?Bfn#o7rIb&yvfe|^Dq4#7nDnQyfXcopnUiA^fAKl-6xcl^=j##1waSPd0yCrd&Jnz)s{BiTaw-`r`d&w!klf`b#I2 zjQFfPdLeJ5hF?VW9O~Um^_-sxDd&J^$8sxsJ0j0+r%jSK^6ccUFS^d5dd`5y{T~6( zjwQ73;EFuEjq7UOK%SlGI)~~xFD#>zfM>^AS05aPJUgR;Yd1AKJ4$P6okM9owQFt6 zKVW?9$wQcbU=q-E1?C^fU0-yaLuq~5UKjHZSb&@D5X?Uav+}-<`3Iuw97^l3!bzBa z!0rrMum|%G-hVUq*W@3F`kLxFGu7u71J8~rx-M(8-Lgt)(OK^(^6ccUFY0Tm=WOiw z&=z=hOp$W3lij^4WqZ2?e#o;E^)=OVT3Y$+0G=IiU>3&gH~Xr5ZaMg=hG!@0Yf9^9 zO&skZ{{Y^>jUKR?Ka?kb7KCB`f!y^)eNAaSDW?|ZAAm;VMvYqjLwPo2UnJ%qi29n+ zdh-1lMv#BNAfv-DH}{9~@w+hzn*0N5CPD4pQ!|OVr|mE=lhv-QJsR^eujodYVqT`) z^{JTzwRca=B&Kdl#k@?GVw3p^^D>qD4?MxVOwn@^)J#IpDhTs3S))T%iI|u9(jYcU zlb17yTE|D^VSX|*HHmk{{N&8Y27fR=S?>Cj*3_OQrFG1lahRXXjOU%7j`_)l zr+UxB{AAH{5|q}@suf^TBv<{m0D@^D^0y#q}BHWzO4D z+79zF<*qO4YicGjW!!hn%Vd3Pgji!7DEBh1Sb^))q<=;#xRd6{gj|MW=A%j{vX z^?)WXQ`Fa#)+de!Vtz7PI$OUP<|n^+vbzoDC(B)5)Yp{O6SlR-{A6bSX~hc6PqzJ- zJrMJgMSV?aJtOf1<|nf|>hTVkpS;Xr$Qez3GSzday?d(X3`-oPGiYTRtJc8v#!ui@ zd}!Lf#ne~9Oz!$r&!P72sh;C-yW*lDaEn&VJ$%W$~x&72i9P_Pf66R%!u5+lF#FA%6FfWr`Ja?cI=4I|Utv;#A%M@Mb zP+EJ>9gX?POl#7~1k6wV)~EA6%uklPzUVrK(z^VkALb`BC))#?Fh99yl=FMcPZnM0 zP+F%2&cOU+)^>Ea5tyHBVrYCzlbOqWw3+^{A8wRRkAN2^UK)?EK0z+2jp19)GRBhwUm~ zR#gwa(d4c#y3V0`PSJp2`70J?u|bz+j$~)|u&<9AUOxrjXrk*Js^=`JF(5tkWfr@; z$s%T}>0ah@-Q%7XT9a! z^g3E6BzRsH%w|niT<-8gxm*9?!i(P`nB4V6eNFY8pO?ouX6M1*KXq66`er|rw+Gb9 zw0RN1M14*5oVb{+4lYx(f&cyDO^EgnrR|)T62t9Z{wJ`<_94 z`}!}BZ?9A~?%@Al^?$wSQsC0(pdDA3o$EQvxFXhcVgD|uUE8cc^fgxmrO$aBQyDvb zx%o2mby}z0*mHPP8lP}&Se$yy!GUV*IVt%@Q{a8T9S^)ofcF9KZQbt=z7OQCU-O|F zwd0EAe^s^{D`J~2xKn`IH8(6oU&Ft*54l!e8C&dfItP7ytzZZC96l=4HPNFoPTkfi z@d5UnarUYBKH#wvXPd$MfJf@qkHzxWpdI8z*?v9yup;)!_2Z?eosE67nc!>4 z?5??2r%P3AYUP7k=<7E#Q?cjpr$rxrdWOZTFDD$2#GbQQHxAziyv~5dH{gB1o1Kn$ zta%?$T2qV2co{tolUy-t`(!FL!gXxiOBMZIPF zZ%>>_IMlP92lFC+`>BNj=0&`Jv;A#wUL@-4YnHEZCczb}u4L(tNL90+{2b0Cp3cd_ zc@eLn-~SQJi@0|CuOXUwk*KdJt$#i2-vss#_@J0_eb_(XyIFV_+&_@JzNoJ$t=;Xf zN#3XqVc;I zGTW}zaMh`Z?VMJYir>9d&+$Lj+8*Xb{LmENy)ZB0=36h$!+DY1^{Jjy&?64Nd*MBI zFQa$UirD+}2A1P@FV%C(lkB2lUc~*bq%4DZ5!bu=^9jz2sGgI$^A&#ga>YwW_2M=a zvF;uEe8BHss^`2+GnowYBEDeFIzyNj@n1IKr#15;N^2^UkaU&KJpDpVFGjWTdp-^ez$47xD5Q%?H8xBL3CSc6*w_a4llh4?SCz1z%xlDa>4N*oa@QC2HPv%U%KhiV zelm|;YGDQY$$Y|}i9>KdS=855&w00RdIjt!^P|S)5wM@k7mWLLRI{Hf>T62t1sCGH z;r;>7Z_yE|}9bJ6+u+L3;@9S+*ITg8Rw* z@~<|Lu%FD^)&A25_mkzWFS^dDX3}uQRd`R#yy0`aS9pdxs`-dcxMx|^bp0jRPv+kg z78hVYnV(Mx^TGXO(REJry$am32_LF&KS=BMvPv$R* zAI#J2CyTCgD6N%Ru5NI?i2En+SOMpY_y+%H`nhnvNbdTg>l{k!+}$tnd=Wo0Z0R&O zU&Qa&WUj~aMWX8*O6%XRoz3BV5%1Wynj&yiQbYa7*H1-uGgS2=kM@+w4q?Xd6-G7)%<(7lB=F4^j*CTm|H6XD{4Iqtx# z;3w`ixdyxnzH!>?N)4|<)Yp{O$uaBgfM3K%-_tP!ei854{x(N`k=*q~eNAcoS?NNB%RPOqs>l~`*blSJ<7o0ccDHm-&!g*6(cX56eo;MX;=TJRo z*v}!kaNd-E*?Knt&YN<_0ao#v^QNNf97^jC?qwcuKbb!=S^5#~C-V-z{rBMgWV!2$ zu5&1@&o~sHhx^I=&F~-V;eIk-XL(wU_mf4}Ih58DExez@{bWA<$MNHEKbiOIQ0A_= zpDgNYs^|P!-G3MG?D)n|UF6yENH_Qz1fHGT^+kP6^&ESJ|3%>0ar099ZNRhRb58G_ zjXXP1UsFA&Tu@a*`OTBovrXUFe0)H|->*@^m^()xD0l03*i;EwB`b%gu_UjKIe zY|KB9yS}KeDXoWQ9Vv$V13oJ-b~xl8@Dt1G-NO6>QD0M9FC5yRLH+?>;1Qh(`3GD% zf5CT6{sA?U=w55(Qn**c72g(4EzBy5U1gN)hxclznZ&Y$O&8!k8b4Xv!U681@wIP^ zx8i*?x$9FiiBUNRCc(WLnEM|!=$TX&`|!NpUc6UB%_Lsr6(qrZH2!B+m*sFDjSooK zEZs+=W)dkKZG7Nf4OiUk>uVQN7VGsXDG~40P&0{R?PfiJ`)J&$pwAV!kH!tYjEL6U zN29c+GTSMwp9H*5g!{?7xqbW#xS!0m^KaC``^j?Gr?jRr+bOLLybU|U{bWAsdlS5$ z%o}ujCEZV^w5BrKDXsl|i_74CGXG`W9`7gfWi4B0Y3?VB`kIzIR&fmgwYjDIo=c@?6*re+eYPo{4KUIkzJ z-K<1i$r}*X+5gx$FaaK;ug1(6~Hgzn~U$;YWPJ|&)HT~`W)`vbA{rP z$MRQsv3ENR>x1|1sh%@)$fOUztKdV=%|l)V&y1^~i@XZC>r*`^?B?X>aPJ=W54yB3 z>5>;4_97@6@7+^9=hn>eKESKs`>VZA0bT{~wsX#9*nLD7zt@a)o04c$Hfo*AC;67@CJb3%Ll*9~}f{8Z=|){5mBH)hwar?#9iX#d}n@qEXfj8#K{ zS0Ul{k!&vz}i0pFA_ zwd#X>Q+`(8>%apsmsIRG+gx?SEu8^0>4{1BtLtZBTpw;*@=4Fcdnwm*eUEVSd z@-n$|pZ_2~ne*Q>Fh7~MPF&m^^ONPSFY0Sb>)x*4 zpF(~zw>x>MALJ+VYvX3wV1BZwuPLql$ED;$elnj>=Q!pk^9NgkbT#?OqU#*0=agRP z{2F-I{BiT`$h+phONY52?^^EqqU#*0=R7VnH~_qB{^!rnFW}bq{_}y^%WRO1eVZIFT zlezA|RbHC>WYKdHRL@DAnVA83nS6TaA1lbqehn1Ou6feo|B-S%>k~b3wfFR zZioB>keA6(@`Cj;RZ9iv5F}l?BLxzYW~B9 zq3}+YyS}Kesh)G`=IS={U}iUZ+R?=JFtc0gysz~_812el+`wmI2 z@NQ1i@|$G#@0~2_Yf9^8gI?dvfq7BqI_pQ;!o0}eu6^McIBzOQ8srZ1 zqCFw6=PreL(X+RM{j1=-si?0ht#tz~rKn+El=We>>r9vzMcysG^KtZf>31Ny&Y8Nl zA{#Q>`O=14=AE~QQ{S2HQ5gxD?OwYN&iP>;!9UNOc2foW$(v`lyk{Td!Q`$ly3Pr1 zz9;}P+qwOTMfJ6R$Ei!N9PHN@GTV3C&0G)%`^n)y=5MY8`^m#Q&)?c``7$QD&e?s; z$Qv@-d6$L@E5f7V)f3Bp?kIxHcAuFB`io&7c|faf4I|+{b7q}z=-IKq^bADTIh58u z9}c@-hP}+dJq>m{z#4UVaN+VTur`&uzUVrK()!ukHP#APlkR&z)pi=JQD-)p=F}M0 zrlRW{N^A2mKT=x4nsjkjzk|JCjcT%TiRqVrS?;2~Rs_XH)-L)Tt9Uuz)8=S}a{iJ$ z2eY6lEN@_>sx=qY5kB_6({b?kr`+G)LtUZ#9e?uJ z>ncaSdXVM8HNWCi2LAT`_&d2brYGx%zYkGgZ**F7ECA}88UFY5RfRI~`O7wjF01)~ zdt)rxN5!i~9qY0Qe z>T62t3$r!{uMSL6U5UGFgWqsvYn46l+fCHhl-73NHjT};NL8&$b#ujUxUpC7CE~Xm zJ_FHp4%Ks&2^FfWq3zUVrK>NywQzbKyt`v*r$hMxTn zb0hcD_tqrCyhwDNL-m}>MlB-}T88tfbuH}r!Q9AbtWj9vzj=}9I)~D_)$+N6{@4cd zCHp@_IKu3Du-^M258y78-1SA*Ih59G{;p};4`$ciL1r=2VRr4__}cqUa2HB+okMBu z`tYG`OPF21Y~(#s2WHo16Qa_3{>$za^)=OVHrn(lTJ$1ab@RvT!8li)*JfN4&Qs;C zFY0Tm=X5*TZba<03{|JWi!sZZXE_wK$9bx#uc@9Bw?AQ4+QCfK%e(C{%bI^!kbhW{ zVJ+%wO6%EU{=6Lo{Muw4GwIIYzQRc^cvn#F`l7z3w4UgEcYS}4EY(gMQ|Zp&*$uTe z;9WscUsGDot)Jw6!8u#?UyqT}oxwd`n_g-1lSR)-P&0{=B<;Cp++w-w@TMOOHl*q?z)Fz2zUVm# zO6zpf($4#@Msk}C0sYs(J>`3Ap7+!TR)Xj`2}T7By;iOu3el0wMua*9n zSHm;7y>tBZ7`P)Y>T7By@niC+!SV16>X!Q2d%!a=ZLVsijtrHafvB%3t=Eilo;{`f zGJjRJWKA@1j0!_1<+O%gBzJvLUsGB?+G>>ha#IM8%vyAJJaCNq4V+;g0li4n*Ob=X z3Jo@o)eq(Rx<=O3`eyLo`jvy5{_91e=On0}Q?Rz@u@=Co&`XS5YYp5ARpp0I2Z3KC zcYV=w5>(Gg{B66?4!A{q9L{MC1#ZQfWvxE01AdX{ISH!gl)6}EM{YgGeaElN=?2`2 zYC{TUo0vsOzXQ>85|q}8c8$}vw{N8W)wJqcJ}^Z_dq=nH_biIZU0?K^1f{j2-Oa9Z z`ZZEtwHOx?2TT#8jnReOfjK35PJ+}rx0ZSyYovZ*e)DD+Fh%UGr|ujL%qjd%7WFmN zbGi*WQFrZz<9yek7RQr-Yv-Wu_;5Nf1?8?U>T9a!bo1KT?3rF5ub+V(lZeCHKn!U@y`X}MSAMiit^sZ zz#M#DSf^tZ@UG>qFY0SbYejdvmIo^J)CLWTw=4$c;ND*EZ{G*rwWzO2t<7$zifT4g zU%y#7tpu2Z#yxXPy|VuI3`EyC)J$SsozVsnSI%x-^)D6Jc%xeS~Hy=bh!tVaM_Q*@m}X&sR38F%DdAkXJxJ2Zq|v^L|yo|*r$SVVnI%_Q{G>Ul2AN>^3c zWneZDKm4QqP~>~bU0>AK)J$ULLSwyUVHv7jomOHt5x+cZLKozFiTawFNp!qh_k3^P zOqJ>5DH7*vV$1$FHJL=BzNWNJsLQ(<6(_447O%xz9DX;}U1HD5U0>AKl-3i|ZO>0S zlA_Y(xtNQ?kJwslK=!PtuPLn`Os>0GH9S?7abc#!p#6SW`@AL(M|7P-^_&VG{X0W< z4^Iy4>2wmf*YAclzTF7E(d4c#y3V0`&e!TBXrk*J zs^|3b@Xa1*vX^(Vp7Zh+aIar#&FZN0?;A~YokMA@NIx}oorkX4YsA4VH-JeVt6J?I z3XF2u>l15AYMn!At$69*b>}f%_0;e4uNDK7{F!bOrnFZ0`^`Vpyrw$e;)KNVk9+>K5E=fWz9zL!oME!5 zdrkE>J!^^OKilfgdrkfU#Y8ZO2+sgs1yh(+?$|x*m2&*F1$~iyKrsLtz^h<2qxws{ ziXi2<5ad-*Oa!k@6TSeif+=R%Zq~~#R0fn>szP3cf?y(?uza5kyb9L$mun~FRa753 zd!&X}L1|58wo_WSJ=-w>_(g1=yK4^ei*ll@>L9;J?)sG0RAxJ+b+e$EhQKdk*7osk z$S?9)=6QP-uo5V(smyjt>qA`)KLWpq9scn!7x_guZt`pmzew2E1QVgx`lt@Tt6<-J z^d(-!k(FUnkyjymeS-BtX-zQ^AgL)7copo%j0F;};!pC*_Q@PK?`H{n|)P2 zw;cQwd3IFK@jILH2zYkvj;~V)^6a8Q@@{K*c9hmsCL^WwAXERLkbl4q=vnQ<{DTRN zSFcxD6OeXMoQ}^HpPt~|9}Nr z6>P)&12*|ifF}Py)YnwciTLfL1fCt!p4(avd3GnY76%~DPVV}mzNUIkZur^^;MuV# z3$G;P*$vHnKLUAnqQ0hj&e=A5J%MM(b_e8!AkS`c>VqH+&ra0Wl-7C?S2sca0h`(K z!9C1B=+`@0%0H01zNoJ$t<&$BjfeaL*8D<%F6JL-_X?2m4@7-UX`R^q!~n=YV0+Cc z24nt#$-7(smw!O9pQv{)#Y9L@1!VxQf;CQGE%7RL_Sv-@c@?tPC)iKayO&}jeA(u= z33wIEcK%X{R}r?^^AGYW6r$@Kf{9=}()T^^Dwt!jLgG~z2l&p_@G2;+sh&e={oD0! z8t{wQ)69Ajzo^ek^UugHlDj^oHPv${trN{nHvqqg)tgsO;uk#}(kdDGMWX8*O6%&q zT{;54h}pGDlK4g1y54g%{32mr6HEjPuUqlJt6ia4X~?Tki29mfBHX=wvH|can8U7IiB~cDL%>K4uR_$9?H2fl}=TPrns^`4^Q2hw-?3nciONnP!qwTgJ^_~b@j{-R=sJh$IdwOm(gvO#i%_+Ycy>*iIyTer>?p0N zbq=L<(mrc<$Uk72QyWP62l+;GcVqs6-1RA~sdWyebwy9}GmwA4oL?Q6@(;F}KZ(Ws z1JQL3rS%)ts|?6LVE39cDgU5TwWj}s^|1t8M*>^cC7n>4-(I=bDO~u&ra_8 zqQ0hjj@}!~j=-~H?tx1rp55%@dJ@l0)Ynwc3Aj1h5_op(Ps%chXIJp->i_ZVM14(Z zose7j9P$s?MIJBZALt*+%fS2tx$FNw+P*Zd=B|7HL_($vcV;5<7-dT5Y{@*&gk;DV z$~@0xx+NuqD3U~jCQ+STArdOh3Q@^eNrjOATHnXH|8IWl#eVU<-}rpISdZu0*I9dg z*IN6!M0rhN{b#vLF8Cj?2HpK7|ASS#GH;>(fheyjtV8sA4hH`NW;@YW@;_L2u-7+@ z{{hvLi1)6EewomJ@78C}%Ad*r{WGcPmr3;`I(XNe1Adw8;6B>2=Mg=$A?LBm#b$ zpkF4;B1dG*Qv6cp+Po-6zf7tpF=#>0Y2cU1+OLXjgMOJ^oBV&jObTntvz@|Pwfo&P z@K0vjww<|+{>j~pt4^SQvfS}0tSQfS3hN<#ObWq2ndztQSdae6J#wZOqkl4mHRah( zVQsAY^B(voGvkLf2cUm)`1T%A8vkTbUQ>6qo6h6$;Frl74bd_~zf9M{7M;*9Q||bp zyr%AIa}%S*;Frmw4xSo=ewofMmkvO`Oi^A_J&BcxrKMI4QpC%SE5SdR)r%Q26aABKJoEC=_$O02r>9%38@`J(px#cPo0Z8FdfD%8)UH%cIi^!h z4|Y|koYSW&HN`JCgPpZW@H}AV$Y<{wb6@-Kei4QB?4EBkYh}a#hkW-V5lt$T1I}4@ z#61ZL>wZJSmIlI^O~0L-bJ{#c?d~<}4(_(d9iPHF&(IY2BbdT0++|!+xzhJ(>I>YH zps+6Q7&84OoY~BE4$kw-V$?*t<7Sgvf`{!Ah6Nl;k3 zz5eaA0nV(qHt)5_KUUp-&W1qEZi^_dshqR-L8SL>xRV3-U0eDc?&P3mKb=!x=U49d zqP(VZPILW#Y+~R}ZkhY_+x$lE{Df6C{UNaPE6Qss=hV$QySn<$bk;#DJ?1vt$wqS~=Z?DN&r0)&4zqVYv}O<{fU_-NcSWuC{4 zmfty@q#j{5X+7?nit?JmdSQM-bY@XHtLHq{>54_NI_`7AInBPQXq`j#BpSP0n1Ej< z%b&Yp82V+_jQ-aG{W9f_FIwkNJ&DX7k(J<=$tJ%t`;2~>ZS2pbqF<(HokR5`X89j_ z4t|+zQokGb(JwQtTDe-|mnmB3P*|sy=B)w$WOhp3#UA~W+YWK4gZ{~K#}}=0D6CJe zn&J-r$t=Y|cNY35dvqVY4*iow>l_Mele}C8{>dySBg_N+lN%kia?|)Hi}IT4NgVpR z_8|CWvX_sNAE006bL-bT&@WT&_@ca~dJ;p=r+I;2CL8v(o#dApkgR8kewm`Yrg{>m zdJc&Pzf5NGW3l9y>AqvNi^eZgl-Cs2JCk*$f`2kAt~sea`X^sI(kTr6ljV*t%4-Vi zUBgG1gMTvXJ#N+o^iN*dsjLP1CyVl$!unC|{YvmpW(ED7x1)dZnH^Il`ZQ z=nk1#djRx{z6aI4`4jp@lMb9~I}6^U$sJ#`&Y^P7hKJRB4543Cx3-Cwd5bIjpz_6* zm+&4p8;QQW#$r>{Z3sJ71DX~lovqlwlb6xQn=w!=42tj(~rEm;|< z>g{}EAAAcXcYM)0hr)Wz#%Mcu1I13g+&Jh}K$<$a^k{2*3ng0TP*`uO;gxs+dVDSm z7Ju<{N>``OXa1VEP@=r1a?aO`ZOtdcj@O?b)Abik^X2)|JFD-4hLqg#MR`r-oQz4P zN5;aASFhNS>%3sc>$USn)ApbtCCY0m=Qxgy-B<}bUXjfw-aHRGUhDc@E9~-DLrRp_ z6xJR+?EjhaGKrb4^{;;GM7sLsm+-}?UnF;YQC?G6^BcA`zL_Mmsk4o41aHew_xLet zJ?a;U@|wcB;_*5y%VWvRNVngm!7i#n>D6FZOC583ieL8B;FJh%1W=Z-*E47|W`bBcbr?94GmK4^n&ug~< z{USC&6)ou(J@;td6ZMNk_arE+|E$?^1@w#9h>W?Ce$n9>um00766Q6biC`E1kQ8S*TwmcYF$Kiq{m@ zVV!m*fPN9{H|w&bU$kDeW)A8XiSnAlI<3d7HK1R_+U!{==@(r)*Wf?>A}Z%lvwJG% zl-5~28}#hh^}$L>&#u_fMho@q(EyT{X%Y^z7K2)2AdoyN!3_^fh{R6xLKvg2MWTQR8Odf57_A2$B2`_8Tm8K>q`| z<5O5uJqZeH|B|H+;D5l@JyQHe|AXlV2YI6Zf#{wDg|)Tmn8o0Kz-)4Aw?zMg_T9R- z)%YKX@|wyyC;wgA4fO07RHTxg-C5@VOVqQIJH9Bdshsor;J6E*XUAfavL!vcy*Z}S zP|r@3*Hq3)IHzX>dUh=KW2&TQ=YFc4zDCbZl-Cs2=GUqM!T*3+->H!N4^HiPnt}cY za>p0tHHCGD&_BJw|A6WEEwVxX1Lt%1j-mg7D6c83rw$5w2>u6bwDUp9|Db5lgfNZ& z0j2#!-My41LiLM}?txwf`)x2?(yQnhzdZu=DrApOXg^VRFQtjl^{Z|t(5qnfJXX@H zXkuJ;5A`Y(qIC|TiSYexmy4iR!P*V}SJJD9_K!Ez=v7cyQ#psi+Qg(o`JWx8)Gv}dK7}=vb11C653GW2+xgIMbIsKB4tNVNGcwI9zRF zIjJH>ojrb*q*u|U_0b17(^rV{`r-eriLkxn)KOCcW7VAuFG+e8u0yOtG&6lsUQ<|) zIlp#I_S-~t@}PT?e$irk{TDb-mODO$HN|TRYn|IB9=>Oi)B|dCwMG4+-M8*;#CftP zuPLn0gntZk9hR&^>CUNJBk4 zx#Lqghq`;IvvIw%=ql*hu}#`0v8ZP^Zq@t4sAnfy=TJFk)vL@~pl8RHK3Xp6+3Ejo zbWfvaM`2CXITY3lC-3YE{s-*#@S(xzf6!COyQBYs-0>-_sXB+kx@hIThTwm|y8b>T z`5#y|9`+Lb4@B!63hP~Yal^p>fW6kwk^B$pjn%%d@jnpdHI;L|vloM6GE>#Y3B4Dg zo?W?~Q%{`X$Q@sl*Hq5Y8}zaNgTOTPsSQbzp54yoj?Hm~Bg$(k=Zvagt&TuXe{^N2 zq-WRVL|0$U3`dmLB-ZmsUUhL!S9>~Nl>86gsQFx+AIKeFl-Cs2*2a!LPTMoox;gHW z|H1u>wQO;IAj)eBYfH|oOA70i z`s^C`C$m}(Qzido>!WR(p?|XIo&<%p>C4@N!9STYt_E zpMqZ|YhWBD`DMmWcdLPZnR3S$iym3FRc!(NBKi7TpqfuB>dlD4Z&rTQ(08cKqIzgvL zMRuy%&GnQQdUMGgpTe5jlc2D^V|24Icyh4;og%M}2}x6Te*U&TdUJ{HNl;jK=)Pnm zcyh6se01Uh*L1a_vDZG0H^GygUCWU}qum-LD*NmPe*@QXpeOu6HW@|wchK;3Kt_++wx<{P)vy^y4y zFh;)?`elmpn!?&9?}a`1WU^H4E4}Q-B&$EAn(xv0Ws250R8OK-#-1YZ%Vayw9?eC+ zOiz`%Z~@#0a>o~~bEuw#>w<9w;Frm~zUE7Qne&#|%|gFS(K?6fNkqRNdma2TnZx{Y z$uF~WpAN@0ewm_m4u$o=?p;@ce=@u9BU18DUS_2=0R5BYjxSp0P*|^+qdgM*lR%Qo+M1pdiva?nUu^iMuj-TqJM-&`chYZIa;VIF0W8k>`< zPO1zXgnpSn9z8pbce32^MR`s2BuuV)IyDVVQ-4@9P#^s=%RN{9XaIMzD6grW#2CwZ zMZHd?tM?8pll(Gk^j^J6b0>@Pn!>sw{^RL^rHSgHv#s;dKY2%$|9QL*Du(oO0W_kCq$?EYo%hsWPa>ME0`)KY1(K?69 zIWzf%xs~QYZ2f>cw^g8bZFk&o-`t38Oz!xibqe&YIy{4Kf@J& z^z2tL(K?69IWI$6*Ixm8*Cnsr8ijz~_0N8}b-HvgmbCXp>l_N}YWw%3fQLJ4w(X|( z=rJ*B?Los|hk%#6-0?;09181BXG>kd!=25jVU!v8J4S7_BJ^Z$@NyTeb11Ao=-%E5 z9`3A3^92*E!eiCFn!oecc)5%6n#wtoy`B#`JvK%iHZvy!&Vu_i?w)~XAa{IGUQ;>8 zDJLlGcU6qKX8zo>a2CA%ksW4u2BN&Ca?a}+`c~I&$Ex?<@4FAqg8$f88mBn}QC?G6 zE4uZr?eT97wIbTUK-v67)bUV9L*_9Pz6<~F97X+ZMB|I{n!;LqV+v;ql}Fg*&Etza4_Bn zvd1T|8!%^ZbNyd2s-?yQ-GJAumYs!PPvCWv1O3k-&%cFkaA&_`RFkZ(H0=R@Uu96a z1nvVKeY;{W+z0GstCl@EGe7*AW@cZ}F z#5Ims=TKO?g{G{5?*krm%m%*?*a`m&P4W9c?)Vhe^=Cc7u8{>*HObq>oc$#{l! z4u$m>pF88=`+(o6<~s$x57>S8k+tyqfWmrc|2ZZ%M<=TqyJxS(@86G6ZjG?cp|G}I zXpscp2YjM&?n(GQVC~wkf1&w45asohxRFxrdA+jj*Ucq~suYK23-D}Ip?wZRzlguw zZ;$;Vrs_Fr8up9ijxWkST+H-2s!iEKKMzh{7@xrs&vV(tzei1kAHWT|r?3;F< zuGlXU<@L}GqovyO$gFuEEyg6PI$3QA?gVEu>F$hl=oj(gDc|-&zle29dYfDF_Zv-= z*A&(+=k~h5{DAv68IJP=#ivr{&?;HT-0{%Dnb5{O)DtTEUC(yO+v2%~xyRgnkk4 zI580WMQnY1;2urCh{C!HzK+v50^Z z8Mt3WVXfsDQT@v5OjV0U2XRk=tzAFt0q#jqShpC`XcFuf@j6o<;C>OyGyfB%*)J01 zHI;L6eKw`QJejwD7lHF+=KF2?uNyE=mOH*Euc@5l{&_$vm?!hO>6H&*p3ElAzx~Aq z=gGufN?}dqoO3l=Pl0(d_n5XD=gG{(c0*^)JXw_26xI&=r-#A&18#A+)`%&fTM?v6 zO~m&Ha>p0tHHG!hC-eHj`vaab+aBK^up#|y4&eI(3Tukj6xNr(iI}zr|>`Z-Z)~MzwVu{)iTg!z#}}=0D6H>V9=;CyMf~90Gq_*G)((IBbw1=G(K?61I{C>H zCF~dRz@Q$PhG#f4%4_^_aG1N)vlHdDeHoK>89(Yz8=3-l_Atk^k$6|{vRv2#^eVVc z!>XpBSHVn&HM@y=6>`TH<#m~*nY7DjHlzQ~$8cwB&s#hl?`oqH?T&(81rG?%T@QK{ zO#4Q)o~TzL%IkWwhfBMRtz4Y`sO&OTpJs*L!n=Ah^W6k`6?|7@W)SFAu*eUOvov}Y zqP(WCe(&>Q8t51C4_o>ef_@SEc5mNh)Gv}dz9_FLtX=flSb%;Julu6!MbIx|dc7ai zMExRBUQ<}RTH32XzlfiIws;5V7qKZRr?WKrMWS^Mm2=jqN2J2ODbLcLiTkE(P>pY8 zHDKRV?)aj04wZ9m|GQfk_D%WLF~%5#J}Xsdwr(*1S&^W?_3Y%1FUo5w=Tvp&-k@j4mj+)(Jv&ys(Q*vx*@^O+$~kA+c1r*~ zJ09~sANA~5^Nxe}YxL|yc}-z`u;=S+@IT`2kk2GYA4q9+j5wTYl&H?=*es=wR)GuOA*IiOkzew)*6xI~4DXfEUeM$rUBEEQ3Eb13A-OmeW zpnj1kuPLnSu1Veo`bGTE+{36}#MU|9?W_M+d!NcVi=%efEisQ#tu{Y^@9tSxqs6!I z-943a)-G845%el}gR)(ySHZL@?uMdXh1~I}oKrf_vqz7z7}a&VJ^1dPo!WFa8sFVh zXLBXumc=BXQ4*DNZSStpq-Zbu75;d=3 z;LUtEv)-1<-qXReg2Gxcq^RzI?6N3Br_*legP{AZ!CUk{AhAA_|Mo}wlBmg1r6~4|;ZdqTNW;vturq9e=EC5h9&|D6c836}bcU z?2ajjnz_@wcR#43>{pw%t#}_I{JRm2FUo5QYsDn5!pFV^QE5B-?q30QRN#ZWTlaHA zgn3P3J?V{p2bY4V&^i{g3!skjb#NW9E8~A>AX?{8J&F78hgpN39dDSq0`=_J?30gw zc7og`cYM)0hw4d88Rz{T^z3-?hC`@l#~QXDmw|eAqIC|{lStpv*#`9N_~!6v)U#tZ z&)v||=-G+ZITY3(R=l_N}E@SfAgT5(W-=zZe zOdJ?=(#5(ZH%};NLIGYpN%)BYpZf z@XO>jL0ZQqL(lH)qU<#1u>YNbD6c83e|D^%2mZ-?^33L;vtgcWP$#Pq`X|dBUzFDr z)>+q9Jp%t^o^)$*-wiNNcCb7?8vT<+c}-!x?ctE;;GfJdul#}j$*gwCF&&M6vS^({ z<(&1YStmg6n$LNuw{<-9i>&VN`-pnia>o~~bEuqC)c2b;=v{N>s;Ds-`b8(6uYbQ4 z{4z!B94hBbUGp{(^sf2H3r|t+n#Cu4o2${g7Og`ltf4OYw89}Zs%qzyA1|TD_iRVU zR(-(ZRPOkqbq=GdSC|G21d}_yD6c836^cE&BR#y7^>i$>YWGWK zM}B{bxVSljiSnAp+BU@MlDD!(ol8-+u$MPG#3lhr>RVm+7M)j@Mk(X%LvE@~M8a z8qO8!cfF2&IJ~)m&Q|F!!g{OK;UVCU#=G0Kxdi@btb2z!hc*6aqP)Jh z!CcbLEa;P%^RysQHD|ba9L{zwu55c1{4)8k2cdT0m&v*e9WME0${k;n*TWuop(i7+ z(9N9Q{!Egp`mEa4INLexe-azzE7UeaCb=Pade}aE9Z>HUK4)`atz;>&= zH2%q=dlFR6xf%Qe{WAH2=ha4lUnWcK^}!7NGUbjhx+g*9ob3nK`~<&De&Kyw5%^`Y z*H`X8LBCAVJqar3EE=H92ER;x{%+p~;Frlx5C3pWa; z&rX!r6xP$tBKJ0feN&V1FXNWMzUk22(Qmp#zew)*qP(WC4t+jyvJvc?9$ytyzb@>X z9;rP$d?)maM0rhN9iG{?AA^0<4{pnE%z}N>n|yfvu{q(=_kn1glQb|3J=^)l;d=M{ zhsLPZj@T81bLKyd=FM;^P3ONpK8gEQ{}_wkHLU7p-#??bFe-otIR44qsm$ zqe}1N7J+l-(D0ux)-225n})O*!n}?#bya1!C)@|3bxyy6Hj=LUp6v?{uM3D(Z4E2A zigRYYiyu<_ax-`{r$3okL-5y0qH#P^e9f<{Y}a8){Rt`c0ju zz${bl_@Z?Vg>{D;V+$jpHeF_6vopV;JBwJow74G3GDYhg3hTn!diSe#r1NIx`9lUm zZJPY4U)1uyvrJK5E4rm*!i=9QjGpOfslO^c+AMC1v;HSjUtITv--p6z%py--uFUe@ zvDyb_`f|q?<#m!zN1X9eT4Bx^IO{)o!@xQdejf@WpT={pmn%=ZC1y{A znZ78mUpz=T12cZE@HjL0oBLPgzIJt9ZWtW?i*MCXq73CI;Ai84D)2Udwf9k6;p?~na4S(m!qIC|Hb8NH+&oP32QB=i8$5YTR zTD$VhqN6ZBkUPF;okQiE$_eYjTS33*#YDr^nc$n*i&usQj8sxPbY+`-Nkde?Ht7p-$B ztmhwiv0@m!h04377}yQoLY4hW4=#Y-wP>9~Vg2}A`^I14E!4s5zmK+tw@{1c9y;aw zw|6bdYbxgy_?OO&g?_3cNoS}b?~S$sJ#G zPlD=6v@-nE!VKP{U5_+1-3{;2?tQwy(RI{qCb}m<^(0Ix2iFOO_h^fn7JvQ;@6qO5 zI#$2Ns@u{xn&_Sch4rpw-7`*tRzhg&l9h8nD`8UF_V@MbM>4tNi|$EKSWnpPx#u8g zB|KCW@9G3v2~%=)ACymsWTJZ#6xO+B$(_%ER)WsH87u9_K>Dlz@oQd+9!g`cr*l2UmFDh;S zuhx3dFY;{AVEoZ=F8vNf_avyC)7&a6Cj<10LSHW}-3|Ihr@MMBpYkV)$sJ#GPlC!h zX?@}kWrBWD(oU;&bq+jYD-uV1TQE({ME4}9oHMAfONlo49^9|`^4jF~+)DqWd+=-{zp(42~W6|~=M#s7W>qP(VZPJie99aBM5 zFz5Fm+a;hW_@T?N>G#(~Gr8l7@|wyys%}#QyMdE6873;m**D6grU zGa`Rm`3TSy^z-a9YyfBquC`2@xaEDc^g9scHHEc8aX4n~-Z#pt+iU9g2OY?vKik}C zk`T@0jxWk<3TwrrE8A~ed853ZVdMM)M0B9v_Zef@K~Jy2Kq%Ot{e&L^!Gjxt#c@> z-_I+L8wA>#%ILN)bU<74<1$v`9^@jqo;11w&wOTQ*O_P zTqIiOP*{ha{#c{~+L{r&cE8#I+M4YqY)sqqHy4TWn(9eZ6(q0l1>XaO{WcdD*GgrF z?Xz6p!#T?xUzFEWPol}*;wed>@1@AS;4;m=Qu*kqPg%R^gVR*B*dS; z+&*Mj^j(hcuxg71M-!#P6+g5LF$$%eYiR>m;7`>%m!2A{BPZjVRFY8t#c@>6?(glx7u2$?6S_mV#V(uRygI&kW$@PCR*o^SSLNT zh`n5>%xE@$VL51&Z|nADzsdAisV)-bHI;K7S@ZO>pa-v*m1Fwbv%2b6=8;}_26D$2 za#eg2XLeJ*&yvXXRu6gIm zjQ^g2D6c836$TyJmel;9RJ`kR&}_Mm%52uQ_c|?OncVS3c}-!hh|Vg{Ves#Fdi{BY zgN`bF)rqAEi({E6uW77xZxnBY__b!((V_4 zdqfIC6XE1Amr&5FfG4X-l3qpA(F@vW^eQN-W<-956YPq+WGoxujoI^75b^?zmG}kMRue z2>M0LtwENgUoj=~@qOk5dGOw`puw?as{Ew1;Q76xWQ@Ggp$}*+Mrj# z?pAY>^eVpeSR(0F2=kiIMA%Sgz%|gTV4XJFNO~3h9b5jVS0T!43Tx}n#Zqtm+JyE7{1_}CdY20c3#xU=7B)U&f| zaZS>*lRG|@bLLw-vmF_jruK=bvj+9-Mmb+?g740#vx$DV)erRSSeqHUBt5(7DQEVe zo*k8Q_RgMUQtL#zy8Vn2Nzbm``G2n9yE7{1YGG&mi!NGtGqGbf55I9bQ*#F2Ul0w{^x%n%4;g;=vEl?0zEs{{pu%4&u$r? zJpuLX{i+K(B%&&*&lP zRk-BLm-H%Rk56bnQFkw;iQqZ$;ylo+U|UmLNqQA+mQ9!RDior14xx!~q>=Me(5qli zygo{L6)Ur{|I@3Wu%>bjg?0XK1Jo~Kwsvune$nxO`+unrm40riVm z#ND1Is9&_ztEdy|7m3z66xPjOyI%wSB35n2BT2ugx4Zp+`bEOLCNvSw7;i?s3Rbzf z?tRp&_}%BXq*oz-{Qq0)gTk89M7X_qF6vdVah0W#UPaOC6iKf_AQ zU|H5PCA|tmgK_`qRfzJM!utBBZ*M@qi0$xpko1d8-}Y~e`bBcbr?94YO=0cS)&}*9 zm{ZMkNx!J?`yxrdNR-zU)~j9)MExRm{U75xs9%($SNxxT5tVbOyO+v2{S(&B13f!- zc-JjS&u&w}#Y)t(lRG|@bEvzQI-6qC^mm|V$G%U@ko4?&Hfo)YdUm394wZ8Pw#-@t zdUovk!-k_#&u+HPwgwtKI|^&6&Y`f@PO6Rm2dsJ5aLNB5?!eNG=zk!0dl_N}1!G2_{{egAt=9_u4+5J1Q(xnMAj)eh=Zw1S zi+Xk}Cc5=h)U!L#@|~n-CwF{NUQ;>8CbHdUm3`rgF|I?-8hH z$6|~xNqTlahClgF&rX!r6xOc}C!qfUYn5276Z#+YdN|?<`X9(0UzFDr)(<-GK>q_a z`Duyde=u)Cg+KZqi1M1k`rWmr=zqYD94wan4#=b3$|P zmFR;;^&|@3boB?nOr~wTdnWp2vNrSgYWy-Otap@@Cz|%tpRPk6sGv0qstX z2LEI>E9dZ3^iOVkyHjiQPo}UgD~h+~W0KWL6li9@= z_g&FHdFv+6QjLGID6grzdPeo@$>5jC1}xiOgnpS_Ej+)YU#8sgMR`r#)uu@+T7X|B zYg;ou9{n<>e1G#E{W3*)P4y(+9Ra7c)137hI_)v~Wp>=MtF6W_Qf)>dW*oWVbt zW#+VTM*rl&?*{p!f3n>1MR`qOZ86+^0Qe`f!Iskk(LZ_G3;)OHpDfC23hO!_OtQg0 znWg!Ky+!|I)}`qrjejzgbB^}2s8cp3M*Y6+@=Sl&RT&qOpA(X7�W!GdOczrFjs` zK79X<3igZQKfSY^8?lYa9iPfM6?aEYoKY2{p3^Mv?0DE!+32y(ruX|XjLJE6?|CnZ zxfsOOxTWu}1^Y$n?@ylS*{@<$&Uw0Pna}U=SambwpnVyztFqI7y*kfqF{g6Q*Bi_9 zR={4-f~s?kL*UPjeJgD`bugB8TPUm*9#J7R^L|Gua+mCw=vJm&V*A2sh{sGuVSQz* zRT`evX_rOUbH~J}X1R_0jAtfyd+GS<> zX*?_Q%RSiesup0tHI;L= zb-DAcKitKEWrO_p!QJe4Cwx#txRXVBP34@k116`{hP!x_U;h5fa5ujTd+?>ic8~Ns z5al(6b@9*%d00Q~U9iS_l%J|FeX{WCtftgItqUWqIC|{lW0~aJq7$SnRiCB3Fw!( zerHY*`en);U$oAldJ@|$^w2MpEj?8EKfg@dh1=0DQ?$;ZdJ^TK`skO*>iVy@LBGrr zt116z{4z!B9181ipTj$Ye=?ibAif9sCx>-=(j5Jh<&G~}=TKP3{c3>z$;@r@fRE^( z{O0DeFX*2vTIWz$ca5Hg{>d!9LwXkaColTo{7&PaEXr%DC(*!tEc#`#fbVtV(J%9S z)vM*`mnnCAQC?F$2}{oa^vh&^ORW#0U*@0o>xQ9UrYNteoocs^ls-SL7&KAR^Mm-@E+~N z0R5};pkE|+e9=0G$~ku)KWLE!eIljSrcUAT9xc06ryv3PMWS^Mm2;MRf8N>v`b3x0 zKd>$E9?kOi%Yu93T%|i%w9cWh4*6O%5PN%CReMcm=cKB}{|GpMJwCbPi`F?5){z*xO5gz0&orYr5*Nt>Z%M@nKyg%4;g; ze6e2HVlvE;x6PZ;V$n2T7Jp*1|6SPek~_XAuc@39GqK0fvEX}PHfZD~FVK*RjGW)C zJ?wai@|wyymqSB0RrXJ2BbSU1Jr5dE8(g|(bZLA^Is;K&Q&<=FStRY_)$=vm+srv# zwZXZ=6x_>`JH9BdDXeF?&X@M_b`CmU*m!$}>Sk5D`M8%S%4-ViMb0aveY_6uzMMH| zm#ONXofnIHdH6e^yd$XDJ*A28A+`HU(5qlq8eW$4DvU-gU4eQP^2h(by(6gEJ*A1z z`%>E@pjW}X`k$5bD#HHM{f>GS3WeaF1fhu#`zsCgDp*aG(FfG4Fg}p2Zc4Ii7@@l#+#s5!N$8kk@PA?#;K8HPDqYksk~=po+28R{1e zv)!GF`bDC=rm(L5mW=`ZA{KI}S~b)!a_s3fL8D(pw|3T8&;gbJ>-0>-_sh$LdwXseW`X8`W zhnww0|AS$n8Qw(f@#H_g72)2mRjI{O5ll%4;g;oSC7+K+lei2(WF8 zdUk%Z>$#$yo!s$7c}?Y6CDcNGKDl!pTe5TITY46 z>*u0=5i|C)+JyQ=ea9KyLH#1pI)}o#{HkX$=ohggN9}S^zv$_R>FFB%B4J(=ng}<( z&&ve83Z~Gr80uPBQ)Rw?ac|VCkU#$at@S}+O=%*Gy7Dvw^eUKQ($Yy8^=qjPrH0%_ zy$Xdeum5jNgpR$wm4jXdQ+PxM@99@dWoNykT%%VZ%4-VioArO7ei3}rRH%}iYN~n+ zeEt#ji{y?^VNLOx!n&>V;vVv3SUk)7_;R8`wM^$_Y8iSnAlIwjj|D(DvhKI{5D zdQ(&NtLu?r8vPu5B*chkAB$$ER`*b@x(dGv}-4 zBha&B`e&Q3M?JfE^Jbx_XD3?cP&udjxb}}h&yJm{>J^B3cFpSr+|uaTQCL%T4u!R) zYdw}s<5O5ubq%D$Ef8a+Ev zUQ<|)i%&uS1Ez2~Qd)6XOQknE?*{rG$Q@sl*A&(doVK9<0aMI1Nx9&ur3yW=_9prt zi1M1kI&|@4^gnNSzs5gVl-JZ<{bTq+^vh(+{?%vb zm+2Ex)&~7D<&H1PYwE7vd*}oDWwN}lH^-x2=8wKLjnFSsl-E>GqG78W=$FZghs+B> zzsx6fterG|nWDU=u>RBOD*7k0cTFbhqknSX4a-*OpDcHLQC?G6N9cB44gSe&*WBy& z=%3tTx}ypDCyVl$!g|>tFZ54l(~1wxMgQa*v(h~@{>fC%p=S3~&MBW~t~&s9ix!3( zhyDcJiWbYSH=6a@m&qNU$~n~Rp2|6^-+%To%uHd)ah5*jEv~SpPN!GC1pOk>Jqar3 zjGFuGP(xKJtL0Vw{59}Lo3b~5MDea`(oCPin%a|~u+BQxSkjtOtu#AloSCYs@iwk4 zYEH==pTe5jlc2DE5)f>UT2oUBf6bEr9mTNR?brgF}YZzr$bc28sLSLNA#2hURb1Cx)ZfTp0_@kMz}!AmFN&1h9O(|1ro=H-jvs$|pbs$A~ zO=0~o_NJu&*f?zRzaDjuZN2zB&DNKKfys82yu3erVHs=%2jz^pRuepDcHL(K?61dNxl*|77-{UiDq*pM3dK^au1$ z7OitAtcP7Xi~h+hdYQ91`X|?%yCGcTpDfC2swXjP{x&P{%Vdf!z8kzPekhMx)qIA2 znR3S$7GC8al{Gu-TA_ck-0?+uO<_IZpB(g01}*8~OK!z| zRF>Usv=aT3MR`qOy~nukHSkYnirkW*RVP0x-OmMBYW$N$>l`ZQ)ON2n)1@?>jjmT7 z^R50d*8OvbihA&kCU<<%I)}@e8l9{<>F{(*nRU9?SMe7_2>z6!1(k5T~srvet{;{e94c%@?8s%^X zqP(VZj@$CfH+&aoFvadY*=w%GDBJ!yw%q4?5R*H;D6grUQ@8DV>&zGM?^pPTe>uaI z$C^D(pO_uQM0ri+oWQuY^Xz72Li{~O*ApU?17|;}rWf_U-+?HvDXb%&yJXeMhW`&9 zhl5=fzK*h3+{Caha4_>PyptYJYQ)qFYKPV8<0k&-Ap^UzHwh7Pm$0o=kml9e53YKW339 zFIQ%H?^unzHvS&HG!Xv2uAcGIZRRm*9o@Q8&WZZMX25;G4}ZRJ6Yc|EbU*nQ-Uo8W zPx9%A*nxk?r>frmUzJ*6&KZc^$r}d9Yxwv3G@f(4TzT3pF&lY(G%0HOOZfZw0roA1 zmc^)V#;c^9Gq?YX7jPf&;xRhi;XdHJ_sIyn4_-V-IRn^ng~yq}-@qSzU%NW55xe?7 zrda{6xgyXz`t7N5Wwd>MP2_d$y*ik4co$FA>7o9yYWuqZ*D>c5hb5%IeZWn!t32R7 z;97C55;XS#g*D~bPGPOP_~SJAKHwgSS8v1j0Y5g#(gnW{`CyZ zubP|)&QDYyRqCZ-Pa><`dvE9$@jBYsi=kh{XH5z}h5aJAKzABUQE6VV@m&!SwF;>_w;G*MntIj5Ru3!EqOg}OTi!#tVqZfs?U^JKZ> zi}ISvIqNI8<2;%BjBhpu=E?ki^HU}`PZs4hm2>j0Dsi67*@b(7Fi++?cKqY2nJ0_# zn!>u)IfuRQ{(v{D|Kl#aKj7PwALMrjy=%GSi}ISndgY)-`2K*m)Jd%a?+zt8;Y!@cNJyCei{)}fvs#Jc)t7G6$$|YMx#Nr0ITY59Vq@H4 zzlfWec9;qKMLcGIo7K2qBwFWCSQ~jx#Qh?E!!qR}>=$u|-c3(x_KQS$z1%S4pD8bs zcz)sZ7c0OQ$@#~jMW~5T`BWSAD)^_xX^lazf}8c3qei_7x#NrSI;Vb0jc+E&yv>?z z-Ii?2P&?XuS%aDg*JJJ;0lf-t^WdI6=vDCLL!(SluR@g9JDsz&ERQAgn#CjHvsP!S zO*a{&p(et|F5^(Ig8#l1?+@|wc>K=GlWpkKsem!+Nu{UWZ~=%$bQMWVc>u=dvOg!)DNp+SRfpkKth&3gS# zqhBOi=TJGv>EmM)*f-_-eoVH2eN&!n)p`i-o5~$uw9cV&PP>pwla1>Ps~JL9H(hWE+* z$%mH?zixs53sJLDb(Cm)(K?61+RO78zE9?Zc-}U6pUg)o4Z=IZzhAV@p|D;a+ymby z^G9C~p3B>Pm-iUN)*cTFmFgle=-tvopl8P~mcESxJv(k{(0McJ*~uMWl-E?w z$ugdadUpJ1@PGrLX9s=XUQe?`o)?M8^Bq z@$epvcU29vh4*Ou`;i7A_#TbwN!$wbJp}J+xMKct`;gI<%3^orBYamw^(2-T%uj&# zXxySG;yS!X<2S>7!!+;FD6A>Zb_(kQ7r(^8`(*B2Pd5wRC-Xj?-L&z2vfS}0tSQfS z3Ts#Y&iFo=4?SbO7Tzawm*YO+_&%A!n(}O?uwKyMT@k!b=6jl~%7*vJJY@O9RL%Qj zQC?Gb>AZbwRiIbF$6R?)AM`4Cb_91py$ZSGi}IT4Nh}(=KOOWcxIyjVk3g@2D{N{g zQLjRj*HllU#dg=7pjW|@Yx-OVy$bHQ%=3&!uR@g96xM^D4@3PTZqwmaCFmFNy-6!| zP`^m-_@ca~u#W0f-VgMP_?OhDUqHWzzy7!;3H6Iac}-!xwtetu&@bW#bR7PygFnBX zJUB+9Uqt1cN7g(Y-raM>tQ^zdp4C;qGLQ7aclT7zS&_TyJ?K^N{tLDZ9teNNdbY2< zY-J3SJ3f_jT=#WPhj;hT-|wv47*Jhx;9hz-zPqP#&gzTH9YC*wU(m1i6!a?iRc`Hr zdKFa8DR-Jw2Jh~S zRaf4!QaCdTYlY(f)hjs>R*8<#2W_MwIL`fiuG$FHv4oIVYlEOGnVN z;}83#FIxh4vvGCnsrTVdmOH*Euc@4~qvz;A(6i&ZGrPJjgu8j&@3NMtXD7;QD(C## zzqTCo?Dz-Y1oHuKH;*0lV!TGrPL$Ua)(Vdatv>I4r`-HzV&ncpV|fkB0pU$RuR`wl zqP(WCR?I)IT`TyVGJJsV@jn-1`RwXn?FWNig($B{tcNsP+9v;<@}k9yC%%7T`N095 z$7+oNO%?c0w9cV=5>~+usAtEcUt8B71Ai{rYH|kk?BtFwTIWzbiHRdpT7aG%A6vL} z;BNTS;acQI)Uy+kV&yJV4_4x#Pc08c`P<@S_ooJmyVO?6cBkG&-dJ8|U z0DV*LmNWY;>YK_PU$oAlumHaa0jxWkYYOW>RYTA} znSVI+y*BtKb0f!9)u;TOCyUlORL)6m6rTfn*L-C_lak%gH`*M0*#`Bl<&G~}=TJGv z_FmK>(7Wa{Zq;5>=fERw@ZG2s^{z$h94hD7p6rNv*ZfeVdm-pu^Da-bXY_>r8{WyH zbq|-6>1m69( z!?VWF;}flONUS~F5^h8M=;`@KilNtMSlEA;8T9z@3`BWN<(!&6a}zehZr9|UiFcR7 zTzSOngTL3oj+fl=MR`r-oIU+V1_eI%<@d}&2Q7xVa=)RY9~Qxmmng5Poa5Ndve`wO z>pZ%sS^s|qB=cp4HwK0L-SHCTHHEce$d1Tqa~>-vMvR}i2xk60w3@bV2>L~G#~0-_ zg*E8R#fR5PP!9L5rl{R7nZH~0GxQ?p7m4zk#`^uS^WpKzxIgnFZDHm=G+b+OY;c71 zI}qKINLgxxo+~^d)O($4cB;Bg?^>78bHy>_@>%f1;UQ+PG9{S!ynt8ta1OZACBmr1ci0m!z{9N8{6NrZ!yw{%E3m5){^F zyVXX2G#(!9p9cPDyzI!u24-RI((gc&*Vi6bM^8pRw?^Bttl~uV*0t7g=*f86_WMEb z%jD7HXWsxQ_Xe=?ujyNN&e zCv&U7?#1YzEXr#N>qZZ{qkl3Ftg}8J{F8a)g&XTM{>h?y5>(C^ntL1lGP%~;gt6e4 z$umsTtI#h~?)aj664cqG9W1X6ewjS)rtgV;@TcF!tbfrjQ*=*)$~i~Z{q_aFOg?8# z`-R|_$$J`(?WFO`6y1}cuwL4}MT_Rpt5|q#Xf+k|Dr!H{-g6ChTjY)}x+g(lovpY2 zkQMYQB2$k2oCdv$+U7pzd||gmbWei9I&;mOrB2YR7yGRrk{jxSMOVtzXLz?qs>+i}ISvImb;ev{(#1yK|?;KWGO%yPAgr zYcGa7S(MjQ&Y4!|FliI??8cu+i(dt$@R8Tv)K`%Wxe3H_p|1GdL1VBb`f*A&)uE5ANgL%(Qe z@%c6LpkGwz^(wSzq_^}t5Uq1oeUC=ZcK*HHQmu)jW7NLS2k3@^XM6XvhSrtlL44(; z2Vp9hCp%Q6-I*J)jmaHfw9aWgW#Ad`Z0AgC`thJ&G3wqwt@Dk*v)!`Wmqjr!PafDX zb$2b8CvUKiUDMEh6%(y<{uFEP0nc{6-nquaIk#ff7w47ELeKVlYlHMwz#RF^AlC*V z@aJv=!&}`u7)xg$TIWz$kL}pQ?mEme?F&8~v4R>k+$e79KB!IQjxSp0P+0ft)M$(X zYSM#-{U%#Njp}oB`iRC*n~K&s6xKxz&pv4ZHR(;G)h<1uMzw9Rx%Jn-wW%nt6}c-e z9BueQshBWjr}L$6$~gz#9Utm3lg;bZZTieHF{*zyRFi%u1{>dc;_pN5_@ca4IK>oh ztd*lwm{^sBtof$w;1@k1zrvcQRps3b{~e>!|FZow{!SM7{d&?5ejmcTHb{&3Uiw_A zsIa}^ullB3c&g`t_I4Y2w?0n2t;1tgJ5$H)#@|U6!))L??RHQ-{Hgw_zkyed6TL5?Iy}=3hNH{ zqAXXANmkW8cV-QK!+pGVyJZ*ncEd9et#hcH^Us4RpdO3t84!1FZER?q%Q7={Zvt2Q#t4O z@}-L9;9dP}YtcyTtIkOaJFD@q7UeaC^;4Ix^PiqbSEYnGN^b@mKJv50w}NuV7v(jD z^&+L#;|JR^RQF?#NN)x|6kBh_w}PU)rm)sIwPgDs@J+T~<|4freAqtcg~mTwbWein zNz5&JKJW6nD1P_VrcW(*CUMp9mA^W{8xXnUi|$EKJ&7liEZ)D zRT&ZR21ImEg6c`!yqd6R0K7vxAD=i<3GdR*X0L5%@K-}hbWei9dS1udvGd?v^2nH@ zFK59!Wqq4P7p8+&g52>%_arE+-u?c7;i0(;HSf5l{ z_&#qL%5C)=b$7u#<;{ut)v`PO?>(9*uc@9yBcr#yM!^|e80xWYI-Ei1&r{P1;ElN4 z@kMz}^&~19KEJdD&Y<_8BHxX02J`L}PKbav;-b8!dJ-Fc-Y|)UGidjwV9qu;gIdF! z_0%DLQf(^AYYOY*S4PgBS$dsY8$@or4LU~uu2-hEhFm0fd{JIgShsL(mX^Ec2Dgj6 z9XknhjHX`GpAiJPNR-zU)}M(FlH@7%3{V!1SJ17f|EOS=>5yB}??7}-#~0-_ zm2<{5oL2jpo;#n`Wnu4OdrgBcnJ?7uS+nvAK-0zJ(=-TDEy{K5* zB1Ad^QC?G6D{_;*O*@|(Rk6mfS3l4kY~ONs%L>rDmOH*EuPLnI-BG3U)z?ui!kX+` z0h)uozcvrK4|><4ye6@pzvs9j{7uw)_l}kYpgFk0DAIUW#{bShw9cV=65bKL{t5E; z;$!E(bR7ow!G&u9&#XYdNbdNebq>{&IJsb7z=vj6`02ZyV*0>+khtNlLq6yiiPkw( zPhxfTrqRRTK6rL?wy^=+2W1tp(dB>ra761I3hS1N=o$*hMb8%KPKbkCwD@>UhYZlx zlsmp?okL-L@ZJ4x(;ycOZsz3p6mrplvnC^lfwrb-okL-*6Hyv^`ieVeBc0ndgj`fB zEqKqozg{e&yry~*pL*D={+5xdYQMdur1RDPsf{t}d&wPNl-E>GBIstfPd@_ERK*5& z>`>=x%H&TIT|MC%+X=L}Lzd7D=E7|*s`vFIY`UKifk7S;&9(d3RVTIWzXr{Yhw zQ8h0e$3ZM!$6ZfcJI8kexOkTmL`W-ZHGJ zwt4^G2BO#<*a{Y6pklE1C1GJ97>E*zfmmP%B4C1|C6l_6J zLBRUYtlRzlzj@}xbsWESzu|c97xQqRYnU}_uDL#EKy(Q($tyx{>(2m2Id7jsY8}1z zP~9HI%IcGo_1%C;zW+jcWNToQNje~nWV8;e$mTzF-zH&@YvLhAUu zc}>-vp-&>7rGWRu1I=Y37k+TCWbB)chjuye60kQCAxF0$v5H>bqa$Rk)RopN_l=$>S5O4@zr_ ziJ&kli~(K+Ya3W}2zeEKX20r)ybA8TCYT5ZinM+JuY%1S;VJSe?6uRTsCgBULqic!iyc9p9aVEaTwHen zcy_Gy^X)T{XE*Xk>^9`tNgbc6Ia|!`S^>|F$?6|i6TU%1k=n!64|#Ufv-y;tb02ti zOlO?!66Dz}UXv^G?5LVEYvm4a;MqZZAJxv|8j7rK+7FOtN7bCy4znHt&yJ1ib~hAx zb}2uXC8~LLl-86dBc-*E-83WcKVbWOLVeNyU_x$f(f>f|_>|U^CnKfxk%QBhg8u=# zzG=cB^gmdjJt-dj4=AlEPew{>o%+Yo|A1|5WWEFa4=mg&kEs0*c=MX7If?O340v{| zO>hU1XSXb6!6D??Ngbazuc?|dv&t=fEb{Dl^O~wT!VSNzz_Vkr z9Zy1$XBX6cWT2X7$D7xb)>Af|-Uj{$Y(a96=znm&O-s@LKuMWIu@-Xgz9>iCq_RL!BZ zet)a}HsBYr!X0}>eo=%~omWEX-(xd zrL|o>(+9vWVmlXg7x_i4Bma&;ei3h8Q(E6DY*zsMBIfL#C-RGOLrZI``9)OCq29ez z&FQ$$F%Wom%=}V2k!LqO!BFJcNgbc6In=wCdNxz8#ApN0jujdX5_xukhZ97e9dDmQ z)tss|bsGWCjtLuHe?*?0Vcqrr$FrlfrurO8>j_?GH-Y~FTep9L=zkDhyJ`peA4naa z(wgdXD6P|ikDdqr1LiR-PxLx9czbH%J|yH?OIhQ#$k|^6c2AB~3-1-QZj8-}Q#y zpEs|mnzM3h5%TO?X<0UP5VBl;igdu;ke?SDYc zB&;1a+k#&vll2Z=d%?6?zQFO!-{-0gZE{W3wH zd9NMps^yh0+LfYTCN+~_nnmcB$$ULtihh}zty%=B{W2-7DbIFF>->YKXMle)Tb=hN z8U2$d<;Lts|75A-Q(9A=?UdFRRUgqmnH{|QY(4rXPnllz3H_5PttroTO6xPVV{d_f zGJEK^ekA%Q8_!o>Q~M|L<~8-MzVJ8+{W969xcj1CW?Pqf-O(>o>iE2QO}(p6{^+s@ z{4$xfr=#eXS(!4x0R1v~^O~ATtX%&J{W96;C5>ayFEhDIkiXh5lQ*v^tsgAPLjPoT zEo5yA^iRI$SJoB%lckQ&o7a@qL*LF=0shJCdsNOc^iMvOW-%+(t6~zO!QA? z#!iXW=%0MVH0-?EKbfjIflJPXhc3*6^Hm%qMCF2b$~Rdc+*w3y}pHJ4T9 zPZfSN@fD&rFRN+{&p_(SB0uMX$i4~u1+KWIz%!s~PW0_Ab}Oy(A-?Cpum?%<;Wm$IH-lXjs^+}9dM-PnD37(= z{v~pU(LQ0(P~Y2i|Lqr1S|>gAe^S2y{u@qv$NA5HsWe^KO5Bs6v^Fo_;}i|g?2Xau z6+tE`s`L9+UcucKspC^xkH~qu6ZRt@haP^}+4Pk%_MnHjCqZf5%GhUSAw08^TFo0T zsYp@nY?7OTyDgN~2j(9Q-2WKjJDfbeL-k4-QaDoFlc2OVx6eKRJ1y+=hD#fkMWw2C z$83pI@3!#fHC1z}*2jljg?F-l%kt$_@J^n6njOl5onNWr^X4^Gb5`wbyf6jc$=3!> z`Qy;yq@a6rQrpq6^UIspRLz;(=+MeK?zt?ix#``j@J=>q@aN==e>=atc};1ZJoiFm4zNysldGngmdWK5A)NXSwn{?L0XLLY@s<}oz7u+}H&1*{Q zQR&xH^2>7B*g6mO^~^F=4N?{!Rqvbf_Bqr{!g=_@q2QOv`ZOsP{W4RN1~T-^lsZ0d zpF_w(765M*rmOmt&Wse=={M zLuq~C@(=V+W+g*(&!B&Dgm!O#wSO{iUQ;uP^8P=)!7r1o*k^eM{W9y?FJFs(nNr8+ z&1-5VG4ZQI2>4~PX5&=PiA8rtxutUa?pF<*65$i zo7a@qQ5}=eKbcwgJLHD`$@3M7Z`Jelhj*cP`y5K^MR5@sr(nj{d$20A`+;0lpYaZt)OVqH^O~wT0d+5S zmOU`bDSP3=qC-megha%Qpe}bYpUkVd%1YeMA-4_nrFW@7g+J;4h$*Y zyryc7^N{lOHL&BQe}7rr3E1)aG$HYQ?|%#_-n^!??)G(i8_U8B7J0FiXWxUlDt7q8 zJmeQi9iKO^DXm}c@~rn`NG2Qd-fD9%w>;IkWg}dWU&Nc&l-3VpeQMA4&1466INdq1 zGG8^d!;~yFzlidVpw8|oCW4LiMC4VlJ%2xmyozalE8ZZlLi+gs+dG0fyQi24n%a)Y zt6+mpJ`;HrKWa}dL0*MS#mElc2P=e`bsP zBG%h~yT~tEb8L%D%`f84Yl4Ze$+w^ycopo|?)f6GqH(Jz0eKaY$0t}Hl-3jzVSL0X z4d7L<(6de=ufpQw%ah2fkn!dtBGY?p0NnFOVE z3;#s)KVSxt{-XcEQm3_>(EmW{_>|VvOoGz-=%1$E;D5jpmKTct2gmk5JB0oRyn7Oq z*4nuh=zqX=zIZMA9}HP%(?ackz?;`p&1sg~b^!3~m~FQ-k!Lr}w`>CP?4*v*o7Ys$ zIsKv5IpEo`GvhjjAB}FY|AEx;dGngm`my>e_ zkHG(cS=4_o`X8*;Z6B)kKcLu8)Vr5rA`Ft(egV7+mbBbMkBmyb6P=xQ=RG1*J7r zb11DdT87sFei4h+A0hIK3f7l+A-_oK_>|UE&7rjRE8E-~_(jZno*?pz%=LeoBEN{Y z&!MzFbi_3d_(iOKt5lI+l+(85W;MTvJFf{QLdBRrS|)Io-}K;akyr6#x>Ew4=}R7; zV0}M!lzaPy|)#*4}NNEX}#U3%V6+7U^j(RqW^)8k8Th2Kae^;r8U*(P+ISNKd1xvA25%m zUq$}|yLxkS(Eot9&!M!gcgAZn_#d#DcZQBZ|APYOyGd&Q1Kzx*YEHfOb*@$BWvlF) zx`;fx_L?U<;~9?B@p~P<5Tj`gp z`hMY?$g|UOoO4WlhQpiJl-9EO(Po$7e5mAcwdj9fZ?VP-&kv-I&zskj)-J-0!kEo@ zs>?eoMgN0dl^f^b`2lZUQ(CW1xR$O1=LR?O{)+wwc`CD;>hl9?CPAIuw;*N`{w)t? zgI^{axP1O3^vkTcoLr55nNr86W)jrdJvEakmOs)0zf5MHXM7O-GMBd7wF~_+dG{o! znM7$wb{P0&vf#T{MZe6(bCg}xewmck)R`rv^^lh?y}>`3UC~?K7yXmpT)EsE{gb7R zPiak^SyEcRn{_t~{F7PCmnZMgKiRocRnvLk)yTUiL1}%+v#|pFlbK%I^?B%@oNQaH zsrFCi&1>piZR;>73;Z&f=jj)sU#8vQBY(o-8Au(UH?OI8we8O@+TfST&NS|R7yUBF zE$m*o0GlJn%{W6oY+UTkMGI{fw)cWgefACLcVVdQB=%2h@ z*Y7&|CrcflH?JwJ+njM74*tokzpYLb`X|@j?)?J&lX>%+(t7dp`p>{WnQc3|^ELV> z+XQZORQo4WHHSL8r)o~qruUCq0KX`^b#{6r@QXa2mE5)ie>AD%Q#FS=yQgZ7(~L#h zPk>*vs#3e>72p?5A9q;eA^4;5?nzKJXV>aIo0|f^sI#V^*a-Zh`euVPZ~gN}qqL^> zBq*)7dM+9Po?Ogqp77fGS+>gYSfexO%_VhwN^5FQg3`K^+56TXfH~##p{`>@e2yx& z)l)6>=HlIxptMdJdwK$Ra}U6uMF@LnkRWl9~NH?JwJ zHNKl#fKMj7=QBEY&Z!L5Y**G4{W5v;n$kKrPir;!WU}I$9-UO?nX2q7UwnK0^ULJz zbEuibHKXM}z%P?s^Yed!ewlrRxM%2>DRq3_K8Kn~M7J7O27Z}rhUSNJ=$E-ZrF;_l zW%Bkp)J)>Y&5S7U%Vb>+{g{D%nO127wyFIxdHWnn>sO~vxqyE%D;uY~7yXlueT*H5 z{>hTZCwfY%=1^K2wDg<^{>kjOjf*AvCm$Q2(;fYjdHWnn>x2m!>ENHtF37_Up?~s; zq^?z8|J5SiyryOnXJ6(RJuS#q-SWCL9Q`twFX^gT zyDk=^f3l{Y>v?=1NFARyuPLnq2e-_9dNM;b>1NSX^iLj;68Yf-ybpNun$mjnjL`N* z@E))WxwH!XlY3mNYo>l5@b)=W%{lz4*rCQamX$XOx~c%)^{Eamz2+uuW>UweY7X`8 zrD~2>o2Z2;z`IWF61}ZH@UFiYw*A%2WjW*RbEul5^=x?yXW(6L^IO{@9(dRD0%pAJ z)zbj>S!FWbK8Mme{*g;wBrwX$->w>*Yo4OoI4I>;Jb1ZF9iP&g>T@WqAJu;y2p;Zi z_3=+NqyD6*j`g_YJP^FxdHWnnYa5$LFYs_@f6JXezPXaBTC%p`=Hh=|?!0+T)tr(+ z$BnWkq^R=d+2XSh;_WVc!e=0LeBQjKYEFkyvqBAjr>K5~=i;*vu56tBX#hL}-n^!2 z&ckr9$x2LB)oXhmpM`Lvd0C+P8Sv&crM1lAcZ1{oi{vsr@A5_3-;|yPpWQ~Ep2eh& z&zskj)-u_NUB;VU$YqDt9f+OsP5E=d>P>Gd?FHVvCbjmECmHQ4l*=yuwN!5XrY!t) zpiPgJ>jbPtveU|V&EnrmS<#{`Qv<)t7cG5mi`qRO7O!YIAw|*X)tem5>yg#(TPHwX zpWkGQ|C=54cMK~sPf=J0jw{7{*B-Jc5#9sLs8u6-crOUQ$J$=U_krZ`qbI&Sn*zUo z#;pVWRNv*lQ-hA7cIyny53T=|qRkB#}0qS&#g ze<|ksus>@(;XT08&*VD7dqGe%D|2fO?*o~~g7wFaL;TS_-lo0$F2DPvs55Fem|GnYnkApG#W){Mc#Y#iUsvBjVva24}6=ku+sb$@U14Cd$8CPztNEMKWQC}ZU z&1*_)BLg?&RIo#FzR0Z*MocxXRP!o$`<(c1W%zzzwf=aGw0V}Tm|PvT1K$&R;mwQS z+<^7!I`2E2BM84cRX@k`1F7Tl_BqLUtHm?RW$yO_;mlIe%e}M==cDlk&cSeQ0JU@E zM>t0ix+?6Z;Q0Y>pHsUwTQRRp|p1WWBUMhTG;Ne2@SG$9` zinRxOTH?Dp%5qSgJ@6{-U8*w>c7KH?o2iE2U z4yAQKul+k?;eN7rtR6Degz%C+bCIRS+viYPmj>nLy@UJ7h53QVR1*enIebLTQsd2Q zs^+v!F|Mm<{E+Dx1>s$B!M0Z1R=guFb$sr;CTdQl_krKjfoC^%R@H8}D=rkSnRyKF zi1X$(Rde>#KBa8|JiELng?Lw7c=oqNJM|rL-n^!??i#yKr+>j6mYn5`J_v%VmP-!$ zA4naaH?JwJM~rqiO9cOeEYr^DgCLBN*(Rd@0dHPYTE`xC&$|r%2T`Z4p$~!(eKELF z?SDYcBwA%P-r|&(t*D%tBi^eyQqyD<-m9Tz6355nM*+ujyZ;`WvSy+T3cmXUYLF~L(xch zCf>6=YB|9j?^#k>?;o3A35;d-JLA>@U@i-JOWt)t)-t8_of!?i6sDPqoDe(Y_X^cV z8eG78mXy|I6{_pNSY~&oO+@ChpuKHRn3}cBo7dD#V*1hp$d+YSEz^)SD@0{lw?M|M z)bV-qnwm)rtyqt2S!OAdBWqSz?J{x`GG=-6nwm*07}*EevMl$=17yt#&f8X;Rx@UK z^P1B7-18ggd%)h!yof#sLhZssy6Ar(b$s5urnFw|q_Z4+57_foebEO&2y|)4(Eorp zuPLqPM%SMWz6UI<+j;as5SF(+8m9I?plZ&>$rdk4MyDw1cOQUEec@8lYyt1xQ#D81 zyj?8#9ssqCF9P5QJ4t`}b`Gno?TJ29B_Nank3VjfSlqYGgfFtqc=vhH&ExS8sq-S3ZRaLXe=z}24n{@qCZSa&JwLZP_?cQ@$${nM=+k$t6Fz@z_ z`si7KYZbhCP1PLx>NYCyB4R~8Zs<)U{48v62tA3Uj?bIdRL!~Du*pX7B4Xja@1Qr4 z(5Tb$a`YtP&1zt~Q zulzf|$E8QmM_pWDVY@aDn1%2UZ=XZWB=&_p88dQ?k~v+pL>~lUN8JZz=zkz}eBM5X znn~nVTrAtsP{sa+JU|}=VZe*r9P~fn?Q^J^M7m+uJK%f3etHF<4}x%iL;ve){{!AW zhtj&G^E3Hl=uK_5UP6y9VerNR9rWgsIzDfoLutK7zG+ee=uLliZj2sX!lWPChUm@3 z+viYPKWN@34m`P7(^l!|(IpIRcFgqdCI=y3Ed6%u}W`k-g2 z)bV-qnwm-6zn!!t0%mrbBhH|Asc>t0UJ`ni^5!))ld$XUc(NE~cB?*`pm(XzPS9GX z_AKSiYf5W*hxvD$;5@l`;l*w3;5^xTy!o4cr>`)nXWM z$%fs^lia~SnK!Q~txvnviLl1=S;OUL)P9+~eF&ws?DV@XbJAbSW8Zk7*Qww%#zpivl{!9epF?Ra^BA%1 z*RVJ8_V>Nf>s0u&DOmJ4l8XwW)6#F9Hyryc->g?f% zORYlL!-DBX=CI?nzrt$e(gW9+)bV-qnyNVsXS^HU7k0cNKFnHb20LC!;Q?{4Q?D`J zyryc-jOwbronXhybNH#2-(bgU%1FhL(8_CKE#l2$3)xx!Lcw;^4@59jFf5y-U@mdlgk)jT`iJqb$d zHG?wTW`I9h$nCDkTo&A?T24gPvefZ;_arE-%Zs{(uLFOyxOq*Hxh%Z!T|Wj{%e;FM zl-9bY+P1)0W=EesL*}wDB)FoUnzhWE*9-ISiJpubp4ANqc%7~olWUNMXFFk;Bc|R4 zzs%t;>uCUMRwzh0=!uM3spIqJwaE^B(UVcOKFdloAVYE7cGY4$+lfd_*`Ee}nNdEK zeSkG9RK$gtAY+y{ukWVDh@OlsFHhV11ZvU z(*`wTmN%~{tv`>gv*9ZEC+{hLZv#FE!YtP@ndpBYb$s5urnFY=t8>!={F7^TM(2PJ zf)LcAr4{-g@a8q8^~tVjlWPskVQqCzZU!F&;rRAd1JwQpyn7N<%_+GaS*q=x!G?Nf zqYr}6_WPO|^goa~KDC=go!wJ4$E|aQW)}EmcGx}g1Na~aiyzm@MgIfdJqfDj1oUe( zK@NVI4d-7+9|XbQ&84Z@|A2Q-g3`Lv&o)*Jc3bMztY6awc3V0&f9m}gzV)S!Pialf zBq*)r^^e)kgx!|%;�jz-|lM;_;z3{QkUq5|q}xOwRaDgWZ-V-`h<-3cD>gMwI0q zyqqBRMZ9@U)tpmji(gp7J9%?a{0;CX61v^CKZ%}1Qpe}bYpUjqZ#8+yCU_^e>6F|9 zyorQ4I`%!$lZZF3shZ<)?_j&<*Ob;=hketrD|d|B0rMi>yr#77*MCk! zN7y$Ful1^+81_wPMDDI@`EOpt+viMT+eFXy6XTm#X`7@d5_+A#jOWa?gtH5RzUH#6 zDyz5PgCI2Y_lraS1F7Tl_BpEqJ^%A;@33fA(f1UE#pOQ&o-@yx(0qxzQy#lEzj!?O zAP8%Q-%m#W1KvKTP*scft+Syw?X+e|L-6Pl zYJMn!(3^|5&!M!|GTK+%26|JMWYcKy=n`r_`&6O!=Hkukq^JIP#t--Jz2p4nzf_to zY}EtL`p;Yb+;({({`+M{zo(Wf>pG{oM8KK8)bV-qdPL4!JmZJ^_Yc49Z2C$Wd(fi* z&-%khrpJ^#hu`1jkXhIGa%H>kgZj^aGkxB?J~00%p7Dc!$jRe7RIijFg(F|%S$}t@ zodevg;olE@@+hj^SLHFQVG*tWo$2%DHKp~vS?dQ3)yjoupmH1Xwo=)sa7xEYI8T;3 zK5t%AT92!@bi4Cr`28)WFP$>+hjQ-zZ}GR_JefDIDXsgAcei<62Ju~fRSt{#p>)(d zGi2+(^JLyWhpIWx${jaNhIx_Ixj+8kn<-3hSM?VCGNq2s+viX<=l8Ff8lzxdq}#Fg zeDKW_bRFKkM!!tnK8LC~C+vF{^nrO%+>#!T!8cQ=X|!~<+Aov0&!M!wbzuChC2$vN zi|SI^2CGoP+spayYM5P19iO+)p|l?0IcU^ExC`Z_b0i@e?m|5oV06C>X4kxZ4yAR( z%mBmIa2INu@{)EC+=X&Tiy9UGZ+6X_*Hq1kX?V1{5In49r+fWrf1+ABDe!o^L5~xe z)bV-qnyNVgW6Hx*z{6TL@OiODWwp|_`>dyZ-zGBNyryc7z1;R{6nI$6zKpy!erk=< zb5fMYxZmhYE0gi&HKn!FqAY7i@K2U?ojF+jl4^9I6sg&K5t%AT5mYNr&fcZ zneZ-MJ;6ApS~*rQxDo{C2fTSrX&qu#s$2^G$+GJ!j-RZqR^HB&KjA*wxnBbsCBS{OdqT?=8q2SWe+RsK5|q|4lgG5I2KK&o z&x5Qpuo8B=A3m(rB$-JapVFE-v!t|+_lZmU4D5aP?HBUi04rhc>{p7fQ<53)o&=?} zecw*en!rl%nRjjLP+%qG?%Hwe=D}p~8Sv&cHIvAk8n}E8+!5bBYRHpAa7Vn?arUk; zJOioY^X4@*lkmTC%0n0Ki0kJ>Z`%cT#A{!DbgKfM0dHPYGl^~{^Ofu2j`)vyZ~GjF zJL0hgZ>;bBdj`CDO=<1e+4Z+G@QZRjSm+u7zv#vG08O(>6_YwXZ(dVcXANBVy&3R} zRI}?C76HGg{g|_vGiC^kH?JwJD^GW=y$1M26V}}@HwJ#u%C*T8_goRgzXRSq399Bq z`?@=r0l&!pzVC&xz%TM%ncy+$uaZd}pQ<_3**#Tr8mgA9o&fwJ<5!>VrvbmHX3y=I zInz~)cTa+8Mls+Qg-lxUurBb6IyNzn*Nsq#&wzJNg3?-+zu#WdwU#3JY0}bA zU@fOU*7kc3{~uDvr?jT_Bq*(AMGIpu9j~Pbjw!Ip0M_!O!MZ+;v;@YxCqZgG(4yar zo3#{I+68*Y0CTFnPsq(7&j0-#@a8pDb58d7`Me%51-)mzzcT`uf{!jfu5)wsbtZLu z-n^!2PEKR%_4dFNY_YwY&JbV<8oEA7bcndlc=MX7IUfU;6>R~g;GS9~pX&irP|x&B z*82C?#lHjIyr#64)ejvPJEg87&guJbE#O_d>pghe;=y$$b$s5urnHt_DRF$Usji~M zMTfFX;9XDZ+-LiU`YDVzuSu<~-~GK6TvyRF+R$n<@UCynPNelUVq? z{x=Wc7x~?YFUkXcQO%M~i!9-NAa#7+K8Kn~tZQF)O~W1cnTvna-hAK}`Rca(J`dgp zynPNelV~&K)CW`G7ftZ!_O&kf;d~CdXWjqb`+&F4p|o!C{_@gAz}5`+S=}oe*qR%C z)>_?yS|oLR-ad!Y+UZ2ofce1IToN`mwlAeZ=XYHecLY4%n;a`!=A4G zF&5aG^QV8AyzyTx;>~MnCNW~f%=lv9d&w+kEvYXqK zHPm_l^DdLdzW$X;x3sp|E#aNao7dD#;&_b4zR%!EBrDqBQE{tM>9Z))!|2~TnK!Q~ ztxY;iZ!-+ov$E0KI``>QuI&5nVPY$IA4naaH?JwJryYtv_!ro-GK-pauT0C8!**?R z83peH-n^!?PCHwg5en>CnM=2sZSBgH2NuPe*Bbxd`+&F4p=!<(@1LWG1Mk|SO6Q{{ z_#gaK?$dNymcpct&)ervH7EJdo~&BHyMC3_b9XxMt{v{HF0_kIVZ40~Rde3&%k3Qj zyz7^FyR$X`@49af7DK3n!nyph7Bj?dfY zP+H3}v}(2eqoHV>yESh(Fv`c>*`83YmCAVg98&9MuR6T!T~9H2`-O@WV3fx_)h{=g z0e=_p4{u&mHK+8t->my5;m)(i(HlN?B}&6P$}oHeQpe}bYpUi1|4KG>{tdr>QK{>P zIVDQH6{+&C@ClH-Ggsf;(TNv-og zJ-sjs;|Dcal6+TrJFNc1J6gD2XEyb1%`Y9;OwQA`Ar>Vo^g zs}MAvRExZdpy&60BCkU7_yiMSXGp8r7L_Tgg2~n*ucCLSJ$G=Ah+-nNubN~Ayb2-v zc%aCu7~I$Y8uBVACW4!e)p*P3RF!Uzvm&qJ@M!yZ+#`|^OoVY)rd$JFg|PRukI1X& zKfn8ZHLrrw+Wf>S*Mc|cs>~5LMSjsjmp0FF$DPu;$avan;1>zgD`G@`(VYIj4iCq_^{x&%9U73KGSujUyY7NV;*E8<<4$Q^8GoV9ycQ<~!@Ug`BEM)zi&pKBUqopg zaOLCW!)BSPBa6aCevyBa{VCjWr?hSv?EeV(MZ#pigCf7^t;?ApHNS{EuL&kX|IZ6k zfL9?{jXEy!Dr67je#onkKK}n^eNb9cOoVf8mTthS5W)lBh`frtu6s(6SHYdv1QVg+ zQn3Pf6~f3C5hAbRQ@2@3YF-6zUQ=3cG74=4{30RM@R`Uj3OoP2GxCe1j!$V#z zXR8bz48Z$mLPayb{&??eb zALtW`_s*!A({j$jlfbhRuJ&)`jXb+=xnY)So*kw25QlLg{s(haXATC5{s*sA!dzr3 zQ(8YxwbTax1K~)`UD5x*;ks)g`X5LgpVHdiVBbYQw>*{Bvy-C#!L6YBOGW)BgWPRriq!rGym?L4oQ!@qU4dsO90saR< zsz;dUf6!*t_Veg}Aa#7IhERD;X}z)}Aqe~rgda`HMgN2IJMtQ!{{e4aQ(EV){k8@C z4}=S=o{9bkWjC)BsQnKp_7nB)rI-lv7MJ$`uR_p|y(97}8ntiifV>Lnnidpe3$jUi@XXMZ=XXj5p*+;bp>99U{mL<$g2o>a`uFpS3zk_ z)f`G|?fg_D;1>yIL+*(DqE5`R6Y`6sj!$V#)f`IeUaoT{0>4Oj^Rl1FFKTmSu0HaM zc>5em>#HL*Vu4>I4F6Cj@{9I<4cn*Y7xCxy|7IeznAdYZ@G69yL$gF)MUnS4XXI5# z9-m-+P+C)G{S^+WgMn8eT>5Dv@+z`}iz1O%A>++!f{8Fsp1Kow6@p3d8If0^8MWwu znpeS_*Ob;PmhQU={30QIyBGKP!%q75~{ zvB)pt&1*{Q;L}}0fnOwCFgzymiym3n)>re3sG38)d#RdJ+Mqu2>;&VWwIa`M;IHq) zk!L4$e5&S9?_TQJY#EzX3OqYOr}DGNvs-mB`X2J^c>5fx=EPkbdmDImLYDz5k!N@3 zLVTQ>XGdvG^*NN*k$Qjcfd7FIXWHEd{SV%@uy;cL1F7RvT2p-vrL|?X1^OQdJBxfo z|AVEMo|~Zm0dJo}X?-;92>KrgZWcd9{{!v9?QUxS1Kzx*YEF*-i&Eg(37f};i99=t z5nk@dvy(bLZ(dV1$87HkH{jU`F)9D??2e{xib9?pZ(dV1=i;IkCxB-s-1adOd3L{T zrYuqO?0EB<(t2W_Ej_{iK-e^Ho9KUFRuXv^{STy$&zskj){B#un1cU-&{}td=zlP> z!Y~y54|wyM(%OI2;#}}Q5H{Ay7yS>uE^PN+?SDYcB#f0^vs0gCt1@b$M}QxiP<{XL ze)K`3W)j=ie0U6gnSyVJ4kyqrGx6eq81&1OIzBa%7!q>EuWdq(>Vxa>HXmSjvEt0~ zYV<*)zTM=HhV}x#OyO>WByIG|4DfmO`2qN&Q8Nki*^SHekL0SnhJP3RGBpOQSdKnu z)J)=tFnKiiWeT5*N<_cRF?E`k-<)(43$(LvQ4MG28O6vvMotN~Rn5mk)z-cvj^$N{rR1HFpUP^1n?jv`A zf3mQC_V!%#PyTT1`=>AeynA``ntE4j4k_ygewl(t_k`Z)m+ATAKqvIelsZ0dUQ_St z+yvVX;Fl?k_|YsJ{W9l&)v0y>k2>DGre+dXuKmsjzf7U$ej5e)Wj?#P{mZD!C&j-5 z-n^!?wiK=(0RLp6LC86O^iOVGwlxy{lckQ&o7a@qx#5*j;GZmbCSMg(Dw|#{E z$-H?@Y5mG~_8ag|7Vd1YeTM$Y5qjq=)c(m-&6yVbbi~n#DXPn}p5;ZtuFA!h{qo{7 zT^Ln!_C$0|Qovr(mw>tRYK&urkTaSK=O%4tQpcxij^DG`%fJ7ms5DFF9-9oiD&2jy z7!7!D&ZwHxHr`-I{nN0MlW4Fg1@?>97amAy=CYhoHRt8bHnx|prmDPd4c?XqyDHV& zzTPf2S}0I8$42L8%Xru;Y83pZg){ujSQ*|_x2J*lcR*<^>pP&q>Ae~%+4Zhn?ty}FFb{wD8kA$^giUTYW+w*Pp;F(cc%XZH!Xp;}|yZT;lbGk@rVBqM| zW92%5(t1~s#&&#GP35UCiW5^6kwH@qm%uZ_{UYAHrfQDQ&rf@8;ayzPpowc3?EKcM zcSUa-ypyGl&zsj&&2gPE+-f$wi?tnvH;lR7?cpF_g=;VV+$_^iQ^X75!lW^rpOh4yEEEKQN9zps{NCB^O~ATcood7 z1Adu8&cU0aUuN&jNv`OZDRq3_yryOn8*ZNa1%8>r+3j7u(JwR5!C?>jW%A}VHIq2E zA)*8LWeS_WHE=?|%oNvM_0@ivym?J&efGd;|P+Cv5STD}@YKI=|X>~YP zQE*q+3}<}U7xCsbRdb50{>Iva??JP>3kLtO3>B*P-R*q`cD$sH&zsj&&8cT-m1G0H z2m8A3?Jx)!Qty(R6zIT?7jIruHOF_?(J@+vnL=5euai3fLyGzOrs%aiD?S6>yr#5v z`Od_BypL_Bn`9lxRhaF|8Ham$Qpe}bYf9_yvyH@kJfoSzem>lsr>JYUa0c$>@#Zz9 zbz_?`;yzwSf4{$qmHCQK)(MHYmxq4`ly?Mmc26-8?mhf43wRYmLFZ*6ui}sF>1O0r zNFJZ?jv)6W2qwa$sW+X0S0OYV_CVxS41Rw%9C;No-aQF|iE!V}Eev=Sf@Sa|kyp_{ zzfY8!S3zk_omo;^FM8FiDe#Mg`S+)Z{GzdqKkh+(k<{@it*J9hO6#C&YwiQTNZ8hC zk;pHa7h5_6`9-{Y5|q|GoU_A#UnI<}*eLRgO2?S_srg0Rc}*}8M(1Sb1Fu3Dn0i{| zRXjNQO@_P*$>a0pHNiv}*;02E@G6Aqee!!Euj166l6d4*$awRbU?P}r8BhSc3SoC~ zpvbFOlv$FY=2h_KHKp}J(=Ro^FA`dJd8UK>BG+L9|C$27Nb2~!c};1(JFKuA_(ek9 z1m9uEFVY)UaTobTym?J&eIf49-_>w#;9hI8$S(@9ylbN77g05bI=iQ8POIRACcv{3 zBA1wpJiFP4dif#GPU`qn&7scjshZ=Yl`$N6cEZjXfg;at-k3!zkY~rcCqdPmuWaWh z;Moc5vpzLMo}K=LHyUc59i=rjlc2Qj6g^=Y_#X)2mHkBjgH?%tJ<$I^>iCq_)J%fX zy3pY69PmF7Y;Hc%MgM~tQGfH$wH zn&Ws$Zz=HXgk6D+^O0vaU4G~V^6aFJ&zsj&&3U*~&mDMnLi)koUdXdM8kd!fJUiaJ zrfSZBJs++C&rVoutn7_EyGV=UQ`I~>-n^!?9vbWu2mS{_*9j|&(f?q$%(4OcA4naa zH?JwJ%R=)9g8zX~cIL%>^gk$#Q!(^E;LU4F>#2DLTHt>mIA6W;6#WmzwRHV8?VtYv z#eSmRy%ZC{Z1lcSz^f2y-*EhZyo%H7$7&$2Lh|^$eGb7yux!yX4R{qoz1|MHkXP}s zg<~-CDrCHU4#7k?KhLTM@G6Ay?^7NluVTfd{&F?1g3_9*Ih59q{A(`(evxo+MPFUy z7s(APry{>d>iE2U4yASAx*3OoUnE49R!1YhXodN*cE~T{?QxpkC7SLc17}z-CB)P^NXmOL%n;cnqxV5(HG#^2_3ZdJV&0LQPf{M-v4!09O0?!V9|HC_%9epMD&Y$uKd3L;cP1T$=_I7)KXD7&t77WzSeLQ-BUA(?c<9-fnTOD zaNcUsFLTnV0Tt+%DRq2mCPAIuQ!|O;J+;v{L<&@YpBPlB3BXdE;e z4}O_K=Z=d+zf51PkOpeMOiF9&%#zajdU`X}GooMM9h$-H|Kl-8Sc*USa~WMR#mN1}i7V6AU;)c(o5c}=~m z-?|qSf?uX^!RAyd`eoL%+8ls>nNr8+&1-5V@nPLK1MtffLV8^2j((YSD@S!fzf9h| zre+cqcZ;>aFH>+_(=!zPGP~Ki|9YyE!0;%OH?JwJKiBTF9{iJqvgr+Gqkr;ukIm`m zpDcBJ-n^!?-mK_#1N@VPeWM;{p?`A6KHaCGe==`gQ(E81`#cH!lZ8U>l#S@0d~%m{ zf3<%yRdcAbd#dJ)aqv1~wkAucJ4CB>Y-F%-ZCZ4?E%1w^j!)Gb>g=AXIc+i@qznR1 z#k*a$Q6Ji!7rgEt+ExPmBHld-s^)~p%7%0WZjsT8&w?5Fqv=M>S5|DjB>o*xT2p%x zl-56vC(cLK)V>k2v7Pd>6-nXS+ltI7V&~=mReVZoYEOdFy2GxSB5TV2){{eBqjMA& z>%`PX<`nOq1f}&oi;*I0suoMWU*wmosBN_RATp=07V+jaRdX(X_nvb5WR9S(?O$64 zo~3oCeczS`OhKvR^X4^Gb54$Gwm%J6f;)74dQAtvYAOF~+x%11WWU-n^!? zHoFt-jQqz+t^66?12PnKtQ*)P2a-3hDXqW8>x%rxX77_$HZ;prG@PBd138fRcfi}{ zP&0|B1|}Q9FH_jqZ<8PTWtJD)6r*3J)bV-y9BL-fVe{a+;Fl?U)d;dfzf8Z+&0eBk zCU2iZ%_Qp2sr3~6GKGolj-;Vq=C9s^NuTY-PF@S(3)UGz_0?EUNn`X}@DIh59KrtEzM{>j4n+tV+hfAYX$ z`)_LhWZt}{W)kww%d5aI6K1o{eP*tBC*S;Perxp0lsZ0dUQ;uPpU-cGfnTN|g921G z>z%w$by+U@W%A}VHIq2+qJIVaGC`kl<#PEux$lD3ThxA;ym?J&J>1Dzf2dZjAUnOa z_MOvZ@-;f`wb4IW>iE2QO=-RT)tqhMpA64nko)O1W%85OC*sgQnK!Q~t!vd;J{$a# zVJ7&YM?8SdonlqcLuYc1D2>nl{1UrW9K-)Ovj4pnpVPSh&8Rs`{Duem&I@dfl0_H=k0SSt%ttuEV9Wx-R66*j!IQ*y8I(z!dol<`#a#xYf5Wb->m%l1@PZ6@Nu(Rs&Des zofmd#1hq)&_`G>dX`OLVXXk+#pi#!BKLaAa$!lFcm=^@Kh&Qh(t)Jh?3fTV`;-3x} z_2c9>`L=EiU)KIti;|xDKLPE4zvvz3KmVoDbYZI=sNH$XpP1M1@0S_zw9- zd2L?42mKJ(8>81Nf=p6W=l89=f;Fej__GV(eZXpNug!(`0ZZ0Rj>q?b)bU5;yhZKc zp2Nd0JDa{z#vb%2K<&aurpG{Dg9p$dv##;w%68ue^~bz!W$c4~2&|-5^TtamQdB#e z8nx^0v~vLD zHSA?Sc@)*|tMZuDun5d+bNg(pIqdX?OBC90(>dZPz7JUAs_Z2Y$iueZVSaC044x4|wxBdGr&UNwDaX zedYDQqvEBGe=g1>bPI+}g?SO%dvfIOzgtc){ka>rNLnTK~3O4DihtFv6 zsA#TH&jn`^6T7K!Uc|1v`PB{PMJ%-cNpW7po7bb#uj5RDjji)gUk^Mg8l)^diZcnl z{1luQu_IGnhr_&xX;r(rtLH_$c};1(u6rq-AFx&zJ!Zi90V^{ckcH<5Qpe}bYf9@| z4xjM+fc5S?FaypHSdmAq8F+rco7a@qjo&3qhVuhfc2Q#;oF6c+9lLC*o!|z_I<``H{+XV9>_Sx>DFU*Trl)=7tI4_bqK2>u9#*~M{ zcQ4fI=fxV8)k@p$v!3F2FI95_1yc=}7csw^NpoOc#0Cd%d4uyJs^-|sZLh+2FOz*4 zd2RgE8l~r?C=dMZrE1Ph_f`ckFJh~EpLqcDBKBzf-`VPU5v4Wd$w+D4`pCdruwTSB zx7{e%FJev$45MJbhy_g9ejoRXD6J_^MoQ~` zFWn7czlg229?=^1i&)^<@1xZFMZ9@U)f}6H$`^2+%$hepbsf%=nY{LaQ+S>%b$s5u zrfSZGOFBc~Jej?He5@;+C$n$KzB+iG%$wI#&55dgqamCpv!lNnUxf2y<~*%Qoq~Vo z$-H?@X+3`WiUhcSz!pr~@e=MIFwYg1&GG($)bV-qn$o&ujA}64KVXwDm)?c@2kh*= zdgJl_0dHPYS{uLb)e!C5!bfVaGH;(lX&rMt2KS5DYJ(uX@lcCSy&g3U_lu;C&)ervTCch(3y1w8 zwq@3~)39H}9L<(a#Qh@PK8Mn}d~G%E7cs$Z?BG#n;c=Nj3*X_teV38MF zdG-a*B6j%0JY*tNju{pQyb4y_*Q7u2Dp&{eU_InjNFARyuV3%-L?!|o^4@B5FYqin zw`_z9G7(;G^hRC<^J_J(HSj9fj%`{mkXON**AHWTkcq$!?r^$uVr9N+YKJLV$V6D) z=n?WNm|Jh3@xZHK$DKoOsd*K=c};1poK$Zt@QavpT9-KB7cuQ&U9yp1Bz1h=yr#7N zCAcgAei6I;sOw~MH>+<})$S-2{owDw$;C)c`vwf_ZU&PzzP&H>) z_U2O9H)T)P7tcMnTFDYFoc)RWrc%e}?Q^J_)3E(<+&5*vpR%{GZ_2v#@12MHro4R) zRdaf1=HR|5GwP8%&+xp81zg|POucW)+viYPH~)6b67DCn%k_V)>Y_?z=U?~r#{0=q z$LH;HD6N~?x3`7+$;_dBpbp$mW?@sVzsLK@ynPO(bCcKeRG z>ifyOc}>-vdf&^CXU8J$DcZM(XRs>U-W_>%Qpe}bYpUj?JCFDcJUdo#z<(F;>{!2U zPXmx=$D7ww%_-b4S{rzFZ1lagi-2dx{&rm2M$NP1&1*_)UyEKk;D5k6zvxy_^z;t< zTQw~M{STy$&zskj)&cFCn}PoUo4ob$%f;|B@?q;(^grOuYf9_1gAMb*|A6&ux6Eh( z{JfJoxLoajK+PmZjF=e@_iC8TBIu>9L!~mnX4wk7S3}JtUUzx?74D<4k0v?=a375s zWQQ!m`)E?fr)Cn_&27ryUJcBm zXUo6ieKcw&aXdz2AKa?}U+@hc6}KvtK8rFv@LmlylV~+@!3VgH#`0_?e}?;L?8V$K z*6RCcl-87IJEgULYqzCvKbb|$nV|#sliBHq?pAm|S?c(d)|6*ErS*ZnSx4c1GV8KH z2k$4dn;n{W#QVvV)|6*ErS({gjb?B^nbqGrwF2%ZvsX?PlhyZ=dGnf@Nu*C1xEy#D z?D*R)xxlMn(NBN9MP7x}@p`qV69Fs;8n2m7x!DMc@?~QO=*34jmB@_7cm#T>Bui)morY(GXs8+)bV-qn$r5q z^vZ9*FJkLsHIQG#*4rJ;KzM5(tTRgU-WshZO}rEwPUDwu!#%bmchU?=s{FCedisyX^Kf%DO(~`S$o#PrP?e)tuhDX7mPL1)H~GEb=N?C!0{F=2cKy%Lpa{rM0Yg-K}>kYbn;X zO3ph5&n|zLc|Uw+QpcyXmJv(@N^4pD=a2QYYAYVU>7O?op53p`8xzXmnNeEH2qprl z^@b38Z|&NO*1^^lDe&x;UF}kC0M86}ym<4PsyT~G8$So09b0g7I`ZsT^JA8EkY^`# zeBQjKYRPO=|7(Fnm~EjeLmE#%lB5RG6zB8U1%0Fje3m-adz#Nf@^5-VAtl%&Mbx zKz;b>b1XIhd3I9A=k0T-nMBNdqmjU~V-5EOcQAsV4_X_nN1h#TpF_)a~}#l>Plmc05#u~H?OIggw2E# zW#E^|9NQk+Sk>qud!wzl6a6x!j?bId)J($7vG*46%Vg%ab_F+rxsm6rxrykP$(z^I zOv13sJ@m_Dvu^ats|$0Z6+@S}s{Jx~^P1AS(S(7w!9STjo%`9V9h@T=d)_Zb|75A- z^X4_B_3H7pn}B~Z(_7{~Z6KT@%ztZp75$TW^P1B7$e_e9@K0v1i;U1enZ2(xX{z>5 z=IwK+n&ZDz-wt@!EY+b6@~+t~GdB<9T}vIGx6h$!&dvDxuE4uy`&+L^-ZguBetQh^ zu6g?$s^%1}JQWGNYi8B2=0RPU8wq)T;?%rr-ad!YTJ|Nk0Bbf2=_Bo{17Rwu#w)rAo=H~7l z1GB#Bs3y0Dz>E+3BHp~FYR>!2l!K+P+vV-o%g7wgl}iTZtXK*=UQ);B&1-SXS8Og64M;&5o?HSrnn z<~60Y%wt_d?4<(vKKHD-W^m>ov%E*BCh&`-j?bIdl-4qj4*ItfpU5XwZMyjm&itF+ z3yeGr{371GrnS~?a3nrY{&{!$U~@S04;kO~dSYCX_zZaWB({1kLeCXuHcxnM4gSiG z#~PhM&lUf7Gtm!+*>=?X3Vt}ucF(zT^uv)lKJT7HC$snHxx$=2)OC!A&r#*JdaBi| z-eET2Mi~0xuv-;B-hm$uTk>YyeDuTN-IGWfdm24gSe{;2&zpyHRT~C3-in?p6V4@~ z9}cTGse1(Y;jnubO9Rw?IJ|ojl-8ZnFQPvh%icdU1N_mLw{w_A1o)##9iMkkg3|g? zm-pz8#{9pQwgZ1OcFJpC-G$(f#=9p$Y5hXO1^v<3#t9FmfIk}B=WuJJ+8>QKuk#}s zqbDP~qKXZ(FHcvceh9pWo{Wiax1e7p>wT|kIrwF=J%>`Q(Jxc#_`G?o@!b?X8QDFb z(YbTL<9fC$YkCqq88`g7X$gLrY=ddec<{?)V@#i$qhBU(UI*uCp(i6NzUk3PWuB?Z zzVgKvy{&g-|2hPIne6NttB2s1$tEeix_$fSm&u#gl-7*}WAsmEO@dz>0RLpx`%~?Q z=$|ZgeBQjKwC>!075XQ$XEVMJ1OH^Utf9_F^iSr^Yf9^9i_%|%e=<9Ke)m)GPi9uF z6)V;L$-H|KRLz;#Sr`2>*+$(aS+(Hj(_U+Tp zeQQndLBCAiJqfDjxJ9f-zf4y6YCLTE((d?%?+EjvC%d*qnZmxQ)bY9Vn$Y_5#2u&l!MrHP*1=>5 z%!>lXG;383`=-2kO=;~uJ>G32%!}&yG*=kIyr|6zy*Ka1T@arEZ=VzY$R!Uv+gbVB zRfBWQQ&bxVrTjwAcGaKNR~2xMY-%>ruEsc)6}?Uio}09pNgbcJ&v{h;bs%`Qv(?8x z)r|U+qB_>&lJh|DY)|dJWm|nXPmb(U%P|Gclb5zunKXA<&UpJA8=FWk@N8#)%bh>I zxss|{vbN#oV(@I&+&!{IJe(t2dphx|9q?==N=uLV197^j}ca685 zgPv5!G&H*b^r#bRjpeS16ZV{A;P)`b}B*=|GzvE7u9`d7FJ^UrAMX zUAOVTzmsS--W>lvaAwDw*Ob;K{w2rqzN9N+LR7Et8*V^9kGuHoCUtz?yr#77HO1?5 zt5X?@&bz)(!EZRXJyuQd+l@D`DXkm)`8s)xd8T4@za&@uhO4`H!(RM$!)L(T=TJ4L z@|r9?63!1YTiPD8fVq)QSnU)$m={SMpSRDUYR;O}0X_b|+Rps1=I#6ACq>DWp&~=Z z%tRWf(=lWoi(DZ^kxD6Ksz{j=3T3!a8Vu1SWa#YbA}JD$icpjkl6y-TzH6Ot=kt^G z!~Ws@53EOC?`6IB*?YgA&ty11=o{A1qZsy$nuZ?vodx?vyz3ll<^)VMG2I5|2RnN9 z+};WHjfP)sqj0~)pUk_?p|Y-iZKP?p=(6O{{#IJ%uy@^XF#7gC@GX?k^{F)_wa%fk zUZvA@QzYzNf6PfcI2-n^^;3%t`op(Syz3k)YbX0wiGRc1_12JCTbf|+x{b@{@*yqm zUf#Z@X3nf0;p;n==9wP5)71#~RiE6kKa2aRLf7Z*Yij093Tuk70e5w$VZq&SUsb8B zrwZ<;^7b_~b7IX;d2IxDb<~eJGjLzk!Nt%^?y%s-fDRs`PUOWAK)<2QrS(rRb?R*<)^sjNNv{n#}NyvZp) zs}=B@!Cp26Pvrh&-g^?%o&@vNo;qPuie$O+%W=)#PbCXal$-X4Z$N~u&wEdT+LH+V zT2{NKom67~&ayckzDxU-myvS=z5(IACqeB=D7KsN;01h#Hs@mhk{<9~+A;TkZ)vqK zq9{#lhB@e#;p>b z!J9m{4JPmmiajPNmBKgTLf7Z*YwBG+#(l#QZFmOh1!bYj;Ti0=%5+YEZ^U`~n%a{v zF7F!d3eUhK+jXKXJcH7kVS(3@uFBS?ynRh&J#w~hm)Z3(lF4zd-OdBY=-R-O{WM_~ z30R`g4DE&^`R zg?-t{y?|ShuyDV{serXB^R$laqJwZ_1l!r#ffeND}roSpAee0KgVu?x1gybfGD z&)ScH7Qhq~x;}4TQ!{64S>)D=&Y_a?{+BB60@qHj-NerwfGNn^*VN4U?AylJ=t8KZ zAl&W#FYf}$vHUJ6j#^2wXTaOnRMsM&Kzp4%87Tp3OXp|bKQYo#lzm@IIK8@v^NUybp4ND|;;gev#1idDl7Ap2VQ;>!K!UMM<`dI(Tsa zybmlKyKi{~{370U4z(wty!rRZOYlBOKIpo?GrSM3+=(&&+~ULGUFT3)f3lqJ+z7L% z)$;KdZow>SZuRQeb6{%X%Q;>Zc%p#Q$A2lX|i-mWcLuGw3 zOXL16m_@T|y02G6#^wp*P|Fq<3vXXjdlKy@NbS0Ti%5ITDRdJ_M(Lj&gnTcd>+|+C zwI{J@_|gHI;35hiH(tj1N>&kfM!pwsUsHP$pZuh01>hp;Qr%m|`3f}h2$wsFc>9{l zdfnh7>q5bUPQ9+LZTc-J{pc5I^K-G+_lBv}48;y6JLuD-zML&0V6e0c=wCKuZV3MbIZR2_o80A9O=UwMeS&Kw! zzSq^g#DV^&ZHjj*JU`T9PkwPzYV0x38(0qj5idl^(b!3cn7NapCX%t$Y`G@VtFZ z&773{_mjqedm=A#y^ITQ*?QYWxr2hYuc@rT-Eg#9yDIV4OH*Ymf1kSUy^-NBbba2w zrm`0Kq=ziNRVj8fI4xuOzrM2|9~u6)o3O40uW;OCELJibgBdm{C9L-GmYB0wxYjwJ0-)k>Z;8noi-^Hk(<5TgI<0i()tDu+&f$ugy0bT{GADg`o zc@-U1M?@m8f?^^J9%OC}yb7p)d-7GoTybw5?VZT05D`oS~)6;#%gvz^L%=*KaRz%OE*PW2g#{364VzVne^By@c$Ys%S9Wv$U??+M@+v97V- z5|CdMXh_p7xse$ne2tuv8d#M{?Y)*tm=|K1po#@<|vn}z(MxBW9G%lSpr%!wa6}MWyrHj3kbJGo}JM3shOko;O0l**}>m` zV(+?;x8egu^ByD5j(Rq`Iykxm&yMx&n(2=`J58^~TIAVLGiQ^{3V-0)G0|&#-BCqv z#pT}jwjs}snmMQIe^vs|jw!0d)gjNW($dL7&aK5~XV=Ex^)>SBgs#uq*VN2O8D_Nucy>&s z=W##e*-c2in1Vbz-oB=0&VRYx?*Pw^wNZXN6nS=I7dS1D^Xz#0n#ww@%l2#FKY(?` z-U{>|Y-$_c5&Z{3*XQkPDr@Co4|TzRz#<4z=s)1?Ybxv3$A_ze|A2k- z$}K_v!Fta()qaUlU9O^SGBsfmgxuPA`%1DiS|Sosd@{czuHP zL1j%b5x$!r9S6J$CY@p}<5kqZE{#QAg^0Va|I9>CD~sF@yb4x5CQ8PuSk`N4pqy91 z+t*aq^YYu?2YwNAJS+7>eo@8NnwQ8g61qN>HPzQt)+L2i?SWs!%97_>AiqfM(ZXxU zFXHWMD(i5Qn9IN~VwwA+W&9%ddF|TD`9;*sq29gJ%<0*uup{v77^^*Z3VC)0Y8!)) zXD4)hYUWVyUh3J{Sk)T=&yFqb{iP@J?CSNeyCKhxcb!AcoVn&PHNdlDD+$O%6eke^=05cVD}bI z7>xddi)Z@0L;nHqI)}>Icg-z-@E@=aj>dD*f3R=Nv|PFWfVZ!ynG-nZRW0!B*a)S^ z+mUB?fVue~&raz2ynRj0oVVlFdjZdmjZnE4g*>}fh*3Zs%B^j{an! z>r+`%&UPwm{hg8B^4b?L?~DiSo`JhJY~7%6^e0nUQ_glO>vO-hz6F0WdwkwM9{tJD zQ_~yd{$$?1rry;b!XCZ?UncYE^6DS-Wd;|YK8LsnKQ{^XJA-b>J*%-h#g)~-)7=72w$sm^}ph5qEY zCvV2d{mInKITha7@5bsvn5S!%*TvR}Ki5B8e--YkP%~##!7aOUUkljro{MVo+XqNa zv`tIu0?$C``qa!RzF8h$@e=<2iQ5LK2G@y0LexiC!d(?==Ct<2?x~YrV!)hj@5q?bY-gJSU|| z+o-m`;eIATB6NK!>!%k7;r$3EdRp^qbZos?;YLs)-jkrR_J5jU`3jzyZJe9d(}lL-jkrRZtL3y@3gQX;ie&zW7DN>|HQ?} z@3!#vH8pcSY_pyB0N%;c$hfv5cqebV8@#s=?)(Z}pSQ26nUnAH{_qocCyR3eM=gSP za+2|l>1J@}m$$E}nWL`wpmSRuzW%}M8=7qjYq3iSZHMJ-4 zJ)-Xq@MW?OldKP*FSF$NzF_oa^7b{gC$Ug6R118WEMTnJTJ&Z5xioc^`!ad^n##KO z&lEH8C$l`Q_#5a??w@r21^SbPuFu=oRMw69L+rqx%yxXerH}sPx%0OlMSn7HUsGAn z@~(0Re=@Ua@9_lv$%cE&ljZ(o-gOQ&bJWr&1e?Kr(Uf*KG!tVZB^5sTh1T#rn$Y!m z*E!V8NjmSBHU{>KlomH9RCkY(oK&p!uYm8-c-J}9%!#PEsn-|wizbcRnQjLAMMIW+ zzW;B__h`K994hPSXn}tC9FEl+g8g*Ev+yqQ`CV8z@#f$hcp0 zT)xzD@wFcKEfnuMhst`RNo~$q*yD3upB(IPv_QIe=zz=eZ=rbmnwmM2uN_UafjeIF zhW8rx%i^*m`OHea3}8qJU7xqFshP9tc=_$6aL4O`-`HMb;Evb8k6+#m0)`ZCUsE$D zzF23fD%|l>P}(@3%vKk*SRCVh_|n)tak*JAh(Dqf1Z|E>{=w9u%&OloL@w_ zBdD``iiyymYKy!I_H1^Fj8|cLzqKv$Dul29XLkg3c26-8z8y|LUIh!<+bH8z{F|fo z5P20M5$`>TKQj@Yta3qK1v89&CF50W_dNW!oL50*O`Ta%Sv$suBEN`fx4t9e7cCh3 ztrzl(gsx9zO`Ta%S(jZgL4Fa7Y}a=U@{81q{u_z>BHnuvRMy|}YZ8E8#4Oj9%lJh; z3Z{qU{38Co{xcK7be9qXUIi;$<|E@(%wKor2=Xd~um5M(2bDF&M6iCVwjFpCOw}S= z#;b5PeN%(H3XzDjuL&l?8P9YR;8n07`(rX*MYUzrZ8@)kx38(JW21K=zleR9o-E@R z#eSaF8~H^-*Qc_k`kKnRF!w$3i`bsg#;(XO3Y_bwkNhIuzNWICkoFk)MJ#9SO&Pyv zS;Mn%IlqXSIn>!bHFHM1S!M=2JLa^(Q6G7BKicjck32h}>r*p_I=iQyP08({$g^Yr zxRlFycHcLcJwl!x?>z}>=CltxcOQ6m?ELfu8P6`F$Lu6I&yLEP+LNHN9-y*11N;X} zUsv@Y`VZ=@)LqbjAas2yYidt|%6g5>jXmH$U|)B-T|obV;iv1<(SN{uPlC$Y=~8n6 z_zzg;6R)z+e-QIy=PtSbfVZ!yne+T*Kjhi5GjA(oJiFb;sv1jR&raz2ynRj09JRK} z$g^YO1&hui&(11mkS+4;c>9`~Iok?vgaOZv?Y>Yhl`ZU&E3BL1N@l0R% z`hRA9P+3zokPu>7ylg!0iGS3 z#tLLSyPC1m`Es5el{K}_p|X|?J=G2T2TUpFkj#HDQ#m;i{RcwVr?RHjIaJm=-F)4_ zf50+U&hbb8!IG-PQRqM5UFT3)dk)_f2mS+=vNm7lKL`mpTq5@$@b)z|b6)liu!@JX zn-7QE+9A&_>PL|#p5X{xpSQ26ne*HzS}zICZcZ(_B;(mxJs78mXE?llP0gGaS6+SJ z1ZOwSUjt-3JIm0a$K+=?ynRh&ooB4C?Fr`wTO-{x(0`CHh|R_G1EK5l_BECDkpro5 z9&m1Oa-*AMo}ym9@o;#=t5#H~2Phb1U>881K@T%Fhp|JqhaUp4yX` z5^O#Pe3{HEA|(%fnH5`J$b6YX*QfR*sIz-&PvXGXF+t$VWT#hsS4CfDg6<8OFO&D4 z1hpsOxbV^{@MW?LTZJp=%RHd;>wmsXDr@S@lFE9@nW!4@C$owPSDvFk+21}=3H`}J z*Qc_k&Mc{{|NW~n1N_PChTXbe=uiHk7}#nxxEguyNl;mv++F_&{K@Rx26GGaCo6Y+ zVI=n_^Y%6MuCCTzGZcK8tZTFN5cFkgn#O3OFH`9HynRi*t5+81RDv&)Wu+#?pf58= zb@NYKcm}+EP3=iE);Jb`FO!`d_g##>Obwm-Z^kXYOy0hxvhF+Cr>`KcdRfpUm6WRMti(Md(jvQ&R?)qCa`S z-wzz*{$y(AP-pkl%(>TfOxb$i7pV+vZnOt}k@t$k#+Bel6S_V%bEva>YUb?J-5c?J z_anA@R^ew`;1~7mcGNKk{Aj%QB&eAay&}kR4DgEzd!LSfpcOAMN_7u=*y2Z{vZn4y zP+32B&_yQ~%eo%D@kB|UwB2Zx^XTRhx;~XPbx(rIy1q6Mom{M;gXfy6Yx&alhU40! zn~V3J1eLW|>WEG*_N=#J^Uvc2(s`F3`N-W|ynRj0oc(^j{?@=04B0%_JqDP9;()&w z`hl}l==!{UP0gHI#Vm6xUe zN0)?#gg%wJE9SRHUnXx~Q&~T~*=hrLGMRUr@4?2Ixzf<9&kxFdnY`;9YEL3et0W72 znao!6(i!w+D&3fI0ezW5*XLd5P-Akrp5vEC%gQ7vIYIgLf7YA=TKRP z)Q*b*e=^IojA%rEa?-qabJ3s7yUwAqHqtHK1^#6A)?s`B`jdOF%WWg~C-e3-wI@;1 z_pV|Byo;9t|i2h`cqJPcg z?*ra-4mET3Kh3{s3cTytE*>kJCM8OSRmUG+e#cWHbba1+4mEQQ8kX-;0p9f}?@4RY zfOnnkuQs5wn~Q{ZokPu>#BNb)Nx-}AD^cj`1iWjghsxTeF%=!| ztg)Yi!u*+O(#tOlRng@xbba1+4wZFp3!{_ZaAy()%Oi::read_checkpoint() setup_dofs(); std::vector x_system(1 + previous_solutions.size()); - VectorType distributed_system(locally_owned_dofs, this->mpi_communicator); + VectorType distributed_system(locally_owned_dofs, + this->locally_relevant_dofs, + this->mpi_communicator); x_system[0] = &(distributed_system); std::vector distributed_previous_solutions; @@ -1584,7 +1586,9 @@ NavierStokesBase::read_checkpoint() for (unsigned int i = 0; i < previous_solutions.size(); ++i) { distributed_previous_solutions.emplace_back( - VectorType(locally_owned_dofs, this->mpi_communicator)); + VectorType(locally_owned_dofs, + this->locally_relevant_dofs, + this->mpi_communicator)); x_system[i + 1] = &distributed_previous_solutions[i]; } parallel::distributed::SolutionTransfer system_trans_vectors( @@ -2336,15 +2340,8 @@ NavierStokesBase::init_temporary_vector() { VectorType tmp; - if constexpr (std::is_same_v || - std::is_same_v) - tmp.reinit(locally_owned_dofs, this->mpi_communicator); + tmp.reinit(locally_owned_dofs, locally_relevant_dofs, this->mpi_communicator); - else if constexpr (std::is_same_v>) - tmp.reinit(locally_owned_dofs, - locally_relevant_dofs, - this->mpi_communicator); return tmp; }