909
909
C also done with J, so that the output I() are always within range.
910
910
C If IOPT=0, nothing is done.
911
911
C *********** USAGE **************************************************
913
913
C with periodic boundary conditions in a space of variable dimension
914
914
C SUBROUTINE LAPLACIAN( ND, N, DX, F, FLAPL )