Commit 3ed184fd authored by Andreas Marek's avatar Andreas Marek

Rename check_correctness function

parent 5748ad54
......@@ -274,7 +274,7 @@ int main(int argc, char** argv) {
}
/* check the results */
status = check_correctness_complex_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_complex_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
if (status !=0){
printf("The computed EVs are not correct !\n");
......
......@@ -278,7 +278,7 @@ int main(int argc, char** argv) {
/* check the results */
status = check_correctness_real_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_real_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
if (status !=0){
printf("The computed EVs are not correct !\n");
......
......@@ -274,7 +274,7 @@ int main(int argc, char** argv) {
}
/* check the results */
status = check_correctness_complex_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_complex_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
if (status !=0){
printf("The computed EVs are not correct !\n");
......
......@@ -277,7 +277,7 @@ int main(int argc, char** argv) {
}
/* check the results */
status = check_correctness_real_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_real_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
if (status !=0){
printf("The computed EVs are not correct !\n");
......
......@@ -241,9 +241,9 @@ int main(int argc, char** argv) {
/* check the results */
#ifdef DOUBLE_PRECISION_COMPLEX
status = check_correctness_complex_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_complex_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#else
status = check_correctness_complex_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_complex_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#endif
if (status !=0){
......
......@@ -233,9 +233,9 @@ int main(int argc, char** argv) {
#ifdef DOUBLE_PRECISION_REAL
/* check the results */
status = check_correctness_real_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_real_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#else
status = check_correctness_real_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_real_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#endif
if (status !=0){
printf("The computed EVs are not correct !\n");
......
......@@ -236,9 +236,9 @@ int main(int argc, char** argv) {
/* check the results */
#ifdef DOUBLE_PRECISION_COMPLEX
status = check_correctness_complex_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_complex_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#else
status = check_correctness_complex_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_complex_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#endif
if (status !=0){
printf("The computed EVs are not correct !\n");
......
......@@ -230,9 +230,9 @@ int main(int argc, char** argv) {
/* check the results */
#ifdef DOUBLE_PRECISION_REAL
status = check_correctness_real_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_real_double_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#else
status = check_correctness_real_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
status = check_correctness_evp_numeric_residuals_real_single_f(na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid);
#endif
if (status !=0){
......
......@@ -323,7 +323,7 @@ program test_complex2
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -315,7 +315,7 @@ program test_real2
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -323,7 +323,7 @@ program test_complex2
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -314,7 +314,7 @@ program test_real2
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -271,7 +271,7 @@ program test_complex_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -270,7 +270,7 @@ program test_complex_gpu_version_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -264,7 +264,7 @@ program test_real_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -264,7 +264,7 @@ program test_real_gpu_version_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -270,7 +270,7 @@ program test_complex_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -270,7 +270,7 @@ program test_complex_gpu_version_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -263,7 +263,7 @@ program test_real_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -264,7 +264,7 @@ program test_real_gpu_version_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -326,7 +326,7 @@ program test_solve_tridi_single
endif
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, a, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, a, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -327,7 +327,7 @@ program test_solve_tridi
endif
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, a, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, a, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
......
......@@ -253,7 +253,7 @@ program test_complex2_double_banded
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -193,14 +193,14 @@ program test_interface
call elpa_deallocate(e2)
call elpa_uninit()
status = check_correctness(na, nev, as1, z1, ev1, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as1, z1, ev1, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a1)
deallocate(as1)
deallocate(z1)
deallocate(ev1)
status = check_correctness(na, nev, as2, z2, ev2, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as2, z2, ev2, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a2)
deallocate(as2)
......
......@@ -265,7 +265,7 @@ program test_complex2_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -397,7 +397,7 @@ program test_complex2_choose_kernel_with_api_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -293,7 +293,7 @@ program test_complex2_default_kernel_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -292,7 +292,7 @@ program test_complex2_gpu_version_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -270,7 +270,7 @@ program test_real2_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -432,7 +432,7 @@ program test_real2_choose_kernel_with_api_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -290,7 +290,7 @@ program test_real2_default_kernel_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -291,7 +291,7 @@ program test_real2_gpu_version_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -334,7 +334,7 @@ program test_real2_default_kernel_qr_decomposition_double_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
ret = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
ret = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -264,7 +264,7 @@ program test_complex2_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -397,7 +397,7 @@ program test_complex2_choose_kernel_with_api_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -290,7 +290,7 @@ program test_complex2_default_kernel_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -289,7 +289,7 @@ program test_complex2_gpu_version_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -269,7 +269,7 @@ program test_real2_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -428,7 +428,7 @@ program test_real2_choose_kernel_with_api_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -289,7 +289,7 @@ program test_real2_default_kernel_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -292,7 +292,7 @@ program test_real2_gpu_version_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -338,7 +338,7 @@ program test_real2_default_kernel_qr_decomposition_single_precision
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
ret = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
ret = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -250,7 +250,7 @@ program test_real2_double_banded
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
......
......@@ -252,7 +252,7 @@ program test_complex2_single_banded
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -248,7 +248,7 @@ program test_real2_single_banded
!-------------------------------------------------------------------------------
! Test correctness of result (using plain scalapack routines)
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
deallocate(a)
deallocate(as)
......
......@@ -1013,10 +1013,10 @@ program test_all_real
if (input_options%datatype .eq. 1) then
status = check_correctness(na, nev, as_real, z_real, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_real, z_real, ev, sc_desc, myid)
endif
if (input_options%datatype .eq. 2) then
status = check_correctness(na, nev, as_complex, z_complex, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_complex, z_complex, ev, sc_desc, myid)
endif
if (status .eq. 1) then
......@@ -1099,7 +1099,7 @@ program test_all_real
elpa_int_value_to_string("real_kernel", ELPA_2STAGE_REAL_DEFAULT),' default kernel:',tEnd - tStart
if (myid == 0) print *," "
status = check_correctness(na, nev, as_real, z_real, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_real, z_real, ev, sc_desc, myid)
if (myid == 0) print *," "
if (status .eq. 1) then
......@@ -1189,7 +1189,7 @@ program test_all_real
trim(elpa_real_kernel_name(this_kernel)),' kernel:',tEnd - tStart
if (myid == 0) print *," "
status = check_correctness(na, nev, as_real, z_real, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_real, z_real, ev, sc_desc, myid)
if (myid == 0) print *," "
if (status .eq. 1) then
......@@ -1283,7 +1283,7 @@ program test_all_real
trim(elpa_real_kernel_name(input_options%this_real_kernel)),' kernel:',tEnd - tStart
if (myid == 0) print *," "
status = check_correctness(na, nev, as_real, z_real, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_real, z_real, ev, sc_desc, myid)
if (myid == 0) print *," "
if (status .eq. 1) then
......@@ -1380,7 +1380,7 @@ program test_all_real
elpa_int_value_to_string("complex_kernel", ELPA_2STAGE_COMPLEX_DEFAULT),' default kernel:',tEnd - tStart
if (myid == 0) print *," "
status = check_correctness(na, nev, as_complex, z_complex, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_complex, z_complex, ev, sc_desc, myid)
if (myid == 0) print *," "
if (status .eq. 1) then
......@@ -1471,7 +1471,7 @@ program test_all_real
trim(elpa_complex_kernel_name(this_kernel)),' kernel:',tEnd - tStart
if (myid == 0) print *," "
status = check_correctness(na, nev, as_complex, z_complex, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_complex, z_complex, ev, sc_desc, myid)
if (myid == 0) print *," "
if (status .eq. 1) then
......@@ -1568,7 +1568,7 @@ program test_all_real
trim(elpa_complex_kernel_name(input_options%this_complex_kernel)),' kernel:',tEnd - tStart
if (myid == 0) print *," "
status = check_correctness(na, nev, as_complex, z_complex, ev, sc_desc, myid)
status = check_correctness_evp_numeric_residuals(na, nev, as_complex, z_complex, ev, sc_desc, myid)
if (myid == 0) print *," "
if (status .eq. 1) then
......
......@@ -503,8 +503,11 @@ program test
#ifdef TEST_MATRIX_ANALYTIC
status = check_correctness_analytic(na, nev, ev, z, nblk, myid, np_rows, np_cols, my_prow, my_pcol, check_all_evals)
#else
!#elif defined(TEST_MATRIX_FRANK)
! status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows,np_cols, my_prow, my_pcol)
!#elif defined(TEST_MATRIX_RANDOM)
if (nev .ge. 1) then
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows,np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows,np_cols, my_prow, my_pcol)
else
! zero eigenvectors and no analytic test => toeplitz
status = check_correctness_eigenvalues_toeplitz(na, diagonalElement, &
......@@ -521,7 +524,7 @@ program test
#ifdef TEST_SOLVE_TRIDIAGONAL
! check eigenvectors
status = check_correctness(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
status = check_correctness_evp_numeric_residuals(na, nev, as, z, ev, sc_desc, nblk, myid, np_rows, np_cols, my_prow, my_pcol)
call check_status(status, myid)
#endif
#endif
......
......@@ -45,14 +45,14 @@
module test_check_correctness
use test_util
interface check_correctness
module procedure check_correctness_complex_double
module procedure check_correctness_real_double
interface check_correctness_evp_numeric_residuals
module procedure check_correctness_evp_numeric_residuals_complex_double
module procedure check_correctness_evp_numeric_residuals_real_double
#ifdef WANT_SINGLE_PRECISION_REAL
module procedure check_correctness_real_single
module procedure check_correctness_evp_numeric_residuals_real_single
#endif
#ifdef WANT_SINGLE_PRECISION_COMPLEX
module procedure check_correctness_complex_single
module procedure check_correctness_evp_numeric_residuals_complex_single
#endif
end interface
......
......@@ -41,7 +41,7 @@
!
! Author: A. Marek, MPCDF
function check_correctness_&
function check_correctness_evp_numeric_residuals_&
&MATH_DATATYPE&
&_&
&PRECISION&
......@@ -323,11 +323,11 @@
#if REALCASE == 1
#ifdef DOUBLE_PRECISION_REAL
!c> int check_correctness_real_double_f(int na, int nev, int na_rows, int na_cols,
!c> int check_correctness_evp_numeric_residuals_real_double_f(int na, int nev, int na_rows, int na_cols,
!c> double *as, double *z, double *ev,
!c> int sc_desc[9], int myid);
#else
!c> int check_correctness_real_single_f(int na, int nev, int na_rows, int na_cols,
!c> int check_correctness_evp_numeric_residuals_real_single_f(int na, int nev, int na_rows, int na_cols,
!c> float *as, float *z, float *ev,
!c> int sc_desc[9], int myid);
#endif
......@@ -336,22 +336,22 @@
#if COMPLEXCASE == 1
#ifdef DOUBLE_PRECISION_COMPLEX
!c> int check_correctness_complex_double_f(int na, int nev, int na_rows, int na_cols,
!c> int check_correctness_evp_numeric_residuals_complex_double_f(int na, int nev, int na_rows, int na_cols,
!c> complex double *as, complex double *z, double *ev,
!c> int sc_desc[9], int myid);
#else
!c> int check_correctness_complex_single_f(int na, int nev, int na_rows, int na_cols,
!c> int check_correctness_evp_numeric_residuals_complex_single_f(int na, int nev, int na_rows, int na_cols,
!c> complex float *as, complex float *z, float *ev,
!c> int sc_desc[9], int myid);
#endif
#endif /* COMPLEXCASE */
function check_correctness_&
function check_correctness_evp_numeric_residuals_&
&MATH_DATATYPE&
&_&
&PRECISION&
&_f (na, nev, na_rows, na_cols, as, z, ev, sc_desc, myid) result(status) &
bind(C,name="check_correctness_&
bind(C,name="check_correctness_evp_numeric_residuals_&
&MATH_DATATYPE&
&_&
&PRECISION&
......@@ -371,7 +371,7 @@ function check_correctness_&
! TODO: I did not want to add all the variables to the C interface as well
! TODO: I think that we should find a better way to pass this information
! TODO: to all the functions anyway (get it from sc_desc, pass elpa_t, etc..)
status = check_correctness_&
status = check_correctness_evp_numeric_residuals_&
&MATH_DATATYPE&
&_&
&PRECISION&
......