Home
last modified time | relevance | path

Searched refs:NDOF (Results 1 – 25 of 117) sorted by relevance

12345

/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/solver/las/
H A Dhecmw_solver_scaling_nn.f9028 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 Dhecmw_solver_las_nn.f90176 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 Dhecmw_jadm_nn.f9034 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 Dhecmw_precond_DIAG_nn.f9035 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 Dhecmw_precond_RIF_nn.f9055 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 Dhecmw_precond_SAINV_nn.f9054 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 Dhecmw_precond_SSOR_nn.f9062 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 Dhecmw_precond_BILU_nn.f9058 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 Dbeam.cc240 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 Ddynamic_mat_ass_bc_ac.f9046 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 Ddynamic_mat_ass_bc_vl.f9056 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 Ddynamic_mat_ass_bc.f9042 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 Dhecmw_mat_ass.f9041 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 Dhecmw_matrix_reorder.f9091 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 Dfstr_EIG_output.f9026 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 Dfstr_EIG_lanczos_util.f9017 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 Dhecmw_solver.f9027 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 Dhecmw_solver_GPBiCG.f9037 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 Dhecmw_solver_CG.f9040 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 Dhecmw_solver_BiCGSTAB.f9035 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 Dhecmw_solver_GMRES.f9037 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 Dhecmw_util_f.F90894 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 DitkFEMElement2DC1Beam.cxx140 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 Dsolve_LINEQ_direct_serial_lag.f9041 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 Dfstr_contact.f90378 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 …]

12345