1!$Id:$ 2 subroutine formfe(pnu,pnb,pna,pnl,aufl,bfl,alfl,dfl, 3 & isw,nl1,nl2,nl3) 4 5! * * F E A P * * A Finite Element Analysis Program 6 7!.... Copyright (c) 1984-2017: Regents of the University of California 8! All rights reserved 9 10!-----[--.----+----.----+----.-----------------------------------------] 11! Purpose: Forms finite element arrays as required 12 13! Inputs: 14! pnu - Pointer to current nodal solution vectors 15! pnb - Pointer for FEM Vector 16! pna - Pointer for FEM diagonal and upper part of array 17! pnl - Pointer for FEM lower part of array 18! aufl - If true assemble 'a' array (which includes 'au') 19! bfl - If true assemble 'b' array 20! alfl - If true assemble 'al' array (feap only) 21! dfl - If true assembel 'b' uncompressed 22! isw - Solution switch controlling action to be taken 23! nl1 - First element to be processed 24! nl2 - Last element to be processed 25! nl3 - Increment to 'nl1' 26 27! Outputs: 28! hr(pnb)- Values for FEM Vector 29! hr(pna)- Values for FEM diagonal and upper array part 30! hr(pnl)- Values for FEM lower array part 31!-----[--.----+----.----+----.-----------------------------------------] 32 implicit none 33 34 include 'cdata.h' 35 include 'cdat1.h' 36 include 'hdatam.h' 37 include 'ndata.h' 38 include 'p_formfe.h' 39 include 'pointer.h' 40 include 'prflag.h' 41 include 'sdata.h' 42 include 'comblk.h' 43 44 logical aufl,bfl,alfl,dfl 45 integer isw, nl1, nl2, nl3 46 47 save 48 49! Form appropriate finite element arrays 50 51 call pform(hr(np(41)),hr(np(44)),hr(np(39)),mr(np(34)),hr(np(35)), 52 & hr(np(36)),mr(np(32)),hr(np(25)),mr(np(31)),hr(np(43)), 53 & mr(np(33)),hr(np(30)),hr(np(38)), 54 & mr(np(21)),hr(pnu),hr(np(42)),hr(pnb),hr(pna), 55 & hr(pnl),ndd,nie,ndf,ndm,nen1,nst,aufl,bfl,dfl, 56 & isw,nl1,nl2,nl3) 57 58! Reset update flag for history variables 59 60 hflgu = .false. 61 h3flgu = .false. 62 63 end 64