Lines Matching refs:amg_la
32 real(kind=8),dimension(amg_nshg(level1)) :: amg_la local
245 amg_la(n) = amg_A_lhs(level1)%p(k,1)
257 rtp = amg_la(aLoc(j))/rtp
267 amg_la(m) = 0
270 amg_la(m) = amg_la(m) - rtp*amg_A_lhs(level1)%p(kj,1)
292 if ( amg_la(k).lt.0 ) then ! a-hat < 0
293 alpha = alpha + amg_la(k)
295 beta = beta + amg_la(k)
300 if (amg_la(k).lt.0) then ! a-hat < 0 in coarse
301 alphac = alphac + amg_la(k)
303 betac = betac + amg_la(k)
306 amg_la(I_nnz) = amg_la(k)
313 if ( amg_la(k).lt.0) then
314 amg_la(k) = amg_la(k)*alpha
316 amg_la(k) = amg_la(k)*beta
318 rtp = max(rtp,abs(amg_la(k)))
329 tmp3 = tmp3 + min(amg_la(k),zero)
330 tmp4 = tmp4 + max(amg_la(k),zero)
331 if (abs(amg_la(k)).gt.rtp) then
332 tmp1 = tmp1 + min(amg_la(k),zero)
333 tmp2 = tmp2 + max(amg_la(k),zero)
336 amg_la(I_nnz) = amg_la(k)
342 if (amg_la(k).lt.0) then
343 amg_la(k) = amg_la(k)*tmp1
345 amg_la(k) = amg_la(k)*tmp2
359 if (amg_la(k).lt.0) then
360 amg_I%pp(i)%p(k) = amg_la(k)
362 amg_I%pp(i)%p(k) = amg_la(k)