36
37
38
39 USE elbufdef_mod
40 use my_alloc_mod
41
42
43
44#include "implicit_f.inc"
45
46
47
48#include "com01_c.inc"
49#include "sphcom.inc"
50#include "param_c.inc"
51#include "task_c.inc"
52#include "spmd_c.inc"
53
54
55
56 INTEGER IPARG(NPARG,*),EL2FA(*),NBF,IOFF(*),
57 . NBPART, IADG(NSPMD,*), SPH2SOL(*), SIOFF
58 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
59
60
61
62
63 INTEGER I, NG, NEL, NFT, IAD, ITY, LFT,
64 . N, LLT, MLW, NB1,
65 . NN1,NN2,RBUF
66 INTEGER, DIMENSION(:),ALLOCATABLE :: IOFFBUF
67 INTEGER, DIMENSION(:),ALLOCATABLE :: IOFF_LOC
68 TYPE(G_BUFEL_) ,POINTER :: GBUF
69
70 CALL my_alloc(ioffbuf,nbf)
71 CALL my_alloc(ioff_loc,sioff)
72 nn1 = 1
73 nn2 = nn1 + (numsph+maxpjet)
74
75
76 DO 490 ng=1,ngroup
77 mlw =iparg(1,ng)
78 nel =iparg(2,ng)
79 ity =iparg(5,ng)
80 nft =iparg(3,ng)
81 iad =iparg(4,ng)
82 lft = 1
83 llt = nel
84
85 IF (ity == 51) THEN
86
87
88
89 IF (mlw == 0) THEN
90 DO i=lft,llt
91 n = nft + i
92 ioff(el2fa(nn1+n)) = 0
93 ENDDO
94 ELSEIF (nsphsol/=0) THEN
95
96 gbuf => elbuf_tab(ng)%GBUF
97 DO i=lft,llt
98 n = nft + i
99 ioff(el2fa(nn1+n)) = nint(
min(gbuf%OFF(i),one))
100 IF (sph2sol(n)==0) THEN
101 ioff_loc(el2fa(nn1+n)) = ioff(el2fa(nn1+n))
102 ELSE
103 ioff_loc(el2fa(nn1+n)) = 0
104 ENDIF
105 ENDDO
106 ELSE
107 gbuf => elbuf_tab(ng)%GBUF
108 DO i=lft,llt
109 n = nft + i
110 ioff(el2fa(nn1+n)) = nint(
min(gbuf%OFF(i),one))
111 ioff_loc(el2fa(nn1+n)) = ioff(el2fa(nn1+n))
112 ENDDO
113 ENDIF
114 ELSE
115 ENDIF
116
117 490 CONTINUE
118
119 IF (nspmd == 1) THEN
121 ELSE
122 DO i = 1, nbf
123 ioffbuf(i) = ioff_loc(i)
124 ENDDO
125
126 IF (ispmd == 0) THEN
127 rbuf = numsphg
128 ELSE
129 rbuf = 1
130 ENDIF
131
133
134 ENDIF
135
136 DEALLOCATE(ioffbuf)
137 DEALLOCATE(ioff_loc)
138 RETURN
subroutine spmd_iget_partn(size, nbf_l, np, nbpart, iadg, srbuf, iflag)
void write_c_c(int *w, int *len)