dualDynamics.cse.h

Go to the documentation of this file.

/*  This file is part of the OpenLB library
 *
 *  Copyright (C) 2022 Adrian Kummerlaender
 *  E-mail contact: info@openlb.net
 *  The most recent release of OpenLB can be downloaded at
 *  <http://www.openlb.net/>
 *
 *  This program is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU General Public License
 *  as published by the Free Software Foundation; either version 2
 *  of the License, or (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public
 *  License along with this program; if not, write to the Free
 *  Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
 *  Boston, MA  02110-1301, USA.
*/
 
#ifndef DUAL_DYNAMICS_CSE_H
#define DUAL_DYNAMICS_CSE_H
 
 
#ifndef DISABLE_CSE
 
#include "dualDynamics.h"
#include "dynamics/latticeDescriptors.h"
 
namespace olb {
 
namespace collision {
 
template <typename... FIELDS>
struct DualPorousBGK::type<descriptors::D3Q19<FIELDS...>,momenta::BulkTuple,equilibria::SecondOrder> {
 
template <CONCEPT(MinimalCell) CELL, typename PARAMETERS, typename V=typename CELL::value_t>
CellStatistic<V> apply(CELL& cell, PARAMETERS& parameters) any_platform
{
auto x22 = cell.template getFieldComponent<olb::descriptors::F>(2);
auto x20 = cell.template getFieldComponent<olb::descriptors::F>(0);
auto x32 = cell.template getFieldComponent<olb::descriptors::F>(12);
auto x49 = cell.template getFieldComponent<olb::descriptors::DJDF>(10);
auto x50 = cell.template getFieldComponent<olb::descriptors::DJDF>(11);
auto x38 = cell.template getFieldComponent<olb::descriptors::F>(18);
auto x27 = cell.template getFieldComponent<olb::descriptors::F>(7);
auto x35 = cell.template getFieldComponent<olb::descriptors::F>(15);
auto x44 = cell.template getFieldComponent<olb::descriptors::DJDF>(5);
auto x53 = cell.template getFieldComponent<olb::descriptors::DJDF>(14);
auto x42 = cell.template getFieldComponent<olb::descriptors::DJDF>(3);
auto x30 = cell.template getFieldComponent<olb::descriptors::F>(10);
auto x37 = cell.template getFieldComponent<olb::descriptors::F>(17);
auto x48 = cell.template getFieldComponent<olb::descriptors::DJDF>(9);
auto x46 = cell.template getFieldComponent<olb::descriptors::DJDF>(7);
auto x51 = cell.template getFieldComponent<olb::descriptors::DJDF>(12);
auto x43 = cell.template getFieldComponent<olb::descriptors::DJDF>(4);
auto x23 = cell.template getFieldComponent<olb::descriptors::F>(3);
auto x54 = cell.template getFieldComponent<olb::descriptors::DJDF>(15);
auto x56 = cell.template getFieldComponent<olb::descriptors::DJDF>(17);
auto x33 = cell.template getFieldComponent<olb::descriptors::F>(13);
auto x25 = cell.template getFieldComponent<olb::descriptors::F>(5);
auto x39 = cell.template getFieldComponent<olb::descriptors::DJDF>(0);
auto x47 = cell.template getFieldComponent<olb::descriptors::DJDF>(8);
auto x28 = cell.template getFieldComponent<olb::descriptors::F>(8);
auto x29 = cell.template getFieldComponent<olb::descriptors::F>(9);
auto x36 = cell.template getFieldComponent<olb::descriptors::F>(16);
auto x21 = cell.template getFieldComponent<olb::descriptors::F>(1);
auto x45 = cell.template getFieldComponent<olb::descriptors::DJDF>(6);
auto x55 = cell.template getFieldComponent<olb::descriptors::DJDF>(16);
auto x26 = cell.template getFieldComponent<olb::descriptors::F>(6);
auto x24 = cell.template getFieldComponent<olb::descriptors::F>(4);
auto x40 = cell.template getFieldComponent<olb::descriptors::DJDF>(1);
auto x31 = cell.template getFieldComponent<olb::descriptors::F>(11);
auto x19 = cell.template getFieldComponent<olb::descriptors::POROSITY>(0);
auto x34 = cell.template getFieldComponent<olb::descriptors::F>(14);
auto x41 = cell.template getFieldComponent<olb::descriptors::DJDF>(2);
auto x57 = cell.template getFieldComponent<olb::descriptors::DJDF>(18);
auto x61 = parameters.template get<olb::descriptors::OMEGA>();
auto x52 = cell.template getFieldComponent<olb::descriptors::DJDF>(13);
auto x58 = x22 + x34;
auto x59 = x23 + x36;
auto x60 = x20 + x21 + x24 + x25 + x26 + x27 + x28 + x29 + x30 + x31 + x32 + x33 + x35 + x37 + x38 + x58 + x59;
auto x62 = x60 + V{1};
auto x63 = V{1} / (x62);
auto x64 = x63*(x60 + V{1});
auto x65 = V{0.0277777777777778}*cell[5];
auto x66 = x19*x63;
auto x67 = V{3}*x66;
auto x68 = x24 - x33;
auto x69 = x26 - x35;
auto x70 = x68 + x69;
auto x71 = x27 - x36;
auto x72 = -x30;
auto x73 = x21 + x25 - x34 + x72;
auto x74 = x70 + x71 + x73;
auto x75 = -x74;
auto x76 = x66*x75;
auto x77 = x76 + V{-1};
auto x78 = -x77;
auto x79 = x75*x78;
auto x80 = -x79;
auto x81 = x28 - x37;
auto x82 = x68 + x81;
auto x83 = -x31;
auto x84 = x29 - x38;
auto x85 = x83 + x84;
auto x86 = -x25 + x58;
auto x87 = x82 + x85 + x86;
auto x88 = -x87;
auto x89 = x66*x88;
auto x90 = x89 + V{1};
auto x91 = x88*x90;
auto x92 = x69 + x81;
auto x93 = -x32;
auto x94 = -x29 + x38;
auto x95 = x93 + x94;
auto x96 = -x27 + x59;
auto x97 = x92 + x95 + x96;
auto x98 = -x97;
auto x99 = x98*x98;
auto x100 = x66*x99;
auto x101 = x100 + x91;
auto x102 = V{3}*x29;
auto x103 = V{3}*x34;
auto x104 = V{3}*x25;
auto x105 = V{3}*x38;
auto x106 = V{3}*x24 - V{3}*x33;
auto x107 = V{3}*x28 - V{3}*x37;
auto x108 = x102 + x103 - x104 - x105 + x106 + x107 + V{3}*x22 - V{3}*x31;
auto x109 = -x108*x63;
auto x110 = V{1} / ((x62)*(x62));
auto x111 = V{4.5}*x110;
auto x112 = x21 + x71 + x72;
auto x113 = -x28 + x37;
auto x114 = x112 + x113 - x22 + V{2}*x25 + x31 - V{2}*x34 + x69 + x94;
auto x115 = -x114;
auto x116 = V{3}*x27;
auto x117 = V{3}*x36;
auto x118 = V{3}*x26 - V{3}*x35;
auto x119 = -x103 + x104 + x106 + x116 - x117 + x118 + V{3}*x21 - V{3}*x30;
auto x120 = -x119*x63;
auto x121 = V{1.5}*x110;
auto x122 = x75*x75;
auto x123 = x121*x122;
auto x124 = x88*x88;
auto x125 = x121*x124;
auto x126 = x121*x99;
auto x127 = x125 + x126 + V{-1};
auto x128 = x123 + x127;
auto x129 = -x120 + x128;
auto x130 = x109 - x111*x115*x115 + x129;
auto x131 = V{0.0277777777777778}*cell[7];
auto x132 = x66*x98;
auto x133 = x132 + V{1};
auto x134 = x133*x98;
auto x135 = x124*x66;
auto x136 = x134 + x135;
auto x137 = -x102 + x105 + x107 - x116 + x117 + x118 + V{3}*x23 - V{3}*x32;
auto x138 = -x137*x63;
auto x139 = -x23 + x32 + x68;
auto x140 = x113 + x139 + V{2}*x27 - V{2}*x36 + x73 + x84;
auto x141 = -x140;
auto x142 = -x111*x141*x141 + x129 + x138;
auto x143 = V{0.0277777777777778}*cell[9];
auto x144 = x89 + V{-1};
auto x145 = -x144;
auto x146 = x145*x88;
auto x147 = -x146;
auto x148 = x122*x66;
auto x149 = x134 + x148;
auto x150 = x83 + x86;
auto x151 = x139 + x150 - x26 + V{2}*x29 + x35 - V{2}*x38 + x71;
auto x152 = -x151;
auto x153 = -x109 + x128;
auto x154 = -x111*x152*x152 + x138 + x153;
auto x155 = V{0.0277777777777778}*cell[13];
auto x156 = x76 + V{1};
auto x157 = x156*x75;
auto x158 = x112 + x22 + V{2}*x24 - V{2}*x33 + x85 + x92;
auto x159 = -x158;
auto x160 = x111*(x159*x159);
auto x161 = x120 + x128;
auto x162 = x109 - x160 + x161;
auto x163 = V{0.0277777777777778}*cell[14];
auto x164 = x100 + x157;
auto x165 = -x111*x114*x114 + x120 + x153;
auto x166 = V{0.0277777777777778}*cell[15];
auto x167 = x23 + V{2}*x26 - V{2}*x35 + x73 + x82 + x95;
auto x168 = -x167;
auto x169 = x111*(x168*x168);
auto x170 = x138 + x161 - x169;
auto x171 = V{0.0277777777777778}*cell[16];
auto x172 = x132 + V{-1};
auto x173 = -x172;
auto x174 = x173*x98;
auto x175 = -x174;
auto x176 = x135 + x157;
auto x177 = -x138;
auto x178 = -x111*x140*x140 + x161 + x177;
auto x179 = V{0.0277777777777778}*cell[17];
auto x180 = x150 + V{2}*x28 - V{2}*x37 + x70 + x93 + x96;
auto x181 = -x180;
auto x182 = x111*(x181*x181);
auto x183 = x109 + x128;
auto x184 = x138 - x182 + x183;
auto x185 = V{0.0277777777777778}*cell[18];
auto x186 = x148 + x91;
auto x187 = -x111*x151*x151 + x177 + x183;
auto x188 = V{0.0555555555555556}*cell[10];
auto x189 = x156*x74;
auto x190 = x87*x87;
auto x191 = x190*x66;
auto x192 = x97*x97;
auto x193 = x192*x66;
auto x194 = x191 + x193;
auto x195 = V{3}*x110;
auto x196 = x122*x195;
auto x197 = x120 + x127 - x196;
auto x198 = V{0.0555555555555556}*cell[11];
auto x199 = x87*x90;
auto x200 = x74*x74;
auto x201 = x200*x66;
auto x202 = x193 + x201;
auto x203 = x124*x195;
auto x204 = x123 + V{-1};
auto x205 = x109 + x126 - x203 + x204;
auto x206 = V{0.0555555555555556}*cell[12];
auto x207 = x133*x97;
auto x208 = x191 + x201;
auto x209 = x195*x99;
auto x210 = x125 + x138 + x204 - x209;
auto x211 = V{0.333333333333333}*cell[0];
auto x212 = x19*x19;
auto x213 = -x100;
auto x214 = x146 + x213;
auto x215 = x121*x190;
auto x216 = x121*x192 + V{-1};
auto x217 = x215 + x216;
auto x218 = x119*x63 + x217;
auto x219 = x121*x200;
auto x220 = x108*x63 + x219;
auto x221 = V{4.5}*x110*(x158*x158) - x218 - x220;
auto x222 = -x135;
auto x223 = x174 + x222;
auto x224 = x137*x63;
auto x225 = x219 + x224;
auto x226 = V{4.5}*x110*(x167*x167) - x218 - x225;
auto x227 = -x148;
auto x228 = x146 + x174;
auto x229 = V{4.5}*x110*(x180*x180) - x217 - x220 - x224;
auto x230 = x74*x77;
auto x231 = V{3}*x110*x200 - x218;
auto x232 = x144*x87;
auto x233 = V{3}*x110*x190 - x216 - x220;
auto x234 = x172*x97;
auto x235 = V{3}*x110*x192 - x215 - x225 + V{1};
auto x236 = x39 + x61*(x64*(-x128*x211*(x195*x212*(x190 + x192 + x200) + V{1}) - x130*x65*(-x67*(-x101 - x80) + V{1}) - x131*x142*(-x67*(-x136 - x80) + V{1}) - x143*x154*(-x67*(-x147 - x149) + V{1}) - x155*x162*(x67*(x101 + x157) + V{1}) - x163*x165*(-x67*(-x147 - x164) + V{1}) - x166*x170*(x67*(x136 + x157) + V{1}) - x171*x178*(-x67*(-x175 - x176) + V{1}) - x179*x184*(x67*(x149 + x91) + V{1}) - x185*x187*(-x67*(-x175 - x186) + V{1}) - x188*x197*(x67*(-x189 + x194) + V{1}) - x198*x205*(x67*(-x199 + x202) + V{1}) - x206*x210*(x67*(-x207 + x208) + V{1}) + V{0.0277777777777778}*x221*(-x67*(x214 + x79) + V{1})*cell[4] + V{0.0277777777777778}*x226*(-x67*(x223 + x79) + V{1})*cell[6] + V{0.0277777777777778}*x229*(-x67*(x227 + x228) + V{1})*cell[8] + V{0.0555555555555556}*x231*(x67*(x194 - x230) + V{1})*cell[1] + V{0.0555555555555556}*x233*(x67*(x202 - x232) + V{1})*cell[2] + V{0.0555555555555556}*x235*(x67*(x208 - x234) + V{1})*cell[3]) - cell[0]) + cell[0];
auto x237 = -cell[1];
auto x238 = x63*x75;
auto x239 = x238 + V{-1};
auto x240 = -x239;
auto x241 = x240*x75;
auto x242 = x19*x241;
auto x243 = -x134;
auto x244 = x242 + x243;
auto x245 = x146 + x244;
auto x246 = -x91;
auto x247 = x174 + x242;
auto x248 = x246 + x247;
auto x249 = -x242;
auto x250 = x134 + x249;
auto x251 = x250 + x91;
auto x252 = V{3}*x19;
auto x253 = -x199*x63;
auto x254 = x239*x77;
auto x255 = x110*x19;
auto x256 = x192*x255;
auto x257 = x254 + x256;
auto x258 = -x207*x63;
auto x259 = x190*x255;
auto x260 = x254 + x259;
auto x261 = x156*x239;
auto x262 = x256 + x261;
auto x263 = -x232*x63;
auto x264 = x259 + x261;
auto x265 = -x234*x63;
auto x266 = x101 + x249;
auto x267 = x136 + x249;
auto x268 = V{3}*x63;
auto x269 = x239*x74;
auto x270 = x190*x63;
auto x271 = x192*x63;
auto x272 = x270 + x271;
auto x273 = -x67*(x228 + x242) + V{1};
auto x274 = x214 + x242;
auto x275 = x223 + x242;
auto x276 = x49 + x61*(x237 + x64*(-x128*x211*(x212*x268*(-x269 + x272) + V{1}) - x130*x65*(x252*(x253 + x257) + V{1}) - x131*x142*(x252*(x258 + x260) + V{1}) - x143*x154*(-x245*x67 + V{1}) - x155*x162*(x252*(x253 + x262) + V{1}) - x163*x165*(x252*(x262 + x263) + V{1}) - x166*x170*(x252*(x258 + x264) + V{1}) - x171*x178*(x252*(x264 + x265) + V{1}) - x179*x184*(x251*x67 + V{1}) - x185*x187*(-x248*x67 + V{1}) - x188*x197*(x252*(x259 + x262) + V{1}) - x198*x205*(x266*x67 + V{1}) - x206*x210*(x267*x67 + V{1}) + V{0.0277777777777778}*x221*(x252*(x257 + x263) + V{1})*cell[4] + V{0.0277777777777778}*x226*(x252*(x260 + x265) + V{1})*cell[6] + V{0.0277777777777778}*x229*x273*cell[8] + V{0.0555555555555556}*x231*(x252*(x257 + x259) + V{1})*cell[1] + V{0.0555555555555556}*x233*(-x274*x67 + V{1})*cell[2] + V{0.0555555555555556}*x235*(-x275*x67 + V{1})*cell[3])) + cell[1];
auto x277 = -cell[2];
auto x278 = -x157;
auto x279 = x63*x88;
auto x280 = x279 + V{-1};
auto x281 = -x280;
auto x282 = x281*x88;
auto x283 = x19*x282;
auto x284 = x174 + x283;
auto x285 = x278 + x284;
auto x286 = x283 + x79;
auto x287 = x243 + x286;
auto x288 = -x283;
auto x289 = x157 + x288;
auto x290 = x134 + x289;
auto x291 = -x230*x63;
auto x292 = x280*x90;
auto x293 = x256 + x292;
auto x294 = x144*x280;
auto x295 = x200*x255;
auto x296 = x294 + x295;
auto x297 = -x189*x63;
auto x298 = x256 + x294;
auto x299 = x292 + x295;
auto x300 = x164 + x288;
auto x301 = x149 + x288;
auto x302 = x280*x87;
auto x303 = x200*x63;
auto x304 = x271 + x303;
auto x305 = -x67*(x174 + x286) + V{1};
auto x306 = x213 + x286;
auto x307 = x227 + x284;
auto x308 = x50 + x61*(x277 + x64*(-x128*x211*(x212*x268*(-x302 + x304) + V{1}) - x130*x65*(x252*(x291 + x293) + V{1}) - x131*x142*(-x287*x67 + V{1}) - x143*x154*(x252*(x258 + x296) + V{1}) - x155*x162*(x252*(x293 + x297) + V{1}) - x163*x165*(x252*(x297 + x298) + V{1}) - x166*x170*(x290*x67 + V{1}) - x171*x178*(-x285*x67 + V{1}) - x179*x184*(x252*(x258 + x299) + V{1}) - x185*x187*(x252*(x265 + x299) + V{1}) - x188*x197*(x300*x67 + V{1}) - x198*x205*(x252*(x293 + x295) + V{1}) - x206*x210*(x301*x67 + V{1}) + V{0.0277777777777778}*x221*(x252*(x291 + x298) + V{1})*cell[4] + V{0.0277777777777778}*x226*x305*cell[6] + V{0.0277777777777778}*x229*(x252*(x265 + x296) + V{1})*cell[8] + V{0.0555555555555556}*x231*(-x306*x67 + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x295 + x298) + V{1})*cell[2] + V{0.0555555555555556}*x235*(-x307*x67 + V{1})*cell[3])) + cell[2];
auto x309 = -cell[3];
auto x310 = x63*x98;
auto x311 = x310 + V{-1};
auto x312 = -x311;
auto x313 = x312*x98;
auto x314 = x19*x313;
auto x315 = x146 + x314;
auto x316 = x278 + x315;
auto x317 = x314 + x79;
auto x318 = x246 + x317;
auto x319 = -x314;
auto x320 = x157 + x91;
auto x321 = x319 + x320;
auto x322 = x133*x311;
auto x323 = x259 + x322;
auto x324 = x263 + x295;
auto x325 = x172*x311;
auto x326 = x259 + x325;
auto x327 = x253 + x295;
auto x328 = x176 + x319;
auto x329 = x186 + x319;
auto x330 = x311*x97;
auto x331 = x270 + x303;
auto x332 = -x67*(x146 + x317) + V{1};
auto x333 = x222 + x317;
auto x334 = x227 + x315;
auto x335 = x51 + x61*(x309 + x64*(-x128*x211*(x212*x268*(-x330 + x331) + V{1}) - x130*x65*(-x318*x67 + V{1}) - x131*x142*(x252*(x291 + x323) + V{1}) - x143*x154*(x252*(x322 + x324) + V{1}) - x155*x162*(x321*x67 + V{1}) - x163*x165*(-x316*x67 + V{1}) - x166*x170*(x252*(x297 + x323) + V{1}) - x171*x178*(x252*(x297 + x326) + V{1}) - x179*x184*(x252*(x322 + x327) + V{1}) - x185*x187*(x252*(x325 + x327) + V{1}) - x188*x197*(x328*x67 + V{1}) - x198*x205*(x329*x67 + V{1}) - x206*x210*(x252*(x295 + x323) + V{1}) + V{0.0277777777777778}*x221*x332*cell[4] + V{0.0277777777777778}*x226*(x252*(x291 + x326) + V{1})*cell[6] + V{0.0277777777777778}*x229*(x252*(x324 + x325) + V{1})*cell[8] + V{0.0555555555555556}*x231*(-x333*x67 + V{1})*cell[1] + V{0.0555555555555556}*x233*(-x334*x67 + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x295 + x326) + V{1})*cell[3])) + cell[3];
auto x336 = -cell[4];
auto x337 = x281*x90;
auto x338 = x174*x63;
auto x339 = x241*x66;
auto x340 = x338 + x339;
auto x341 = x185*(-x252*(x337 + x340) + V{1});
auto x342 = x134*x63;
auto x343 = -x342;
auto x344 = x339 + x343;
auto x345 = x337 + x344;
auto x346 = x156*x240;
auto x347 = x282*x66;
auto x348 = x338 + x347;
auto x349 = x171*(-x252*(x346 + x348) + V{1});
auto x350 = x343 + x347;
auto x351 = x346 + x350;
auto x352 = x145*x281;
auto x353 = -x339;
auto x354 = x342 + x353;
auto x355 = x352 + x354;
auto x356 = -x347;
auto x357 = x240*x78;
auto x358 = x342 + x357;
auto x359 = x356 + x358;
auto x360 = x244 + x283;
auto x361 = -x302*x66;
auto x362 = -x269*x66;
auto x363 = x63*x99;
auto x364 = x241 + x282 - x363;
auto x365 = -x352;
auto x366 = x340 + x365;
auto x367 = -x357;
auto x368 = x348 + x367;
auto x369 = -x67*(x247 + x283) + V{1};
auto x370 = x52 + x61*(x336 + x64*(-x128*x211*(-x212*x268*x364 + V{1}) - x130*x65*(x252*(x257 + x292) + V{1}) - x131*x142*(x252*x359 + V{1}) - x143*x154*(x252*x355 + V{1}) - x155*x162*(x252*(x262 + x292) + V{1}) - x163*x165*(x252*(x262 + x294) + V{1}) - x166*x170*(-x252*x351 + V{1}) - x178*x349 - x179*x184*(-x252*x345 + V{1}) - x187*x341 - x188*x197*(x252*(x262 + x361) + V{1}) - x198*x205*(x252*(x293 + x362) + V{1}) - x206*x210*(-x360*x67 + V{1}) + V{0.0277777777777778}*x221*(x252*(x257 + x294) + V{1})*cell[4] + V{0.0277777777777778}*x226*(-x252*x368 + V{1})*cell[6] + V{0.0277777777777778}*x229*(-x252*x366 + V{1})*cell[8] + V{0.0555555555555556}*x231*(x252*(x257 + x361) + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x298 + x362) + V{1})*cell[2] + V{0.0555555555555556}*x235*x369*cell[3])) + cell[4];
auto x371 = -cell[5];
auto x372 = x279 + V{1};
auto x373 = x372*x88;
auto x374 = x373*x66;
auto x375 = x131*(x252*(x358 + x374) + V{1});
auto x376 = x145*x372;
auto x377 = x344 + x376;
auto x378 = -x374;
auto x379 = x338 + x378;
auto x380 = x346 + x379;
auto x381 = -x346;
auto x382 = x342 + x374;
auto x383 = x381 + x382;
auto x384 = x372*x90;
auto x385 = -x384;
auto x386 = x340 + x385;
auto x387 = x354 + x384;
auto x388 = x144*x372;
auto x389 = x19*x373;
auto x390 = x250 + x389;
auto x391 = x372*x87;
auto x392 = -x391*x66;
auto x393 = x256 + x362;
auto x394 = -x241;
auto x395 = x363 + x373;
auto x396 = x394 + x395;
auto x397 = -x252*(x340 + x376) + V{1};
auto x398 = -x338;
auto x399 = x374 + x398;
auto x400 = x357 + x399;
auto x401 = -x389;
auto x402 = x247 + x401;
auto x403 = x53 + cell[5];
auto x404 = x403 + x61*(x371 + x64*(-x128*x211*(x212*x268*x396 + V{1}) - x130*x65*(x252*(x257 + x384) + V{1}) - x142*x375 - x143*x154*(-x252*x377 + V{1}) - x155*x162*(x252*(x262 + x384) + V{1}) - x163*x165*(x252*(x262 + x388) + V{1}) - x166*x170*(x252*x383 + V{1}) - x171*x178*(-x252*x380 + V{1}) - x179*x184*(x252*x387 + V{1}) - x185*x187*(-x252*x386 + V{1}) - x188*x197*(x252*(x262 + x392) + V{1}) - x198*x205*(x252*(x384 + x393) + V{1}) - x206*x210*(x390*x67 + V{1}) + V{0.0277777777777778}*x221*(x252*(x257 + x388) + V{1})*cell[4] + V{0.0277777777777778}*x226*(x252*x400 + V{1})*cell[6] + V{0.0277777777777778}*x229*x397*cell[8] + V{0.0555555555555556}*x231*(x252*(x257 + x392) + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x388 + x393) + V{1})*cell[2] + V{0.0555555555555556}*x235*(-x402*x67 + V{1})*cell[3]));
auto x405 = -cell[6];
auto x406 = x133*x312;
auto x407 = x146*x63;
auto x408 = x339 + x407;
auto x409 = x143*(-x252*(x406 + x408) + V{1});
auto x410 = x63*x91;
auto x411 = -x410;
auto x412 = x339 + x411;
auto x413 = x406 + x412;
auto x414 = x313*x66;
auto x415 = x407 + x414;
auto x416 = x163*(-x252*(x346 + x415) + V{1});
auto x417 = x411 + x414;
auto x418 = x346 + x417;
auto x419 = x173*x312;
auto x420 = x353 + x410;
auto x421 = x419 + x420;
auto x422 = -x414;
auto x423 = x357 + x410;
auto x424 = x422 + x423;
auto x425 = x242 + x246 + x314;
auto x426 = -x330*x66;
auto x427 = x124*x63;
auto x428 = x241 + x313 - x427;
auto x429 = -x419;
auto x430 = x408 + x429;
auto x431 = x367 + x415;
auto x432 = -x67*(x242 + x315) + V{1};
auto x433 = x54 + x61*(x405 + x64*(-x128*x211*(-x212*x268*x428 + V{1}) - x130*x65*(x252*x424 + V{1}) - x131*x142*(x252*(x260 + x322) + V{1}) - x154*x409 - x155*x162*(-x252*x418 + V{1}) - x165*x416 - x166*x170*(x252*(x264 + x322) + V{1}) - x171*x178*(x252*(x264 + x325) + V{1}) - x179*x184*(-x252*x413 + V{1}) - x185*x187*(x252*x421 + V{1}) - x188*x197*(x252*(x264 + x426) + V{1}) - x198*x205*(-x425*x67 + V{1}) - x206*x210*(x252*(x323 + x362) + V{1}) + V{0.0277777777777778}*x221*(-x252*x431 + V{1})*cell[4] + V{0.0277777777777778}*x226*(x252*(x260 + x325) + V{1})*cell[6] + V{0.0277777777777778}*x229*(-x252*x430 + V{1})*cell[8] + V{0.0555555555555556}*x231*(x252*(x260 + x426) + V{1})*cell[1] + V{0.0555555555555556}*x233*x432*cell[2] + V{0.0555555555555556}*x235*(x252*(x326 + x362) + V{1})*cell[3])) + cell[6];
auto x434 = -cell[7];
auto x435 = x310 + V{1};
auto x436 = x435*x98;
auto x437 = x436*x66;
auto x438 = x65*(x252*(x423 + x437) + V{1});
auto x439 = x173*x435;
auto x440 = x412 + x439;
auto x441 = -x437;
auto x442 = x407 + x441;
auto x443 = x346 + x442;
auto x444 = x410 + x437;
auto x445 = x381 + x444;
auto x446 = x133*x435;
auto x447 = -x446;
auto x448 = x408 + x447;
auto x449 = x420 + x446;
auto x450 = x172*x435;
auto x451 = x19*x436;
auto x452 = x451 + x91;
auto x453 = x249 + x452;
auto x454 = x435*x97;
auto x455 = -x454*x66;
auto x456 = x259 + x362;
auto x457 = x427 + x436;
auto x458 = x394 + x457;
auto x459 = -x252*(x408 + x439) + V{1};
auto x460 = -x407;
auto x461 = x437 + x460;
auto x462 = x357 + x461;
auto x463 = -x451;
auto x464 = x146 + x463;
auto x465 = x242 + x464;
auto x466 = x55 + cell[7];
auto x467 = x466 + x61*(x434 + x64*(-x128*x211*(x212*x268*x458 + V{1}) - x130*x438 - x131*x142*(x252*(x260 + x446) + V{1}) - x143*x154*(-x252*x448 + V{1}) - x155*x162*(x252*x445 + V{1}) - x163*x165*(-x252*x443 + V{1}) - x166*x170*(x252*(x264 + x446) + V{1}) - x171*x178*(x252*(x264 + x450) + V{1}) - x179*x184*(x252*x449 + V{1}) - x185*x187*(-x252*x440 + V{1}) - x188*x197*(x252*(x264 + x455) + V{1}) - x198*x205*(x453*x67 + V{1}) - x206*x210*(x252*(x446 + x456) + V{1}) + V{0.0277777777777778}*x221*(x252*x462 + V{1})*cell[4] + V{0.0277777777777778}*x226*(x252*(x260 + x450) + V{1})*cell[6] + V{0.0277777777777778}*x229*x459*cell[8] + V{0.0555555555555556}*x231*(x252*(x260 + x455) + V{1})*cell[1] + V{0.0555555555555556}*x233*(-x465*x67 + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x450 + x456) + V{1})*cell[3]));
auto x468 = -cell[8];
auto x469 = x63*x79;
auto x470 = x414 + x469;
auto x471 = x65*(-x252*(x337 + x470) + V{1});
auto x472 = x157*x63;
auto x473 = -x472;
auto x474 = x414 + x473;
auto x475 = x337 + x474;
auto x476 = x347 + x469;
auto x477 = x131*(-x252*(x406 + x476) + V{1});
auto x478 = x347 + x473;
auto x479 = x406 + x478;
auto x480 = x352 + x472;
auto x481 = x422 + x480;
auto x482 = x356 + x472;
auto x483 = x419 + x482;
auto x484 = x278 + x283 + x314;
auto x485 = x295 + x361;
auto x486 = x122*x63;
auto x487 = x282 + x313 - x486;
auto x488 = x365 + x470;
auto x489 = x429 + x476;
auto x490 = -x67*(x286 + x314) + V{1};
auto x491 = x56 + x61*(x468 + x64*(-x128*x211*(-x212*x268*x487 + V{1}) - x130*x471 - x142*x477 - x143*x154*(x252*(x296 + x322) + V{1}) - x155*x162*(-x252*x475 + V{1}) - x163*x165*(x252*x481 + V{1}) - x166*x170*(-x252*x479 + V{1}) - x171*x178*(x252*x483 + V{1}) - x179*x184*(x252*(x299 + x322) + V{1}) - x185*x187*(x252*(x299 + x325) + V{1}) - x188*x197*(-x484*x67 + V{1}) - x198*x205*(x252*(x299 + x426) + V{1}) - x206*x210*(x252*(x322 + x485) + V{1}) + V{0.0277777777777778}*x221*(-x252*x488 + V{1})*cell[4] + V{0.0277777777777778}*x226*(-x252*x489 + V{1})*cell[6] + V{0.0277777777777778}*x229*(x252*(x296 + x325) + V{1})*cell[8] + V{0.0555555555555556}*x231*x490*cell[1] + V{0.0555555555555556}*x233*(x252*(x296 + x426) + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x325 + x485) + V{1})*cell[3])) + cell[8];
auto x492 = -cell[9];
auto x493 = x163*(x252*(x437 + x480) + V{1});
auto x494 = x439 + x478;
auto x495 = x441 + x469;
auto x496 = x337 + x495;
auto x497 = -x337;
auto x498 = x437 + x472;
auto x499 = x497 + x498;
auto x500 = x447 + x476;
auto x501 = x446 + x482;
auto x502 = x289 + x451;
auto x503 = -x282;
auto x504 = x436 + x486;
auto x505 = x503 + x504;
auto x506 = -x252*(x439 + x476) + V{1};
auto x507 = -x469;
auto x508 = x437 + x507;
auto x509 = x352 + x508;
auto x510 = x286 + x463;
auto x511 = x57 + cell[9];
auto x512 = x511 + x61*(x492 + x64*(-x128*x211*(x212*x268*x505 + V{1}) - x130*x65*(-x252*x496 + V{1}) - x131*x142*(-x252*x500 + V{1}) - x143*x154*(x252*(x296 + x446) + V{1}) - x155*x162*(x252*x499 + V{1}) - x165*x493 - x166*x170*(x252*x501 + V{1}) - x171*x178*(-x252*x494 + V{1}) - x179*x184*(x252*(x299 + x446) + V{1}) - x185*x187*(x252*(x299 + x450) + V{1}) - x188*x197*(x502*x67 + V{1}) - x198*x205*(x252*(x299 + x455) + V{1}) - x206*x210*(x252*(x446 + x485) + V{1}) + V{0.0277777777777778}*x221*(x252*x509 + V{1})*cell[4] + V{0.0277777777777778}*x226*x506*cell[6] + V{0.0277777777777778}*x229*(x252*(x296 + x450) + V{1})*cell[8] + V{0.0555555555555556}*x231*(-x510*x67 + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x296 + x455) + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x450 + x485) + V{1})*cell[3]));
auto x513 = -cell[10];
auto x514 = x238 + V{1};
auto x515 = x514*x75;
auto x516 = x19*x515;
auto x517 = x134 + x516;
auto x518 = x179*(x67*(x517 + x91) + V{1});
auto x519 = x147 + x516;
auto x520 = x134 + x519;
auto x521 = x175 + x516;
auto x522 = x521 + x91;
auto x523 = x514*x77;
auto x524 = x256 + x523;
auto x525 = x259 + x523;
auto x526 = x156*x514;
auto x527 = x256 + x526;
auto x528 = x259 + x526;
auto x529 = x67*(x101 + x516) + V{1};
auto x530 = x67*(x136 + x516) + V{1};
auto x531 = x514*x74;
auto x532 = -x516;
auto x533 = x228 + x532;
auto x534 = x100 + x519;
auto x535 = x135 + x521;
auto x536 = x40 + cell[10];
auto x537 = x536 + x61*(x513 + x64*(-x128*x211*(x212*x268*(x272 - x531) + V{1}) - x130*x65*(x252*(x253 + x524) + V{1}) - x131*x142*(x252*(x258 + x525) + V{1}) - x143*x154*(x520*x67 + V{1}) - x155*x162*(x252*(x253 + x527) + V{1}) - x163*x165*(x252*(x263 + x527) + V{1}) - x166*x170*(x252*(x258 + x528) + V{1}) - x171*x178*(x252*(x265 + x528) + V{1}) - x184*x518 - x185*x187*(x522*x67 + V{1}) - x188*x197*(x252*(x259 + x527) + V{1}) - x198*x205*x529 - x206*x210*x530 + V{0.0277777777777778}*x221*(x252*(x263 + x524) + V{1})*cell[4] + V{0.0277777777777778}*x226*(x252*(x265 + x525) + V{1})*cell[6] + V{0.0277777777777778}*x229*(-x533*x67 + V{1})*cell[8] + V{0.0555555555555556}*x231*(x252*(x259 + x524) + V{1})*cell[1] + V{0.0555555555555556}*x233*(x534*x67 + V{1})*cell[2] + V{0.0555555555555556}*x235*(x535*x67 + V{1})*cell[3]));
auto x538 = -cell[11];
auto x539 = x157 + x389;
auto x540 = x166*(x67*(x134 + x539) + V{1});
auto x541 = x175 + x389;
auto x542 = x157 + x541;
auto x543 = x389 + x80;
auto x544 = x134 + x543;
auto x545 = x256 + x384;
auto x546 = x295 + x388;
auto x547 = x256 + x388;
auto x548 = x295 + x384;
auto x549 = x67*(x164 + x389) + V{1};
auto x550 = x67*(x149 + x389) + V{1};
auto x551 = x174 + x401 + x79;
auto x552 = x100 + x543;
auto x553 = x148 + x541;
auto x554 = x41 + cell[11];
auto x555 = x554 + x61*(x538 + x64*(-x128*x211*(x212*x268*(x304 - x391) + V{1}) - x130*x65*(x252*(x291 + x545) + V{1}) - x131*x142*(x544*x67 + V{1}) - x143*x154*(x252*(x258 + x546) + V{1}) - x155*x162*(x252*(x297 + x545) + V{1}) - x163*x165*(x252*(x297 + x547) + V{1}) - x170*x540 - x171*x178*(x542*x67 + V{1}) - x179*x184*(x252*(x258 + x548) + V{1}) - x185*x187*(x252*(x265 + x548) + V{1}) - x188*x197*x549 - x198*x205*(x252*(x295 + x545) + V{1}) - x206*x210*x550 + V{0.0277777777777778}*x221*(x252*(x291 + x547) + V{1})*cell[4] + V{0.0277777777777778}*x226*(-x551*x67 + V{1})*cell[6] + V{0.0277777777777778}*x229*(x252*(x265 + x546) + V{1})*cell[8] + V{0.0555555555555556}*x231*(x552*x67 + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x295 + x547) + V{1})*cell[2] + V{0.0555555555555556}*x235*(x553*x67 + V{1})*cell[3]));
auto x556 = -cell[12];
auto x557 = x155*(x67*(x320 + x451) + V{1});
auto x558 = x147 + x451;
auto x559 = x157 + x558;
auto x560 = x452 + x80;
auto x561 = x259 + x446;
auto x562 = x259 + x450;
auto x563 = x67*(x176 + x451) + V{1};
auto x564 = x67*(x186 + x451) + V{1};
auto x565 = x464 + x79;
auto x566 = x135 + x451 + x80;
auto x567 = x148 + x558;
auto x568 = x42 + cell[12];
auto x569 = x568 + x61*(x556 + x64*(-x128*x211*(x212*x268*(x331 - x454) + V{1}) - x130*x65*(x560*x67 + V{1}) - x131*x142*(x252*(x291 + x561) + V{1}) - x143*x154*(x252*(x324 + x446) + V{1}) - x162*x557 - x163*x165*(x559*x67 + V{1}) - x166*x170*(x252*(x297 + x561) + V{1}) - x171*x178*(x252*(x297 + x562) + V{1}) - x179*x184*(x252*(x327 + x446) + V{1}) - x185*x187*(x252*(x327 + x450) + V{1}) - x188*x197*x563 - x198*x205*x564 - x206*x210*(x252*(x295 + x561) + V{1}) + V{0.0277777777777778}*x221*(-x565*x67 + V{1})*cell[4] + V{0.0277777777777778}*x226*(x252*(x291 + x562) + V{1})*cell[6] + V{0.0277777777777778}*x229*(x252*(x324 + x450) + V{1})*cell[8] + V{0.0555555555555556}*x231*(x566*x67 + V{1})*cell[1] + V{0.0555555555555556}*x233*(x567*x67 + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x295 + x562) + V{1})*cell[3]));
auto x570 = -cell[13];
auto x571 = x515*x66;
auto x572 = x342 + x571;
auto x573 = x179*(x252*(x384 + x572) + V{1});
auto x574 = x166*(x252*(x382 + x526) + V{1});
auto x575 = -x376;
auto x576 = x572 + x575;
auto x577 = x514*x78;
auto x578 = -x577;
auto x579 = x382 + x578;
auto x580 = x398 + x571;
auto x581 = x384 + x580;
auto x582 = x399 + x526;
auto x583 = x206*(x67*(x389 + x517) + V{1});
auto x584 = -x531*x66;
auto x585 = x212*x268*(x395 + x515) + V{1};
auto x586 = -x571;
auto x587 = x338 + x586;
auto x588 = x376 + x587;
auto x589 = x379 + x577;
auto x590 = x389 + x521;
auto x591 = x43 + cell[13];
auto x592 = x591 + x61*(x570 + x64*(-x128*x211*x585 - x130*x65*(x252*(x384 + x524) + V{1}) - x131*x142*(x252*x579 + V{1}) - x143*x154*(x252*x576 + V{1}) - x155*x162*(x252*(x384 + x527) + V{1}) - x163*x165*(x252*(x388 + x527) + V{1}) - x170*x574 - x171*x178*(x252*x582 + V{1}) - x184*x573 - x185*x187*(x252*x581 + V{1}) - x188*x197*(x252*(x392 + x527) + V{1}) - x198*x205*(x252*(x545 + x584) + V{1}) - x210*x583 + V{0.0277777777777778}*x221*(x252*(x388 + x524) + V{1})*cell[4] + V{0.0277777777777778}*x226*(-x252*x589 + V{1})*cell[6] + V{0.0277777777777778}*x229*(-x252*x588 + V{1})*cell[8] + V{0.0555555555555556}*x231*(x252*(x392 + x524) + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x547 + x584) + V{1})*cell[2] + V{0.0555555555555556}*x235*(x590*x67 + V{1})*cell[3]));
auto x593 = -cell[14];
auto x594 = x143*(x252*(x352 + x572) + V{1});
auto x595 = x350 + x577;
auto x596 = x337 + x587;
auto x597 = x497 + x572;
auto x598 = -x526;
auto x599 = x348 + x598;
auto x600 = x342 + x526;
auto x601 = x356 + x600;
auto x602 = x288 + x517;
auto x603 = x363 + x515;
auto x604 = x503 + x603;
auto x605 = -x252*(x348 + x577) + V{1};
auto x606 = x352 + x580;
auto x607 = x284 + x532;
auto x608 = x44 + cell[14];
auto x609 = x608 + x61*(x593 + x64*(-x128*x211*(x212*x268*x604 + V{1}) - x130*x65*(x252*(x293 + x523) + V{1}) - x131*x142*(-x252*x595 + V{1}) - x154*x594 - x155*x162*(x252*(x293 + x526) + V{1}) - x163*x165*(x252*(x298 + x526) + V{1}) - x166*x170*(x252*x601 + V{1}) - x171*x178*(-x252*x599 + V{1}) - x179*x184*(x252*x597 + V{1}) - x185*x187*(-x252*x596 + V{1}) - x188*x197*(x252*(x361 + x527) + V{1}) - x198*x205*(x252*(x293 + x584) + V{1}) - x206*x210*(x602*x67 + V{1}) + V{0.0277777777777778}*x221*(x252*(x298 + x523) + V{1})*cell[4] + V{0.0277777777777778}*x226*x605*cell[6] + V{0.0277777777777778}*x229*(x252*x606 + V{1})*cell[8] + V{0.0555555555555556}*x231*(x252*(x361 + x524) + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x298 + x584) + V{1})*cell[2] + V{0.0555555555555556}*x235*(-x607*x67 + V{1})*cell[3]));
auto x610 = -cell[15];
auto x611 = x410 + x571;
auto x612 = x179*(x252*(x446 + x611) + V{1});
auto x613 = x155*(x252*(x444 + x526) + V{1});
auto x614 = -x439;
auto x615 = x611 + x614;
auto x616 = x444 + x578;
auto x617 = x460 + x571;
auto x618 = x446 + x617;
auto x619 = x461 + x526;
auto x620 = x198*(x67*(x452 + x516) + V{1});
auto x621 = x212*x268*(x457 + x515) + V{1};
auto x622 = x407 + x586;
auto x623 = x439 + x622;
auto x624 = x442 + x577;
auto x625 = x451 + x519;
auto x626 = x45 + cell[15];
auto x627 = x61*(x610 + x64*(-x128*x211*x621 - x130*x65*(x252*x616 + V{1}) - x131*x142*(x252*(x446 + x525) + V{1}) - x143*x154*(x252*x618 + V{1}) - x162*x613 - x163*x165*(x252*x619 + V{1}) - x166*x170*(x252*(x446 + x528) + V{1}) - x171*x178*(x252*(x450 + x528) + V{1}) - x184*x612 - x185*x187*(x252*x615 + V{1}) - x188*x197*(x252*(x455 + x528) + V{1}) - x205*x620 - x206*x210*(x252*(x561 + x584) + V{1}) + V{0.0277777777777778}*x221*(-x252*x624 + V{1})*cell[4] + V{0.0277777777777778}*x226*(x252*(x450 + x525) + V{1})*cell[6] + V{0.0277777777777778}*x229*(-x252*x623 + V{1})*cell[8] + V{0.0555555555555556}*x231*(x252*(x455 + x525) + V{1})*cell[1] + V{0.0555555555555556}*x233*(x625*x67 + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x562 + x584) + V{1})*cell[3])) + x626;
auto x628 = -cell[16];
auto x629 = x185*(x252*(x419 + x611) + V{1});
auto x630 = x417 + x577;
auto x631 = x406 + x622;
auto x632 = -x406;
auto x633 = x611 + x632;
auto x634 = x415 + x598;
auto x635 = x410 + x526;
auto x636 = x422 + x635;
auto x637 = x319 + x516 + x91;
auto x638 = -x313;
auto x639 = x427 + x638;
auto x640 = x515 + x639;
auto x641 = -x252*(x415 + x577) + V{1};
auto x642 = x419 + x617;
auto x643 = x315 + x532;
auto x644 = x46 + cell[16];
auto x645 = x61*(x628 + x64*(-x128*x211*(x212*x268*x640 + V{1}) - x130*x65*(-x252*x630 + V{1}) - x131*x142*(x252*(x323 + x523) + V{1}) - x143*x154*(-x252*x631 + V{1}) - x155*x162*(x252*x636 + V{1}) - x163*x165*(-x252*x634 + V{1}) - x166*x170*(x252*(x323 + x526) + V{1}) - x171*x178*(x252*(x326 + x526) + V{1}) - x179*x184*(x252*x633 + V{1}) - x187*x629 - x188*x197*(x252*(x426 + x528) + V{1}) - x198*x205*(x637*x67 + V{1}) - x206*x210*(x252*(x323 + x584) + V{1}) + V{0.0277777777777778}*x221*x641*cell[4] + V{0.0277777777777778}*x226*(x252*(x326 + x523) + V{1})*cell[6] + V{0.0277777777777778}*x229*(x252*x642 + V{1})*cell[8] + V{0.0555555555555556}*x231*(x252*(x426 + x525) + V{1})*cell[1] + V{0.0555555555555556}*x233*(-x643*x67 + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x326 + x584) + V{1})*cell[3])) + x644;
auto x646 = -cell[17];
auto x647 = x155*(x252*(x384 + x498) + V{1});
auto x648 = x374 + x472;
auto x649 = x166*(x252*(x446 + x648) + V{1});
auto x650 = x498 + x575;
auto x651 = x614 + x648;
auto x652 = x384 + x508;
auto x653 = x374 + x507;
auto x654 = x446 + x653;
auto x655 = x188*(x67*(x451 + x539) + V{1});
auto x656 = x295 + x392;
auto x657 = x212*x268*(x373 + x504) + V{1};
auto x658 = x376 + x495;
auto x659 = x378 + x469;
auto x660 = x439 + x659;
auto x661 = x451 + x543;
auto x662 = x47 + cell[17];
auto x663 = x61*(x64*(-x128*x211*x657 - x130*x65*(x252*x652 + V{1}) - x131*x142*(x252*x654 + V{1}) - x143*x154*(x252*(x446 + x546) + V{1}) - x162*x647 - x163*x165*(x252*x650 + V{1}) - x170*x649 - x171*x178*(x252*x651 + V{1}) - x179*x184*(x252*(x446 + x548) + V{1}) - x185*x187*(x252*(x450 + x548) + V{1}) - x197*x655 - x198*x205*(x252*(x455 + x548) + V{1}) - x206*x210*(x252*(x446 + x656) + V{1}) + V{0.0277777777777778}*x221*(-x252*x658 + V{1})*cell[4] + V{0.0277777777777778}*x226*(-x252*x660 + V{1})*cell[6] + V{0.0277777777777778}*x229*(x252*(x450 + x546) + V{1})*cell[8] + V{0.0555555555555556}*x231*(x661*x67 + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x455 + x546) + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x450 + x656) + V{1})*cell[3]) + x646) + x662;
auto x664 = -cell[18];
auto x665 = x171*(x252*(x419 + x648) + V{1});
auto x666 = x376 + x474;
auto x667 = x406 + x659;
auto x668 = x632 + x648;
auto x669 = x385 + x470;
auto x670 = x384 + x472;
auto x671 = x422 + x670;
auto x672 = x319 + x539;
auto x673 = x373 + x486 + x638;
auto x674 = -x252*(x376 + x470) + V{1};
auto x675 = x419 + x653;
auto x676 = x317 + x401;
auto x677 = x48 + cell[18];
auto x678 = x61*(x64*(-x128*x211*(x212*x268*x673 + V{1}) - x130*x65*(-x252*x669 + V{1}) - x131*x142*(-x252*x667 + V{1}) - x143*x154*(x252*(x322 + x546) + V{1}) - x155*x162*(x252*x671 + V{1}) - x163*x165*(-x252*x666 + V{1}) - x166*x170*(x252*x668 + V{1}) - x178*x665 - x179*x184*(x252*(x322 + x548) + V{1}) - x185*x187*(x252*(x325 + x548) + V{1}) - x188*x197*(x67*x672 + V{1}) - x198*x205*(x252*(x426 + x548) + V{1}) - x206*x210*(x252*(x322 + x656) + V{1}) + V{0.0277777777777778}*x221*x674*cell[4] + V{0.0277777777777778}*x226*(x252*x675 + V{1})*cell[6] + V{0.0277777777777778}*x229*(x252*(x325 + x546) + V{1})*cell[8] + V{0.0555555555555556}*x231*(-x67*x676 + V{1})*cell[1] + V{0.0555555555555556}*x233*(x252*(x426 + x546) + V{1})*cell[2] + V{0.0555555555555556}*x235*(x252*(x325 + x656) + V{1})*cell[3]) + x664) + x677;
auto x679 = x236 + x276 + x308 + x335 + x370 + x404 + x433 + x467 + x491 + x512 + x537 + x555 + x569 + x592 + x609 + x627 + x645 + x663 + x678 + V{1};
auto x680 = -x184;
auto x681 = -x187;
auto x682 = -x154;
auto x683 = V{0.0277777777777778}*cell[8];
auto x684 = -x125;
auto x685 = V{1} - x126;
auto x686 = x684 + x685;
auto x687 = -x123;
auto x688 = x109 + x687;
auto x689 = x138 + x182 + x686 + x688;
auto x690 = V{0.0277777777777778}*cell[4];
auto x691 = -x110*x19*x99;
auto x692 = x577 + x691;
auto x693 = x120 + x686;
auto x694 = x160 + x688 + x693;
auto x695 = x255*x99;
auto x696 = x578 + x695;
auto x697 = -x130;
auto x698 = V{0.0277777777777778}*cell[6];
auto x699 = -x110*x124*x19;
auto x700 = x577 + x699;
auto x701 = x138 + x687;
auto x702 = x169 + x693 + x701;
auto x703 = x124*x255;
auto x704 = x578 + x703;
auto x705 = -x142;
auto x706 = -x162;
auto x707 = x526 + x695;
auto x708 = -x165;
auto x709 = -x170;
auto x710 = x526 + x703;
auto x711 = -x178;
auto x712 = V{0.0555555555555556}*cell[1];
auto x713 = x196 + x693;
auto x714 = V{0.0555555555555556}*cell[2];
auto x715 = x203 + x685 + x688;
auto x716 = V{0.0555555555555556}*cell[3];
auto x717 = x209 + x684 + x701 + V{1};
auto x718 = -x197;
auto x719 = -x205;
auto x720 = -x210;
auto x721 = -x128;
auto x722 = x571 + x695;
auto x723 = x61*(x593 + x64*(x131*x705*(-x252*x595 + V{1}) + x155*x706*(x252*(x497 + x707) + V{1}) + x163*x708*(x252*(x352 + x707) + V{1}) + x166*x709*(x252*x601 + V{1}) + x171*x711*(-x252*x599 + V{1}) + x179*x680*(x252*x597 + V{1}) + x185*x681*(-x252*x596 + V{1}) + x188*x718*(x252*(x356 + x707) + V{1}) + x198*x719*(x252*(x497 + x722) + V{1}) + x206*x720*(x602*x67 + V{1}) + x211*x721*(x212*x268*x604 + V{1}) + x594*x682 + x605*x698*x702 + x65*x697*(x252*(-x337 - x692) + V{1}) + x683*x689*(x252*x606 + V{1}) + x690*x694*(x252*(x352 + x696) + V{1}) + x712*x713*(x252*(-x347 - x692) + V{1}) + x714*x715*(x252*(x352 + x722) + V{1}) + x716*x717*(-x607*x67 + V{1})));
auto x724 = x571 + x703;
auto x725 = x61*(x628 + x64*(x131*x705*(x252*(-x406 - x700) + V{1}) + x143*x682*(-x252*x631 + V{1}) + x155*x706*(x252*x636 + V{1}) + x163*x708*(-x252*x634 + V{1}) + x166*x709*(x252*(x632 + x710) + V{1}) + x171*x711*(x252*(x419 + x710) + V{1}) + x179*x680*(x252*x633 + V{1}) + x188*x718*(x252*(x422 + x710) + V{1}) + x198*x719*(x637*x67 + V{1}) + x206*x720*(x252*(x632 + x724) + V{1}) + x211*x721*(x212*x268*x640 + V{1}) + x629*x681 + x641*x690*x694 + x65*x697*(-x252*x630 + V{1}) + x683*x689*(x252*x642 + V{1}) + x698*x702*(x252*(x419 + x704) + V{1}) + x712*x713*(x252*(-x414 - x700) + V{1}) + x714*x715*(-x643*x67 + V{1}) + x716*x717*(x252*(x419 + x724) + V{1})));
auto x726 = x357 + x695;
auto x727 = x357 + x703;
auto x728 = x381 + x695;
auto x729 = x346 + x691;
auto x730 = x381 + x703;
auto x731 = x346 + x699;
auto x732 = x339 + x691;
auto x733 = x353 + x695;
auto x734 = x61*(x371 + x64*(x143*x682*(-x252*x377 + V{1}) + x155*x706*(x252*(x384 + x728) + V{1}) + x163*x708*(x252*(-x376 - x729) + V{1}) + x166*x709*(x252*x383 + V{1}) + x171*x711*(-x252*x380 + V{1}) + x179*x680*(x252*x387 + V{1}) + x185*x681*(-x252*x386 + V{1}) + x188*x718*(x252*(x374 + x728) + V{1}) + x198*x719*(x252*(x384 + x733) + V{1}) + x206*x720*(x390*x67 + V{1}) + x211*x721*(x212*x268*x396 + V{1}) + x375*x705 + x397*x683*x689 + x65*x697*(x252*(x384 + x726) + V{1}) + x690*x694*(x252*(x575 + x726) + V{1}) + x698*x702*(x252*x400 + V{1}) + x712*x713*(x252*(x374 + x726) + V{1}) + x714*x715*(x252*(-x376 - x732) + V{1}) + x716*x717*(-x402*x67 + V{1})));
auto x735 = x339 + x699;
auto x736 = x353 + x703;
auto x737 = x61*(x434 + x64*(x131*x705*(x252*(x446 + x727) + V{1}) + x143*x682*(-x252*x448 + V{1}) + x155*x706*(x252*x445 + V{1}) + x163*x708*(-x252*x443 + V{1}) + x166*x709*(x252*(x446 + x730) + V{1}) + x171*x711*(x252*(-x439 - x731) + V{1}) + x179*x680*(x252*x449 + V{1}) + x185*x681*(-x252*x440 + V{1}) + x188*x718*(x252*(x437 + x730) + V{1}) + x198*x719*(x453*x67 + V{1}) + x206*x720*(x252*(x446 + x736) + V{1}) + x211*x721*(x212*x268*x458 + V{1}) + x438*x697 + x459*x683*x689 + x690*x694*(x252*x462 + V{1}) + x698*x702*(x252*(x614 + x727) + V{1}) + x712*x713*(x252*(x437 + x727) + V{1}) + x714*x715*(-x465*x67 + V{1}) + x716*x717*(x252*(-x439 - x735) + V{1})));
auto x738 = x336 - x52 + x591 - x61*(x336 + x64*(x131*x705*(x252*x359 + V{1}) + x143*x682*(x252*x355 + V{1}) + x155*x706*(x252*(-x337 - x729) + V{1}) + x163*x708*(x252*(x352 + x728) + V{1}) + x166*x709*(-x252*x351 + V{1}) + x179*x680*(-x252*x345 + V{1}) + x188*x718*(x252*(-x347 - x729) + V{1}) + x198*x719*(x252*(-x337 - x732) + V{1}) + x206*x720*(-x360*x67 + V{1}) + x211*x721*(-x212*x268*x364 + V{1}) + x341*x681 + x349*x711 + x369*x716*x717 + x65*x697*(x252*(x497 + x726) + V{1}) + x683*x689*(-x252*x366 + V{1}) + x690*x694*(x252*(x352 + x726) + V{1}) + x698*x702*(-x252*x368 + V{1}) + x712*x713*(x252*(x356 + x726) + V{1}) + x714*x715*(x252*(x352 + x733) + V{1}))) + x61*(x570 + x64*(x131*x705*(x252*x579 + V{1}) + x143*x682*(x252*x576 + V{1}) + x155*x706*(x252*(x384 + x707) + V{1}) + x163*x708*(x252*(x575 + x707) + V{1}) + x171*x711*(x252*x582 + V{1}) + x185*x681*(x252*x581 + V{1}) + x188*x718*(x252*(x374 + x707) + V{1}) + x198*x719*(x252*(x384 + x722) + V{1}) + x211*x585*x721 + x573*x680 + x574*x709 + x583*x720 + x65*x697*(x252*(x384 + x696) + V{1}) + x683*x689*(-x252*x588 + V{1}) + x690*x694*(x252*(-x376 - x692) + V{1}) + x698*x702*(-x252*x589 + V{1}) + x712*x713*(x252*(x374 + x696) + V{1}) + x714*x715*(x252*(x575 + x722) + V{1}) + x716*x717*(x590*x67 + V{1})));
auto x739 = x405 - x54 - x61*(x405 + x64*(x131*x705*(x252*(x632 + x727) + V{1}) + x155*x706*(-x252*x418 + V{1}) + x166*x709*(x252*(-x406 - x731) + V{1}) + x171*x711*(x252*(x419 + x730) + V{1}) + x179*x680*(-x252*x413 + V{1}) + x185*x681*(x252*x421 + V{1}) + x188*x718*(x252*(-x414 - x731) + V{1}) + x198*x719*(-x425*x67 + V{1}) + x206*x720*(x252*(-x406 - x735) + V{1}) + x211*x721*(-x212*x268*x428 + V{1}) + x409*x682 + x416*x708 + x432*x714*x715 + x65*x697*(x252*x424 + V{1}) + x683*x689*(-x252*x430 + V{1}) + x690*x694*(-x252*x431 + V{1}) + x698*x702*(x252*(x419 + x727) + V{1}) + x712*x713*(x252*(x422 + x727) + V{1}) + x716*x717*(x252*(x419 + x736) + V{1}))) + x61*(x610 + x64*(x131*x705*(x252*(x446 + x704) + V{1}) + x143*x682*(x252*x618 + V{1}) + x163*x708*(x252*x619 + V{1}) + x166*x709*(x252*(x446 + x710) + V{1}) + x171*x711*(x252*(x614 + x710) + V{1}) + x185*x681*(x252*x615 + V{1}) + x188*x718*(x252*(x437 + x710) + V{1}) + x206*x720*(x252*(x446 + x724) + V{1}) + x211*x621*x721 + x612*x680 + x613*x706 + x620*x719 + x65*x697*(x252*x616 + V{1}) + x683*x689*(-x252*x623 + V{1}) + x690*x694*(-x252*x624 + V{1}) + x698*x702*(x252*(-x439 - x700) + V{1}) + x712*x713*(x252*(x437 + x704) + V{1}) + x714*x715*(x625*x67 + V{1}) + x716*x717*(x252*(x614 + x724) + V{1}))) + x626;
auto x740 = x237 + x371 + x434 - x49 - x53 + x536 - x55 + x608 - x61*(x237 + x64*(x131*x705*(x252*(x358 + x703) + V{1}) + x143*x682*(-x245*x67 + V{1}) + x155*x706*(x252*(x410 + x728) + V{1}) + x163*x708*(x252*(-x407 - x729) + V{1}) + x166*x709*(x252*(x342 + x730) + V{1}) + x171*x711*(x252*(-x338 - x731) + V{1}) + x179*x680*(x251*x67 + V{1}) + x185*x681*(-x248*x67 + V{1}) + x188*x718*(x252*(x703 + x728) + V{1}) + x198*x719*(x266*x67 + V{1}) + x206*x720*(x267*x67 + V{1}) + x211*x721*(x212*x268*(x363 + x394 + x427) + V{1}) + x273*x683*x689 + x65*x697*(x252*(x423 + x695) + V{1}) + x690*x694*(x252*(x460 + x726) + V{1}) + x698*x702*(x252*(x398 + x727) + V{1}) + x712*x713*(x252*(x703 + x726) + V{1}) + x714*x715*(-x274*x67 + V{1}) + x716*x717*(-x275*x67 + V{1}))) + x61*(x513 + x64*(x131*x705*(x252*(x342 + x704) + V{1}) + x143*x682*(x520*x67 + V{1}) + x155*x706*(x252*(x635 + x695) + V{1}) + x163*x708*(x252*(x460 + x707) + V{1}) + x166*x709*(x252*(x600 + x703) + V{1}) + x171*x711*(x252*(x398 + x710) + V{1}) + x185*x681*(x522*x67 + V{1}) + x188*x718*(x252*(x703 + x707) + V{1}) + x198*x529*x719 + x206*x530*x720 + x211*x721*(x212*x268*(x427 + x603) + V{1}) + x518*x680 + x65*x697*(x252*(x410 + x696) + V{1}) + x683*x689*(-x533*x67 + V{1}) + x690*x694*(x252*(-x407 - x692) + V{1}) + x698*x702*(x252*(-x338 - x700) + V{1}) + x712*x713*(x252*(x696 + x703) + V{1}) + x714*x715*(x534*x67 + V{1}) + x716*x717*(x535*x67 + V{1}))) + x644 + x723 + x725 - x734 - x737 + x738 + x739;
auto x741 = x469 + x691;
auto x742 = x384 + x695;
auto x743 = -x110*x122*x19;
auto x744 = x376 + x743;
auto x745 = x122*x255;
auto x746 = x342 + x745;
auto x747 = x575 + x695;
auto x748 = x384 + x745;
auto x749 = x575 + x745;
auto x750 = x374 + x745;
auto x751 = x61*(x64*(x131*x705*(-x252*x667 + V{1}) + x143*x682*(x252*(-x406 - x744) + V{1}) + x155*x706*(x252*x671 + V{1}) + x163*x708*(-x252*x666 + V{1}) + x166*x709*(x252*x668 + V{1}) + x179*x680*(x252*(x632 + x748) + V{1}) + x185*x681*(x252*(x419 + x748) + V{1}) + x188*x718*(x67*x672 + V{1}) + x198*x719*(x252*(x422 + x748) + V{1}) + x206*x720*(x252*(x632 + x750) + V{1}) + x211*x721*(x212*x268*x673 + V{1}) + x65*x697*(-x252*x669 + V{1}) + x665*x711 + x674*x690*x694 + x683*x689*(x252*(x419 + x749) + V{1}) + x698*x702*(x252*x675 + V{1}) + x712*x713*(-x67*x676 + V{1}) + x714*x715*(x252*(-x414 - x744) + V{1}) + x716*x717*(x252*(x419 + x750) + V{1})) + x664);
auto x752 = x352 + x695;
auto x753 = x352 + x745;
auto x754 = x497 + x695;
auto x755 = x337 + x743;
auto x756 = x497 + x745;
auto x757 = x347 + x743;
auto x758 = x446 + x745;
auto x759 = x61*(x492 + x64*(x131*x705*(-x252*x500 + V{1}) + x143*x682*(x252*(x446 + x753) + V{1}) + x155*x706*(x252*x499 + V{1}) + x166*x709*(x252*x501 + V{1}) + x171*x711*(-x252*x494 + V{1}) + x179*x680*(x252*(x446 + x756) + V{1}) + x185*x681*(x252*(-x439 - x755) + V{1}) + x188*x718*(x502*x67 + V{1}) + x198*x719*(x252*(x437 + x756) + V{1}) + x206*x720*(x252*(x356 + x758) + V{1}) + x211*x721*(x212*x268*x505 + V{1}) + x493*x708 + x506*x698*x702 + x65*x697*(-x252*x496 + V{1}) + x683*x689*(x252*(x614 + x753) + V{1}) + x690*x694*(x252*x509 + V{1}) + x712*x713*(-x510*x67 + V{1}) + x714*x715*(x252*(x437 + x753) + V{1}) + x716*x717*(x252*(-x439 - x757) + V{1})));
auto x760 = x419 + x745;
auto x761 = x468 - x56 - x61*(x468 + x64*(x143*x682*(x252*(x632 + x753) + V{1}) + x155*x706*(-x252*x475 + V{1}) + x163*x708*(x252*x481 + V{1}) + x166*x709*(-x252*x479 + V{1}) + x171*x711*(x252*x483 + V{1}) + x179*x680*(x252*(-x406 - x755) + V{1}) + x185*x681*(x252*(x497 + x760) + V{1}) + x188*x718*(-x484*x67 + V{1}) + x198*x719*(x252*(-x414 - x755) + V{1}) + x206*x720*(x252*(-x406 - x757) + V{1}) + x211*x721*(-x212*x268*x487 + V{1}) + x471*x697 + x477*x705 + x490*x712*x713 + x683*x689*(x252*(x419 + x753) + V{1}) + x690*x694*(-x252*x488 + V{1}) + x698*x702*(-x252*x489 + V{1}) + x714*x715*(x252*(x422 + x753) + V{1}) + x716*x717*(x252*(x356 + x760) + V{1}))) + x61*(x64*(x131*x705*(x252*x654 + V{1}) + x143*x682*(x252*(x446 + x749) + V{1}) + x163*x708*(x252*x650 + V{1}) + x171*x711*(x252*x651 + V{1}) + x179*x680*(x252*(x446 + x748) + V{1}) + x185*x681*(x252*(x614 + x748) + V{1}) + x198*x719*(x252*(x437 + x748) + V{1}) + x206*x720*(x252*(x446 + x750) + V{1}) + x211*x657*x721 + x647*x706 + x649*x709 + x65*x697*(x252*x652 + V{1}) + x655*x718 + x683*x689*(x252*(-x439 - x744) + V{1}) + x690*x694*(-x252*x658 + V{1}) + x698*x702*(-x252*x660 + V{1}) + x712*x713*(x661*x67 + V{1}) + x714*x715*(x252*(x437 + x749) + V{1}) + x716*x717*(x252*(x614 + x750) + V{1})) + x646) + x662;
auto x762 = x277 + x403 - x44 + x492 - x50 + x554 - x57 + x593 - x61*(x277 + x64*(x131*x705*(-x287*x67 + V{1}) + x143*x682*(x252*(x352 + x746) + V{1}) + x155*x706*(x252*(x472 + x754) + V{1}) + x163*x708*(x252*(x480 + x695) + V{1}) + x166*x709*(x290*x67 + V{1}) + x171*x711*(-x285*x67 + V{1}) + x179*x680*(x252*(x497 + x746) + V{1}) + x185*x681*(x252*(-x338 - x755) + V{1}) + x188*x718*(x300*x67 + V{1}) + x198*x719*(x252*(x745 + x754) + V{1}) + x206*x720*(x301*x67 + V{1}) + x211*x721*(x212*x268*(x363 + x486 + x503) + V{1}) + x305*x698*x702 + x65*x697*(x252*(-x337 - x741) + V{1}) + x683*x689*(x252*(x398 + x753) + V{1}) + x690*x694*(x252*(x507 + x752) + V{1}) + x712*x713*(-x306*x67 + V{1}) + x714*x715*(x252*(x745 + x752) + V{1}) + x716*x717*(-x307*x67 + V{1}))) + x61*(x538 + x64*(x131*x705*(x544*x67 + V{1}) + x143*x682*(x252*(x575 + x746) + V{1}) + x155*x706*(x252*(x670 + x695) + V{1}) + x163*x708*(x252*(x472 + x747) + V{1}) + x171*x711*(x542*x67 + V{1}) + x179*x680*(x252*(x384 + x746) + V{1}) + x185*x681*(x252*(x398 + x748) + V{1}) + x188*x549*x718 + x198*x719*(x252*(x742 + x745) + V{1}) + x206*x550*x720 + x211*x721*(x212*x268*(x395 + x486) + V{1}) + x540*x709 + x65*x697*(x252*(x507 + x742) + V{1}) + x683*x689*(x252*(-x338 - x744) + V{1}) + x690*x694*(x252*(-x376 - x741) + V{1}) + x698*x702*(-x551*x67 + V{1}) + x712*x713*(x552*x67 + V{1}) + x714*x715*(x252*(x745 + x747) + V{1}) + x716*x717*(x553*x67 + V{1}))) + x677 - x723 + x734 + x738 + x751 - x759 + x761;
auto x763 = x469 + x699;
auto x764 = x446 + x703;
auto x765 = x407 + x743;
auto x766 = x614 + x703;
auto x767 = x410 + x745;
auto x768 = x419 + x703;
auto x769 = x632 + x703;
auto x770 = x309 - x46 + x466 - x48 - x51 + x511 + x568 - x61*(x309 + x64*(x131*x705*(x252*(-x406 - x763) + V{1}) + x143*x682*(x252*(-x406 - x765) + V{1}) + x155*x706*(x321*x67 + V{1}) + x163*x708*(-x316*x67 + V{1}) + x166*x709*(x252*(x472 + x769) + V{1}) + x171*x711*(x252*(x472 + x768) + V{1}) + x179*x680*(x252*(x632 + x767) + V{1}) + x185*x681*(x252*(x410 + x760) + V{1}) + x188*x718*(x328*x67 + V{1}) + x198*x719*(x329*x67 + V{1}) + x206*x720*(x252*(x745 + x769) + V{1}) + x211*x721*(x212*x268*(x486 + x639) + V{1}) + x332*x690*x694 + x65*x697*(-x318*x67 + V{1}) + x683*x689*(x252*(x460 + x760) + V{1}) + x698*x702*(x252*(x507 + x768) + V{1}) + x712*x713*(-x333*x67 + V{1}) + x714*x715*(-x334*x67 + V{1}) + x716*x717*(x252*(x703 + x760) + V{1}))) + x61*(x556 + x64*(x131*x705*(x252*(x507 + x764) + V{1}) + x143*x682*(x252*(x460 + x758) + V{1}) + x163*x708*(x559*x67 + V{1}) + x166*x709*(x252*(x472 + x764) + V{1}) + x171*x711*(x252*(x472 + x766) + V{1}) + x179*x680*(x252*(x410 + x758) + V{1}) + x185*x681*(x252*(x614 + x767) + V{1}) + x188*x563*x718 + x198*x564*x719 + x206*x720*(x252*(x703 + x758) + V{1}) + x211*x721*(x212*x268*(x457 + x486) + V{1}) + x557*x706 + x65*x697*(x560*x67 + V{1}) + x683*x689*(x252*(-x439 - x765) + V{1}) + x690*x694*(-x565*x67 + V{1}) + x698*x702*(x252*(-x439 - x763) + V{1}) + x712*x713*(x566*x67 + V{1}) + x714*x715*(x567*x67 + V{1}) + x716*x717*(x252*(x745 + x766) + V{1}))) + x628 + x664 - x725 + x737 + x739 - x751 + x759 + x761;
cell[0] = x236;
cell[1] = x276;
cell[2] = x308;
cell[3] = x335;
cell[4] = x370;
cell[5] = x404;
cell[6] = x433;
cell[7] = x467;
cell[8] = x491;
cell[9] = x512;
cell[10] = x537;
cell[11] = x555;
cell[12] = x569;
cell[13] = x592;
cell[14] = x609;
cell[15] = x627;
cell[16] = x645;
cell[17] = x663;
cell[18] = x678;
return { x679, (x740*x740 + x762*x762 + x770*x770)/((x679)*(x679)) };
}
 
};
 
}
 
}
 
#endif
 
#endif