111
|
1 ! PR debug/71642
|
|
2 ! { dg-do compile }
|
|
3 ! { dg-options "-g" }
|
|
4
|
|
5 MODULE gauss_colloc
|
|
6 INTEGER, PARAMETER :: dp=8
|
|
7 CONTAINS
|
|
8 SUBROUTINE collocGauss(h,h_inv,grid,poly,alphai,posi,max_r2,&
|
|
9 periodic,gdim,local_bounds,local_shift,poly_shift,scale,lgrid,error)
|
|
10 REAL(dp), DIMENSION(0:, 0:, 0:), &
|
|
11 INTENT(inout) :: grid
|
|
12 INTEGER, INTENT(inout), OPTIONAL :: lgrid
|
|
13 CONTAINS
|
|
14 SUBROUTINE kloop6
|
|
15 IF (kJump/=1 .AND. (ikstart+kmax-kstart>=ndim(2)+l_shift(2) .OR.&
|
|
16 ikstart2+kmin-kstart2<=l_ub(2)-ndim(2))) THEN
|
|
17 DO
|
|
18 DO k=kstart2,kend2,-1
|
|
19 IF ( PRESENT ( lgrid ) ) THEN
|
|
20 grid(ik,ij,ii) = grid(ik,ij,ii) + p_v*res_k
|
|
21 END IF
|
|
22 END DO
|
|
23 END DO
|
|
24 END IF
|
|
25 END SUBROUTINE
|
|
26 END SUBROUTINE
|
|
27 END MODULE gauss_colloc
|