/dports/science/getdp/getdp-3.4.0-source/Kernel/ |
H A D | Cal_AnalyticIntegration.cpp | 19 double DetJ ; in Cal_AnalyticIntegration() local 28 DetJ = (E->x[2] - E->x[0]) * (E->y[1] - E->y[0]) - in Cal_AnalyticIntegration() 38 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() 42 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() 46 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() 57 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() 61 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() 65 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() 76 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() 80 / (2. * DetJ * DetJ) ) ; in Cal_AnalyticIntegration() [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/ModelPDEevol/ |
H A D | ModelPDEevol.F90 | 169 REAL(KIND=dp), ALLOCATABLE, SAVE :: Basis(:,:),dBasisdx(:,:,:), DetJ(:) local 213 DEALLOCATE(Basis,dBasisdx, DetJ, MASS, STIFF, FORCE, VeloCoeff ) 218 ALLOCATE(Basis(ngp,nd), dBasisdx(ngp,nd,3), DetJ(ngp), & 241 DetJ(t) = IP % s(t) * Detj(t) 248 …CALL LinearForms_GradUdotGradU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, DetJ, STIFF, DiffCoe… 255 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, STIFF, ReactCoeff ) 262 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, MASS, TimeCoeff ) 268 CALL LinearForms_UdotF(ngp, nd, Basis, DetJ, SourceCoeff, FORCE) 289 … CALL LinearForms_GradUdotU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, Basis, DetJ, STIFF, & 313 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ local [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/Lua/ |
H A D | ModelPDEevol.F90 | 140 REAL(KIND=dp), ALLOCATABLE, SAVE :: Basis(:,:),dBasisdx(:,:,:), DetJ(:) local 185 DEALLOCATE(Basis,dBasisdx, DetJ, MASS, STIFF, FORCE, VeloCoeff ) 190 ALLOCATE(Basis(ngp,nd), dBasisdx(ngp,nd,3), DetJ(ngp), & 213 DetJ(t) = IP % s(t) * Detj(t) 220 …CALL LinearForms_GradUdotGradU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, DetJ, STIFF, DiffCoe… 227 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, STIFF, ReactCoeff ) 234 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, MASS, TimeCoeff ) 240 CALL LinearForms_UdotF(ngp, nd, Basis, DetJ, SourceCoeff, FORCE) 261 … CALL LinearForms_GradUdotU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, Basis, DetJ, STIFF, & 285 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local [all …]
|
H A D | ModelPDEthreaded.F90 | 141 REAL(KIND=dp), ALLOCATABLE, SAVE :: Basis(:,:),dBasisdx(:,:,:), DetJ(:) local 165 DEALLOCATE(Basis,dBasisdx, DetJ, Coeffs, CoeffsInBasis, & 171 ALLOCATE(Basis(ngp,nd), dBasisdx(ngp,nd,3), DetJ(ngp), & 257 DetJ(t) = IP % s(t)*Detj(t) 262 …CALL LinearForms_GradUdotGradU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, DetJ, STIFF, diff_co… 265 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, STIFF, ReactInBasis) 268 CALL LinearForms_GradUdotU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, Basis, DetJ, STIFF, & 272 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, MASS, rho) 275 CALL LinearForms_UdotF(ngp, nd, Basis, DetJ, LoadAtIP, FORCE) 328 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/ModelPDEthreaded/ |
H A D | ModelPDEthreaded.F90 | 141 REAL(KIND=dp), ALLOCATABLE, SAVE :: Basis(:,:),dBasisdx(:,:,:), DetJ(:) local 165 DEALLOCATE(Basis,dBasisdx, DetJ, Coeffs, CoeffsInBasis, & 171 ALLOCATE(Basis(ngp,nd), dBasisdx(ngp,nd,3), DetJ(ngp), & 257 DetJ(t) = IP % s(t)*Detj(t) 262 …CALL LinearForms_GradUdotGradU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, DetJ, STIFF, DiffInB… 265 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, STIFF, ReactInBasis) 268 CALL LinearForms_GradUdotU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, Basis, DetJ, STIFF, & 272 CALL LinearForms_UdotU(ngp, nd, Element % TYPE % DIMENSION, Basis, DetJ, MASS, rho) 275 CALL LinearForms_UdotF(ngp, nd, Basis, DetJ, LoadAtIP, FORCE) 328 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/PoissonThreaded3/ |
H A D | Poisson.F90 | 100 REAL(KIND=dp), ALLOCATABLE, SAVE :: Basis(:,:), dBasisdx(:,:,:), DetJ(:), LoadAtIPs(:) 117 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 122 DEALLOCATE(Basis, dBasisdx, DetJ, LoadAtIPs) 124 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 132 …stat = ElementInfoVec( Element, Nodes, ngp, IP % U, IP % V, IP % W, DetJ, SIZE(Basis,2), Basis, dB… 134 DetJ(1:ngp) = IP % s(1:ngp)*Detj(1:ngp) 137 CALL LinearForms_GradUdotGradU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, DetJ, STIFF) 145 CALL LinearForms_UdotF(ngp, nd, Basis, DetJ, LoadAtIPs, FORCE)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/PoissonThreaded2/ |
H A D | Poisson.F90 | 127 REAL(KIND=dp), ALLOCATABLE, SAVE :: Basis(:,:), dBasisdx(:,:,:), DetJ(:), LoadAtIPs(:) 144 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 150 DEALLOCATE(Basis, dBasisdx, DetJ, LoadAtIPs) 152 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 162 IP % U, IP % V, IP % W, DetJ, SIZE(Basis,2), Basis, dBasisdx ) 164 DetJ(1:ngp) = IP % s(1:ngp)*Detj(1:ngp) 167 CALL LinearForms_GradUdotGradU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, DetJ, STIFF) 175 CALL LinearForms_UdotF(ngp, nd, Basis, DetJ, LoadAtIPs, FORCE)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/PoissonThreaded/ |
H A D | Poisson.F90 | 157 REAL(KIND=dp), ALLOCATABLE, SAVE :: Basis(:,:), dBasisdx(:,:,:), DetJ(:), LoadAtIPs(:) 174 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 179 DEALLOCATE(Basis, dBasisdx, DetJ, LoadAtIPs) 181 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 189 …stat = ElementInfoVec( Element, Nodes, ngp, IP % U, IP % V, IP % W, DetJ, SIZE(Basis,2), Basis, dB… 191 DetJ(1:ngp) = IP % s(1:ngp)*Detj(1:ngp) 194 CALL LinearForms_GradUdotGradU(ngp, nd, Element % TYPE % DIMENSION, dBasisdx, DetJ, STIFF) 202 CALL LinearForms_UdotF(ngp, nd, Basis, DetJ, LoadAtIPs, FORCE)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/ |
H A D | StreamSolver.F90 | 262 REAL(KIND=dp) :: DetJ, U, V, W, S, Radius local 287 stat = ElementInfo( Element, Nodes, U, V, W, DetJ, Basis, dBasisdx ) 288 s = s * DetJ 333 REAL(KIND=dp) :: DetJ, U, V, W, S, L, Normal(3) local 358 stat = ElementInfo( Element, Nodes, U, V, W, DetJ, & 360 s = s * DetJ
|
H A D | Poisson.F90 | 125 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 152 STIFF(1:nd,1:nd) = STIFF(1:nd,1:nd) + IP % s(t) * DetJ * & 154 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd) 195 FORCE(j) = FORCE(j) + IP % s(t)*DetJ*LoadAtIP*Basis(j)
|
H A D | ElementSizeSolver.F90 | 133 REAL(KIND=dp) :: Basis(n),DetJ,LoadAtIP,Weight local 157 LoadAtIP = DetJ ** Power 161 Weight = IP % s(t) * DetJ
|
H A D | ModelPDE.F90 | 116 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 174 Weight = IP % s(t) * DetJ 217 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 251 Weight = IP % s(t) * DetJ
|
H A D | WaveSolver.F90 | 170 REAL(KIND=dp) :: Weight, Basis(nd), dBasisdx(nd,3), DetJ local 211 Weight = IP % s(t) * DetJ 273 REAL(KIND=dp) :: Weight, Basis(nd), dBasisdx(nd,3), DetJ local 309 Weight = IP % s(t) * DetJ
|
H A D | DirectionSolver.F90 | 333 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 360 STIFF(1:nd,1:nd) = STIFF(1:nd,1:nd) + IP % s(t) * DetJ * & 362 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd) 403 FORCE(j) = FORCE(j) + IP % s(t)*DetJ*LoadAtIP*Basis(j)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/1dtests/ |
H A D | Poisson.f90 | 97 REAL(KIND=dp) :: Basis(n),dBasisdx(n,3),DetJ,LoadAtIP local 124 STIFF(1:n,1:n) = STIFF(1:n,1:n) + IP % s(t) * DetJ * & 126 FORCE(1:n) = FORCE(1:n) + IP % s(t) * DetJ * LoadAtIP * Basis(1:n)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/LinearFormsAssembly/ |
H A D | LinearFormsAssembly.F90 | 292 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP,Weight 323 Weight = IP % s(t) * DetJ 339 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd) 357 DetJ(:), LoadAtIPs(:) 376 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 382 DEALLOCATE(Basis, dBasisdx, DetJ, LoadAtIPs) 384 DetJ(ngp), LoadAtIPs(ngp), STAT=allocstat) 393 IP % U, IP % V, IP % W, DetJ, SIZE(Basis,2), Basis, dBasisdx ) 398 DetJ(i) = Ip % s(i)*Detj(i) 403 dBasisdx, DetJ, STIFF) [all …]
|
/dports/math/mfem/mfem-4.3/fem/tmop/ |
H A D | tmop_pa_jp2.cpp | 26 Vector &DetJ, in MFEM_REGISTER_TMOP_KERNELS() argument 40 auto E = Reshape(DetJ.Write(), Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 72 return DetJ.Min(); in MFEM_REGISTER_TMOP_KERNELS()
|
H A D | tmop_pa_jp3.cpp | 26 Vector &DetJ, in MFEM_REGISTER_TMOP_KERNELS() argument 38 auto E = Reshape(DetJ.Write(), Q1D, Q1D, Q1D, NE); in MFEM_REGISTER_TMOP_KERNELS() 74 return DetJ.Min(); in MFEM_REGISTER_TMOP_KERNELS()
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/pmultg/ |
H A D | Poisson.f90 | 97 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 124 STIFF(1:nd,1:nd) = STIFF(1:nd,1:nd) + IP % s(t) * DetJ * & 126 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/PoissonPFEM/ |
H A D | Poisson.f90 | 97 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 124 STIFF(1:nd,1:nd) = STIFF(1:nd,1:nd) + IP % s(t) * DetJ * & 126 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/CurvedBndryPFEM2/ |
H A D | Poisson.f90 | 108 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),ddBasisddx(nd,3,3),DetJ,LoadAtIP local 135 STIFF(1:nd,1:nd) = STIFF(1:nd,1:nd) + IP % s(t) * DetJ * & 137 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/ExtrudedMesh/ |
H A D | Poisson.f90 | 99 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 126 STIFF(1:nd,1:nd) = STIFF(1:nd,1:nd) + IP % s(t) * DetJ * & 128 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/ExtrudedMeshBaseline1/ |
H A D | Poisson.f90 | 99 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 126 STIFF(1:nd,1:nd) = STIFF(1:nd,1:nd) + IP % s(t) * DetJ * & 128 FORCE(1:nd) = FORCE(1:nd) + IP % s(t) * DetJ * LoadAtIP * Basis(1:nd)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/ModelPDEhandle/ |
H A D | ModelPDEhandle.F90 | 105 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 146 Weight = IP % s(t) * DetJ 189 REAL(KIND=dp) :: Basis(nd),DetJ,LoadAtIP local 216 Weight = IP % s(t) * DetJ
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/ModelPDEipadapt/ |
H A D | ModelPDEhandle.F90 | 105 REAL(KIND=dp) :: Basis(nd),dBasisdx(nd,3),DetJ,LoadAtIP local 147 Weight = IP % s(t) * DetJ 190 REAL(KIND=dp) :: Basis(nd),DetJ,LoadAtIP local 217 Weight = IP % s(t) * DetJ
|