/dports/math/mfem/mfem-4.3/fem/ |
H A D | kernels.hpp | 793 DeviceCube QQQ(sQQQ, Q1D, Q1D, Q1D); in EvalZ() 811 const DeviceCube QQQ(sQQQ, Q1D, Q1D, Q1D); in PullEval() 927 DeviceCube XxBBB(sQQQ[0], Q1D, Q1D, Q1D); in EvalZ() 928 DeviceCube XyBBB(sQQQ[1], Q1D, Q1D, Q1D); in EvalZ() 929 DeviceCube XzBBB(sQQQ[2], Q1D, Q1D, Q1D); in EvalZ() 977 DeviceCube XxBBB(sQQQ[0], Q1D, Q1D, Q1D); in PushEval() 978 DeviceCube XyBBB(sQQQ[1], Q1D, Q1D, Q1D); in PushEval() 979 DeviceCube XzBBB(sQQQ[2], Q1D, Q1D, Q1D); in PushEval() 1239 DeviceCube XxBBG(sQQQ[0], Q1D, Q1D, Q1D); in GradZ() 1240 DeviceCube XxBGB(sQQQ[1], Q1D, Q1D, Q1D); in GradZ() [all …]
|
H A D | bilininteg_hdiv.cpp | 33 const int NQ = Q1D*Q1D; in PAHdivSetup2D() 65 const int NQ = Q1D*Q1D*Q1D; in PAHdivSetup3D() 297 auto op = Reshape(op_.Read(), Q1D, Q1D, Q1D, 6, NE); in PAHdivMassAssembleDiagonal3D() 369 auto op = Reshape(op_.Read(), Q1D, Q1D, Q1D, 6, NE); in PAHdivMassApply3D() 549 const int NQ = Q1D*Q1D; in PADivDivSetup2D() 575 const int NQ = Q1D*Q1D*Q1D; in PADivDivSetup3D() 740 auto op = Reshape(op_.Read(), Q1D, Q1D, Q1D, NE); in PADivDivApply3D() 1040 auto op = Reshape(op_.Read(), Q1D, Q1D, Q1D, NE); in PADivDivAssembleDiagonal3D() 1132 const int NQ = Q1D*Q1D*Q1D; in PADivL2Setup3D() 1246 auto op = Reshape(op_.Read(), Q1D, Q1D, Q1D, NE); in PAHdivL2Apply3D() [all …]
|
H A D | bilininteg_mass_pa.cpp | 101 MFEM_FORALL_2D(e, NE, Q1D,Q1D,1, in AssemblePA() 123 const auto W = Reshape(ir->GetWeights().Read(), Q1D,Q1D,Q1D); in AssemblePA() 124 const auto J = Reshape(geom->J.Read(), Q1D,Q1D,Q1D,3,3,NE); in AssemblePA() 126 Reshape(coeff.Read(), Q1D,Q1D,Q1D,NE); in AssemblePA() 127 auto v = Reshape(pa_data.Write(), Q1D,Q1D,Q1D,NE); in AssemblePA() 128 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in AssemblePA() 282 auto D = Reshape(d.Read(), Q1D, Q1D, Q1D, NE); in PAMassAssembleDiagonal3D() 353 auto D = Reshape(d_.Read(), Q1D, Q1D, Q1D, NE); in SmemPAMassAssembleDiagonal3D() 355 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in SmemPAMassAssembleDiagonal3D() 802 auto D = Reshape(d_.Read(), Q1D, Q1D, Q1D, NE); in PAMassApply3D() [all …]
|
H A D | bilininteg_divergence.cpp | 31 const int NQ = Q1D*Q1D; in PADivergenceSetup2D() 61 const int NQ = Q1D*Q1D*Q1D; in PADivergenceSetup3D() 183 auto op = Reshape(op_.Read(), Q1D*Q1D, 2,2, NE); in PADivergenceApply2D() 321 auto op = Reshape(op_.Read(), Q1D*Q1D, 2,2, NE); in PADivergenceApplyTranspose2D() 437 auto op = Reshape(op_.Read(), Q1D*Q1D*Q1D, 3,3, NE); in PADivergenceApply3D() 542 const int q = qx + (qy + qz * Q1D) * Q1D; in PADivergenceApply3D() 620 auto op = Reshape(op_.Read(), Q1D*Q1D*Q1D, 3,3, NE); in PADivergenceApplyTranspose3D() 700 const int q = qx + (qy + qz * Q1D) * Q1D; in PADivergenceApplyTranspose3D() 800 auto Q = Reshape(q_.Read(), Q1D*Q1D*Q1D, 3,3, NE); in SmemPADivergenceApply3D() 804 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in SmemPADivergenceApply3D() [all …]
|
H A D | bilininteg_vecmass.cpp | 117 const int Q1D = T_Q1D ? T_Q1D : q1d; in PAVectorMassApply2D() local 120 MFEM_VERIFY(Q1D <= MAX_Q1D, ""); in PAVectorMassApply2D() 123 auto op = Reshape(op_.Read(), Q1D, Q1D, NE); in PAVectorMassApply2D() 214 const int Q1D = T_Q1D ? T_Q1D : q1d; in PAVectorMassApply3D() local 217 MFEM_VERIFY(Q1D <= MAX_Q1D, ""); in PAVectorMassApply3D() 220 auto op = Reshape(op_.Read(), Q1D, Q1D, Q1D, NE); in PAVectorMassApply3D() 394 MFEM_VERIFY(Q1D <= MAX_Q1D, ""); in PAVectorMassAssembleDiagonal2D() 396 auto op = Reshape(op_.Read(), Q1D, Q1D, NE); in PAVectorMassAssembleDiagonal2D() 406 for (int qx = 0; qx < Q1D; ++qx) in PAVectorMassAssembleDiagonal2D() 446 MFEM_VERIFY(Q1D <= MAX_Q1D, ""); in PAVectorMassAssembleDiagonal3D() [all …]
|
H A D | bilininteg_diffusion_pa.cpp | 101 MFEM_FORALL_2D(e, NE, Q1D,Q1D,1, in PADiffusionSetup2D() 168 MFEM_FORALL_2D(e, NE, Q1D,Q1D,1, in PADiffusionSetup2D() 208 const auto W = Reshape(w.Read(), Q1D,Q1D,Q1D); in PADiffusionSetup3D() 209 const auto J = Reshape(j.Read(), Q1D,Q1D,Q1D,3,3,NE); in PADiffusionSetup3D() 211 Reshape(c.Read(), coeffDim,Q1D,Q1D,Q1D,NE); in PADiffusionSetup3D() 212 auto D = Reshape(d.Write(), Q1D,Q1D,Q1D, symmetric ? 6 : 9, NE); in PADiffusionSetup3D() 213 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in PADiffusionSetup3D() 780 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in SmemPADiffusionDiagonal3D() 1676 double u[Q1D], v[Q1D], w[Q1D]; in SmemPADiffusionApply3D() 1734 double u[Q1D], v[Q1D], w[Q1D]; in SmemPADiffusionApply3D() [all …]
|
H A D | bilininteg_gradient.cpp | 80 const int NQ = Q1D*Q1D; in PAGradientSetup2D() 115 const int NQ = Q1D*Q1D*Q1D; in PAGradientSetup3D() 273 MFEM_VERIFY(Q1D <= MAX_Q1D, ""); in PAGradientApply2D() 277 auto op = Reshape(op_.Read(), Q1D*Q1D, 2,2, NE); in PAGradientApply2D() 407 auto op = Reshape(op_.Read(), Q1D*Q1D*Q1D, 3,3, NE); in PAGradientApply3D() 498 const int q = qx + (qy + qz * Q1D) * Q1D; in PAGradientApply3D() 599 MFEM_VERIFY(TR_D1D <= Q1D, ""); in SmemPAGradientApply3D() 600 MFEM_VERIFY(TE_D1D <= Q1D, ""); in SmemPAGradientApply3D() 606 auto D = Reshape(d_.Read(), Q1D*Q1D*Q1D, 3, 3, NE); in SmemPAGradientApply3D() 610 MFEM_FORALL_3D(e, NE, (Q1D>8)?8:Q1D, (Q1D>8)?8:Q1D, (Q1D>8)?8:Q1D, in SmemPAGradientApply3D() [all …]
|
H A D | bilininteg_vecdiffusion.cpp | 32 const int NQ = Q1D*Q1D; in PAVectorDiffusionSetup2D() 69 const int NQ = Q1D*Q1D*Q1D; in PAVectorDiffusionSetup3D() 271 MFEM_VERIFY(Q1D <= MAX_Q1D, ""); in PAVectorDiffusionApply2D() 276 auto D = Reshape(d_.Read(), Q1D*Q1D, 3, NE); in PAVectorDiffusionApply2D() 397 auto op = Reshape(op_.Read(), Q1D*Q1D*Q1D, 6, NE); in PAVectorDiffusionApply3D() 486 const int q = qx + (qy + qz * Q1D) * Q1D; in PAVectorDiffusionApply3D() 579 const int Q1D = quad1D; in AddMultPA() local 625 auto D = Reshape(d.Read(), Q1D*Q1D, 3, NE); in PAVectorDiffusionDiagonal2D() 690 MFEM_VERIFY(Q1D <= MQ1, ""); in PAVectorDiffusionDiagonal3D() 693 auto Q = Reshape(d.Read(), Q1D*Q1D*Q1D, 6, NE); in PAVectorDiffusionDiagonal3D() [all …]
|
/dports/math/mfem/mfem-4.3/fem/qinterp/ |
H A D | grad.hpp | 51 const auto j = Reshape(j_, Q1D, Q1D, 2, 2, NE); in Derivatives2D() 54 Reshape(y_, Q1D, Q1D, VDIM, 2, NE): in Derivatives2D() 55 Reshape(y_, VDIM, 2, Q1D, Q1D, NE); in Derivatives2D() 57 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in Derivatives2D() 158 const auto j = Reshape(j_, Q1D, Q1D, Q1D, 3, 3, NE); in Derivatives3D() 161 Reshape(y_, Q1D, Q1D, Q1D, VDIM, 3, NE): in Derivatives3D() 162 Reshape(y_, VDIM, 3, Q1D, Q1D, Q1D, NE); in Derivatives3D() 164 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in Derivatives3D() 182 DeviceTensor<3> DQQ0(sm1[0], D1D, Q1D, Q1D); in Derivatives3D() 183 DeviceTensor<3> DQQ1(sm1[1], D1D, Q1D, Q1D); in Derivatives3D() [all …]
|
H A D | eval.hpp | 45 const int Q1D = T_Q1D ? T_Q1D : q1d; in Values2D() local 48 const auto b = Reshape(b_, Q1D, D1D); in Values2D() 51 Reshape(y_, Q1D, Q1D, VDIM, NE): in Values2D() 52 Reshape(y_, VDIM, Q1D, Q1D, NE); in Values2D() 54 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in Values2D() 68 ConstDeviceMatrix B(sB, D1D,Q1D); in Values2D() 80 MFEM_FOREACH_THREAD(qy,y,Q1D) in Values2D() 113 Reshape(y_, Q1D, Q1D, Q1D, VDIM, NE): in Values3D() 114 Reshape(y_, VDIM, Q1D, Q1D, Q1D, NE); in Values3D() 116 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in Values3D() [all …]
|
H A D | det.cpp | 43 const int Q1D = T_Q1D ? T_Q1D : q1d; in Det2D() local 45 const auto B = Reshape(b, Q1D, D1D); in Det2D() 46 const auto G = Reshape(g, Q1D, D1D); in Det2D() 48 auto Y = Reshape(y, Q1D, Q1D, NE); in Det2D() 50 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in Det2D() 68 MFEM_FOREACH_THREAD(qy,y,Q1D) in Det2D() 70 MFEM_FOREACH_THREAD(qx,x,Q1D) in Det2D() 105 auto Y = Reshape(y, Q1D, Q1D, Q1D, NE); in Det3D() 114 MFEM_FORALL_3D_GRID(e, NE, Q1D, Q1D, Q1D, GRID, in Det3D() 134 MFEM_FOREACH_THREAD(qz,z,Q1D) in Det3D() [all …]
|
/dports/math/mfem/mfem-4.3/fem/tmop/ |
H A D | tmop_pa_p3_c0.cpp | 40 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 44 Reshape(c0_.Read(), Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 46 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 49 const auto W = Reshape(w_.Read(), Q1D, Q1D, Q1D); in MFEM_REGISTER_TMOP_KERNELS() 55 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 58 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 103 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 105 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 107 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 144 const int Q1D = PA.maps->nqpt; in AddMultPA_C0_3D() local [all …]
|
H A D | tmop_pa_w3_c0.cpp | 41 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 45 Reshape(c0_.Read(), Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 47 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 50 const auto W = Reshape(w_.Read(), Q1D, Q1D, Q1D); in MFEM_REGISTER_TMOP_KERNELS() 54 auto E = Reshape(energy.Write(), Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 56 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 59 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 104 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 106 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 136 const int Q1D = PA.maps->nqpt; in GetLocalStateEnergyPA_C0_3D() local [all …]
|
H A D | tmop_pa_h3s_c0.cpp | 35 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 40 Reshape(c0_.Read(), Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 42 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 43 const auto W = Reshape(w_.Read(), Q1D, Q1D, Q1D); in MFEM_REGISTER_TMOP_KERNELS() 46 auto H0 = Reshape(h0_.Write(), DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 48 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 75 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 77 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 79 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 115 const int Q1D = PA.maps_lim->nqpt; in AssembleGradPA_C0_3D() local [all …]
|
H A D | tmop_pa_tc3.cpp | 32 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 35 auto J = Reshape(j_.Write(), DIM,DIM, Q1D,Q1D,Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 37 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 40 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 42 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 74 auto J = Reshape(j_.Write(), DIM,DIM, Q1D,Q1D,Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 76 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 97 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 99 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 139 const int Q1D = maps.nqpt; in ComputeAllElementTargets() local [all …]
|
H A D | tmop_pa_h3m.cpp | 33 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 35 const auto b = Reshape(b_.Read(), Q1D, D1D); in MFEM_REGISTER_TMOP_KERNELS() 37 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 39 const auto H = Reshape(h_.Read(), DIM, DIM, DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 42 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 46 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 63 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 65 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 67 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 121 const int Q1D = PA.maps->nqpt; in AddMultGradPA_3D() local [all …]
|
H A D | tmop_pa_h3m_c0.cpp | 33 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 35 const auto H0 = Reshape(h0_.Read(), DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 36 const auto b = Reshape(b_.Read(), Q1D, D1D); in MFEM_REGISTER_TMOP_KERNELS() 41 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 45 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 57 kernels::internal::LoadB<MD1,MQ1>(D1D,Q1D,b,B); in MFEM_REGISTER_TMOP_KERNELS() 63 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 65 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 67 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 102 const int Q1D = PA.maps->nqpt; in AddMultGradPA_C0_3D() local [all …]
|
H A D | tmop_pa_p2_c0.cpp | 42 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 46 Reshape(c0_.Read(), Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 48 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 49 const auto b = Reshape(b_.Read(), Q1D, D1D); in MFEM_REGISTER_TMOP_KERNELS() 51 const auto W = Reshape(w_.Read(), Q1D, Q1D); in MFEM_REGISTER_TMOP_KERNELS() 57 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in MFEM_REGISTER_TMOP_KERNELS() 60 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 96 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 98 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 134 const int Q1D = PA.maps->nqpt; in AddMultPA_C0_2D() local [all …]
|
H A D | tmop_pa_w2_c0.cpp | 43 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 47 Reshape(c0_.Read(), Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 49 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 52 const auto W = Reshape(w_.Read(), Q1D, Q1D); in MFEM_REGISTER_TMOP_KERNELS() 56 auto E = Reshape(energy.Write(), Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 58 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in MFEM_REGISTER_TMOP_KERNELS() 61 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 97 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 99 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 123 const int Q1D = PA.maps->nqpt; in GetLocalStateEnergyPA_C0_2D() local [all …]
|
H A D | tmop_pa_tc2.cpp | 34 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 37 auto J = Reshape(j_.Write(), DIM,DIM, Q1D,Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 39 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in MFEM_REGISTER_TMOP_KERNELS() 42 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 44 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 74 auto J = Reshape(j_.Write(), DIM,DIM, Q1D,Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 76 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in MFEM_REGISTER_TMOP_KERNELS() 96 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 98 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 133 const int Q1D = maps.nqpt; in ComputeAllElementTargets() local [all …]
|
H A D | tmop_pa_jp3.cpp | 32 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 34 const auto b = Reshape(b_.Read(), Q1D, D1D); in MFEM_REGISTER_TMOP_KERNELS() 35 const auto g = Reshape(g_.Read(), Q1D, D1D); in MFEM_REGISTER_TMOP_KERNELS() 38 auto E = Reshape(DetJ.Write(), Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 40 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 43 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 60 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 62 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 64 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 90 const int Q1D = maps.nqpt; in MinDetJpr_3D() local [all …]
|
H A D | tmop_pa_h2s_c0.cpp | 36 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 41 Reshape(c0_.Read(), Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 43 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 44 const auto W = Reshape(w_.Read(), Q1D, Q1D); in MFEM_REGISTER_TMOP_KERNELS() 47 auto H0 = Reshape(h0_.Write(), DIM, DIM, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 49 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in MFEM_REGISTER_TMOP_KERNELS() 52 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 70 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 72 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 108 const int Q1D = PA.maps_lim->nqpt; in AssembleGradPA_C0_2D() local [all …]
|
H A D | tmop_pa_w3.cpp | 87 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 89 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 92 const auto W = Reshape(w_.Read(), Q1D, Q1D, Q1D); in MFEM_REGISTER_TMOP_KERNELS() 95 auto E = Reshape(energy.Write(), Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 97 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in MFEM_REGISTER_TMOP_KERNELS() 100 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 117 MFEM_FOREACH_THREAD(qz,z,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 119 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 121 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 160 const int Q1D = PA.maps->nqpt; in GetLocalStateEnergyPA_3D() local [all …]
|
H A D | tmop_pa_h2m.cpp | 35 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() local 37 const auto b = Reshape(b_.Read(), Q1D, D1D); in MFEM_REGISTER_TMOP_KERNELS() 38 const auto g = Reshape(g_.Read(), Q1D, D1D); in MFEM_REGISTER_TMOP_KERNELS() 39 const auto J = Reshape(j_.Read(), DIM, DIM, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 41 const auto H = Reshape(h_.Read(), DIM, DIM, DIM, DIM, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 44 MFEM_FORALL_2D(e, NE, Q1D, Q1D, NBZ, in MFEM_REGISTER_TMOP_KERNELS() 47 const int Q1D = T_Q1D ? T_Q1D : q1d; in MFEM_REGISTER_TMOP_KERNELS() 63 MFEM_FOREACH_THREAD(qy,y,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 65 MFEM_FOREACH_THREAD(qx,x,Q1D) in MFEM_REGISTER_TMOP_KERNELS() 118 const int Q1D = PA.maps->nqpt; in AddMultGradPA_2D() local [all …]
|
/dports/math/mfem/mfem-4.3/tests/unit/miniapps/ |
H A D | test_sedov.cpp | 260 Q1D,Q1D,Q1D,NE,3,3); in kSmemForceMult3D() 265 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in kSmemForceMult3D() 277 double (*MQQ0)[Q1D][Q1D] = (double (*)[Q1D][Q1D]) (sm1+0); in kSmemForceMult3D() 278 double (*MQQ1)[Q1D][Q1D] = (double (*)[Q1D][Q1D]) (sm1+1); in kSmemForceMult3D() 279 double (*MQQ2)[Q1D][Q1D] = (double (*)[Q1D][Q1D]) (sm1+2); in kSmemForceMult3D() 281 double (*QQQ0)[Q1D][Q1D] = (double (*)[Q1D][Q1D]) (sm0+0); in kSmemForceMult3D() 282 double (*QQQ1)[Q1D][Q1D] = (double (*)[Q1D][Q1D]) (sm0+1); in kSmemForceMult3D() 283 double (*QQQ2)[Q1D][Q1D] = (double (*)[Q1D][Q1D]) (sm0+2); in kSmemForceMult3D() 655 MFEM_FORALL_3D(e, NE, Q1D, Q1D, Q1D, in kSmemForceMultTranspose3D() 666 double (*MQQ0)[Q1D][Q1D] = (double (*)[Q1D][Q1D]) (sm1+0); in kSmemForceMultTranspose3D() [all …]
|