impbufdef_mod.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    impbufdef_mod    ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    dealloc_impbuf   ../engine/source/implicit/imp_solv.F
!||    imp_buck         ../engine/source/implicit/imp_buck.F
!||    imp_chkm         ../engine/source/implicit/imp_solv.F
!||    imp_fout         ../engine/source/implicit/imp_solv.F
!||    imp_restarcp     ../engine/source/implicit/imp_sol_init.F
!||    imp_sol_init     ../engine/source/implicit/imp_sol_init.F
!||    imp_solv         ../engine/source/implicit/imp_solv.F
!||    resol            ../engine/source/engine/resol.F
!||    wrrestp          ../engine/source/output/restart/wrrestp.F
!||====================================================================
      MODULE impbufdef_mod
C----------------------------------------------------------------------- 
#include "my_real.inc"
C=======================================================================
c DEFINE SIZES (integers arrays)
C=======================================================================
        integer ::   s_iadk               
        integer ::   s_jdik     
          integer ::   s_isij 
          integer ::   s_imij 
          integer ::   s_nss  
          integer ::   s_iss  
          integer ::   s_isij2
          integer ::   s_nss2 
          integer ::   s_iss2 
          integer ::   s_nsc2 
          integer ::   s_nrowk
          integer ::   s_icok 
          integer ::   s_icokm
          integer ::   s_nmij2
          integer ::   s_nss3 
          integer ::   s_isb2 
          integer ::   s_nsrb2
              integer ::   s_ikc  
              integer ::   s_ikud 
              integer ::   s_w_ddl
              integer ::   s_iadm 
        integer ::   s_jdim     
              integer ::   s_ndofi
              integer ::   s_iddli
              integer ::   s_ikinw
          integer ::   s_inbuf_c
C-
        integer ::   s_i_stok     !  replace NUM_IMP            
        integer ::   s_cand_n     !  replace NS_IMP  
        integer ::   s_cand_e     !  replace NE_IMP   
        integer ::   s_indsubt    !  replace IND_IMP  
C=======================================================================
c DEFINE SIZES (float arrays)
C=======================================================================
        integer ::   s_diag_k           
        integer ::   s_lt_k           
        integer ::   s_diag_m           
        integer ::   s_lt_m           
              integer ::   s_lb
              integer ::   s_lb0
              integer ::   s_bkud
              integer ::   s_d_imp
              integer ::   s_dr_imp
        integer ::   s_elbuf_c
          integer ::   s_bufmat_c
          integer ::   s_x_c
          integer ::   s_dd
          integer ::   s_ddr
          integer ::   s_x_a
          integer ::   s_fext
          integer ::   s_dg
          integer ::   s_dgr
          integer ::   s_dg0
          integer ::   s_dgr0
          integer ::   s_bufin_c
          integer ::   s_ac
          integer ::   s_acr
c=======================================================================      
c define type IMPBUF_STRUCT_ for implicit structure array
c=======================================================================      
C=======================================================================
      TYPE impbuf_struct_
C=======================================================================
c integers & small integer arrays
C=======================================================================
        integer ::   nddl               
        integer ::   nnzk               
C----NS_IMP... will be changed to POINTER and pt to %I_STOK... in resol
        integer ::   nrbyac               
        integer ::   nint2               
        integer ::   NMC               
        integer ::   nmc2               
        integer ::   nmonv               
        integer, DIMENSION(20) ::  lsize      
        integer, DIMENSION(20) ::  i_imp      
C=======================================================================
c DEFINE ARRAYS (integers arrays)
C=======================================================================
        integer, DIMENSION(:) , ALLOCATABLE ::  irbyac      
        integer, DIMENSION(:) , ALLOCATABLE ::  nsc      
        integer, DIMENSION(:) , ALLOCATABLE ::  iint2      
        integer, DIMENSION(:) , ALLOCATABLE ::  nkud      
        integer, DIMENSION(:) , ALLOCATABLE ::  imonv      
        integer, DIMENSION(:) , ALLOCATABLE ::  iddl      
        integer, DIMENSION(:) , ALLOCATABLE ::  ndof      
        integer, DIMENSION(:) , ALLOCATABLE ::  inloc    
        integer, DIMENSION(:) , ALLOCATABLE ::  iadk    
        integer, DIMENSION(:) , ALLOCATABLE ::  jdik    
          integer, DIMENSION(:) , ALLOCATABLE ::  isij 
          integer, DIMENSION(:) , ALLOCATABLE ::  imij 
          integer, DIMENSION(:) , ALLOCATABLE ::  nss  
          integer, DIMENSION(:) , ALLOCATABLE ::  iss  
          integer, DIMENSION(:) , ALLOCATABLE ::  isij2
          integer, DIMENSION(:) , ALLOCATABLE ::  nss2 
          integer, DIMENSION(:) , ALLOCATABLE ::  iss2 
          integer, DIMENSION(:) , ALLOCATABLE ::  nsc2 
          integer, DIMENSION(:) , ALLOCATABLE ::  nrowk
          integer, DIMENSION(:) , ALLOCATABLE ::  icok 
          integer, DIMENSION(:) , ALLOCATABLE ::  icokm
          integer, DIMENSION(:) , ALLOCATABLE ::  nmij2
          integer, DIMENSION(:) , ALLOCATABLE ::  nss3 
          integer, DIMENSION(:) , ALLOCATABLE ::  isb2 
          integer, DIMENSION(:) , ALLOCATABLE ::  nsrb2
              integer, DIMENSION(:) , ALLOCATABLE ::  ikc  
              integer, DIMENSION(:) , ALLOCATABLE ::  ikud 
              integer, DIMENSION(:) , ALLOCATABLE ::  w_ddl
              integer, DIMENSION(:) , ALLOCATABLE ::  iadm 
              integer, DIMENSION(:) , ALLOCATABLE ::  jdim 
              integer, DIMENSION(:) , ALLOCATABLE ::  ndofi
              integer, DIMENSION(:) , ALLOCATABLE ::  iddli
              integer, DIMENSION(:) , ALLOCATABLE ::  ikinw
              integer, DIMENSION(:) , ALLOCATABLE ::  inbuf_c
C---for contact
        integer, DIMENSION(:) , ALLOCATABLE ::  i_stok      
        integer, DIMENSION(:) , ALLOCATABLE ::  cand_n      
        integer, DIMENSION(:) , ALLOCATABLE ::  cand_e    
        integer, DIMENSION(:) , ALLOCATABLE ::  indsubt   
C=======================================================================
c DEFINE ARRAYS (float arrays)
C=======================================================================
        my_real, DIMENSION(25) ::  r_imp      
        my_real, DIMENSION(:) , ALLOCATABLE ::   diag_k           
        my_real, DIMENSION(:) , ALLOCATABLE ::   lt_k           
        my_real, DIMENSION(:) , ALLOCATABLE ::   diag_m           
        my_real, DIMENSION(:) , ALLOCATABLE ::   lt_m           
              my_real, DIMENSION(:) , ALLOCATABLE ::   lb
              my_real, DIMENSION(:) , ALLOCATABLE ::   lb0
              my_real, DIMENSION(:) , ALLOCATABLE ::   bkud
              my_real, DIMENSION(:) , ALLOCATABLE ::   d_imp
              my_real, DIMENSION(:) , ALLOCATABLE ::   dr_imp
        my_real, DIMENSION(:) , ALLOCATABLE ::   elbuf_c
          my_real, DIMENSION(:) , ALLOCATABLE ::   bufmat_c
          my_real, DIMENSION(:) , ALLOCATABLE ::   x_c
          my_real, DIMENSION(:) , ALLOCATABLE ::   dd
          my_real, DIMENSION(:) , ALLOCATABLE ::   ddr
          my_real, DIMENSION(:) , ALLOCATABLE ::   x_a
          my_real, DIMENSION(:) , ALLOCATABLE ::   fext
          my_real, DIMENSION(:) , ALLOCATABLE ::   dg
          my_real, DIMENSION(:) , ALLOCATABLE ::   dgr
          my_real, DIMENSION(:) , ALLOCATABLE ::   dg0
          my_real, DIMENSION(:) , ALLOCATABLE ::   dgr0
          my_real, DIMENSION(:) , ALLOCATABLE ::   bufin_c
          my_real, DIMENSION(:) , ALLOCATABLE ::   ac
          my_real, DIMENSION(:) , ALLOCATABLE ::   acr
C=======================================================================
      END TYPE impbuf_struct_
C=======================================================================
c---------------
      END MODULE impbufdef_mod
!||====================================================================
!||    imp_ktan_def   ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    cmatch3        ../engine/source/elements/shell/coqueba/cmatc3.F
!||    cmatip3        ../engine/source/elements/shell/coqueba/cmatc3.F
!||    etfac_ini      ../engine/source/implicit/imp_init.F
!||    get_etfac_s    ../engine/source/elements/solid/solide8z/get_etfac_s.F
!||    gethkt3        ../engine/source/elements/solid/solide8z/gethkt3.F
!||    iktmat_ini     ../engine/source/implicit/imp_init.F
!||    imp_ktan       ../engine/share/modules/impbufdef_mod.F
!||    ktbuf_ini      ../engine/source/implicit/imp_init.F
!||    mmats          ../engine/source/elements/solid/solide8z/mmats.F
!||    put_etfac      ../engine/source/elements/solid/solide8z/put_etfac.F
!||    putsignor3     ../engine/source/elements/solid/solide8z/putsignor3.F
!||    putsignorc3    ../engine/source/elements/shell/coqueba/cmatc3.F
!||    sktcons2       ../engine/source/elements/solid/solide8z/sktcons2.F
!||====================================================================
      MODULE imp_ktan_def
#include "my_real.inc"
C------
      TYPE mlaw_tag_    ! element buffer variables depending on material law
        INTEGER ::   l_etfac      ! = 1
        INTEGER ::   l_sige       ! = 6 ; = 5 for shell
        INTEGER ::   l_a_kt       
        INTEGER ::   l_s ubkt     ! =6
      END TYPE mlaw_tag_  
      
C      TYPE L_KTBUFG_      ! general law (hyper-elastic, elasto-plastic)
C        my_real, DIMENSION(:) , POINTER ::   ETFAC    !  relative tangent modulus ratio per ele
C      END TYPE L_KTBUFG_
      
      TYPE l_ktbufep_      ! elastoplastic material per integration point
        my_real, DIMENSION(:) , POINTER ::   a_kt     !  alpha   dgama for shell
        my_real, DIMENSION(:) , POINTER ::   sige     !  trail stress                              
      END TYPE l_ktbufep_ 
      
      TYPE l_ktbufgm_      !  modulus KT general matrix per integration point
C----------------|DD(L_SUBKT) | DG(L_SUBKT) |
C----------------|  sym       | GG(L_SUBKT) |
        my_real, DIMENSION(:) , POINTER ::   kt_dd
        my_real, DIMENSION(:) , POINTER ::   kt_dg
        my_real, DIMENSION(:) , POINTER ::   kt_gg
      END TYPE l_ktbufgm_
 
      TYPE ktbuf_struct_
C------the following integer arguments are not necessary for the moment
c        integer     :: MLAW           !  1
c        integer     :: NEL            !  2
c        integer     :: NPTR           !  3
c        integer     :: NPTS           !  4
c        integer     :: NPTT           !  5
C ------------add MLAW_TAG here to avoid pass everywhere the arguments
        my_real  ,          DIMENSION(:)       , ALLOCATABLE :: etfac   ! relative tangent modulus ratio - per element(global)
        TYPE(mlaw_tag_)   , DIMENSION(:)       , ALLOCATABLE :: mlaw_tag ! variable dimensions - per element(global)
        TYPE (l_ktbufep_) , DIMENSION(:,:,:)   , ALLOCATABLE :: ktbufep  ! EP mat - per integration point
        TYPE (l_ktbufgm_) , DIMENSION(:,:,:)   , ALLOCATABLE :: ktbufmg  ! ET mat - per integration point
      END TYPE ktbuf_struct_
C      
      END MODULE imp_ktan_def
!||====================================================================
!||    imp_bfgs       ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    bfgs_0         ../engine/source/implicit/imp_bfgs.F
!||    bfgs_1         ../engine/source/implicit/imp_bfgs.F
!||    bfgs_1p        ../engine/source/implicit/imp_bfgs.F
!||    bfgs_2         ../engine/source/implicit/imp_bfgs.F
!||    bfgs_2p        ../engine/source/implicit/imp_bfgs.F
!||    bfgs_h1        ../engine/source/implicit/imp_bfgs.F
!||    bfgs_h1p       ../engine/source/implicit/imp_bfgs.F
!||    bfgs_h2        ../engine/source/implicit/imp_bfgs.F
!||    bfgs_h2p       ../engine/source/implicit/imp_bfgs.F
!||    bfgs_ini       ../engine/source/implicit/imp_bfgs.F
!||    bfgs_ls        ../engine/source/implicit/imp_bfgs.F
!||    deallocm_imp   ../engine/source/implicit/imp_solv.F
!||    get_slin       ../engine/source/implicit/imp_bfgs.F
!||    nsloan_0       ../engine/source/implicit/imp_bfgs.F
!||    nsloan_5       ../engine/source/implicit/imp_bfgs.F
!||====================================================================
      MODULE imp_bfgs
C
#include "my_real.inc"
      INTEGER n_bfgs,iactb
      my_real s_lin
C
C tableau de int de : V,W pour BFGS algorithm----
      my_real, DIMENSION(:,:), ALLOCATABLE :: bfgs_v,bfgs_w
C 
      END MODULE imp_bfgs
!||====================================================================
!||    imp_dyna       ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    deallocm_imp   ../engine/source/implicit/imp_solv.F
!||    du_ini         ../engine/source/implicit/imp_solv.F
!||    du_ini_hp      ../engine/source/implicit/imp_solv.F
!||    dyna_cpk0      ../engine/source/implicit/imp_dyna.F
!||    dyna_cpr0      ../engine/source/implicit/imp_dyna.F
!||    dyna_in0       ../engine/source/implicit/imp_dyna.F
!||    dyna_ina       ../engine/source/implicit/imp_dyna.F
!||    dyna_ini       ../engine/source/implicit/imp_dyna.F
!||    dyna_iniv      ../engine/source/implicit/imp_dyna.F
!||    dyna_ivfac     ../engine/source/implicit/imp_dyna.F
!||    dyna_wex       ../engine/source/implicit/imp_dyna.F
!||    ecrit          ../engine/source/output/ecrit.F
!||    freimpl        ../engine/source/input/freimpl.F
!||    fv_fint0       ../engine/source/constraints/general/impvel/fv_imp0.F
!||    getdyna_a      ../engine/source/implicit/imp_dyna.F
!||    imp_dycrb      ../engine/source/implicit/imp_dyna.F
!||    imp_dykv       ../engine/source/implicit/imp_dyna.F
!||    imp_dykv0      ../engine/source/implicit/imp_dyna.F
!||    imp_dynam      ../engine/source/implicit/imp_dyna.F
!||    imp_dynar      ../engine/source/implicit/imp_dyna.F
!||    imp_fhht       ../engine/source/implicit/imp_dyna.F
!||    imp_fhht1      ../engine/source/implicit/imp_dyna.F
!||    imprrest       ../engine/source/output/restart/rdresb.F
!||    impwrest       ../engine/source/output/restart/wrrest.F
!||    inte_dyna      ../engine/source/implicit/imp_dyna.F
!||    rbe2cor        ../engine/source/constraints/general/rbody/rgbcor.F
!||    rbe2cor0       ../engine/source/constraints/general/rbody/rgbcor.F
!||    rbycor         ../engine/source/constraints/general/rbody/rbycor.F
!||    rgbcor         ../engine/source/constraints/general/rbody/rgbcor.F
!||    stop_sensor    ../engine/source/tools/sensor/stop_sensor.F
!||    wfv_imp        ../engine/source/constraints/general/impvel/fv_imp0.F
!||====================================================================
      MODULE imp_dyna
C
#include "my_real.inc"
C
C tableau de travail :  implicit dynamic ----
      my_real  dy_g,dy_b,dy_edamp
      my_real, DIMENSION(:,:), ALLOCATABLE :: dy_d,dy_dr,dy_v
      my_real, DIMENSION(:,:), ALLOCATABLE :: dy_vr,dy_a,dy_ar
      my_real, DIMENSION(:,:), ALLOCATABLE :: dy_in
      my_real, DIMENSION(:,:), ALLOCATABLE :: dy_dam,dy_damr
      my_real, DIMENSION(:), ALLOCATABLE :: dy_diak0,dy_ltk0
      my_real, DIMENSION(:), ALLOCATABLE :: dy_r0,dy_r1
      my_real, DIMENSION(:,:), ALLOCATABLE :: dy_dam0,dy_damr0
      INTEGER, DIMENSION(:), ALLOCATABLE :: dy_iadk0,dy_jdik0
C 
      END MODULE imp_dyna
!||====================================================================
!||    imp_monv    ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    id_mvini    ../engine/source/airbag/monv_imp0.F
!||    monv_diag   ../engine/source/airbag/monv_imp0.F
!||    monv_fvl    ../engine/source/airbag/monv_imp0.F
!||    monv_m3     ../engine/source/airbag/monv_imp0.F
!||    monv_prem   ../engine/source/airbag/monv_imp0.F
!||    mv_matv     ../engine/source/airbag/monv_imp0.F
!||    recu_kdis   ../engine/source/airbag/monv_imp0.F
!||====================================================================
      MODULE imp_monv
#include "my_real.inc"
      INTEGER numn_mv,ni2_mv,nrb_mv,nnmax_mv,nbc_mv,nfx_mv,nrw_mv,
     .        nrbe3_mv,r3m_max,nspc_mv,nrbe2_mv
C  int de taille NUMN_MV :noeuds locales de mv 
      INTEGER, DIMENSION(:), ALLOCATABLE :: in_mv
C  int IDDL de taille 3*NUMN_MV,6*(NRB_FR ,4*NI2_FR)
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: id_mv,id_mvm,id_mvm4
      INTEGER, DIMENSION(:,:,:), ALLOCATABLE :: id_mvm2,id_mvm3
C int de taille NI2_MV,NRB_MV :------
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: ii2_mv,irb_mv,irbe2_mv
      INTEGER, DIMENSION(:), ALLOCATABLE ::   irbe3_mv
C  K_DIAG reels de taille (3,NUMN_MV) 
      my_real, DIMENSION(:,:),ALLOCATABLE :: diag_mv
C K_DIAG reels de taille (6,NRB_MV) 
      my_real, DIMENSION(:,:),ALLOCATABLE :: diag_mvm,diag_mvm4
C K_DIAG reels de taille (6,4,NI2_MV) 
      my_real, DIMENSION(:,:,:),ALLOCATABLE :: diag_mvm2
      my_real, DIMENSION(:,:,:),ALLOCATABLE :: diag_mvm3
C int (3,NBC_MV),(2,NFX_MV),(NRW_MV)
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: ibc_mv,ifx_mv
      INTEGER, DIMENSION(:), ALLOCATABLE :: irw_mv,ispc_mv
      my_real, DIMENSION(:), ALLOCATABLE :: fcdi_mv,mcdi_mv
C      
      END MODULE imp_monv
!||====================================================================
!||    imp_rwl      ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    diag_int     ../engine/source/mpi/implicit/imp_fri.F
!||    dim_kinkn    ../engine/source/implicit/imp_int_k.F
!||    fr_u2dd      ../engine/source/mpi/implicit/imp_fri.F
!||    fv_rwl       ../engine/source/constraints/general/rwall/srw_imp.F
!||    fv_rwl0      ../engine/source/constraints/general/rwall/srw_imp.F
!||    fv_rwlr0     ../engine/source/constraints/general/rwall/srw_imp.F
!||    imp3_a2b     ../engine/source/airbag/monv_imp0.F
!||    imp3_u2x     ../engine/source/airbag/monv_imp0.F
!||    ind_kinefr   ../engine/source/mpi/implicit/imp_fri.F
!||    ind_kinfrk   ../engine/source/mpi/implicit/imp_fri.F
!||    ini_kinkn    ../engine/source/implicit/imp_int_k.F
!||    monv_fvl     ../engine/source/airbag/monv_imp0.F
!||    rwl_impd     ../engine/source/constraints/general/rwall/srw_imp.F
!||    upd_fr       ../engine/source/mpi/implicit/imp_fri.F
!||    upd_kml      ../engine/source/mpi/implicit/imp_fri.F
!||    upd_ksl      ../engine/source/mpi/implicit/imp_fri.F
!||    updk_mv      ../engine/source/airbag/monv_imp0.F
!||====================================================================
      MODULE imp_rwl
#include "my_real.inc"
      INTEGER n_rwl
C tableau de int de taille N_RWL :noeuds secnds de sliding RW 
      INTEGER, DIMENSION(:), ALLOCATABLE :: in_rwl
C tableau de Direction normale reels de taille (3,N_RWL) 
      my_real, DIMENSION(:,:),ALLOCATABLE :: nor_rwl
C      
      END MODULE imp_rwl
!||====================================================================
!||    imp_lintf   ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    diag_kif    ../engine/source/implicit/imp_solv.F
!||    imp_solv    ../engine/source/implicit/imp_solv.F
!||    ini_kif     ../engine/source/implicit/imp_solv.F
!||    matv_kif    ../engine/source/implicit/imp_solv.F
!||    save_kif    ../engine/source/implicit/imp_solv.F
!||====================================================================
      MODULE imp_lintf
#include "my_real.inc"
      INTEGER nddlif,nzif
C tableau de int de la [k] interface 
      INTEGER, DIMENSION(:), ALLOCATABLE :: iadif,jdiif,iftok
C tableau de reels de la [k] interface 
      my_real, DIMENSION(:),ALLOCATABLE :: diag_if,lt_if
C      
      END MODULE imp_lintf
!||====================================================================
!||    imp_rest     ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    imp_trans    ../engine/source/output/restart/wrrest.F
!||    imp_trans0   ../engine/source/output/restart/wrrest.F
!||    imprrest     ../engine/source/output/restart/rdresb.F
!||    impwrest     ../engine/source/output/restart/wrrest.F
!||====================================================================
      MODULE imp_rest
#include "my_real.inc"
C tableau de reels de R-file 
      my_real, DIMENSION(:),ALLOCATABLE :: imp_rr
C      
      END MODULE imp_rest
C
!||====================================================================
!||    imp_ppat    ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    dim_span    ../engine/source/implicit/ind_glob_k.F
!||    fil_span0   ../engine/source/implicit/ind_glob_k.F
!||    fil_span1   ../engine/source/implicit/ind_glob_k.F
!||    ind_span    ../engine/source/implicit/ind_glob_k.F
!||====================================================================
      MODULE imp_ppat
C tableau de int de pre-filtrage 
      INTEGER, DIMENSION(:),ALLOCATABLE :: pre_fpat
C      
      END MODULE imp_ppat
!||====================================================================
!||    imp_qstat      ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    deallocm_imp   ../engine/source/implicit/imp_solv.F
!||    dis_cp         ../engine/source/implicit/imp_solv.F
!||    qstat_end      ../engine/source/implicit/imp_dyna.F
!||    qstat_ini      ../engine/source/implicit/imp_dyna.F
!||    qstat_it       ../engine/source/implicit/imp_dyna.F
!||====================================================================
      MODULE imp_qstat
C
#include "my_real.inc"
C
C tableau de travail :  linear quasi-static (iterative)----
      my_real, DIMENSION(:), ALLOCATABLE :: qs_u,qs_d
C tableau de travail : D_n-1 used for Gravity case
      my_real, DIMENSION(:), ALLOCATABLE :: d_n_1,dr_n_1
C 
      END MODULE imp_qstat
!||====================================================================
!||    imp_kbcs       ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    deallocm_imp   ../engine/source/implicit/imp_solv.F
!||    freimpl        ../engine/source/input/freimpl.F
!||    imp_init       ../engine/source/implicit/imp_init.F
!||    imp_mumps1     ../engine/source/implicit/imp_mumps.F
!||    lecimpl        ../engine/source/input/lectur.F
!||    prout_buck     ../engine/source/input/lectur.F
!||====================================================================
      MODULE imp_kbcs
C tableau de travail de BCS-EXT 
      INTEGER  mem_bcs
      my_real, DIMENSION(:),ALLOCATABLE :: hold
      INTEGER  msgl_b,maxset_b,msg_lvl,b_order,b_mcore,lmemn
      my_real  emin_b,emax_b,shift_b
C      
      END MODULE imp_kbcs
!||====================================================================
!||    imp_knon     ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    deallocm     ../engine/source/implicit/imp_solv.F
!||    int_matv     ../engine/source/implicit/imp_int_k.F
!||    int_matvp    ../engine/source/implicit/imp_int_k.F
!||    kin_knl      ../engine/source/implicit/imp_int_k.F
!||    recu_kdis0   ../engine/source/implicit/lin_solv.F
!||====================================================================
      MODULE imp_knon
#include "my_real.inc"
      INTEGER numn_kn,ni2_kn,nrb_kn,nbc_kn,nfx_kn,nrw_kn,nrbe3_kn,
     .        rkn_max,nspc_kn,nrbe2_kn
      INTEGER, DIMENSION(:), ALLOCATABLE :: in_kn
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: id_kn,id_knm,id_knm4
      INTEGER, DIMENSION(:,:,:), ALLOCATABLE :: id_knm2,id_knm3
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: ii2_kn,irb_kn,irbe2_kn
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: ibc_kn,ifx_kn
      INTEGER, DIMENSION(:), ALLOCATABLE :: irw_kn,irbe3_kn,ispc_kn
      my_real, DIMENSION(:), ALLOCATABLE :: fcdi_kn,mcdi_kn
C      
      END MODULE imp_knon
!||====================================================================
!||    imp_aspc     ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    bc_imp2      ../engine/source/constraints/general/bcs/bc_imp0.F
!||    bc_impa      ../engine/source/constraints/general/bcs/bc_imp0.F
!||    bc_impr1     ../engine/source/constraints/general/bcs/bc_imp0.F
!||    diag_int     ../engine/source/mpi/implicit/imp_fri.F
!||    dim_kinkn    ../engine/source/implicit/imp_int_k.F
!||    fr_u2dd      ../engine/source/mpi/implicit/imp_fri.F
!||    get_nspc     ../engine/source/constraints/general/bcs/bc_imp0.F
!||    imp3_a2b     ../engine/source/airbag/monv_imp0.F
!||    imp3_u2x     ../engine/source/airbag/monv_imp0.F
!||    imp_init     ../engine/source/implicit/imp_init.F
!||    ind_kinefr   ../engine/source/mpi/implicit/imp_fri.F
!||    ind_kinfrk   ../engine/source/mpi/implicit/imp_fri.F
!||    ini_kinkn    ../engine/source/implicit/imp_int_k.F
!||    kin_knl      ../engine/source/implicit/imp_int_k.F
!||    monv_fvl     ../engine/source/airbag/monv_imp0.F
!||    put_nspc     ../engine/source/constraints/general/bcs/bc_imp0.F
!||    upd_aspc0    ../engine/source/constraints/general/bcs/bc_imp0.F
!||    upd_fr       ../engine/source/mpi/implicit/imp_fri.F
!||    upd_kml      ../engine/source/mpi/implicit/imp_fri.F
!||    upd_ksl      ../engine/source/mpi/implicit/imp_fri.F
!||    updk_mv      ../engine/source/airbag/monv_imp0.F
!||====================================================================
      MODULE imp_aspc
      INTEGER nspcl,nspcnt
#include "my_real.inc"
      INTEGER, DIMENSION(:),ALLOCATABLE :: in_spc,ic_spc
      my_real, DIMENSION(:), ALLOCATABLE :: skew_spc
C      
      END MODULE imp_aspc
C
!||====================================================================
!||    imp_worki   ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    imp_chkm    ../engine/source/implicit/imp_solv.F
!||    imp_solv    ../engine/source/implicit/imp_solv.F
!||====================================================================
      MODULE imp_worki
#include      "my_real.inc"
      INTEGER, DIMENSION(:),ALLOCATABLE :: iadi,jdii,itok
      INTEGER  nddli ,nnzi,iconta,nsrem ,nsl
      my_real,
     .         DIMENSION(:),ALLOCATABLE :: diag_i,lt_i,xi_c
 
      END MODULE imp_worki
C
!||====================================================================
!||    imp_workh      ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    deallocm_imp   ../engine/source/implicit/imp_solv.F
!||    imp_pcgh       ../engine/source/implicit/imp_pcg.F
!||    imp_ppcgh      ../engine/source/implicit/imp_pcg.F
!||    ini_k0h        ../engine/source/implicit/imp_solv.F
!||    ini_kisc       ../engine/source/implicit/lin_solv.F
!||    lin_solv       ../engine/source/implicit/lin_solv.F
!||    lin_solvh0     ../engine/source/implicit/lin_solv.F
!||    lin_solvh1     ../engine/source/implicit/lin_solv.F
!||    lin_solvhm     ../engine/source/implicit/lin_solv.F
!||    lin_solvih2    ../engine/source/implicit/lin_solv.F
!||    mav_liuh       ../engine/source/implicit/produt_v.F
!||    mav_lu_h       ../engine/source/implicit/produt_v.F
!||    mav_lui_h      ../engine/source/implicit/produt_v.F
!||    mmav_lth       ../engine/source/implicit/produt_v.F
!||    prec5hc_solv   ../engine/source/implicit/prec_solv.F
!||    produt_uh      ../engine/source/implicit/produt_v.F
!||    produt_uh2     ../engine/source/implicit/produt_v.F
!||    produt_vmh     ../engine/source/implicit/produt_v.F
!||====================================================================
      MODULE imp_workh
#include      "my_real.inc"
      INTEGER  w_maxl,istop_h
      my_real k_lamda0,k_lamda1
      my_real,
     .         DIMENSION(:),ALLOCATABLE :: diag_t,l_u,l_f0
      my_real,
     .         DIMENSION(:),ALLOCATABLE :: pcg_w1,pcg_w2,pcg_w3
      my_real,
     .         DIMENSION(:),ALLOCATABLE :: tmp_w1,tmp_w2
#if defined knf
!dec$ attributes offload:mic :: IADK0,JDIK0,IADM0,JDIM0
#endif
      INTEGER, DIMENSION(:),ALLOCATABLE :: iadk0,jdik0,iadm0,jdim0
#if defined knf
!dec$ attributes offload:mic :: IADI0,JDII0
#endif
      INTEGER, DIMENSION(:),ALLOCATABLE :: iadi0,jdii0
#if defined knf
!dec$ attributes offload:mic :: LT_K0,LT_M0,LT_I0
#endif
      my_real,
     .        DIMENSION(:),ALLOCATABLE :: lt_k0,lt_m0,lt_i0
 
      END MODULE imp_workh
C
!||====================================================================
!||    imp_workg   ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    imp_buck    ../engine/source/implicit/imp_buck.F
!||====================================================================
      MODULE imp_workg
#include      "my_real.inc"
      INTEGER  iprmes_el(30)
      my_real,
     .         DIMENSION(:),ALLOCATABLE :: diag_kg,lt_kg
 
      END MODULE imp_workg
C
!||====================================================================
!||    imp_fvbcl     ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    bc_impr1      ../engine/source/constraints/general/bcs/bc_imp0.F
!||    fv_imp1       ../engine/source/constraints/general/impvel/fv_imp0.F
!||    fvbc_allo     ../engine/source/constraints/general/impvel/fv_imp0.F
!||    fvbc_deallo   ../engine/source/constraints/general/impvel/fv_imp0.F
!||    recu_ul       ../engine/source/constraints/general/impvel/fv_imp0.F
!||    recukin       ../engine/source/implicit/recudis.F
!||    upd_glob_k    ../engine/source/implicit/upd_glob_k.F
!||====================================================================
      MODULE imp_fvbcl
#include      "my_real.inc"
      INTEGER  nfvbcl,nkud_l,nkud_1
      INTEGER, DIMENSION(:),ALLOCATABLE :: ict_1,icr_1,ikud_1
      my_real,
     .         DIMENSION(:),ALLOCATABLE :: bkud_1,fvbcudl
 
      END MODULE imp_fvbcl
C
!||====================================================================
!||    imp_pcg_proj   ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    deallocm_imp   ../engine/source/implicit/imp_solv.F
!||    freimpl        ../engine/source/input/freimpl.F
!||    imp_inisi      ../engine/source/implicit/imp_pcg.F
!||    imp_inist      ../engine/source/implicit/imp_pcg.F
!||    imp_inix       ../engine/source/implicit/imp_pcg.F
!||    imp_ppcgh      ../engine/source/implicit/imp_pcg.F
!||    imp_pro_p      ../engine/source/implicit/imp_pcg.F
!||    imp_solv       ../engine/source/implicit/imp_solv.F
!||    imp_updst      ../engine/source/implicit/imp_pcg.F
!||    imp_updv2      ../engine/source/implicit/imp_pcg.F
!||    ini_k0h        ../engine/source/implicit/imp_solv.F
!||    lecimpl        ../engine/source/input/lectur.F
!||====================================================================
      MODULE imp_pcg_proj
#include      "my_real.inc"
      INTEGER  m_vs,ncg_run
      my_real
     .         k_tmp(2,2),m_tmp(2,2)
      my_real,
     .         DIMENSION(:,:),ALLOCATABLE :: proj_s,proj_t,proj_k
      my_real,
     .         DIMENSION(:),ALLOCATABLE :: proj_la_1,proj_v,proj_w
 
      END MODULE imp_pcg_proj
C
!||====================================================================
!||    imp_ktan       ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    cmatch3        ../engine/source/elements/shell/coqueba/cmatc3.F
!||    cmatip3        ../engine/source/elements/shell/coqueba/cmatc3.F
!||    etfac_ini      ../engine/source/implicit/imp_init.F
!||    get_etfac_s    ../engine/source/elements/solid/solide8z/get_etfac_s.F
!||    gethkt3        ../engine/source/elements/solid/solide8z/gethkt3.F
!||    iktmat_ini     ../engine/source/implicit/imp_init.F
!||    ktbuf_ini      ../engine/source/implicit/imp_init.F
!||    mmats          ../engine/source/elements/solid/solide8z/mmats.F
!||    put_etfac      ../engine/source/elements/solid/solide8z/put_etfac.F
!||    putsignor3     ../engine/source/elements/solid/solide8z/putsignor3.F
!||    putsignorc3    ../engine/source/elements/shell/coqueba/cmatc3.F
!||    sktcons2       ../engine/source/elements/solid/solide8z/sktcons2.F
!||--- uses       -----------------------------------------------------
!||    imp_ktan_def   ../engine/share/modules/impbufdef_mod.F
!||====================================================================
      MODULE imp_ktan
 
      USE imp_ktan_def
 
      TYPE(ktbuf_struct_),TARGET,DIMENSION(:),ALLOCATABLE :: ktbuf_str
      END MODULE imp_ktan
!||====================================================================
!||    imp_spbrm      ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    deallocm_imp   ../engine/source/implicit/imp_solv.F
!||    freimpl        ../engine/source/input/freimpl.F
!||    lecimpl        ../engine/source/input/lectur.F
!||    spb_ieref3     ../engine/source/implicit/imp_solv.F
!||    spb_ieref_bc   ../engine/source/implicit/imp_solv.F
!||    spb_ref_nds    ../engine/source/implicit/imp_solv.F
!||    spb_rm_rig     ../engine/source/implicit/imp_solv.F
!||    spbrm_pre      ../engine/source/implicit/imp_solv.F
!||    spmd_e_ref     ../engine/source/mpi/implicit/imp_spmd.F
!||    spmd_n_ref     ../engine/source/mpi/implicit/imp_spmd.F
!||====================================================================
      MODULE imp_spbrm
 
      INTEGER  e_ref(4),nbc_b,ne_bc4,ne_bc3,ikce(6),n_seg,ilskew,lskew_g
      INTEGER, DIMENSION(:),ALLOCATABLE :: ibc_b
      INTEGER, DIMENSION(:),ALLOCATABLE :: ie_bc4,ie_bc3
C------introduced due to SPMD-----      
      my_real
     .        x_ref(3,4),d_ref(3,4),rlskew(9)
C-----E_REF(4)4N of elem_ref, NBC_B ,IBC_B(NBC_B): node_id w/ BC;
C------IE_BC(NE_BC) Elem_id which has the BC nodes 
C-------N_SEG=3 3N seg, 0 : no ref;  ----   
 
      END MODULE imp_spbrm
C
!||====================================================================
!||    imp_frk        ../engine/share/modules/impbufdef_mod.F
!||--- called by ------------------------------------------------------
!||    dim_fr_k       ../engine/source/mpi/implicit/imp_fri.F
!||    dim_nrmax      ../engine/source/mpi/implicit/imp_fri.F
!||    fr_dlft        ../engine/source/mpi/implicit/imp_fri.F
!||    get_ikin2g     ../engine/source/mpi/implicit/imp_fri.F
!||    getnddli_g     ../engine/source/mpi/implicit/imp_fri.F
!||    imp_pcgh       ../engine/source/implicit/imp_pcg.F
!||    ind_fr_k       ../engine/source/mpi/implicit/imp_fri.F
!||    ind_fr_k0      ../engine/source/mpi/implicit/imp_fri.F
!||    ind_kine_kp    ../engine/source/mpi/implicit/imp_fri.F
!||    ind_nrfr       ../engine/source/mpi/implicit/imp_fri.F
!||    ini_fr_k       ../engine/source/mpi/implicit/imp_fri.F
!||    kin_nrmax      ../engine/source/mpi/implicit/imp_fri.F
!||    kin_nrmax0     ../engine/source/mpi/implicit/imp_fri.F
!||    nddli_frb      ../engine/source/mpi/implicit/imp_fri.F
!||    pr_deb         ../engine/source/implicit/imp_solv.F
!||    pr_matrix      ../engine/source/implicit/imp_solv.F
!||    pr_solnfo      ../engine/source/implicit/imp_solv.F
!||    set_ikin2g     ../engine/source/mpi/implicit/imp_fri.F
!||    spc_fr_k       ../engine/source/mpi/implicit/imp_fri.F
!||    spmd_inf_g     ../engine/source/mpi/implicit/imp_spmd.F
!||    spmd_max_iv    ../engine/source/mpi/implicit/imp_spmd.F
!||    spmd_sumf_k    ../engine/source/mpi/implicit/imp_spmd.F
!||    spmd_sumf_v    ../engine/source/mpi/implicit/imp_spmd.F
!||    spmd_sumfc_v   ../engine/source/mpi/implicit/imp_spmd.F
!||    upd_fr_k       ../engine/source/mpi/implicit/imp_fri.F
!||    zero_ikin2g    ../engine/source/mpi/implicit/imp_fri.F
!||====================================================================
      MODULE imp_frk
C
      INTEGER nrmax,n_frnn,nddlfr,nzkfr,ndfrmax,len_v,len_k
C     NDDLFRB :nb de noeuds fr. partages par au moins 3 proc,----
C  50d1m     DDLP0 :dernier ddl de noeuds fr.avec procs precedents ----
C  50d1m     DDLP1 :premier ddl de noeuds fr.avec procs suivants ----
      INTEGER nddlfrb,nddlfrb1,ddlp0,ddlp1
C
C tableau de int de : connectivite des noeuds FR. locaux----
      INTEGER, DIMENSION(:), ALLOCATABLE :: ifrloc
      INTEGER, DIMENSION(:), ALLOCATABLE :: iad_rl,fr_icol
C tableau de int de taille NKINE indice loc->global: ----
      INTEGER, DIMENSION(:), ALLOCATABLE :: ikin2g
C tableau de int de : matrice locale des noeuds FR. locaux par proc----
      INTEGER, DIMENSION(:), ALLOCATABLE :: nd_fr,iadfr,jdifr,iddlfr
#if defined knf
!dec$ attributes offload:mic :: IFR2K
#endif
      INTEGER, DIMENSION(:), ALLOCATABLE :: ifr2k,jfr2k
C 
      END MODULE imp_frk
C