/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/MagnetoDynamics/ |
H A D | WhitneyAVSolver.F90 | 1217 Bx = SUM(POT(n+1:nd) * RotWBasis(1:nd-n,1)) 1218 By = SUM(POT(n+1:nd) * RotWBasis(1:nd-n,2)) 1279 Bx = SUM(POT(n+1:nd) * RotWBasis(1:nd-n,1)) 1280 By = SUM(POT(n+1:nd) * RotWBasis(1:nd-n,2)) 1281 Bz = SUM(POT(n+1:nd) * RotWBasis(1:nd-n,3)) 1781 RotWJ(1:3) = RotWBasis(j,1:3) 1818 SUM(M*RotWBasis(i,:)))*detJ*IP%s(t) 1824 + SUM(RotWBasis(i,:) * MATMUL(A_t, RotWBasis(j,:)))*detJ*IP%s(t) 1826 … STIFF(p,q) = STIFF(p,q) + mu * SUM(RotWBasis(i,:)*RotWBasis(j,:))*detJ*IP%s(t) 1849 SUM(RotWBasis(i,:)*RotWBasis(j,:)) * detJ*IP % s(t) [all …]
|
H A D | WhitneyAVHarmonicSolver.F90 | 841 Bx = SUM(POTC(n+1:nd) * RotWBasis(1:nd-n,1)) 842 By = SUM(POTC(n+1:nd) * RotWBasis(1:nd-n,2)) 906 Bx = SUM(POTC(n+1:nd) * RotWBasis(1:nd-n,1)) 907 By = SUM(POTC(n+1:nd) * RotWBasis(1:nd-n,2)) 908 Bz = SUM(POTC(n+1:nd) * RotWBasis(1:nd-n,3)) 1059 REAL(KIND=dp) :: WBasis(nd,3), RotWBasis(nd,3) local 1352 RotWJ(1:3) = RotWBasis(j,1:3) 1372 SUM(M*RotWBasis(i,:)))*detJ*IP%s(t) 1388 SUM(MATMUL(Nu, RotWBasis(j,:))*RotWBasis(i,:))*detJ*IP%s(t) 1420 REAL(KIND=dp) :: WBasis(nd,3), RotWBasis(nd,3) local [all …]
|
H A D | CalcFields.F90 | 1062 EdgeBasis=WBasis, RotBasis=RotWBasis, dBasisdx=dBasisdx, & 1067 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 1097 B(k,:) = MATMUL( SOL(k,np+1:nd), RotWBasis(1:nd-np,:) ) 1242 MATMUL(SOL(1,np+1:nd), RotWBasis(1:nd-np,:))) 1244 MATMUL(SOL(2,np+1:nd), RotWBasis(1:nd-np,:))) 1336 MATMUL(SOL(1,np+1:nd), RotWBasis(1:nd-np,:))) 2207 CALL GetEdgeBasis(Element, WBasis, RotWBasis, Basis, dBasisdx) 2609 B(k,:) = normal*sum( SOL(k,np+1:nd)* RotWBasis(1:nd-np,3) ) 3078 REAL(KIND=dp) :: WBasis(nd,3), RotWBasis(nd,3), B(2, 3), & local 3098 CALL GetEdgeBasis(Element, WBasis, RotWBasis, Basis, dBasisdx) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/ |
H A D | VectorHelmholtz.F90 | 473 REAL(KIND=dp), ALLOCATABLE :: Basis(:),dBasisdx(:,:),WBasis(:,:),RotWBasis(:,:) local 481 SAVE AllocationsDone, WBasis, RotWBasis, Basis, dBasisdx, & 486 ALLOCATE( WBasis(m,3), RotWBasis(m,3), Basis(m), dBasisdx(m,3), & 518 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 535 SUM(RotWBasis(i,:) * RotWBasis(j,:)) * weight 623 ALLOCATE( WBasis(m,3), RotWBasis(m,3), Basis(m), dBasisdx(m,3),& 655 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 963 REAL(KIND=dp) :: RotWBasis(35,3), Basis(35), dBasisdx(35,3), E(2,3) local 1168 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 1170 B = CMPLX(MATMUL( SOL(2,np+1:nd), RotWBasis(1:nd-np,:) ) / (Omega), & [all …]
|
H A D | EMWaveSolver.F90 | 405 REAL(KIND=dp) :: Basis(n), dBasisdx(n,3),WBasis(nd,3),RotWBasis(nd,3) local 441 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 463 SUM(RotWBasis(i,:) * RotWBasis(j,:)) * weight 491 REAL(KIND=dp) :: Basis(n), WBasis(nd,3), RotWBasis(nd,3), dBasisdx(n,3) local 542 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver ) 933 REAL(KIND=dp), ALLOCATABLE :: RotWBasis(:,:), Basis(:), dBasisdx(:,:) local 942 Indexes, WBasis, RotWBasis, Basis, dBasisdx 947 Indexes(n), WBasis(n,3), RotWBasis(n,3), Basis(n), dBasisdx(n,3)) 995 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = pSolver )
|
H A D | DistributeSource.F90 | 111 REAL(KIND=dp), ALLOCATABLE :: Basis(:),BasisB(:),RotWBasis(:,:), & local 173 RhsInd(n), WBasis(n,3), RotWBasis(n,3), dBasisdx(n,3), dBasisdxb(n,3) ) 358 EdgeBasis = Wbasis, RotBasis = RotWBasis, USolver = GSolver )
|
H A D | FourierLoss.F90 | 858 REAL(KIND=dp), ALLOCATABLE :: WBasis(:,:), RotWBasis(:,:) local 871 ALLOCATE( WBasis(nt,3), RotWBasis(nt,3), Indeces(nt), ElemField(nt) ) 973 CALL GetEdgeBasis(Element, WBasis, RotWBasis, Basis, dBasisdx) 1036 CurlAtIp(1) = SUM( ElemField(n+1:nt) * RotWBasis(1:(nt-n),1) ) 1037 CurlAtIp(2) = SUM( ElemField(n+1:nt) * RotWBasis(1:(nt-n),2) ) 1038 CurlAtIp(3) = SUM( ElemField(n+1:nt) * RotWBasis(1:(nt-n),3) ) 1141 DEALLOCATE( WBasis, RotWBasis )
|
H A D | CircuitsAndDynamics.F90 | 419 REAL(KIND=dp) :: WBasis(nd,3), RotWBasis(nd,3) local 479 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 549 REAL(KIND=dp) :: wBase(nn), gradv(3), WBasis(nd,3), RotWBasis(nd,3) local 609 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 678 REAL(KIND=dp) :: wBase(nn), gradv(3), WBasis(nd,3), RotWBasis(nd,3), & local 748 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 1298 REAL(KIND=dp) :: WBasis(nd,3), RotWBasis(nd,3) local 1353 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 1428 REAL(KIND=dp) :: wBase(nn), gradv(3), WBasis(nd,3), RotWBasis(nd,3) local 1482 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/ |
H A D | InterpolateMeshToMesh.F90 | 577 RotWBasis(:,:), WBasis(:,:) local 1029 ALLOCATE(dVals(k,3),WBasis(k,3),RotWBasis(k,3)) 1031 DEALLOCATE(dVals,WBasis,RotWBasis) 1032 ALLOCATE(dVals(k,3),WBasis(k,3),RotWBasis(k,3)) 1039 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Vals,dVals) 1324 Basis(n), dBasisdx(n,3), WBasis(n,3), RotWBasis(n,3) ) 1450 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 1562 IF(EdgeBasis) DEALLOCATE( dBasisdx, WBasis, RotWBasis )
|
H A D | ElementDescription.F90 | 9991 RotWBasis(i,:)=rBase 9995 RotWBasis(i,1)=(dBase(2)*du(3) - dBase(3)*du(2)) 9996 RotWBasis(i,2)=(dBase(3)*du(1) - dBase(1)*du(3)) 9997 RotWBasis(i,3)=(dBase(1)*du(2) - dBase(2)*du(1)) 10016 RotWBasis(i,1)=(dBase(2)*du(3) - dBase(3)*du(2))/n 10017 RotWBasis(i,2)=(dBase(3)*du(1) - dBase(1)*du(3))/n 10018 RotWBasis(i,3) = (dBase(1)*du(2) - dBase(2)*du(1))/n 10060 RotWBasis(i,1)=(dBase(2)*du(3) - dBase(3)*du(2))/n 10061 RotWBasis(i,2)=(dBase(3)*du(1) - dBase(1)*du(3))/n 10062 RotWBasis(i,3)=(dBase(1)*du(2) - dBase(2)*du(1))/n [all …]
|
H A D | BlockSolve.F90 | 994 REAL(KIND=dp), ALLOCATABLE :: WBasis(:,:), RotWBasis(:,:) local 998 ALLOCATE( WBasis(n,3), RotWBasis(n,3), Basis(n), dBasisDx(n,3), Indexes(n) ) 1137 RotBasis = RotWBasis, dBasisdx = dBasisdx, & 1142 CALL GetEdgeBasis(Element, WBasis, RotWBasis, Basis, dBasisdx)
|
H A D | MeshUtils.F90 | 8369 ALLOCATE( dBasisdx(n,3), WBasis(n,3), WBasisM(n,3), RotWBasis(n,3) ) 8641 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 8656 CALL GetEdgeBasis(ElementM,WBasisM,RotWBasis,BasisM,dBasisdx) 8751 DEALLOCATE( dBasisdx, WBasis, WBasisM, RotWBasis ) 8775 REAL(KIND=dp), POINTER, SAVE :: Basis(:),WBasis(:,:),RotWBasis(:,:), & local 8928 ALLOCATE( WBasis(n,3), WBasisM(n,3), RotWBasis(n,3), STAT=AllocStat ) 9686 CALL GetEdgeBasis(Element,WBasis,RotWBasis,Basis,dBasisdx) 9723 CALL GetEdgeBasis(ElementM,WBasisM,RotWBasis,BasisM,dBasisdx) 9997 DEALLOCATE( WBasis, WBasisM, RotWBasis ) 10545 Wbasis2(:,:),RotWBasis(:,:),dBasisdx(:,:) local [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/SaveData/ |
H A D | SaveLine.F90 | 642 REAL(KIND=dp) :: WBasis(54,3),RotWBasis(54,3), NodedBasisdx(54,3) local 747 … detJ, NodeBasis, NodedBasisdx, EdgeBasis = WBasis, RotBasis = RotWBasis, USolver = Var % Solver)
|