/dports/science/elmerfem/elmerfem-release-9.0/fem/src/ |
H A D | ElementDescription.F90 | 6027 EdgeBasis(1,:) = -EdgeBasis(1,:) 6057 EdgeBasis(k,:) = -EdgeBasis(k,:) 6082 EdgeBasis(k,:) = -EdgeBasis(k,:) 6142 EdgeBasis(1,:) = -EdgeBasis(1,:) 6157 EdgeBasis(3,:) = -EdgeBasis(3,:) 6172 EdgeBasis(5,:) = -EdgeBasis(5,:) 6187 EdgeBasis(7,:) = -EdgeBasis(7,:) 6249 EdgeBasis(1,:) = -EdgeBasis(1,:) 6263 EdgeBasis(2,:) = -EdgeBasis(2,:) 6277 EdgeBasis(3,:) = -EdgeBasis(3,:) [all …]
|
H A D | InterpolateMeshToMesh.F90 | 707 EdgeBasis = .FALSE. 712 IF (EdgeBasis) THEN 1051 IF (.NOT.EdgeBasis.OR.(EdgeBasis.AND.j<=np)) THEN 1066 IF (.NOT.EdgeBasis.OR.(EdgeBasis.AND.j<=np)) THEN 1258 EdgeBasis = .FALSE. 1264 IF( EdgeBasis ) THEN 1290 IF(EdgeBasis) THEN 1321 IF(EdgeBasis) THEN 1334 IF(EdgeBasis) THEN 1422 IF(EdgeBasis) THEN [all …]
|
H A D | Interpolation.F90 | 137 GlobalDistance, LocalDistance, EdgeBasis ) RESULT(IsInElement) 149 LOGICAL, OPTIONAL :: EdgeBasis local 287 trans = PRESENT(EdgeBasis) 288 IF(trans) trans=EdgeBasis
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/ |
H A D | DistributeSource.F90 | 115 LOGICAL :: UseQTree, TryQTree, Stat, EdgeBasis, PiolaT local 195 EdgeBasis = GetLogical( Params,'Edge Basis', Stat ) 212 EdgeBasis = .FALSE. 216 IF( m > n ) EdgeBasis = ANY( RhsSol % Perm(n+1:m) > 0 ) 217 IF( EdgeBasis ) THEN 258 IF( EdgeBasis ) THEN 259 IntegStuff = GaussPoints( Element, EdgeBasis = EdgeBasis, PReferenceElement=PiolaT,& 358 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = GSolver )
|
H A D | HelmholtzProjection.F90 | 213 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 217 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) 228 Basis=Basis, EdgeBasis=WBasis, dBasisdx=DBasis, & 519 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 523 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) 535 Basis=Basis, EdgeBasis=WBasis, dBasisdx=DBasis, &
|
H A D | DCRComplexSolve.F90 | 805 REAL(KIND=dp), ALLOCATABLE :: EdgeBasis(:), Basis(:) local 930 EdgeBasis, dBasisdx ) 951 u = SUM( EdgeBasis(1:En) * x(1:En) ) 952 v = SUM( EdgeBasis(1:En) * y(1:En) ) 953 w = SUM( EdgeBasis(1:En) * z(1:En) ) 1056 DEALLOCATE( x, y, z, EdgeBasis, Basis, dBasisdx, & 1142 REAL(KIND=dp), ALLOCATABLE :: EdgeBasis(:) local 1215 EdgeBasis, dBasisdx ) 1262 u = SUM( EdgeBasis(1:En) * x(1:En) ) 1263 v = SUM( EdgeBasis(1:En) * y(1:En) ) [all …]
|
H A D | StatElecSolve.F90 | 1607 REAL(KIND=dp), ALLOCATABLE :: EdgeBasis(:), Basis(:) local 1745 EdgeBasis, dBasisdx ) 1766 u = SUM( EdgeBasis(1:En) * x(1:En) ) 1767 v = SUM( EdgeBasis(1:En) * y(1:En) ) 1768 w = SUM( EdgeBasis(1:En) * z(1:En) ) 1779 Residual = -SUM( Flux(1:En) * EdgeBasis(1:En) ) 1859 REAL(KIND=dp), ALLOCATABLE :: EdgeBasis(:) local 1929 EdgeBasis, dBasisdx ) 1977 u = SUM( EdgeBasis(1:En) * x(1:En) ) 1978 v = SUM( EdgeBasis(1:En) * y(1:En) ) [all …]
|
H A D | FlowSolve.F90 | 1702 EdgeBasis, dEdgeBasisdx ) 1719 u = SUM( EdgeBasis(1:En) * x(1:En) ) 1720 v = SUM( EdgeBasis(1:En) * y(1:En) ) 1721 w = SUM( EdgeBasis(1:En) * z(1:En) ) 1868 EdgeBasis(:), Basis(:), dBasisdx(:,:) local 1930 EdgeBasis, dBasisdx ) 1937 u = SUM( EdgeBasis(1:En) * EdgeNodes % x(1:En) ) 1972 u = SUM( EdgeBasis(1:En) * x(1:En) ) 1973 v = SUM( EdgeBasis(1:En) * y(1:En) ) 1974 w = SUM( EdgeBasis(1:En) * z(1:En) ) [all …]
|
H A D | EMWaveSolver.F90 | 225 LOGICAL :: PiolaVersion, SecondOrder, EdgeBasis local 276 EdgeBasis = .NOT. ListCheckPresent( SolverParams,'Linear System Refactorize' ) .AND. & 290 IF( EdgeBasis ) CALL ListAddLogical( SolverParams,'Linear System Refactorize',.FALSE.) 292 IF ( EdgeBasis ) CALL ListRemove( SolverParams, 'Linear System Refactorize' ) 436 IP = GaussPoints(Element, EdgeBasis = .TRUE., PReferenceElement = PiolaVersion) 441 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 530 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) 539 Basis = Basis, EdgeBasis = Wbasis, dBasisdx = dBasisdx, ApplyPiolaTransform = .TRUE.) 542 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 984 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) [all …]
|
H A D | HeatSolve.F90 | 2325 EdgeBasis, dBasisdx ) 2346 u = SUM( EdgeBasis(1:En) * x(1:En) ) 2347 v = SUM( EdgeBasis(1:En) * y(1:En) ) 2348 w = SUM( EdgeBasis(1:En) * z(1:En) ) 2359 Residual = -SUM( Flux(1:En) * EdgeBasis(1:En) ) 2517 EdgeBasis, dBasisdx ) 2524 u = SUM( EdgeBasis(1:En) * EdgeNodes % x(1:En) ) 2525 v = SUM( EdgeBasis(1:En) * EdgeNodes % y(1:En) ) 2565 u = SUM( EdgeBasis(1:En) * x(1:En) ) 2566 v = SUM( EdgeBasis(1:En) * y(1:En) ) [all …]
|
H A D | VectorHelmholtz.F90 | 209 LOGICAL :: PiolaVersion, EdgeBasis, LowFrequencyModel local 259 EdgeBasis = .NOT. ListCheckPresent( SolverParams,'Linear System Refactorize' ) .AND. & 266 IF( EdgeBasis ) CALL ListAddLogical( SolverParams,'Linear System Refactorize',.FALSE.) 268 IF ( EdgeBasis ) CALL ListRemove( SolverParams, 'Linear System Refactorize' ) 512 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) 518 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 649 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) 655 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 1154 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) 1168 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver )
|
H A D | ElasticSolve.F90 | 4690 EdgeBasis, dEdgeBasisdx ) 4707 u = SUM( EdgeBasis(1:En) * x(1:En) ) 4708 v = SUM( EdgeBasis(1:En) * y(1:En) ) 4709 w = SUM( EdgeBasis(1:En) * z(1:En) ) 4822 REAL(KIND=dp), ALLOCATABLE :: EdgeBasis(:), Basis(:) local 4896 EdgeBasis, dBasisdx ) 4903 u = SUM( EdgeBasis(1:En) * EdgeNodes % x(1:En) ) 4938 u = SUM( EdgeBasis(1:En) * x(1:En) ) 4939 v = SUM( EdgeBasis(1:En) * y(1:En) ) 4940 w = SUM( EdgeBasis(1:En) * z(1:En) ) [all …]
|
H A D | StressSolve.F90 | 2936 EdgeBasis, dEdgeBasisdx ) 2941 u = SUM( EdgeBasis(1:En) * EdgeNodes % x(1:En) ) 2942 v = SUM( EdgeBasis(1:En) * EdgeNodes % y(1:En) ) 2943 w = SUM( EdgeBasis(1:En) * EdgeNodes % z(1:En) ) 2953 u = SUM( EdgeBasis(1:En) * x(1:En) ) 2954 v = SUM( EdgeBasis(1:En) * y(1:En) ) 2955 w = SUM( EdgeBasis(1:En) * z(1:En) ) 3095 EdgeBasis, dBasisdx ) 3135 u = SUM( EdgeBasis(1:En) * x(1:En) ) 3136 v = SUM( EdgeBasis(1:En) * y(1:En) ) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/ |
H A D | EnthalpySolver.F90 | 2308 EdgeBasis, dBasisdx ) 2329 u = SUM( EdgeBasis(1:En) * x(1:En) ) 2330 v = SUM( EdgeBasis(1:En) * y(1:En) ) 2331 w = SUM( EdgeBasis(1:En) * z(1:En) ) 2342 Residual = -SUM( Flux(1:En) * EdgeBasis(1:En) ) 2500 EdgeBasis, dBasisdx ) 2507 u = SUM( EdgeBasis(1:En) * EdgeNodes % x(1:En) ) 2508 v = SUM( EdgeBasis(1:En) * EdgeNodes % y(1:En) ) 2548 u = SUM( EdgeBasis(1:En) * x(1:En) ) 2549 v = SUM( EdgeBasis(1:En) * y(1:En) ) [all …]
|
H A D | FabricSolve.F90 | 1143 REAL(KIND=dp) :: EdgeBasis(n), EdgedBasisdx(n,3), EdgeddBasisddx(n,3,3) local 1184 EdgeBasis, EdgedBasisdx, EdgeddBasisddx, .FALSE. ) 1193 CALL FindParentUVW( Edge,n,LeftParent,n1,U,V,W,EdgeBasis ) 1197 CALL FindParentUVW( Edge,n,RightParent,n2,U,V,W,EdgeBasis ) 1211 cu(i) = SUM( (Velo(i,1:n)-MeshVelo(i,1:n)) * EdgeBasis(1:n) )
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/MagnetoDynamics/ |
H A D | WhitneyAVSolver.F90 | 251 PiolaVersion, SecondOrder, LFact, LFactFound, EdgeBasis, & local 501 IF ( EdgeBasis ) CALL ListRemove( SolverParams, 'Linear System Refactorize' ) 1197 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 1258 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 1448 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 1455 IP % W(t), DetF = DetJ, Basis = Basis, EdgeBasis = WBasis, & 1600 IP % W(t), DetF = DetJ, Basis = Basis, EdgeBasis = WBasis, & 1950 IP % W(t), DetF = DetJ, Basis = Basis, EdgeBasis = WBasis, & 2019 IP % W(t), DetF = DetJ, Basis = Basis, EdgeBasis = WBasis, & 2202 IP % W(t), DetF = DetJ, Basis = Basis, EdgeBasis = WBasis, & [all …]
|
H A D | WhitneyAVHarmonicSolver.F90 | 275 EdgeBasis = .NOT. LFactFound .AND. GetLogical( SolverParams, 'Edge Basis', Found ) 289 IF ( EdgeBasis ) CALL ListRemove( SolverParams, 'Linear System Refactorize' ) 821 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 886 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 946 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 1020 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 1179 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 1187 IP % W(t), DetF = DetJ, Basis = Basis, EdgeBasis = WBasis, & 1445 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 1452 IP % W(t), DetF = DetJ, Basis = Basis, EdgeBasis = WBasis, & [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/SecondOrderEdgeElement2D/ |
H A D | EdgeElementTest.F90 | 230 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 253 IP % W(t), DetF=detJ, Basis=Basis, EdgeBasis=EBasis, RotBasis=CurlEBasis, & 353 IP % W(t), DetF=detJ, Basis=Basis, EdgeBasis=EBasis, & 359 IP % W(t), DetF=detJ, Basis=Basis, EdgeBasis=EBasis, & 446 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, EdgeBasisDegree=2)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/SaveData/ |
H A D | SaveLine.F90 | 112 SkipBoundaryInfo, GotDivisions, EdgeBasis, DG local 234 EdgeBasis = .FALSE. 238 IF( EdgeBasis ) THEN 633 LOGICAL :: UseGivenNode, PiolaVersion, EdgeBasis local 690 EdgeBasis = .FALSE. 710 EdgeBasis = .FALSE. 732 IF( EdgeBasis ) THEN 745 IF( EdgeBasis ) THEN 776 IF( EdgeBasis ) THEN 1893 EdgeBasis = .FALSE. [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/SecondOrderEdgeElement3D/ |
H A D | EdgeElementTest.F90 | 206 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 229 IP % W(t), DetF=detJ, Basis=Basis, EdgeBasis=EBasis, RotBasis=CurlEBasis, & 329 IP % W(t), DetF=detJ, Basis=Basis, EdgeBasis=EBasis, & 335 IP % W(t), detF=detJ, Basis=Basis, EdgeBasis=EBasis, & 425 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, EdgeBasisDegree=2)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/EdgeElementInterpolation/ |
H A D | EdgeElementSolve.f90 | 226 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion) 248 IP % W(t), DetF=detJ, Basis=Basis, EdgeBasis=EBasis, RotBasis=CurlEBasis, & 348 IP % W(t), DetF=detJ, Basis=Basis, EdgeBasis=EBasis, & 354 IP % W(t), detF=detJ, Basis=Basis, EdgeBasis=EBasis, & 442 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion)
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/HelmholtzEdge/ |
H A D | Hrot.f90 | 237 CALL EdgeBasis( WhitneyBasis, RotWhitneyBasis, & 325 CALL EdgeBasis( WhitneyBasis, RotWhitneyBasis, u, v, w, Element, n ) 503 SUBROUTINE EdgeBasis( WhitneyBasis, RotWhitneyBasis, u, v, w, Element, n ) subroutine 547 END SUBROUTINE EdgeBasis
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/SecondOrderEdgeElement2D_BCs/ |
H A D | EdgeFEInterpolationTest.F90 | 211 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, & 219 IP % W(t), detF=detJ, Basis=Basis, EdgeBasis=EBasis, & 327 …IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=PiolaVersion, EdgeBasisDegree=Elemen…
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/FaceElement2D_BCs/ |
H A D | FaceFEInterpolationTest.F90 | 208 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=.TRUE., & 303 …IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=.TRUE., EdgeBasisDegree=ElementOrder…
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/FaceElement3D_BCs/ |
H A D | FaceFEInterpolationTest.F90 | 208 IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=.TRUE., & 303 …IP = GaussPoints(Element, EdgeBasis=.TRUE., PReferenceElement=.TRUE., EdgeBasisDegree=ElementOrder…
|