Changeset 1013 for palm/trunk/SOURCE/poisfft_hybrid.f90
- Timestamp:
- Sep 21, 2012 7:03:55 AM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/poisfft_hybrid.f90
r810 r1013 2 2 3 3 !------------------------------------------------------------------------------ 4 !5 4 ! Current revisions: 6 5 ! ----------------- 7 ! 6 ! FLOAT type conversion replaced by REAL 8 7 ! 9 8 ! Former revisions: … … 914 913 IF ( i >= 0 .AND. i < nnxh ) THEN 915 914 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / & 916 FLOAT( nx+1 ) ) ) / ( dx * dx ) + &915 REAL( nx+1 ) ) ) / ( dx * dx ) + & 917 916 2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / & 918 FLOAT( ny+1 ) ) ) / ( dy * dy )917 REAL( ny+1 ) ) ) / ( dy * dy ) 919 918 ELSEIF ( i == nnxh ) THEN 920 919 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / & 921 FLOAT( nx+1 ) ) ) / ( dx * dx ) + &920 REAL( nx+1 ) ) ) / ( dx * dx ) + & 922 921 2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / & 923 FLOAT(ny+1) ) ) / ( dy * dy )922 REAL(ny+1) ) ) / ( dy * dy ) 924 923 ELSE 925 924 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / & 926 FLOAT( nx+1 ) ) ) / ( dx * dx ) + &925 REAL( nx+1 ) ) ) / ( dx * dx ) + & 927 926 2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / & 928 FLOAT( ny+1 ) ) ) / ( dy * dy )927 REAL( ny+1 ) ) ) / ( dy * dy ) 929 928 ENDIF 930 929 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.