Home
last modified time | relevance | path

Searched refs:conMAT (Results 1 – 18 of 18) sorted by relevance

/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/analysis/static/
H A Dfstr_mat_con_contact.f90267 … if(associated(conMAT%indexL).and.associated(conMAT%indexU))deallocate(conMAT%indexL,conMAT%indexU)
268 allocate(conMAT%indexL(0:conMAT%NP), conMAT%indexU(0:conMAT%NP), stat=ierr)
270 conMAT%indexL = 0 ; conMAT%indexU = 0
271 if(associated(conMAT%itemL).and.associated(conMAT%itemU))deallocate(conMAT%itemL,conMAT%itemU)
274 conMAT%itemL = 0 ; conMAT%itemU = 0
407 if(associated(conMAT%AL)) deallocate(conMAT%AL)
408 allocate(conMAT%AL(nn*conMAT%NPL), stat=ierr)
413 allocate(conMAT%AU(nn*conMAT%NPU), stat=ierr)
417 if(associated(conMAT%B)) deallocate(conMAT%B)
421 if(associated(conMAT%X)) deallocate(conMAT%X)
[all …]
H A Dsolve_LINEQ_contact.f9046 subroutine solve_LINEQ_contact(hecMESH,hecMAT,fstrMAT,istat,rf,conMAT) argument
53 type (hecmwST_matrix),optional :: conMAT local
67 if(paraContactFlag.and.present(conMAT)) then
68 call solve_LINEQ_iter_contact(hecMESH,hecMAT,fstrMAT,istat,conMAT)
75 if(paraContactFlag.and.present(conMAT)) then
76 call solve_LINEQ_MKL_contact(hecMESH,hecMAT,fstrMAT,istat,conMAT)
82 if(paraContactFlag.and.present(conMAT)) then
83 call solve_LINEQ_mumps_contact(hecMESH,hecMAT,fstrMAT,istat,conMAT)
105 if(paraContactFlag.and.present(conMAT)) then
H A Dfstr_solve_NonLinear.f90429 type (hecmwST_matrix), optional :: conMAT local
476 if(paraContactFlag.and.present(conMAT)) then
499 if( paraContactFlag.and.present(conMAT) ) then
500 call hecmw_mat_clear_b(conMAT)
505 if(paraContactFlag.and.present(conMAT)) then
534 call hecmw_mat_clear( conMAT )
535 conMAT%X = 0.0d0
619 call hecmw_mat_clear_b( conMAT )
695 call hecmw_mat_clear_b( conMAT )
712 if(paraContactFlag.and.present(conMAT)) then
[all …]
H A Dfstr_AddBC.f9013 subroutine fstr_AddBC(cstep,hecMESH,hecMAT,fstrSOLID,fstrPARAM,fstrMAT,iter,conMAT) argument
28 type(hecmwST_matrix), optional :: conMAT !< hecmw matrix for contact only local
91 if(present(conMAT)) then
92 call hecmw_mat_ass_bc(hecMAT, in, idof, RHS, conMAT)
98 if(present(conMAT)) then
99 call fstr_mat_ass_bc_contact(conMAT,fstrMAT,in,idof,RHS)
138 if(present(conMAT)) then
139 call hecmw_mat_ass_bc(hecMAT, in, idof, RHS, conMAT)
145 if(present(conMAT)) then
146 call fstr_mat_ass_bc_contact(conMAT,fstrMAT,in,idof,RHS)
H A Dsolve_LINEQ_mkl.F9054 subroutine solve_LINEQ_MKL_contact(hecMESH,hecMAT,fstrMAT,istat,conMAT) argument
59 type (hecmwST_matrix), intent(in),optional :: conMAT local
73 if(paraContactFlag.and.present(conMAT)) then
74 call sparse_matrix_para_contact_set_vals(spMAT, hecMAT, fstrMAT, conMAT)
75 call sparse_matrix_para_contact_set_rhs(spMAT, hecMAT, fstrMAT, conMAT)
H A Dsolve_LINEQ_MUMPS_contact.f9065 subroutine solve_LINEQ_MUMPS_contact(hecMESH,hecMAT,fstrMAT,istat,conMAT) argument
71 type (hecmwST_matrix), intent(in),optional :: conMAT local
92 if(paraContactFlag.and.present(conMAT)) then
93 call sparse_matrix_para_contact_set_vals(spMAT, hecMAT, fstrMAT, conMAT)
94 call sparse_matrix_para_contact_set_rhs(spMAT, hecMAT, fstrMAT, conMAT)
H A Dfstr_Residual.f9023 subroutine fstr_Update_NDForce(cstep,hecMESH,hecMAT,fstrSOLID,conMAT) argument
37 type(hecmwST_matrix), intent(inout), optional :: conMAT local
60 if(present(conMAT)) call fstr_Update_NDForce_SPC( cstep, hecMESH, fstrSOLID, conMAT%B )
190 function fstr_get_norm_para_contact(hecMAT,fstrMAT,conMAT,hecMESH) result(rhsB)
196 type(hecmwST_matrix), intent(in) :: conMAT local
229 if(conMAT%b((i-1)*ndof+j) /= 0.0D0) then
230 print *,myrank,'i0',i,'conMAT%b',conMAT%b((i-1)*ndof+j)
235 rhs_con_all(i0+1:i0+ndof) = conMAT%b((i-1)*ndof+1:i*ndof)
242 rhs_con(i) = rhs_con_all(offset+i) + hecMAT%B(i) + conMAT%B(i)
247 rhs_con(nndof+i) = conMAT%B(i0+i)
H A Dsolve_LINEQ_iter_contact.f90331 write(1000+myrank,*) conMAT%B(1:conMAT%N*ndof)
332 write(1000+myrank,*) 'conMAT%B(external): ',conMAT%N*ndof+1,'-',conMAT%NP*ndof
333 write(1000+myrank,*) conMAT%B(conMAT%N*ndof+1:conMAT%NP*ndof)
336 write(1000+myrank,*) conMAT%B(conMAT%NP*ndof+1:conMAT%NP*ndof+fstrMAT%num_lagrange)
354 write(1000+myrank,*) Btot(conMAT%NP*ndof+1:conMAT%NP*ndof+fstrMAT%num_lagrange)
373 write(1000+myrank,*) Btot(conMAT%NP*ndof+1:conMAT%NP*ndof+fstrMAT%num_lagrange)
928 Btot(i) = conMAT%B(i)
1180 je = conMAT%indexL(irow)
1182 jcol = conMAT%itemL(j)
1193 je = conMAT%indexU(irow)
[all …]
H A Dfstr_solve_NLGEOM.f9028 subroutine FSTR_SOLVE_NLGEOM(hecMESH,hecMAT,fstrSOLID,fstrMAT,fstrPARAM,conMAT) argument
35 type (hecmwST_matrix ),optional :: conMAT local
154 if(paraContactFlag.and.present(conMAT)) then
156 …tep_num, restart_substep_num, sub_step, fstr_get_time(), fstr_get_timeinc(), infoCTChange, conMAT )
191 if(paraContactFlag.and.present(conMAT)) then
192 … call fstr_mat_con_contact( tot_step, hecMAT, fstrSOLID, fstrMAT, infoCTChange, conMAT)
193 conMAT%B(:) = 0.0d0
H A Dsparse_matrix_contact.f90503 type(hecmwST_matrix), intent(in) :: conMAT local
551 spMAT%A(m)=hecMAT%D(offset_d+jdof) + conMAT%D(offset_d+jdof)
617 spMAT%A(m) = conMAT%AL(offset_l+jdof)
633 spMAT%A(m) = conMAT%D(offset_d+jdof)
654 spMAT%A(m) = conMAT%AU(offset_u+jdof)
751 type (hecmwST_matrix), intent(in) :: conMAT local
763 if(conMAT%b((i-1)*ndof+j) /= 0.0D0) then
764 print *,myrank,'i0',i,spMAT%N,'conMAT%b',conMAT%b((i-1)*ndof+j)
773 if(conMAT%b((i-1)*ndof+j) /= 0.0D0) then
774 rhs_con(i0+j) = conMAT%b((i-1)*ndof+j)
[all …]
H A Dfstr_AddContactStiff.f90260 subroutine fstr_Update_NDForce_contact(cstep,hecMESH,hecMAT,fstrMAT,fstrSOLID,conMAT) argument
266 type(hecmwST_matrix), optional :: conMAT !< type hecmwST_matrix for contact part only local
310 if(present(conMAT)) then
311 … call update_NDForce_contact(nnode,ndLocal,id_lagrange,lagrange,ctNForce,ctTForce,fstrSOLID,conMAT)
322 if(present(conMAT)) call fstr_Update_NDForce_SPC(cstep, hecMESH, fstrSOLID, conMAT%B)
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/hecmw1/src/solver/matrix/
H A Dhecmw_mat_ass.f90344 if(present(conMAT)) conMAT%B(NDOF*inode-(NDOF-idof)) = 0.0D0
354 if(present(conMAT)) then
357 conMAT%B(idx) = conMAT%B(idx) - val
368 if(present(conMAT)) conMAT%D(ndof2*inode-ii+i)=0.d0
376 if(present(conMAT)) conMAT%D(ndof2*inode-ii) = 0.d0
379 if(present(conMAT)) conMAT%D(ndof2*inode-ii) = 0.d0
395 if(present(conMAT)) conMAT%AL(ndof2*k-ii+i) = 0.d0
411 if(present(conMAT)) then
414 conMAT%B(idx) = conMAT%B(idx) - val
434 if(present(conMAT)) conMAT%AU(ndof2*k-ii+i) = 0.d0
[all …]
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/analysis/dynamic/transit/
H A Dfstr_dynamic_nlimplicit.f90414 ,conMAT ) argument
436 type(hecmwST_matrix), optional :: conMAT local
557 if(paraContactFlag.and.present(conMAT)) then
563 if(paraContactFlag.and.present(conMAT)) then
652 call hecmw_mat_clear( conMAT )
653 call hecmw_mat_clear_b( conMAT )
654 conMAT%X = 0.0d0
682 …t_ass_bc (hecMESH, hecMATmpc, fstrSOLID, fstrDYNAMIC, fstrPARAM, fstrMAT, stepcnt, conMAT=conMAT)
683 …t_ass_bc_vl(hecMESH, hecMATmpc, fstrSOLID, fstrDYNAMIC, fstrPARAM, fstrMAT, stepcnt, conMAT=conMAT)
684 …t_ass_bc_ac(hecMESH, hecMATmpc, fstrSOLID, fstrDYNAMIC, fstrPARAM, fstrMAT, stepcnt, conMAT=conMAT)
[all …]
H A Ddynamic_mat_ass_bc.f9013 …utine DYNAMIC_MAT_ASS_BC(hecMESH, hecMAT, fstrSOLID ,fstrDYNAMIC, fstrPARAM, fstrMAT, iter, conMAT) argument
29 type(hecmwST_matrix), optional :: conMAT local
102 if(present(conMAT)) then
103 call hecmw_mat_ass_bc(hecMAT, in, idof, RHS, conMAT)
109 if(present(conMAT)) then
110 call fstr_mat_ass_bc_contact(conMAT,fstrMAT,in,idof,RHS)
152 if(present(conMAT)) then
153 call hecmw_mat_ass_bc(hecMAT, in, idof, RHS, conMAT)
159 if(present(conMAT)) then
160 call fstr_mat_ass_bc_contact(conMAT,fstrMAT,in,idof,RHS)
H A Ddynamic_mat_ass_bc_vl.f9013 …ne DYNAMIC_MAT_ASS_BC_VL(hecMESH, hecMAT, fstrSOLID, fstrDYNAMIC, fstrPARAM, fstrMAT, iter, conMAT) argument
27 type(hecmwST_matrix), optional :: conMAT local
98 if(present(conMAT)) then
99 call hecmw_mat_ass_bc(hecMAT, in, idof, RHS, conMAT)
105 if(present(conMAT)) then
106 call fstr_mat_ass_bc_contact(conMAT,fstrMAT,in,idof,RHS)
H A Ddynamic_mat_ass_bc_ac.f9014 …ne DYNAMIC_MAT_ASS_BC_AC(hecMESH, hecMAT, fstrSOLID ,fstrDYNAMIC, fstrPARAM, fstrMAT, iter, conMAT) argument
29 type(hecmwST_matrix), optional :: conMAT local
88 if(present(conMAT)) then
89 call hecmw_mat_ass_bc(hecMAT, in, idof, RHS, conMAT)
95 if(present(conMAT)) then
96 call fstr_mat_ass_bc_contact(conMAT,fstrMAT,in,idof,RHS)
H A Dfstr_solve_dynamic.f9023 ,conMAT ) argument
37 type(hecmwST_matrix), optional :: conMAT local
212 if(paraContactFlag.and.present(conMAT)) then
216 ,conMAT )
/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/main/
H A Dfistr_main.f9025 type(hecmwST_matrix), save :: conMAT variable
126 call hecmw_nullify_matrix ( conMAT )
307 call fstr_solve_NLGEOM( hecMESH, hecMAT, fstrSOLID, fstrMAT, fstrPR, conMAT )
379 conMAT )