diff options
Diffstat (limited to 'tangle/LLBM/kernel/collect_q_criterion.h')
| -rw-r--r-- | tangle/LLBM/kernel/collect_q_criterion.h | 72 |
1 files changed, 39 insertions, 33 deletions
diff --git a/tangle/LLBM/kernel/collect_q_criterion.h b/tangle/LLBM/kernel/collect_q_criterion.h index 19b7f68..fa19dc7 100644 --- a/tangle/LLBM/kernel/collect_q_criterion.h +++ b/tangle/LLBM/kernel/collect_q_criterion.h @@ -56,39 +56,45 @@ __device__ static void apply( T x28 = -x13; T x29 = x2 + x28; T x30 = -rho*(x14 + x21 - x24) + rho*(x14 + x24 + x27 + x29) - T{72.0000000000000}*f_curr[11] - T{72.0000000000000}*f_curr[7]; - T x31 = T{72.0000000000000}*f_curr[1]; - T x32 = T{72.0000000000000}*f_curr[17]; - T x33 = x4 + u_2; - T x34 = T{6.00000000000000}*u_2; - T x35 = x11 - x34; - T x36 = rho*(x15 + x35 - T{9.00000000000000}*x33*x33); - T x37 = -u_2; - T x38 = x37 + u_0; - T x39 = x11 + x34; - T x40 = x13 + x39; - T x41 = rho*(x17 + x40 - T{9.00000000000000}*x38*x38); - T x42 = u_0 + u_2; - T x43 = T{9.00000000000000}*(x42*x42); - T x44 = x27 + x34; - T x45 = x14 + x28; - T x46 = -rho*(x15 + x39 - x43) + rho*(x43 + x44 + x45) - T{72.0000000000000}*f_curr[15] - T{72.0000000000000}*f_curr[3]; - T x47 = T{72.0000000000000}*f_curr[4]; - T x48 = T{72.0000000000000}*f_curr[14]; - T x49 = x18 + u_2; - T x50 = rho*(x20 + x35 - T{9.00000000000000}*x49*x49); - T x51 = x37 + u_1; - T x52 = rho*(x3 + x40 - T{9.00000000000000}*x51*x51); - T x53 = u_1 + u_2; - T x54 = T{9.00000000000000}*(x53*x53); - T x55 = -rho*(x20 + x39 - x54) + rho*(x29 + x44 + x54) - T{72.0000000000000}*f_curr[0] - T{72.0000000000000}*f_curr[18]; - T x56 = T{2.00000000000000}*rho; - T x57 = T{6.00000000000000}*x8; - T x58 = -x31 - x32 - x36 - x41 + x46; - T x59 = -x0 - x1 - x16 - x22 + x30; - T x60 = T{6.00000000000000}*x6; - T x61 = -x47 - x48 - x50 - x52 + x55; - T x62 = T{6.00000000000000}*x12; - T strain = T{0.0277777777777778}*sqrt((x0 + x1 + x16 + x22 + x30)*(x0 + x1 + x16 + x22 + x30) + (x31 + x32 + x36 + x41 + x46)*(x31 + x32 + x36 + x41 + x46) + (x47 + x48 + x50 + x52 + x55)*(x47 + x48 + x50 + x52 + x55) + T{0.500000000000000}*((-x56*(x39 - x62) + x56*(x44 + x62) + x58 + x61 - 72*f_curr[16] - 72*f_curr[2])*(-x56*(x39 - x62) + x56*(x44 + x62) + x58 + x61 - 72*f_curr[16] - 72*f_curr[2])) + T{0.500000000000000}*((-x56*(x10 + x20 - x60) + x56*(x26 + x29 + x60) + x59 + x61 - 72*f_curr[12] - 72*f_curr[6])*(-x56*(x10 + x20 - x60) + x56*(x26 + x29 + x60) + x59 + x61 - 72*f_curr[12] - 72*f_curr[6])) + T{0.500000000000000}*((-x56*(x15 - x57 + x7 - 2) + x56*(x25 + x45 + x57 + 2) + x58 + x59 - 72*f_curr[10] - 72*f_curr[8])*(-x56*(x15 - x57 + x7 - 2) + x56*(x25 + x45 + x57 + 2) + x58 + x59 - 72*f_curr[10] - 72*f_curr[8]))); + T x31 = x0 + x1 + x16 + x22 + x30; + T x32 = T{72.0000000000000}*f_curr[1]; + T x33 = T{72.0000000000000}*f_curr[17]; + T x34 = x4 + u_2; + T x35 = T{6.00000000000000}*u_2; + T x36 = x11 - x35; + T x37 = rho*(x15 + x36 - T{9.00000000000000}*x34*x34); + T x38 = -u_2; + T x39 = x38 + u_0; + T x40 = x11 + x35; + T x41 = x13 + x40; + T x42 = rho*(x17 + x41 - T{9.00000000000000}*x39*x39); + T x43 = u_0 + u_2; + T x44 = T{9.00000000000000}*(x43*x43); + T x45 = x27 + x35; + T x46 = x14 + x28; + T x47 = -rho*(x15 + x40 - x44) + rho*(x44 + x45 + x46) - T{72.0000000000000}*f_curr[15] - T{72.0000000000000}*f_curr[3]; + T x48 = x32 + x33 + x37 + x42 + x47; + T x49 = T{72.0000000000000}*f_curr[4]; + T x50 = T{72.0000000000000}*f_curr[14]; + T x51 = x18 + u_2; + T x52 = rho*(x20 + x36 - T{9.00000000000000}*x51*x51); + T x53 = x38 + u_1; + T x54 = rho*(x3 + x41 - T{9.00000000000000}*x53*x53); + T x55 = u_1 + u_2; + T x56 = T{9.00000000000000}*(x55*x55); + T x57 = -rho*(x20 + x40 - x56) + rho*(x29 + x45 + x56) - T{72.0000000000000}*f_curr[0] - T{72.0000000000000}*f_curr[18]; + T x58 = x49 + x50 + x52 + x54 + x57; + T x59 = T{2.00000000000000}*rho; + T x60 = T{6.00000000000000}*x8; + T x61 = -x32 - x33 - x37 - x42 + x47; + T x62 = -x0 - x1 - x16 - x22 + x30; + T x63 = -x59*(x15 - x60 + x7 + T{-2.00000000000000}) + x59*(x25 + x46 + x60 + T{2.00000000000000}) + x61 + x62 - T{72.0000000000000}*f_curr[10] - T{72.0000000000000}*f_curr[8]; + T x64 = T{6.00000000000000}*x6; + T x65 = -x49 - x50 - x52 - x54 + x57; + T x66 = -x59*(x10 + x20 - x64) + x59*(x26 + x29 + x64) + x62 + x65 - T{72.0000000000000}*f_curr[12] - T{72.0000000000000}*f_curr[6]; + T x67 = T{6.00000000000000}*x12; + T x68 = -x59*(x40 - x67) + x59*(x45 + x67) + x61 + x65 - T{72.0000000000000}*f_curr[16] - T{72.0000000000000}*f_curr[2]; + T strain = T{0.0277777777777778}*sqrt(x31*x31 + x48*x48 + x58*x58 + T{0.500000000000000}*(x63*x63) + T{0.500000000000000}*(x66*x66) + T{0.500000000000000}*(x68*x68)); float vorticity = cell_curl_norm[gid]; float q = vorticity*vorticity - strain*strain; |