Home
last modified time | relevance | path

Searched refs:StiffMatrix (Results 1 – 25 of 74) sorted by relevance

123

/dports/science/elmerfem/elmerfem-release-9.0/fem/src/
H A DTimeIntegrate.F9067 NB1 = SIZE( StiffMatrix,1 )
68 NB2 = SIZE( StiffMatrix,2 )
103 NB1 = SIZE( StiffMatrix,1 )
104 NB2 = SIZE( StiffMatrix,2 )
114 StiffMatrix(i,j) = Beta * StiffMatrix(i,j) + (1/dt)*MassMatrix(i,j)
289 StiffMatrix(i,j) = (1.0d0/dt)*MassMatrix(i,j) + StiffMatrix(i,j)
303 StiffMatrix(i,j) = (1.5d0/dt)*MassMatrix(i,j) + StiffMatrix(i,j)
419 StiffMatrix(i,j) = a(1)*MassMatrix(i,j) + StiffMatrix(i,j)
464 StiffMatrix(i,j) = StiffMatrix(i,j) + &
561 StiffMatrix(i,j) = StiffMatrix(i,j) / 3 + &
[all …]
H A DStressGeneral.F9056 SUBROUTINE StressGeneralCompose( MassMatrix, StiffMatrix,ForceVector, & argument
60 REAL(KIND=dp) :: MassMatrix(:,:), StiffMatrix(:,:),NodalHeatExpansion(:,:,:) local
115 StiffMatrix = 0.0D0
253 StiffMatrix( 2*(p-1)+i,2*(q-1)+j ) = &
254 StiffMatrix( 2*(p-1)+i,2*(q-1)+j ) + s * A(i,j)
263 StiffMatrix( dim*(p-1)+i,dim*(q-1)+j ) = &
264 StiffMatrix( dim*(p-1)+i,dim*(q-1)+j ) + s*A(i,j)
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/
H A DStatMagSolve.F90587 StiffMatrix = 0.0D0
649 StiffMatrix(p,q) = StiffMatrix(p,q) + mat
715 StiffMatrix = 0.0D0
777 StiffMatrix(2*p-1,2*q-1) = StiffMatrix(2*p-1,2*q-1) + a11
778 StiffMatrix(2*p,2*q) = StiffMatrix(2*p,2*q) + a22
780 StiffMatrix(2*p,2*q-1) = StiffMatrix(2*p,2*q-1) + a21
781 StiffMatrix(2*p-1,2*q) = StiffMatrix(2*p-1,2*q) + a12
952 StiffMatrix = 0.0D0
1006 StiffMatrix(3*p-2,3*q-2) = StiffMatrix(3*p-2,3*q-2) + mat
1007 StiffMatrix(3*p-1,3*q-1) = StiffMatrix(3*p-1,3*q-1) + mat
[all …]
H A DRateOfChange.F9074 TYPE(Matrix_t), POINTER :: StiffMatrix local
111 StiffMatrix => Solver % Matrix
112 ForceVector => StiffMatrix % RHS
203 CALL InitializeToZero( StiffMatrix, ForceVector )
290 CALL UpdateGlobalEquations( StiffMatrix, LocalStiffMatrix, &
319 CALL SolveSystem( StiffMatrix, ParMatrix, ForceVector, &
334 SUBROUTINE LocalMatrix( StiffMatrix, Force, Element, n, &
337 REAL(KIND=dp) :: StiffMatrix(:,:), Force(:), LocalU0(:), & local
356 StiffMatrix = 0.0d0
380 StiffMatrix( (p-1)*VariableComponents+i, &
[all …]
H A DTransportEquation.F9072 TYPE(Matrix_t), POINTER :: StiffMatrix local
117 StiffMatrix => Solver % Matrix
118 ForceVector => StiffMatrix % RHS
487 REAL(KIND=dp) :: StiffMatrix(:,:), Force(:), & local
504 StiffMatrix = 0.0d0
574 StiffMatrix = 0.0d0
701 StiffMatrix(:,:),Force(:),dt local
740 StiffMatrix, Force, X, V, A )
748 StiffMatrix, Force, X, V, A )
770 s = s - StiffMatrix(i,j) * X(j)
[all …]
H A DMovingElstatSolver.F9067 TYPE(Matrix_t), POINTER :: StiffMatrix local
116 StiffMatrix => Solver % Matrix
117 ForceVector => StiffMatrix % RHS
155 MatrixValues(SIZE(StiffMatrix % Values)), &
291 CALL InitializeToZero( StiffMatrix, ForceVector )
316 MatrixValues = StiffMatrix % Values
482 StiffMatrix % Values = MatrixValues
612 StiffMatrix % Values = MatrixValues
881 REAL(KIND=dp) :: StiffMatrix(:,:) local
894 StiffMatrix = 0.0d0
[all …]
H A DFDiffusion.F9071 TYPE(Matrix_t),POINTER :: StiffMatrix local
120 StiffMatrix => Solver % Matrix
121 ForceVector => StiffMatrix % RHS
304 CALL UpdateGlobalEquations( StiffMatrix, LocalStiffMatrix, &
424 SUBROUTINE LocalMatrix( StiffMatrix, Force, AngularFrequency, &
427 REAL(KIND=dp) :: StiffMatrix(:,:), Force(:), AngularFrequency, & local
527 StiffMatrix( 2*(i-1)+1, 2*(j-1)+1 ) = REAL( LSTIFF(i,j) )
528 StiffMatrix( 2*(i-1)+1, 2*(j-1)+2 ) = -AIMAG( LSTIFF(i,j) )
530 StiffMatrix( 2*(i-1)+2, 2*(j-1)+2 ) = REAL( LSTIFF(i,j) )
611 StiffMatrix( 2*(i-1)+1, 2*(j-1)+1 ) = REAL( LSTIFF(i,j) )
[all …]
H A DTransientPhaseChange.F9068 TYPE(Matrix_t),POINTER :: StiffMatrix local
370 StiffMatrix => Solver % Matrix
536 SUBROUTINE VelocityLocalMatrix( StiffMatrix, MassMatrix, ForceVector,&
540 REAL(KIND=dp) :: StiffMatrix(:,:), MassMatrix(:,:), ForceVector(:) local
559 StiffMatrix = 0.0d0
688 StiffMatrix(p,q) = StiffMatrix(p,q) + s * Basis(p) * Basis(q) * &
696 StiffMatrix(p,q) = StiffMatrix(p,q) + &
718 SUBROUTINE SurfaceLocalMatrix( StiffMatrix, MassMatrix, ForceVector,&
722 REAL(KIND=dp) :: StiffMatrix(:,:), MassMatrix(:,:), ForceVector(:) local
738 StiffMatrix = 0.0d0
[all …]
H A DArteryOutlet.F90110 TYPE(Matrix_t),POINTER :: StiffMatrix local
169 StiffMatrix => Solver % Matrix
170 IF(.NOT. ASSOCIATED(StiffMatrix)) THEN
173 ForceVector => StiffMatrix % RHS
406 CALL InitializeToZero( StiffMatrix, ForceVector )
461 CALL InitializeToZero( StiffMatrix, ForceVector )
569 ELSE IF ( StiffMatrix % FORMAT == MATRIX_CRS .AND. StiffMatrix % Symmetric ) THEN
573 CALL ZeroRow( StiffMatrix, k )
658 SUBROUTINE LocalMatrix( StiffMatrix,MassMatrix,&
677 StiffMatrix = 0.0d0
[all …]
H A DFDiffusion3D.F9071 TYPE(Matrix_t),POINTER :: StiffMatrix local
120 StiffMatrix => Solver % Matrix
121 ForceVector => StiffMatrix % RHS
320 CALL UpdateGlobalEquations( StiffMatrix, LocalStiffMatrix, &
401 SUBROUTINE LocalMatrix( StiffMatrix, Force, AngularFrequency, &
404 REAL(KIND=dp) :: StiffMatrix(:,:), Force(:), AngularFrequency, & local
539 StiffMatrix( 2*dim*(p-1)+2*i, 2*dim*(q-1)+2*j ) = &
625 StiffMatrix( 2*(i-1)+1, 2*(j-1)+1 ) = REAL( LSTIFF(i,j) )
626 StiffMatrix( 2*(i-1)+1, 2*(j-1)+2 ) = -AIMAG( LSTIFF(i,j) )
627 StiffMatrix( 2*(i-1)+2, 2*(j-1)+1 ) = AIMAG( LSTIFF(i,j) )
[all …]
H A DElasticSolve.F90410 StiffMatrix => Solver % Matrix
660 ALLOCATE(StiffMatrix % BulkRHS(SIZE(StiffMatrix % RHS)))
909 StiffMatrix % RHS => StiffMatrix % BulkRHS
1120 StiffMatrix % RHS => StiffMatrix % BulkRHS
1155 StiffMatrix % DValues(i) = StiffMatrix % DValues(i) - Displacement(i)
1178 CALL ParallelUpdateRHS(StiffMatrix, StiffMatrix % BulkRHS)
1518 StiffMatrix = 0.0D0
2139 StiffMatrix = 0.0D0
2529 StiffMatrix = 0.0D0
2766 StiffMatrix(DOFs*p,DOFs*(q-1)+i) = StiffMatrix(DOFs*p,DOFs*(q-1)+i) + &
[all …]
H A DAcoustics.F90207 StiffMatrix => Solver % Matrix
208 ForceVector => StiffMatrix % RHS
4507 StiffMatrix = 0.0d0
4607 StiffMatrix = 0.0d0
4680 StiffMatrix = 0.0d0
4750 StiffMatrix = 0.0d0
4867 StiffMatrix = 0.0d0
4953 StiffMatrix = 0.0d0
5038 StiffMatrix = 0.0d0
5142 StiffMatrix = 0.0d0
[all …]
H A DPhaseChangeSolve.F9093 TYPE(Matrix_t),POINTER :: StiffMatrix local
353 StiffMatrix => Solver % Matrix
368 CALL InitializeToZero( StiffMatrix, ForceVector )
393 CALL UpdateGlobalEquations( StiffMatrix, LocalStiffMatrix, &
451 CALL InitializeToZero( StiffMatrix, ForceVector )
955 StiffMatrix = 0.0d0
1077 StiffMatrix(p,q) = StiffMatrix(p,q) + s * Basis(p) * Basis(q) * &
1080 StiffMatrix(p,q) = StiffMatrix(p,q) + &
1128 StiffMatrix = 0.0d0
1177 StiffMatrix(p,q) = StiffMatrix(p,q) + &
[all …]
H A DStatCurrentSolve.F90100 TYPE(Matrix_t), POINTER :: StiffMatrix local
151 StiffMatrix => Solver % Matrix
152 ForceVector => StiffMatrix % RHS
251 IF ( .NOT.ASSOCIATED( StiffMatrix % MassValues ) ) THEN
252 ALLOCATE( StiffMatrix % Massvalues( LocalNodes ) )
253 StiffMatrix % MassValues = 0.0d0
815 SUBROUTINE StatCurrentCompose( StiffMatrix,Force,Conductivity, &
818 REAL(KIND=dp) :: StiffMatrix(:,:),Force(:),Load(:), Conductivity(:,:,:) local
837 StiffMatrix = 0.0d0
895 StiffMatrix(p,q) = StiffMatrix(p,q) + S*A
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/
H A DAdjointSolver.F9069 TYPE(Matrix_t),POINTER :: InitMat,TransMat,StiffMatrix local
104 StiffMatrix => Solver % Matrix
105 ForceVector => StiffMatrix % RHS
167 IF ( SIZE(StiffMatrix % Values) .NE. SIZE(TransMat % Values) ) THEN
171 StiffMatrix % Values = TransMat % Values
172 StiffMatrix % Rows = TransMat % Rows
173 StiffMatrix % Cols = TransMat % Cols
174 IF(ASSOCIATED(TransMat % Diag)) StiffMatrix % Diag = TransMat % Diag
H A DForceToStress.F9072 TYPE(Matrix_t), POINTER :: StiffMatrix local
153 StiffMatrix => Solver % Matrix
154 ForceVector => StiffMatrix % RHS
309 SUBROUTINE LocalMatrix( MassMatrix, StiffMatrix, &
314 REAL(KIND=dp) :: StiffMatrix(:,:), MassMatrix(:,:) local
339 StiffMatrix = 0.0D0
375 StiffMatrix(p,q) = &
376 StiffMatrix(p,q) + s*Basis(q)*Basis(p)
H A DComputeStrainRate.F9074 TYPE(Matrix_t),POINTER :: StiffMatrix local
161 StiffMatrix => Solver % Matrix
162 ForceVector => StiffMatrix % RHS
322 SUBROUTINE LocalMatrix(COMP, MassMatrix, StiffMatrix, ForceVector, &
327 REAL(KIND=dp) :: StiffMatrix(:,:), MassMatrix(:,:) local
356 StiffMatrix = 0.0D0
418 StiffMatrix(p,q) = &
419 StiffMatrix(p,q) + s*Basis(q)*Basis(p)
H A DTwoMeshes.F90105 TYPE(Matrix_t), POINTER :: StiffMatrix local
627 StiffMatrix => Solver % Matrix
628 ForceVector => StiffMatrix % RHS
633 CALL SetDirichtletPoint( StiffMatrix, ForceVector,1,1, &
638 CALL SetDirichtletPoint( StiffMatrix, ForceVector,1,1, &
643 CALL SetDirichtletPoint( StiffMatrix, ForceVector,1,1, &
820 TYPE(Matrix_t), POINTER :: StiffMatrix local
837 ELSE IF ( StiffMatrix % FORMAT == MATRIX_CRS .AND. &
838 StiffMatrix % Symmetric ) THEN
841 s = StiffMatrix % Values(StiffMatrix % Diag(PermIndex))
[all …]
H A DComputeDevStressNS.F9075 TYPE(Matrix_t),POINTER :: StiffMatrix local
182 StiffMatrix => Solver % Matrix
183 ForceVector => StiffMatrix % RHS
372 SUBROUTINE LocalNSMatrix(COMP, MassMatrix, StiffMatrix, ForceVector, &
379 REAL(KIND=dp) :: StiffMatrix(:,:), MassMatrix(:,:) local
418 StiffMatrix = 0.0_dp
493 StiffMatrix(p,q) = &
494 StiffMatrix(p,q) + s*Basis(q)*Basis(p)
H A DExportVertically.F9076 TYPE(Matrix_t),POINTER :: StiffMatrix local
119 StiffMatrix => Solver % Matrix
120 ForceVector => StiffMatrix % RHS
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/Adjoint/
H A DAdjoint_LinearSolver.F9077 TYPE(Matrix_t),POINTER :: InitMat,TransMat,StiffMatrix local
117 StiffMatrix => Solver % Matrix
118 ForceVector => StiffMatrix % RHS
230 StiffMatrix % Values = TransMat % Values
231 StiffMatrix % Rows = TransMat % Rows
232 StiffMatrix % Cols = TransMat % Cols
233 IF(ASSOCIATED(TransMat % Diag)) StiffMatrix % Diag = TransMat % Diag
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/AdjointSSA/
H A DAdjointSSA_AdjointSolver.F9074 TYPE(Matrix_t),POINTER :: InitMat,TransMat,StiffMatrix local
118 StiffMatrix => Solver % Matrix
119 ForceVector => StiffMatrix % RHS
235 StiffMatrix % Values = TransMat % Values
236 StiffMatrix % Rows = TransMat % Rows
237 StiffMatrix % Cols = TransMat % Cols
238 IF(ASSOCIATED(TransMat % Diag)) StiffMatrix % Diag = TransMat % Diag
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/LevelSet/
H A DLevelSetCurvature.F9056 TYPE(Matrix_t),POINTER :: StiffMatrix local
93 StiffMatrix => Solver % Matrix
94 ForceVector => StiffMatrix % RHS
239 CALL SolveSystem( StiffMatrix, ParMatrix, ForceVector, &
277 SUBROUTINE LocalMatrix( StiffMatrix,ForceVector, Surf, &
282 REAL(KIND=dp), DIMENSION(:,:) :: StiffMatrix local
306 StiffMatrix = 0.0D0
361 StiffMatrix(p,q) = StiffMatrix(p,q) + s * A
H A DLevelSetSolver.F9061 TYPE(Matrix_t),POINTER :: StiffMatrix local
101 StiffMatrix => Solver % Matrix
102 ForceVector => StiffMatrix % RHS
292 SUBROUTINE LocalMatrix( MassMatrix,StiffMatrix,ForceVector, &
297 REAL(KIND=dp), DIMENSION(:,:) :: MassMatrix,StiffMatrix,NodalVelo local
329 StiffMatrix = 0.0D0
474 StiffMatrix(p,q) = StiffMatrix(p,q) + s * A
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/MeshAdaptation_2D/
H A DMMG2D_MetricAniso.F9093 TYPE(Matrix_t),POINTER :: StiffMatrix local
147 StiffMatrix => Solver % Matrix
148 ForceVector => StiffMatrix % RHS
288 SUBROUTINE LocalMatrix(COMP, MassMatrix, StiffMatrix, ForceVector, &
291 REAL(KIND=dp) :: StiffMatrix(:,:), MassMatrix(:,:) local
314 StiffMatrix = 0.0D0
352 StiffMatrix(p,q) = &
353 StiffMatrix(p,q) + s*Basis(q)*Basis(p) + &

123