/dports/science/elmerfem/elmerfem-release-9.0/fem/src/ |
H A D | TimeIntegrate.F90 | 67 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 D | StressGeneral.F90 | 56 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 D | StatMagSolve.F90 | 587 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 D | RateOfChange.F90 | 74 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 D | TransportEquation.F90 | 72 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 D | MovingElstatSolver.F90 | 67 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 D | FDiffusion.F90 | 71 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 D | TransientPhaseChange.F90 | 68 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 D | ArteryOutlet.F90 | 110 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 D | FDiffusion3D.F90 | 71 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 D | ElasticSolve.F90 | 410 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 D | Acoustics.F90 | 207 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 D | PhaseChangeSolve.F90 | 93 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 D | StatCurrentSolve.F90 | 100 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 D | AdjointSolver.F90 | 69 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 D | ForceToStress.F90 | 72 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 D | ComputeStrainRate.F90 | 74 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 D | TwoMeshes.F90 | 105 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 D | ComputeDevStressNS.F90 | 75 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 D | ExportVertically.F90 | 76 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 D | Adjoint_LinearSolver.F90 | 77 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 D | AdjointSSA_AdjointSolver.F90 | 74 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 D | LevelSetCurvature.F90 | 56 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 D | LevelSetSolver.F90 | 61 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 D | MMG2D_MetricAniso.F90 | 93 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) + &
|