/dports/science/elmerfem/elmerfem-release-9.0/fem/src/ |
H A D | SolidMechanicsUtils.F90 | 363 Stiff(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(R(1:DOFs,1:DOFs)), & 364 MATMUL(Stiff(1:DOFs,1:DOFs),R(1:DOFs,1:DOFs))) 365 Force(1:DOFs) = MATMUL(TRANSPOSE(R(1:DOFs,1:DOFs)),Force(1:DOFs)) 368 Mass(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(R(1:DOFs,1:DOFs)), & 369 MATMUL(Mass(1:DOFs,1:DOFs),R(1:DOFs,1:DOFs))) 403 DOFs = (nd-nb)*DOFs 404 Stiff(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(R(1:DOFs,1:DOFs)), & 405 MATMUL(Stiff(1:DOFs,1:DOFs),R(1:DOFs,1:DOFs))) 406 Force(1:DOFs) = MATMUL(TRANSPOSE(R(1:DOFs,1:DOFs)),Force(1:DOFs)) 413 Mass(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(R(1:DOFs,1:DOFs)), & [all …]
|
H A D | ElementUtils.F90 | 2019 REAL(KIND=dp) :: s,R(n*DOFs,n*DOFs),Q(n*DOFs,n*DOFs),N1(3),T1(3),T2(3) 2034 R(DOFs*(i-1)+1,DOFs*(i-1)+1) = N1(1) 2043 R(DOFs*(i-1)+1,DOFs*(i-1)+1) = N1(1) 2044 R(DOFs*(i-1)+1,DOFs*(i-1)+2) = N1(2) 2045 R(DOFs*(i-1)+1,DOFs*(i-1)+3) = N1(3) 2047 R(DOFs*(i-1)+2,DOFs*(i-1)+1) = T1(1) 2048 R(DOFs*(i-1)+2,DOFs*(i-1)+2) = T1(2) 2049 R(DOFs*(i-1)+2,DOFs*(i-1)+3) = T1(3) 2051 R(DOFs*(i-1)+3,DOFs*(i-1)+1) = T2(1) 2052 R(DOFs*(i-1)+3,DOFs*(i-1)+2) = T2(2) [all …]
|
H A D | DefUtils.F90 | 2579 DOFs = x % DOFs 2583 ALLOCATE( M(DOFs*n,DOFs*n), A(DOFs*n,DOFs*n), F(DOFs*n) ) 2897 DOFs = x % DOFs 2901 ALLOCATE( M(DOFs*n,DOFs*n), A(DOFs*n,DOFs*n), B(DOFs*n,DOFs*n), F(DOFs*n) ) 3682 DOFs = x % DOFs 3708 ALLOCATE( G(DOFs*n,DOFs*n), F(DOFs*n) ) 3822 DOFs = x % DOFs 4094 DOFs = x % DOFs 4263 DOFs = x % DOFs 4419 DOFs = x % DOFs [all …]
|
H A D | MGPrec.F90 | 36 INTEGER :: i,j,k,me,n, DOFs local 69 DOFs = PSolver % Variable % DOFs 72 DOFs, PSolver, PSolver % MultiGridLevel, FirstCall(stack_pos))
|
H A D | Multigrid.F90 | 559 DO i=1,DOFs 583 DO i=1,DOFs 616 INTEGER :: DOFs local 727 DO i=1,NB/DOFs 4775 INTEGER :: DOFs local 4825 v(DOFs*(Cols(j)-1)+k) + u(DOFs*(i-1)+k) * Values(j) 4834 v(DOFs*(i-1)+k) + u(DOFs*(Cols(j)-1)+k) * Values(j) 4850 INTEGER :: DOFs local 5608 INTEGER :: DOFs local 5645 v(DOFs*(i-1)+k) = u(DOFs*(j-1)+k) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/ |
H A D | SpringAssembly.F90 | 65 INTEGER :: BC, DOFs, i, istat, j, n local 96 DOFs = SolverPtr % Variable % DOFs 97 ALLOCATE(K(DOFs,DOFs), M(DOFs,DOFs), D(DOFs,DOFs), f(DOFs), Assemble(DOFs), & 98 AssembleMass(DOFs)) 121 ALLOCATE(Work(DOFs,n), MassVals(DOFs,n), STAT=istat) 125 ALLOCATE(Work(DOFs,n), MassVals(DOFs,n), STAT=istat) 137 DO i=1,DOFs 161 DO j=1,DOFs 169 DO j=1,DOFs
|
H A D | ShellSolver.F90 | 4257 Stiff(1:DOFs,1:DOFs) = Stiff(1:DOFs,1:DOFs) + & 4265 Stiff(1:DOFs,1:DOFs) = Stiff(1:DOFs,1:DOFs) + & 4273 Stiff(1:DOFs,1:DOFs) = Stiff(1:DOFs,1:DOFs) + ( & 4281 Stiff(1:DOFs,1:DOFs) = Stiff(1:DOFs,1:DOFs) + & 4414 Stiff(1:DOFs,1:DOFs) = Stiff(1:DOFs,1:DOFs) - ( & 4422 Stiff(1:DOFs,1:DOFs) = Stiff(1:DOFs,1:DOFs) - ( & 4581 …Stiff(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(TMat(1:DOFs,1:DOFs)),MATMUL(Stiff(1:DOFs,1:DOFs),TMat(1:DO… 4582 …Stiff(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(Q(1:DOFs,1:DOFs)),MATMUL(Stiff(1:DOFs,1:DOFs),Q(1:DOFs,1:D… 4593 …Mass(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(TMat(1:DOFs,1:DOFs)),MATMUL(Mass(1:DOFs,1:DOFs),TMat(1:DOFs… 4594 …Mass(1:DOFs,1:DOFs) = MATMUL(TRANSPOSE(Q(1:DOFs,1:DOFs)),MATMUL(Mass(1:DOFs,1:DOFs),Q(1:DOFs,1:DOF… [all …]
|
H A D | RigidBodyReduction.F90 | 83 INTEGER :: DOFs, n local 562 DOFs*DOFs*( RigidNodes ) 618 = Permutation( DOFs * Perm(t) - DOFs + i ) 643 Permutation( DOFs * Perm(t) - DOFs + 1 ) 650 Permutation( DOFs * Perm(t) - DOFs + 2 ) 661 Permutation( DOFs * Perm(t) - DOFs + 2 ) 667 Permutation( DOFs * Perm(t) - DOFs + 1 ) 673 Permutation( DOFs * Perm(t) - DOFs + 1 ) 681 Permutation( DOFs * Perm(t) - DOFs + 3 ) 687 Permutation( DOFs * Perm(t) - DOFs + 3 ) [all …]
|
H A D | Smitc.F90 | 104 DOFs = Solver % Variable % DOFs 108 STIFF( DOFs*N, DOFs*N ), & 109 MASS( DOFs*N, DOFs*N ), & 110 DAMP( DOFs*N, DOFs*N ), & 299 INTEGER :: n, DOFs local 410 i = DOFs*(p-1)+1 412 j = DOFs*(q-1)+1 422 i = DOFs*(p-1)+1 429 i = DOFs*(p-1)+1 431 j = DOFs*(q-1)+1 [all …]
|
H A D | ElasticSolve.F90 | 107 DOFs = dim 2733 MassMatrix(DOFs*(p-1)+i,DOFs*(q-1)+i) & 2766 StiffMatrix(DOFs*p,DOFs*(q-1)+i) = StiffMatrix(DOFs*p,DOFs*(q-1)+i) + & 2770 StiffMatrix(DOFs*p,DOFs*(q-1)+i) = StiffMatrix(DOFs*p,DOFs*(q-1)+i) + & 2776 … StiffMatrix(DOFs*p,DOFs*(q-1)+i) = StiffMatrix(DOFs*p,DOFs*(q-1)+i) + DetDefG**2 * & 2782 … StiffMatrix(DOFs*p,DOFs*q) = StiffMatrix(DOFs*p,DOFs*q) + PressurePar * Basis(p) * Basis(q) * s 2837 StiffMatrix(DOFs*(p-1)+i,DOFs*(q-1)+j) & 2844 StiffMatrix(DOFs*(p-1)+i,DOFs*q) & 3127 … BoundaryMatrix((p-1)*DOFs+i,(q-1)*DOFs+j) = BoundaryMatrix((p-1)*DOFs+i,(q-1)*DOFs+j) + & 3144 … BoundaryMatrix((p-1)*DOFs+i,(q-1)*DOFs+j) = BoundaryMatrix((p-1)*DOFs+i,(q-1)*DOFs+j) + & [all …]
|
H A D | Spalart-Allmaras.F90 | 55 INTEGER :: i,j,k,n,iter,t,body_id,eq_id,istat,LocalNodes,bf_id,DOFs local 78 DOFs = KE % DOFs 92 ALLOCATE( MASS( 2*DOFs*N,2*DOFs*N ), & 93 STIFF( 2*DOFs*N,2*DOFs*N ),LOAD( DOFs,N ), & 94 FORCE( 2*DOFs*N ), TimeForce( 2*DOFs*N ), STAT=istat ) 170 CALL Condensate( DOFs*N, STIFF, FORCE, TimeForce )
|
H A D | CoordinateTransform.F90 | 372 DO k=1,RotMvar % DOFs 373 RotMvar % Values( RotMvar % DOFs*(RotMvar % Perm( & 379 DO k=1,alphavecvar % DOFs 380 alphavecvar % Values( alphavecvar % DOFs*(alphavecvar % Perm( & 386 DO k=1,betavecvar % DOFs 387 betavecvar % Values( betavecvar % DOFs*(betavecvar % Perm( & 393 DO k=1,gammavecvar % DOFs 394 gammavecvar % Values( gammavecvar % DOFs*(gammavecvar % Perm( & 468 DO k=1,RotMvar % DOFs 470 RotMvar % DOFs*(RotMvar % Perm(Element % DGIndexes(j))-1)+k)
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/Adjoint/ |
H A D | Adjoint_GradientValidation.F90 | 107 integer,SAVE :: DOFs local 151 DOFs = Var % DOFs 157 IF (GradVar%DOFs.NE.DOFs) & 165 IF (PVar%DOFs.NE.DOFs) & 220 DO c=1,DOFs 221 x(DOFs*(i-1)+c)=Values(DOFs*(Perm(ActiveNodes(i))-1)+c) 224 xp(DOFs*(i-1)+c)=PValues(DOFs*(PPerm(ActiveNodes(i))-1)+c) 226 xp(DOFs*(i-1)+c)=-g(DOFs*(i-1)+c) 249 Values(DOFs*(Perm(ActiveNodes(i))-1)+c)=x(DOFs*(i-1)+c)+h*xp(DOFs*(i-1)+c) 280 Do c=1,DOFs [all …]
|
H A D | Adjoint_LinearSolver.F90 | 75 INTEGER :: DOFs local 121 DOFs = Sol % DOFs 127 IF (DIM.eq.(DOFs+1)) CurrentModel % Dimension = DOFs 147 IF (DOFs>1) THEN 216 IF (VbSol % DOFs.NE.DOFs) then 218 'Variable Vb has ',VbSol % DOFs,' DOFs, should be',DOFs 245 CALL RotateNTSystemAll( Vb, VbPerm, DOFs ) 247 c = DOFs 293 IF (DOFs>1) THEN 359 IF (DOFs>1) THEN [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/ |
H A D | CalvingFrontAdvance3D.F90 | 92 DOFs = Var % DOFs 308 Advance((Perm(i)-1)*DOFs + 1) = 0.0_dp 309 Advance((Perm(i)-1)*DOFs + 2) = 0.0_dp 558 Advance((Perm(j)-1)*DOFs + 1) = & 560 Advance((Perm(j)-1)*DOFs + 2) = & 683 Advance((Perm(k)-1)*DOFs + 1) = Advance((Perm(k)-1)*DOFs + 1) + Displace(1) 684 Advance((Perm(k)-1)*DOFs + 2) = Advance((Perm(k)-1)*DOFs + 2) + Displace(2) 808 Advance((Perm(j)-1)*DOFs + 1) = Advance((Perm(j)-1)*DOFs + 1) + Displace(1) 809 Advance((Perm(j)-1)*DOFs + 2) = Advance((Perm(j)-1)*DOFs + 2) + Displace(2) 1039 Advance((Perm(k)-1)*DOFs + 1) = Advance((Perm(k)-1)*DOFs + 1) + Displace(1) [all …]
|
H A D | CalvingGlacierAdvance3D.F90 | 52 INTEGER :: i, j, k, m, n, DOFs, TotalNodes,& local 93 DOFs = Var % DOFs 94 IF(Var % DOFs /= 3) CALL Fatal(SolverName, "Variable should have 3 DOFs...") 171 NodeNormal(1) = NormalVar % Values(((NormalVar % Perm(i)-1)*NormalVar % DOFs) + 1) 172 NodeNormal(2) = NormalVar % Values(((NormalVar % Perm(i)-1)*NormalVar % DOFs) + 2) 181 NodeVelo(1) = VeloVar % Values(((VeloVar % Perm(i)-1)*VeloVar % DOFs) + 1) 182 NodeVelo(2) = VeloVar % Values(((VeloVar % Perm(i)-1)*VeloVar % DOFs) + 2) 183 NodeVelo(3) = VeloVar % Values(((VeloVar % Perm(i)-1)*VeloVar % DOFs) + 3) 216 Advance((Perm(i)-1)*DOFs + 1) = Displace(1) 217 Advance((Perm(i)-1)*DOFs + 2) = Displace(2) [all …]
|
H A D | Calving3D.F90 | 168 IF(CalvingVar % DOFs /= 3) & 173 DOFs = CalvingVar % DOFs 1800 MyOrderedCalvingValues((county*DOFs)-2) = CalvingValues((CalvingPerm(i)*DOFs)-2) 1802 MyOrderedCalvingValues(county*DOFs) = CalvingValues(CalvingPerm(i)*DOFs) 1837 AllCalvingValues((county*DOFs)-2) = WorkReal((i*DOFs)-2) 1838 AllCalvingValues((county*DOFs)-1) = WorkReal((i*DOFs)-1) 1839 AllCalvingValues((county*DOFs)) = WorkReal((i*DOFs)) 1852 k = i * DOFs 1865 k = i * DOFs 2058 k = j*DOFs [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/bentonite/ |
H A D | BentoniteSolver.F90 | 106 DOFs = Solver % Variable % DOFs 132 STIFF(DOFs*N,DOFs*N ), & 133 MASS( DOFs*N,DOFs*N ), Xi(N), Temp(N), Eta(N), STAT=istat ) 526 DO i=1,DOFs 528 DO j=1,DOFs 545 DO i=1,DOFs 546 k = DOFs*(p-1) + i 644 DO i=1,DOFs 646 DO j=1,DOFs 659 DO i=1,DOFs [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/ResultOutputSolve/ |
H A D | VtkOutputSolver.F90 | 100 DO j=1,Var % DOFs 104 IF ( Var % DOFs < 3 ) THEN 118 DO j=1,Var1 % DOFs 122 IF ( Var1 % DOFs < 3 ) THEN 146 Var % Values(Var % DOFs*(k-1)+Var % DOFs) 159 IF ( Var % DOFs == 1 ) THEN 193 n = Var % DOFs - (MaxDOF - 3) 330 l = MeshUpdate % DOFs * (l-1) 603 OrigDOFs = Var % DOFs 606 Var % DOFs = Var % DOFs*2 [all …]
|
H A D | DXOutputSolver.F90 | 66 IF ( SIZE( Var % Values ) == Var % DOFs ) THEN 83 Model % NumberOfNodes, Var % DOFs, 0, & 97 Var % DOFs-1, nTime, MasterUnit, Prefix ) 102 Var % DOFs, 0, nTime, MasterUnit, Prefix ) 107 Var % DOFs, 0, nTime, MasterUnit, Prefix ) 114 Var % DOFs, 0, nTime, MasterUnit, Prefix ) 266 DO j=1, Var % DOFs 268 Var % Values(Var % DOFs*(k-1)+j) 281 DO j=1,Var1 % DOFs 392 k = Displacement % DOFs * (k-1) [all …]
|
/dports/math/alberta/alberta3-920315ae1bbfd1b1fb6672d916619ac37a411e95/doc/ |
H A D | dof-concept.tex | 3 \idx{DOFs!concepts|(} 7 Degrees of freedom\idx{DOFS@{degree of freedom (DOFs)}} (DOFs) connect 20 For the definition of general finite element spaces DOFs located at 22 the interior of elements are needed. DOFs at a vertex are shared by 23 all elements which meet at this vertex, DOFs at an edge or face are 24 shared by all elements which contain this edge or face, and DOFs 67 hierarchical type basis can share the same abstract set of DOFs. 71 freedom. On the other hand, an application uses DOFs only in 73 while the administration of DOFs is handled by the mesh, definition 74 and access to DOFs is mainly done via finite element spaces. [all …]
|
H A D | rc-impl.tex | 96 \item allocation of new DOFs; 162 \idx{DOFs!adding and removing of DOFs|(} 250 3d) the following DOFs are passed on 279 \caption[Mesh-refinement, maintenance of DOFs]{Edge DOFs that are 287 DOFs are adjusted for 288 \idx{refinement!DOFs!newly created} 310 In 3d, there may be also DOFs at faces; the face DOFs in the boundary of the 312 for DOFs at faces): 363 DOFs are removed on the parent for all patch elements (some DOFs are 375 are removed and in 3d the DOFs in the refinement edge and the DOFs in the [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/AdjointSSA/ |
H A D | AdjointSSA_CostTaubSolver.F90 | 101 INTEGER :: DOFs local 217 DOFs = VVar % DOFs 223 IF (VbVar%DOFs.NE.DOFs) CALL FATAL(SolverName,'Dimension error') 239 IF (DOFs.EQ.2 ) THEN 263 Do j=1,DOFs 264 Vnode(i,j)=VValues(VVar%DOFs*(VPerm(NodeIndexes(i))-1)+j) 292 DO j=1,DOFs 318 DO j=1,DOFs 320 VbValues(VbVar%DOFs*(VbPerm(NodeIndexes(i))-1)+j)=& 321 VbValues(VbVar%DOFs*(VbPerm(NodeIndexes(i))-1)+j)+Vnodeb(i,j)
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/examples/Test_SurfaceBoundaryEnth/ |
H A D | SteadyEnthalpyField.sif | 177 Variable DOFs = 1 194 Variable DOFs = 1 198 Exported Variable 1 DOFs = 1 201 Exported Variable 2 DOFs = 1 204 Exported Variable 3 DOFs = 1 207 Exported Variable 4 DOFs = 1 210 Exported Variable 5 DOFs = 1 213 Exported Variable 6 DOFs = 1 216 Exported Variable 7 DOFs = 1 219 Exported Variable 8 DOFs = 1 [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/MagnetoDynamics/ |
H A D | CalcFields.F90 | 731 IF ( ASSOCIATED(MFD) ) DOFs=DOFs+3 732 IF ( ASSOCIATED(MFS) ) DOFs=DOFs+3 733 IF ( ASSOCIATED(VP) ) DOFs=DOFs+3 734 IF ( ASSOCIATED(CD) ) DOFs=DOFs+3 735 IF ( ASSOCIATED(FWP) ) DOFs=DOFs+1 736 IF ( ASSOCIATED(EF) ) DOFs=DOFs+3 737 IF ( ASSOCIATED(JXB) ) DOFs=DOFs+3 738 IF ( ASSOCIATED(MST) ) DOFs=DOFs+6 740 DOFs = DOFs*vDOFs 742 IF ( ASSOCIATED(ML) ) DOFs=DOFs+1 [all …]
|