@@ -100,11 +100,12 @@ SUBROUTINE tce_mrcc_ccsdpt_subg(d_t1,k_t1_offset,
100
100
c
101
101
c MIC
102
102
c
103
+ #ifdef USE_OFFLOAD
103
104
triplesx_alloced= .false.
104
- triplesx1_alloced= .false.
105
105
v2sub_alloced= .false.
106
106
t1sub_alloced= .false.
107
107
t2sub_alloced= .false.
108
+ #endif
108
109
c estimate triplesx size
109
110
range_p4= 0
110
111
do t_p4b = noab+1 ,noab+ nvab
@@ -115,22 +116,16 @@ SUBROUTINE tce_mrcc_ccsdpt_subg(d_t1,k_t1_offset,
115
116
range_h1 = max (range_h1,int_mb(k_range+ t_h1b-1 ))
116
117
enddo
117
118
size= (range_p4** 3 )* (range_h1** 3 )
118
- triplesx_copyback= .true.
119
- c call util_align64(size)
120
- triplesx_mxlgth= size
121
119
if (.not. MA_PUSH_GET(mt_dbl,size,' (T) singles' ,l_singles,
122
120
1 k_singles)) call errquit(' ccsd_t: MA error sgl' ,
123
121
2 size,MA_ERR)
124
122
if (.not. MA_PUSH_GET(mt_dbl,size,' (T) doubles' ,l_doubles,
125
123
1 k_doubles)) call errquit(' ccsd_t: MA error dbl' ,
126
124
2 size,MA_ERR)
127
125
#ifdef USE_OFFLOAD
128
- #if 0
129
- call offl_alloc(dbl_mb(k_singles),size)
130
- triplesx1_alloced=.true.
131
- #else
132
- triplesx1_alloced= .false.
133
- #endif
126
+ triplesx_copyback= .true.
127
+ c call util_align64(size)
128
+ triplesx_mxlgth= size
134
129
call offl_alloc(dbl_mb(k_doubles),size)
135
130
triplesx_alloced= .true.
136
131
#endif
@@ -187,12 +182,11 @@ SUBROUTINE tce_mrcc_ccsdpt_subg(d_t1,k_t1_offset,
187
182
c c call dcopy(size, 0.0d0, 0, dbl_mb(k_doubles), 1)
188
183
c c call offl_zerofill(dbl_mb(k_doubles),size)
189
184
call dcopy(triplesx_mxlgth, 0.0d0 , 0 , dbl_mb(k_singles), 1 )
190
- #if 0
191
- call offl_zerofill(dbl_mb(k_singles),triplesx_mxlgth)
192
- #endif
193
185
194
186
call dcopy(triplesx_mxlgth, 0.0d0 , 0 , dbl_mb(k_doubles), 1 )
187
+ #ifdef USE_OFFLOAD
195
188
call offl_zerofill(dbl_mb(k_doubles),triplesx_mxlgth)
189
+ #endif
196
190
c -----------
197
191
c
198
192
c cx if (.not.MA_PUSH_GET(mt_dbl,size,'(T) singles',l_singles,
@@ -385,7 +379,6 @@ SUBROUTINE tce_mrcc_ccsdpt_subg(d_t1,k_t1_offset,
385
379
triplesx_alloced= .false.
386
380
#if 0
387
381
call offl_free(dbl_mb(k_singles),size)
388
- triplesx1_alloced=.false.
389
382
#endif
390
383
#endif
391
384
0 commit comments