Skip to content

Commit

Permalink
Merge pull request #282 from AlexanderRichert-NOAA/fftpack_tests_feb25
Browse files Browse the repository at this point in the history 
Extended fftpack testing
  • Loading branch information
@AlexanderRichert-NOAA
AlexanderRichert-NOAA authored Feb 19, 2025
2 parents 611836d + 8bb4af7 commit 86fe9a4
Show file tree
Hide file tree
Showing 8 changed files with 681 additions and 33 deletions.
89 changes: 58 additions & 31 deletions tests/CMakeLists.txt
View file
Original file line number Diff line number Diff line change
@@ -1,4 +1,7 @@
# This is the CMake file for the test directory of NCEPLIBS-ip.
# This is the CMake file for the tests/ directory of NCEPLIBS-ip.
# It define dozens of CTest unit tests. All new and modified code
# should have one or more unit tests to verify correct operation,
# including validation of outputs against known cases.
#
# Alex Richert, Mark Potts, Kyle Gerheiser, Eric Engle

Expand Down Expand Up @@ -53,7 +56,7 @@ function(create_sp_test name kind timeout)
set_target_properties(${name}_${kind} PROPERTIES COMPILE_FLAGS "-fconvert=big-endian ${fortran_${kind}_flags}")
endif()
add_test(NAME ${name}_${kind} COMMAND ${name}_${kind})
target_compile_definitions(${name}_${kind} PRIVATE KIND_${kind})
target_compile_definitions(${name}_${kind} PRIVATE "LSIZE=${kind}")
if(TEST_TIME_LIMIT)
set_tests_properties(${name}_${kind} PROPERTIES TIMEOUT ${timeout})
endif()
Expand All @@ -69,29 +72,31 @@ endif()
foreach(kind ${kinds})
string(TOUPPER ${kind} kind_definition)

# Test ipxwafs routines
### Test ipxwafs routines
add_executable(test_ipxwafs_${kind} test_ipxwafs.F90)
target_link_libraries(test_ipxwafs_${kind} PUBLIC ip::ip_${kind})
target_compile_definitions(test_ipxwafs_${kind} PRIVATE "LSIZE=${kind_definition}")
set_target_properties(test_ipxwafs_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
add_test(test_ipxwafs_${kind} test_ipxwafs_${kind})

# Test earth_radius_mod.
### Test earth_radius_mod.
add_executable(test_earth_radius_${kind} test_earth_radius.F90)
target_link_libraries(test_earth_radius_${kind} PUBLIC ip::ip_${kind})
target_compile_definitions(test_earth_radius_${kind} PRIVATE "LSIZE=${kind_definition}")
set_target_properties(test_earth_radius_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
add_test(test_earth_radius_${kind} test_earth_radius_${kind})

# grib-2 tests

### GRIB2 tests
# Create mod files by kind, GRIB1
add_library(test_input_library_grib2_${kind} input_data_mod_grib2_${kind}.F90)
target_link_libraries(test_input_library_grib2_${kind} PUBLIC ip::ip_${kind})
target_compile_definitions(test_input_library_grib2_${kind} PUBLIC "LSIZE=${kind_definition}")
set_target_properties(test_input_library_grib2_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
add_library(test_library_grib2_${kind} interp_mod_grib2_${kind}.F90)
target_link_libraries(test_library_grib2_${kind} PUBLIC test_input_library_grib2_${kind})
set_target_properties(test_library_grib2_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")


# Set up executables, GRIB2
add_executable(tst_gdswzd_grib2_${kind} tst_gdswzd_grib2.c)
set_target_properties(tst_gdswzd_grib2_${kind} PROPERTIES LINKER_LANGUAGE C)
target_compile_definitions(tst_gdswzd_grib2_${kind} PRIVATE "LSIZE=${kind_definition}")
Expand All @@ -103,10 +108,11 @@ foreach(kind ${kinds})
target_link_libraries(test_vector_grib2_${kind} PRIVATE test_library_grib2_${kind})
set_target_properties(test_scalar_grib2_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
set_target_properties(test_vector_grib2_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")


# GDSWZD test, GRIB2
add_test(tst_gdswzd_c_grib2_${kind} tst_gdswzd_grib2_${kind})
# scalar tests

# Scalar tests, GRIB2
add_test(test_lambert_bilinear_scalar_grib2_${kind} test_scalar_grib2_${kind} 218 0)
add_test(test_gaussian_neighbor_scalar_grib2_${kind} test_scalar_grib2_${kind} 127 2)
add_test(test_latlon_bilinear_scalar_grib2_${kind} test_scalar_grib2_${kind} 3 0)
Expand All @@ -115,15 +121,15 @@ foreach(kind ${kinds})
add_test(test_rotatedB_spectral_scalar_grib2_${kind} test_scalar_grib2_${kind} 205 4)
add_test(test_rotatedE_budget_scalar_grib2_${kind} test_scalar_grib2_${kind} 203 3)

# scalar station point tests
# Scalar station point tests, GRIB2
add_test(test_station_points_bilinear_scalar_grib2_${kind} test_scalar_grib2_${kind} -1 0)
add_test(test_station_points_bicubic_scalar_grib2_${kind} test_scalar_grib2_${kind} -1 1)
add_test(test_station_points_neighbor_scalar_grib2_${kind} test_scalar_grib2_${kind} -1 2)
add_test(test_station_points_budget_scalar_grib2_${kind} test_scalar_grib2_${kind} -1 3)
add_test(test_station_points_spectral_scalar_grib2_${kind} test_scalar_grib2_${kind} -1 4)
add_test(test_station_points_neighbor_budget_scalar_grib2_${kind} test_scalar_grib2_${kind} -1 6)
# # vector tests

# Vector tests, GRIB2
add_test(test_lambert_bilinear_vector_grib2_${kind} test_vector_grib2_${kind} 218 0)
add_test(test_gaussian_neighbor_vector_grib2_${kind} test_vector_grib2_${kind} 127 2)
add_test(test_latlon_bilinear_vector_grib2_${kind} test_vector_grib2_${kind} 3 0)
Expand All @@ -135,23 +141,25 @@ foreach(kind ${kinds})
add_test(test_rotatedB_direct_ncep_post_spectral_vector_grib2_${kind} test_vector_grib2_${kind} 32769b 4)
add_test(test_rotatedE_direct_budget_vector_grib2_${kind} test_vector_grib2_${kind} 32768 3)

# vector station point tests
# Vector station point tests, GRIB2
add_test(test_station_points_bilinear_vector_grib2_${kind} test_vector_grib2_${kind} -1 0)
add_test(test_station_points_bicubic_vector_grib2_${kind} test_vector_grib2_${kind} -1 1)
add_test(test_station_points_neighbor_vector_grib2_${kind} test_vector_grib2_${kind} -1 2)
add_test(test_station_points_budget_vector_grib2_${kind} test_vector_grib2_${kind} -1 3)
add_test(test_station_points_spectral_vector_grib2_${kind} test_vector_grib2_${kind} -1 4)
add_test(test_station_points_neighbor_budget_vector_grib2_${kind} test_vector_grib2_${kind} -1 6)

# grib-1 tests

### GRIB1 tests
# Create mod files by kind, GRIB1
add_library(test_input_library_grib1_${kind} input_data_mod_grib1_${kind}.F90)
target_link_libraries(test_input_library_grib1_${kind} PUBLIC ip::ip_${kind})
target_compile_definitions(test_input_library_grib1_${kind} PUBLIC "LSIZE=${kind_definition}")
set_target_properties(test_input_library_grib1_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
add_library(test_library_grib1_${kind} interp_mod_grib1_${kind}.F90)
target_link_libraries(test_library_grib1_${kind} PUBLIC test_input_library_grib1_${kind})
set_target_properties(test_library_grib1_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")


# Set up executables, GRIB1
add_executable(tst_gdswzd_grib1_${kind} tst_gdswzd_grib1.c)
set_target_properties(tst_gdswzd_grib1_${kind} PROPERTIES LINKER_LANGUAGE C)
target_compile_definitions(tst_gdswzd_grib1_${kind} PRIVATE "LSIZE=${kind_definition}")
Expand All @@ -163,17 +171,20 @@ foreach(kind ${kinds})
target_link_libraries(test_vector_grib1_${kind} PRIVATE test_library_grib1_${kind})
set_target_properties(test_scalar_grib1_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
set_target_properties(test_vector_grib1_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")


# GDSWZD test, GRIB1
add_test(tst_gdswzd_c_grib1_${kind} tst_gdswzd_grib1_${kind})

# Scalar tests, GRIB1
add_test(test_lambert_bilinear_scalar_grib1_${kind} test_scalar_grib1_${kind} 218 0)
add_test(test_gaussian_neighbor_scalar_grib1_${kind} test_scalar_grib1_${kind} 127 2)
add_test(test_latlon_bilinear_scalar_grib1_${kind} test_scalar_grib1_${kind} 3 0)
add_test(test_mercator_bicubic_scalar_grib1_${kind} test_scalar_grib1_${kind} 8 1)
add_test(test_polar_stereo_neighbor_budget_scalar_grib1_${kind} test_scalar_grib1_${kind} 212 6)
add_test(test_rotatedB_spectral_scalar_grib1_${kind} test_scalar_grib1_${kind} 205 4)
add_test(test_rotatedE_budget_scalar_grib1_${kind} test_scalar_grib1_${kind} 203 3)
# vector tests

# Vector tests, GRIB1
add_test(test_lambert_bilinear_vector_grib1_${kind} test_vector_grib1_${kind} 218 0)
add_test(test_gaussian_neighbor_vector_grib1_${kind} test_vector_grib1_${kind} 127 2)
add_test(test_latlon_bilinear_vector_grib1_${kind} test_vector_grib1_${kind} 3 0)
Expand All @@ -182,6 +193,31 @@ foreach(kind ${kinds})
add_test(test_rotatedB_spectral_vector_grib1_${kind} test_vector_grib1_${kind} 205 4)
add_test(test_rotatedE_budget_vector_grib1_${kind} test_vector_grib1_${kind} 203 3)

### Various sp_mod tests
create_sp_test(test_ncpus ${kind} 0.3)
create_sp_test(test_splaplac ${kind} 0.3)
create_sp_test(test_splat ${kind} 0.3)
create_sp_test(test_sppad ${kind} 0.3)
create_sp_test(test_sptezv ${kind} 0.3)
create_sp_test(test_fft ${kind} 0.3)
create_sp_test(test_sptrung ${kind} 0.3)
create_sp_test(test_sptrungv ${kind} 2)

### Tests for fftpack subroutines
add_library(test_fftpack_mod_${kind} test_fftpack_mod_${kind}.F90)
target_link_libraries(test_fftpack_mod_${kind} PRIVATE ip::ip_${kind})
set_target_properties(test_fftpack_mod_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
add_executable(test_fftpack_${kind} test_fftpack.F90)
target_compile_definitions(test_fftpack_${kind} PRIVATE "LSIZE=${kind}")
set_target_properties(test_fftpack_${kind} PROPERTIES COMPILE_FLAGS "${fortran_${kind}_flags}")
target_link_libraries(test_fftpack_${kind} PRIVATE test_fftpack_mod_${kind})
target_include_directories(test_fftpack_${kind} PRIVATE ${CMAKE_CURRENT_BINARY_DIR})
target_include_directories(test_fftpack_${kind} PRIVATE ${CMAKE_CURRENT_BINARY_DIR}/../src/include_${kind})
add_test(test_fftpack_${kind} test_fftpack_${kind})

### Add test labels:

# Label slow-running tests so they can be skipped using 'ctest -LE/--label-exclude'
set(slowtests
"test_rotatedB_spectral_vector_grib1"
"test_rotatedB_spectral_scalar_grib1"
Expand All @@ -195,17 +231,8 @@ foreach(kind ${kinds})
set_tests_properties("${slowtest}_${kind}" PROPERTIES LABELS "SLOW_TEST")
endforeach()

# sp tests
create_sp_test(test_ncpus ${kind} 0.3)
create_sp_test(test_splaplac ${kind} 0.3)
create_sp_test(test_splat ${kind} 0.3)
create_sp_test(test_sppad ${kind} 0.3)
create_sp_test(test_sptezv ${kind} 0.3)
create_sp_test(test_fft ${kind} 0.3)
create_sp_test(test_sptrung ${kind} 0.3)
create_sp_test(test_sptrungv ${kind} 2)

# label tests without input data to allow trimmed down unit test runs
# Label tests without input data to allow trimmed down unit test runs
# using 'ctest -L NO_INPUT_DATA'.
set(noinputdatatests
"test_earth_radius"
"test_fft"
Expand Down
2 changes: 1 addition & 1 deletion tests/test_fft.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ program test_fft
use iso_fortran_env, only: real32, real64
implicit none

#ifdef KIND_d
#if(LSIZE==d)
real, parameter :: PI = 3.14159265358979
integer, parameter :: precision = real64
integer, parameter :: imax = 256
Expand Down
13 changes: 13 additions & 0 deletions tests/test_fftpack.F90
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
PROGRAM test_fftpack
#if(LSIZE==4)
use test_fftpack_mod_4
#elif(LSIZE==D)
use test_fftpack_mod_d
#elif(LSIZE==8)
use test_fftpack_mod_8
#endif
implicit none

call run_fftpack_tests()

END PROGRAM test_fftpack
Loading

0 comments on commit 86fe9a4

Please sign in to comment.