/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/solver/las/ |
H A D | hecmw_solver_scaling_nn.f90 | 28 integer(kind=kint) :: i,j,k,ii,ij, ip(hecMAT%NDOF),iq(hecMAT%NDOF) 36 NDOF = hecMAT%NDOF 37 NDOF2 = NDOF*NDOF 50 do k=1, NDOF 51 scale (NDOF*(i-1)+k)= 1.d0/dsqrt(dabs(D(NDOF*NDOF*(i-1)+(k-1)*(NDOF+1)+1))) 103 B(NDOF*(i-1)+k)=B(NDOF*(i-1)+k)*scale(NDOF*(i-1)+k) 122 NDOF = hecMAT%NDOF 123 NDOF2 = NDOF*NDOF 136 do k=1,NDOF 137 X(NDOF*(i-1)+k)=X(NDOF*(i-1)+k)*scale(NDOF*(i-1)+k) [all …]
|
H A D | hecmw_solver_las_nn.f90 | 176 NDOF = hecMAT%NDOF 177 NDOF2 = NDOF*NDOF 512 NDOF = hecMAT%NDOF 513 allocate(W(NDOF*hecMAT%NP,NDOF)) 518 W(NDOF*(ip-1)+i,j) = D(NDOF*NDOF*(ip-1)+(i-1)*NDOF+j) 531 D(NDOF*NDOF*(ip-1)+(i-1)*NDOF+j) = W(NDOF*(ip-1)+i,j) 544 do i = 1, hecMAT1%NP*hecMAT1%NDOF*hecMAT1%NDOF 548 do i = 1, hecMAT1%NPU*hecMAT1%NDOF*hecMAT1%NDOF 564 do i = 1, hecMAT%NP*hecMAT%NDOF*hecMAT%NDOF 568 do i = 1, hecMAT%NPU*hecMAT%NDOF*hecMAT%NDOF [all …]
|
H A D | hecmw_jadm_nn.f90 | 34 allocate(AJAD((hecMAT%NPL+hecMAT%NPU)*hecMAT%NDOF*hecMAT%NDOF)) 72 NDOF = hecMAT%NDOF 73 NDOF2 = NDOF*NDOF 79 do idof=1,hecMAT%NDOF 81 … Y(NDOF*(i-1)+idof) = Y(NDOF*(i-1)+idof) + D(NDOF2*(i-1)+NDOF*(idof-1)+jdof)*X(NDOF*(i-1)+jdof) 104 NDOF = hecMAT%NDOF;NDOF2=NDOF*NDOF 175 NDOF2=NDOF*NDOF 184 do idof = 1, NDOF 185 do jdof = 1, NDOF 197 do idof = 1, NDOF [all …]
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/solver/precond/nn/ |
H A D | hecmw_precond_DIAG_nn.f90 | 35 real (kind=kreal):: ALUtmp(hecMAT%NDOF,hecMAT%NDOF), PW(hecMAT%NDOF) 44 NDOF = hecMAT%NDOF 45 NDOF2 = NDOF*NDOF 75 do i = 1, NDOF 81 do k= 1, NDOF 118 do j=1,NDOF 121 do j=2,NDOF 123 X(j)=X(j)-ALU(NDOF*NDOF*(i-1)+NDOF*(j-1)+k )*X(k) 128 X(j)=X(j)-ALU(NDOF*NDOF*(i-1)+NDOF*(j-1)+k )*X(k) 130 X(j)=ALU(NDOF*NDOF*(i-1)+(NDOF+1)*(j-1)+1 )*X(j) [all …]
|
H A D | hecmw_precond_RIF_nn.f90 | 55 NDOF = hecMAT%NDOF 56 NDOF2 = NDOF*NDOF 95 real(kind=kreal) :: SW(NDOF),X(NDOF) 96 NDOF2=NDOF*NDOF 121 ZP(NDOF*(i-1)+1:NDOF*(i-1)+NDOF) = X(1:NDOF) 126 ZP(NDOF*(i-1)+idof)= ZP(NDOF*(i-1)+idof)*RIFD(NDOF2*(i-1)+(idof-1)*NDOF+idof) 150 X(idof) = X(idof) - RIFD(NDOF*NDOF*(i-1)+NDOF*(jdof-1)+idof)*X(jdof) 174 NDOF = hecMAT%NDOF 175 NDOF2= NDOF*NDOF 333 NDOF=hecMAT%NDOF [all …]
|
H A D | hecmw_precond_SAINV_nn.f90 | 54 NDOF = hecmat%NDOF 55 NDOF2 = NDOF*NDOF 93 real(kind=kreal) :: SW(NDOF),X(NDOF) 119 T(NDOF*(i-1)+idof)=T(NDOF*(i-1)+idof)+SAINVD(NDOF2*(i-1)+NDOF*(jdof-1)+idof)*X(jdof) 121 T(NDOF*(i-1)+idof)=T(NDOF*(i-1)+idof)*SAINVD(NDOF2*(i-1)+NDOF*(idof-1)+idof) 151 ZP(NDOF*(i-1)+idof) = ZP(NDOF*(i-1)+idof)+SAINVD(NDOF2*(i-1)+NDOF*(idof-1)+jdof)*X(jdof) 167 real(kind=krealp) :: dd, dtmp(hecMAT%NDOF), X(hecMAT%NDOF) 201 vv(NDOF*(i-1)+idof) = vv(NDOF*(i-1)+idof) + D(NDOF2*(i-1)+NDOF*(idof-1)+jdof)*X(jdof) 211 … vv(NDOF*(in-1)+idof) = vv(NDOF*(in-1)+idof) + AL(NDOF2*(j-1)+NDOF*(jdof-1)+idof)*X(jdof) 221 … vv(NDOF*(in-1)+idof) = vv(NDOF*(in-1)+idof) + AU(NDOF2*(j-1)+NDOF*(jdof-1)+idof)*X(jdof) [all …]
|
H A D | hecmw_precond_SSOR_nn.f90 | 62 real (kind=kreal) :: ALUtmp(hecMAT%NDOF,hecMAT%NDOF), PW(hecMAT%NDOF) 84 NDOF=hecMAT%NDOF 85 NDOF2=NDOF*NDOF 168 do i = 1, NDOF 179 do i = 1, NDOF 185 do k= 1, NDOF 197 do i = 1, NDOF 221 real(kind=kreal) :: SW(NDOF), X(NDOF) 233 NDOF2=NDOF*NDOF 343 ZP(NDOF*(iold-1)+1:NDOF*(iold-1)+NDOF) = X(1:NDOF) [all …]
|
H A D | hecmw_precond_BILU_nn.f90 | 58 NDOF = hecMAT%NDOF 119 X(ii)=X(ii)-Dlu0(NDOF*NDOF*(i-1)+NDOF*(ii-1)+ij )*X(ij) 124 X(ii)=X(ii)-Dlu0(NDOF*NDOF*(i-1)+NDOF*(ii-1)+ij )*X(ij) 126 X(ii)=Dlu0(NDOF*NDOF*(i-1)+(NDOF+1)*(ii-1)+1 )*X(ii) 155 X(ii)=X(ii)-Dlu0(NDOF*NDOF*(i-1)+NDOF*(ii-1)+ij )*X(ij) 162 X(ii)=Dlu0(NDOF*NDOF*(i-1)+(NDOF+1)*(ii-1)+1 )*X(ii) 215 NDOF2=NDOF*NDOF 369 NDOF2=NDOF*NDOF 750 NDOF2=NDOF*NDOF 1503 real(kind=kreal) :: ALU(NDOF,NDOF), D(NDOF*NDOF), PW(NDOF) [all …]
|
/dports/cad/tochnog/Sources/Tochnog-Latest-jan-2014/ |
H A D | beam.cc | 240 old_beam_dof[NNOL*NDOF], new_beam_dof[NNOL*NDOF], in beam_2d() 245 old_rotation_matrix[NNOL*NDOF*NNOL*NDOF], new_rotation_matrix[NNOL*NDOF*NNOL*NDOF], in beam_2d() 246 local_beam_matrix[NNOL*NDOF*NNOL*NDOF], new_beam_matrix[NNOL*NDOF*NNOL*NDOF], in beam_2d() 247 work[NNOL*NDOF*NNOL*NDOF]; in beam_2d() 428 new_beam_matrix, work, NNOL*NDOF, NNOL*NDOF ); in beam_2d() 430 NNOL*NDOF, NNOL*NDOF ); in beam_2d() 440 backrotated_beam_moment, NNOL*NDOF, NNOL*NDOF, 1 ); in beam_2d() 471 new_beam_matrix[(inol*NDOF+idof)*NNOL*NDOF+jnol*NDOF+jdof]*dtime; in beam_2d() 481 new_beam_matrix[(inol*NDOF+idof)*NNOL*NDOF+jnol*NDOF+jdof]; in beam_2d() 498 new_beam_matrix[(inol*NDOF+idof)*NNOL*NDOF+jnol*NDOF+jdof]*dtime; in beam_2d() [all …]
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/analysis/dynamic/transit/ |
H A D | dynamic_mat_ass_bc_ac.f90 | 46 NDOF = hecMAT%NDOF 78 + b2*fstrDYNAMIC%VEL (NDOF*in-(NDOF-idof),1) & 84 + b2*fstrDYNAMIC%VEL (NDOF*in-(NDOF-idof),1) & 85 + b3*fstrDYNAMIC%ACC (NDOF*in-(NDOF-idof),1) & 137 - fstrDYNAMIC%DISP(NDOF*in-(NDOF-idof),3) & 139 hecMAT%B (NDOF*in-(NDOF-idof)) = RHS 140 fstrDYNAMIC%VEC1(NDOF*in-(NDOF-idof)) = 1.0d0 180 NDOF = hecMAT%NDOF 200 fstrDYNAMIC%ACC (NDOF*in-(NDOF-idof),1) = RHS 236 NDOF = hecMAT%NDOF [all …]
|
H A D | dynamic_mat_ass_bc_vl.f90 | 56 NDOF = hecMAT%NDOF 93 RHS = fstrDYNAMIC%DISP(NDOF*in-(NDOF-idof),1) & 94 + b2*fstrDYNAMIC%VEL (NDOF*in-(NDOF-idof),1) & 144 RHS = fstrDYNAMIC%DISP(NDOF*in-(NDOF-idof),3) & 146 hecMAT%B (NDOF*in-(NDOF-idof)) = RHS 147 fstrDYNAMIC%VEC1(NDOF*in-(NDOF-idof)) = 1.0d0 184 NDOF = hecMAT%NDOF 205 fstrDYNAMIC%VEL (NDOF*in-(NDOF-idof),1) = RHS 252 NDOF = hecMAT%NDOF 273 RHS = fstrDYNAMIC%DISP(NDOF*in-(NDOF-idof),3) & [all …]
|
H A D | dynamic_mat_ass_bc.f90 | 42 ndof = hecMAT%NDOF 180 NDOF = hecMAT%NDOF 199 hecMAT%B (NDOF*in-(NDOF-idof)) = RHS 200 fstrDYNAMIC%VEC1(NDOF*in-(NDOF-idof)) = 1.0d0 203 fstrSOLID%REACTION(NDOF*(in-1)+idof) = fstrSOLID%QFORCE(NDOF*(in-1)+idof) 235 NDOF = hecMAT%NDOF 256 fstrDYNAMIC%DISP(NDOF*in-(NDOF-idof),1) = RHS 291 ndof = hecMAT%NDOF 299 NDOF = hecMAT%NDOF 318 hecMAT%B(NDOF*in-(NDOF-idof)) = RHS*fstrDYNAMIC%VEC1(NDOF*in-(NDOF-idof)) [all …]
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/solver/matrix/ |
H A D | hecmw_mat_ass.f90 | 41 ndof = hecMAT%NDOF 80 NDOF = hecMAT%NDOF 285 NDOF = hecMAT%NDOF 298 idx = NDOF**2 * (k-1) + NDOF * (idof-1) + jdof 314 idx = NDOF**2 * (k-1) + NDOF * (idof-1) + jdof 319 idx = NDOF**2 * (inod - 1) + NDOF * (idof - 1) + jdof 338 NDOF = hecMAT%NDOF 343 hecMAT%B(NDOF*inode-(NDOF-idof)) = RHS 344 if(present(conMAT)) conMAT%B(NDOF*inode-(NDOF-idof)) = 0.0D0 345 ndof2 = NDOF*NDOF [all …]
|
H A D | hecmw_matrix_reorder.f90 | 91 call reorder_off_diag2(N, NDOF, iperm, & 111 j0new = (inew-1)*NDOF 112 j0old = (iold-1)*NDOF 113 do j=1,NDOF 133 j0new = (inew-1)*NDOF 134 j0old = (iold-1)*NDOF 135 do j=1,NDOF 167 NDOF2 = NDOF*NDOF 186 NDOF2 = NDOF*NDOF 218 NDOF2 = NDOF*NDOF [all …]
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/analysis/dynamic/mode/ |
H A D | fstr_EIG_output.f90 | 26 NDOF = hecMESH%n_dof 27 NNDOF = N *NDOF 28 NPNDOF = NP*NDOF 40 allocate(r(NDOF)) 48 do k = 1, NDOF 49 in = NDOF*(j-1) + k 59 do j = 1, NDOF 60 in = NDOF*(i-1) + j 132 NDOF = hecMAT%NDOF 215 integer(kind=kint) :: NDOF local [all …]
|
H A D | fstr_EIG_lanczos_util.f90 | 17 integer(kind=kint) :: N, NP, NDOF, NDOF2, NNDOF, NPNDOF local 26 NDOF = hecMESH%n_dof 27 NDOF2 = NDOF*NDOF 28 NNDOF = N *NDOF 29 NPNDOF = NP*NDOF 41 do j = 1, NDOF 42 …hecMAT%D(NDOF2*(i-1) + (NDOF+1)*(j-1) + 1) = hecMAT%D(NDOF2*(i-1) + (NDOF+1)*(j-1) + 1) + sigma * …
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/solver/ |
H A D | hecmw_solver.f90 | 27 integer(kind=kint) :: i, myrank, NDOF local 29 NDOF=hecMAT%NDOF 88 subroutine hecmw_substitute_solver(hecMESH, hecMATorig, NDOF) argument 100 integer(kind=kint) NDOF local 101 if (NDOF == hecMATorig%NDOF) then 103 else if (NDOF < hecMATorig%NDOF) then 106 call hecmw_blockmatrix_expand(hecMATorig,hecMAT,NDOF) 127 if (NDOF /= hecMATorig%NDOF) then 128 call hecmw_vector_contract(hecMATorig,hecMAT,NDOF)
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/solver/iterative/ |
H A D | hecmw_solver_GPBiCG.f90 | 37 integer(kind=kint ) :: N, NP, NDOF, NNDOF local 81 NDOF = hecMAT%NDOF 82 NNDOF = N * NDOF 95 allocate (WW(NDOF*NP,14)) 120 call hecmw_InnerProduct_R(hecMESH, NDOF, B, B, BNRM2, Tcomm) 128 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,RT), WW(:,R), RHO, Tcomm) 193 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,RT), WW(:,PT), RHO1, Tcomm) 291 call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,R), WW(:,R), RR(1)) 292 call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,R), WW(:,RT), RR(2)) 318 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,R), WW(:,R), DNRM2, Tcomm) [all …]
|
H A D | hecmw_solver_CG.f90 | 40 integer(kind=kint ) :: N, NP, NDOF, NNDOF local 79 NDOF = hecMAT%NDOF 80 NNDOF = N * NDOF 97 allocate (WW(NDOF*NP, 4)) 126 call hecmw_InnerProduct_R(hecMESH, NDOF, B, B, BNRM2, Tcomm) 195 call hecmw_copy_R(hecMESH, NDOF, WW(:,Z), WW(:,P)) 198 call hecmw_xpay_R(hecMESH, NDOF, BETA, WW(:,Z), WW(:,P)) 213 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,P), WW(:,Q), C1, Tcomm) 232 call hecmw_axpy_R(hecMESH, NDOF, ALPHA, WW(:,P), X) 237 call hecmw_axpy_R(hecMESH, NDOF, -ALPHA, WW(:,Q), WW(:,R)) [all …]
|
H A D | hecmw_solver_BiCGSTAB.f90 | 35 integer(kind=kint ) :: N, NP, NDOF, NNDOF local 76 NDOF = hecMAT%NDOF 77 NNDOF = N * NDOF 93 allocate (WW(NDOF*NP, 8)) 120 call hecmw_InnerProduct_R(hecMESH, NDOF, B, B, BNRM2, Tcomm) 157 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,R), WW(:,RT), RHO, Tcomm) 192 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,RT), WW(:,V), C2, Tcomm) 221 call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,T), WW(:,S), CG(1)) 222 call hecmw_InnerProduct_R_nocomm(hecMESH, NDOF, WW(:,T), WW(:,T), CG(2)) 248 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,R), WW(:,R), DNRM2, Tcomm) [all …]
|
H A D | hecmw_solver_GMRES.f90 | 37 integer(kind=kint ) :: N, NP, NDOF, NNDOF local 79 NDOF = hecMAT%NDOF 80 NNDOF = N * NDOF 92 if (NREST >= NDOF*NP-1) NREST = NDOF*NP-2 98 allocate (WW(NDOF*NP,NRK)) 132 call hecmw_InnerProduct_R(hecMESH, NDOF, B, B, BNRM2, Tcomm) 172 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,R), WW(:,R), DNRM2, Tcomm) 215 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,W), WW(:,V+K-1), val, Tcomm) 223 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,W), WW(:,W), val, Tcomm) 374 call hecmw_InnerProduct_R(hecMESH, NDOF, WW(:,R), WW(:,R), DNRM2, Tcomm) [all …]
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/common/ |
H A D | hecmw_util_f.F90 | 894 nBlock = P%NDOF * P%NDOF 949 nBlock = P%NDOF * P%NDOF 1016 hecMAT%NDOF = hecMATorig%NDOF 1029 NDOF = hecMATorig%NDOF 1030 NDOF2 = NDOF*NDOF 1038 hecMAT%NDOF = NDOF 1063 NDOF2 = NDOF*NDOF 1082 hecMAT%NDOF = NDOF 1087 do j = 1, NDOF 1107 do j = 1, NDOF [all …]
|
/dports/science/InsightToolkit/ITK-5.0.1/Modules/Numerics/FEM/src/ |
H A D | itkFEMElement2DC1Beam.cxx | 140 const unsigned int NDOF = this->GetNumberOfDegreesOfFreedom(); in GetStiffnessMatrix() local 142 MatrixType k(NDOF, NDOF); in GetStiffnessMatrix() 143 MatrixType kb(NDOF, NDOF); in GetStiffnessMatrix() 184 const unsigned int NDOF = this->GetNumberOfDegreesOfFreedom(); in GetMassMatrix() local 185 MatrixType m(NDOF, NDOF, 0.0); in GetMassMatrix() 186 MatrixType mb(NDOF, NDOF, 0.0); in GetMassMatrix() 187 MatrixType k(NDOF, NDOF, 0.0); in GetMassMatrix()
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/analysis/static/ |
H A D | solve_LINEQ_direct_serial_lag.f90 | 41 ntdf = hecMAT%NP*hecMAT%NDOF + fstrMAT%num_lagrange 42 ilag_sta = hecMAT%NP*hecMAT%NDOF + 1 43 numNon0 = hecMAT%NPU*hecMAT%NDOF**2+hecMAT%NP*hecMAT%NDOF*(ntdf+1)/2 & 44 + (fstrMAT%numU_lagrange)*hecMAT%NDOF+fstrMAT%num_lagrange
|
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/common/ |
H A D | fstr_contact.f90 | 378 integer(kind=kint), parameter :: NDOF=3 local 418 force(1:nnode*NDOF+NDOF) = force(1:nnode*NDOF+NDOF)*nrlforce 421 do k=1,NDOF 422 hecMAT%B(NDOF*(nd-1)+k)=hecMAT%B(NDOF*(nd-1)+k)-force((m-1)*NDOF+k) 536 integer(kind=kint), parameter :: NDOF=3 local 546 allocate(coord(NDOF*hecMESH%n_node)) 547 do i=1,NDOF*hecMESH%n_node 551 do i=1,NDOF*hecMESH%n_node 566 idx = NDOF*(ndlocal(i)-1) 567 ecoord(1:NDOF,i) = coord(idx+1:idx+NDOF) [all …]
|