"Fossies" - the Fresh Open Source Software Archive  

Source code changes of the file "src/ptests/pcgns_ftest.F90" between
CGNS-4.1.2.tar.gz and CGNS-4.2.0.tar.gz

About: The CGNS (CFD General Notation System) provides a standard for recording and recovering computer data associated with the numerical solution of fluid dynamics equations.

pcgns_ftest.F90  (CGNS-4.1.2):pcgns_ftest.F90  (CGNS-4.2.0)
skipping to change at line 13 skipping to change at line 13
USE mpi USE mpi
USE ISO_C_BINDING USE ISO_C_BINDING
USE CGNS USE CGNS
IMPLICIT NONE IMPLICIT NONE
#ifdef WINNT #ifdef WINNT
INCLUDE 'cgnswin_f.h' INCLUDE 'cgnswin_f.h'
#endif #endif
INTEGER(cgsize_t), PARAMETER :: totcnt = 40320 * 10 INTEGER(cgsize_t), PARAMETER :: totcnt = 40320 * 10
INTEGER, PARAMETER :: NLOOPS = 5000
INTEGER(cgsize_t) npp INTEGER(cgsize_t) npp
INTEGER(C_INT) commsize, commrank, mpi_err INTEGER(C_INT) commsize, commrank, mpi_err
INTEGER i, nb, nz, nerrs INTEGER i, nb, nz, nerrs
INTEGER ierr, F, B, Z, E, S INTEGER ierr, F, B, Z, E, S
INTEGER Cx, Cy, Cz, Fx, Fy, Fz, Ax, Ay, Az INTEGER Cx, Cy, Cz, Fx, Fy, Fz, Ax, Ay, Az
INTEGER(cgsize_t) sizes(3), start, END, n, j INTEGER(cgsize_t) sizes(3), start, END, n, j
INTEGER(cgsize_t), PARAMETER :: start_1 = 1 INTEGER(cgsize_t), PARAMETER :: start_1 = 1
REAL*8 ts, te, tt, dsize REAL*8 ts, te, tt, dsize
REAL*8 dx(totcnt), dy(totcnt), dz(totcnt) REAL*8 dx(totcnt), dy(totcnt), dz(totcnt)
INTEGER(cgsize_t), ALLOCATABLE, DIMENSION(:) :: ie INTEGER(cgsize_t), ALLOCATABLE, DIMENSION(:) :: ie
CHARACTER*32 name CHARACTER*32 name
CHARACTER*11 piomode(2) CHARACTER*11 piomode(2)
INTEGER :: istat INTEGER :: istat
INTEGER :: precision INTEGER :: precision
INTEGER, TARGET :: value
DATA piomode /'independent','collective'/ DATA piomode /'independent','collective'/
CALL MPI_INIT(mpi_err) CALL MPI_INIT(mpi_err)
CALL MPI_COMM_SIZE(MPI_COMM_WORLD,commsize,mpi_err) CALL MPI_COMM_SIZE(MPI_COMM_WORLD,commsize,mpi_err)
CALL MPI_COMM_RANK(MPI_COMM_WORLD,commrank,mpi_err) CALL MPI_COMM_RANK(MPI_COMM_WORLD,commrank,mpi_err)
IF (commsize .GT. 8) THEN IF (commsize .GT. 8) THEN
IF (commrank .EQ. 0) & IF (commrank .EQ. 0) &
PRINT *, 'number of processes must be 8 or less' PRINT *, 'number of processes must be 8 or less'
skipping to change at line 78 skipping to change at line 79
IF (commrank .EQ. 0) THEN IF (commrank .EQ. 0) THEN
PRINT *,'number processes =', commsize PRINT *,'number processes =', commsize
PRINT *,'array size per process =', npp PRINT *,'array size per process =', npp
PRINT *,'total array size =', totcnt PRINT *,'total array size =', totcnt
PRINT *,'total Mb for all data =', dsize PRINT *,'total Mb for all data =', dsize
ENDIF ENDIF
! default is MPI_COMM_WORLD, but can set another communicator with this ! default is MPI_COMM_WORLD, but can set another communicator with this
! call cgp_mpi_comm_f(MPI_COMM_WORLD,ierr) ! call cgp_mpi_comm_f(MPI_COMM_WORLD,ierr)
! Check repeated opening and closing of a file to detect issues with
! missed closed HDF5 objects, CGNS-109.
DO n = 1, NLOOPS
CALL cgp_open_f('pcgns_ftest.cgns',CG_MODE_WRITE,F,ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
CALL cgp_close_f(F,ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
ENDDO
CALL cgp_open_f('pcgns_ftest.cgns',CG_MODE_WRITE,F,ierr) CALL cgp_open_f('pcgns_ftest.cgns',CG_MODE_WRITE,F,ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
DO nb=1,2 DO nb=1,2
WRITE(name,'(a4,i2)') 'Base',nb WRITE(name,'(a4,i2)') 'Base',nb
CALL cg_base_write_f(F,name,3,3,B,ierr) CALL cg_base_write_f(F,name,3,3,B,ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
CALL cgp_pio_mode_f(INT(nb-1,C_INT), ierr) CALL cgp_pio_mode_f(INT(nb-1,C_INT), ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
skipping to change at line 248 skipping to change at line 258
IF (commrank .EQ. 0) THEN IF (commrank .EQ. 0) THEN
PRINT *,'read :',tt,' secs,',dsize/tt,' Mb/sec (', & PRINT *,'read :',tt,' secs,',dsize/tt,' Mb/sec (', &
piomode(B),') errors =',nerrs piomode(B),') errors =',nerrs
ENDIF ENDIF
ENDDO ENDDO
CALL cgp_close_f(F,ierr) CALL cgp_close_f(F,ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
! Disable with gfortran, GCC Bugzilla - Bug 99982
#ifndef __GFORTRAN__
! test cg_configure_f
value = MPI_COMM_SELF
CALL cg_configure_f(CG_CONFIG_HDF5_MPI_COMM, C_LOC(value), ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
IF (commrank .EQ. 0) THEN
CALL cgp_open_f('pcgns_ftest.cgns',CG_MODE_READ,F,ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
CALL cgp_close_f(F,ierr)
IF (ierr .NE. CG_OK) CALL cgp_error_exit_f
ENDIF
#endif
CALL MPI_FINALIZE(mpi_err) CALL MPI_FINALIZE(mpi_err)
END PROGRAM pcgns_ftest END PROGRAM pcgns_ftest
 End of changes. 4 change blocks. 
1 lines changed or deleted 25 lines changed or added

Home  |  About  |  Features  |  All  |  Newest  |  Dox  |  Diffs  |  RSS Feeds  |  Screenshots  |  Comments  |  Imprint  |  Privacy  |  HTTP(S)