/dports/science/frontistr/FrontISTR-c66bdc397de319ca59a0565b3f3b1a3b33f0c50c/fistr1/src/analysis/static/ |
H A D | fstr_mat_con_contact.f90 | 267 … 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 D | solve_LINEQ_contact.f90 | 46 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 D | fstr_solve_NonLinear.f90 | 429 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 D | fstr_AddBC.f90 | 13 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 D | solve_LINEQ_mkl.F90 | 54 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 D | solve_LINEQ_MUMPS_contact.f90 | 65 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 D | fstr_Residual.f90 | 23 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 D | solve_LINEQ_iter_contact.f90 | 331 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 D | fstr_solve_NLGEOM.f90 | 28 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 D | sparse_matrix_contact.f90 | 503 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 D | fstr_AddContactStiff.f90 | 260 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 D | hecmw_mat_ass.f90 | 344 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 D | fstr_dynamic_nlimplicit.f90 | 414 ,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 D | dynamic_mat_ass_bc.f90 | 13 …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 D | dynamic_mat_ass_bc_vl.f90 | 13 …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 D | dynamic_mat_ass_bc_ac.f90 | 14 …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 D | fstr_solve_dynamic.f90 | 23 ,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 D | fistr_main.f90 | 25 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 )
|