6 INTEGER, PARAMETER :: dp
=8
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:), &
12 INTEGER, INTENT(inout
), OPTIONAL
:: lgrid
15 IF (kJump
/=1 .AND
. (ikstart
+kmax
-kstart
>=ndim(2)+l_shift(2) .OR
.&
16 ikstart2
+kmin
-kstart2
<=l_ub(2)-ndim(2))) THEN
19 IF ( PRESENT ( lgrid
) ) THEN
20 grid(ik
,ij
,ii
) = grid(ik
,ij
,ii
) + p_v
*res_k
27 END MODULE gauss_colloc